similar to: basic graphics

Dear friends. I've forgotten how to simulate data with a known correlation, e.g. only two columns. I'd be most pleased to be told. Thanks in advance Troels Ring, M.D Department of Nephrology Aalborg, Denmark tring at mail1.stofanet.dk -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read

Thank you to Peter Dalgaard, who provided help for some of this solution. It is still not ideal since the dots which are equidistant in the x-direction ideally should respect the reciprocal scale in the y-direction. op <- par("cex"= 0.7) x<-y<-NA plot(x,y,xlim=c(0,10),ylim=1/c(2000,70), axes=F) box() axis(1)

Dear Douglas Bates. I just downloaded the compiled version (I'm a poor Windows devil, not yet having found the time to move to a more advanced platform...) from NT- the files are dated 30.5-1999 so they are not old - and the problem persisted....wonder what I did wrong ? R : Copyright 1999, The R Development Core Team Version 0.64.0 Patched (May 3, 1999) R is free software and comes with

Dear friends. When using in windows version the menu function to change directory, what actions are actually performed, i.e. what commands could I write to a file to make a CD ? I couldn't find it mentioned anywhere - undoubtedly an oversigth on my part. Troels Ring, M.D Department of Nephrology Aalborg, Denmark tring at mail1.stofanet.dk

Dear friends. I've made a very simple procedure to make Bland-Altman plots and it works OK except for the inability to take a variate as argument in annotation. I'm sure it is a very simple error, so if you have 5 seconds please tell me. Best wishes Troels bland <- function (x) #accepts two columns { mn <- 0.5*(x[,1]+x[,2]) diff <- x[,1]-x[,2]

move or copy the directories mass, nnet and class to the library directory - then execute link.html.help() Now execute library(MASS) and data(petrol) Should work. "Troels Ring"

Dear friends. I wanted to reproduce a sheet of paper useful for plotting the reciprocal of some biochemical variable (creatinine) over time (Bleyer Am J Kidney Dis 34:576-578,1999). Time on X axis with units years and months. On Y axis, 1/creatinine, but with the original creatinine value indicated - to make it easier for people to use it. When using the graph, the linear decrease in 1/crea over

Dear friends. This is again a very simple question, and I hope you will accept it. I have a dataset with variables in columns all with their names. Now I want to add a column indicating simply identity of patients, numbered according to row. There are 34 rows. So I did: >ID <- c(1:34) but was not allowed to >Ptt <- cbind(patients,ID) and was informed that >Error: names attribute

Dear friends. I have trobpe with path. Below I have tried to paste a path and file name but it will not work > path <- "D:/rw0641/own/own/procedures/GFR_TIME" > path [1] "D:/rw0641/own/own/procedures/GFR_TIME" > path1 <- paste(path,"gfr.txt",sep="/") > path1 [1] "D:/rw0641/own/own/procedures/GFR_TIME/gfr.txt" -- this is good

Dear Peter. I guess you know that Jim Lindseys code include nordr and ordglm in library gnlm - I attach the htmls which do various linear and nonlinear ordinal regressions - exemplified with just the data mentioned, McCullagh (1980) JRSS B42, 109-142. I had it work very fine. -----Oprindelig meddelelse----- Fra: Peter Malewski <p.malewski at tu-bs.de> Til: Troels Ring <tring at

Dear friends. I appreciate the help that plot.date is not yet fully developed but beautiful and efficient as is anyway. I attach again an example from a real patient and now the problem is how to represent a chemical C-reactive protein (CRP) on the plot together with concentration of vancomycin. CRP was 25-200 and vancomycin below about 20. I only could get it about right by dividing CRP by

Dear friends. In the documentation for PAR it is specified that "las" allowing turning axis labels away from default=parallel to axis is not implemented. I need to make graphs with horizontal Y axis labels - how to ? I'm using win98. Thank you in advance. Sincerely yours Troels Ring, MD Department of Nephrology Aalborg Hospital Aalborg Denmark

Dear friends. I would like to link an object file made in C. Dyn.load rather wants to link a whole library and I only wanted to link this small object file. I followed instructions given for S - by Clive Loader in his software for changepoint handling and it is his object file I made from http://cm.bell-labs.com/cm/ms/departments/sia/clive/smcp.shar Any dirty tricks or easy oversights ? Yours

Dear friends. I just downloaded tseries dated 27.10.99 (from Aalborg, Denmark) and it seems it lacked something, library ts was not available ? What to do ? > library(tseries) Loading required package: ts Warning: There is no package called `ts' > > Best wishes Troels Ring, M.D Department of Nephrology Aalborg, Denmark tring at mail1.stofanet.dk

Dear friends. I'm in windows version R1.0.0. The code below produced a very beautiful and efficient graph in R Version 0.99.0 Patched (February 9, 2000) but in the newest, the x-axis is unreadable with dates written unformatted. I tried to set the numbers using par(xaxp) or xaxs but to no avail. Since it was OK just in the last version I guess there might perhaps be something wrong

Dear friends. Can anyone explain what is happening here ? I am told that no data entry editor is available in this version of R ? The example is not far from that shown in help(data.entry) ? Besides, could an empty spreadsheet entry be an option to start inputting data ? R : Copyright 2000, The R Development Core Team Version 1.0.0 (February 29, 2000) R is free software and comes with

Dear friends. I tried the session below with 10 MB in both vsize and nsize but didn't get the example work. Is this a problem in LME or in me or both or somewhere else or undefined ? R : Copyright 1999, The R Development Core Team Version 0.64.0 Patched (May 3, 1999) R is free software and comes with ABSOLUTELY NO WARRANTY. You are welcome to redistribute it under certain conditions. Type

Dear friends. On the Bayesian averaging homepage http://www.research.att.com/~volinsky/bma.html I found some S code some of which perhaps may run in R. There was a call to an algorithm possibly within S but not supported by R 64.1: "leaps and bounds". I guess it is a minimization step. Can anyone clarify the algorithm and perhaps even give a pointer to some code ? I guess this may be

Dear friends. Do Polr in Mass change the sign of the coefficients ? Example (McCullagh 1980) options(contrasts=c("contr.treatment","contr.poly")) library(Mass) freq <- c(19,29,24,497,560,269) yy <- ordered(gl(3,1,6)) z4 <- polr(yy~x,weights=freq) > z4 Call: polr(formula = yy ~ x, weights = freq) Coefficients: x2 -0.6026492 Intercepts: 1|2

Dear friends. As promised ver 1.01 solved the former problem and now dates plot OK in windows. But. In the following code, "cex.axis"=0.5 works OK but "xaxp"=c(14620,14680,4) has seemingly no effect. The limits for xaxp was found from the default plot and inspecting par(). I need to put the date marks at specific points. Could I suppress them altogether and make a new xaxis