[R] adding statistical output to a plot

I have written a function to emulate minitab's QQ plotting output (with SW test and AD test results on the graph): mtab.norm<-function(x) { library(nortest) library(lattice) x<-as.numeric(x) x<-as.vector(x) plot.ht<-4.6 plot.wd<-4.6 pt.ht=plot.ht/5 txt.sz<-(plot.ht/7.5) X11(width=plot.wd, height=

Re: [R] adding statistical output to a plot

Thanks to both of you, I used the below code instead, which has been cleaned and pruned a little. It produces a trellis of qq plots with the proper annotations as per Minitab from a dataframe passed to mtab.matrix(). However, I can't get individual variable names (these are column names in the pa

[R] Indexing from two variables

Hello, thank you all for your patience and time I am essentially trying to get disorganised data into long form for linear modelling. I have 2 dataframes "rec" and "book" Each row in "book" needs to be pasted onto the end of several of the rows of "rec" according to two variables in the row:" M

[R] Importing dates from SPSS file

Hello all, kind regards, I have imported a data.frame from SPSS using "foreign":read.spss but unfortunately it is reading dates in a way neither R nor myself can understand. > book$DATE [1] 13502246400 13443321600 13477795200 13472956800 13501728000 13445395200 13501382400 13502851200 134441856

Re: [R] Importing dates from SPSS file

as.POSIXlt(book$DATE, origin="1582/10/14") works. The Gregorian calendar was the kicker, thanks to ggrothendieck at gmail.com Ross -- View this message in context: http://r.789695.n4.nabble.com/Importing-dates-from-SPSS-file-tp3260293p3260329.html Sent from the R help mailing list archive at

[R] Just another pattern matching / indexing question

Hi all, thank you for your patience. I am dealing with a large dataset detailing patients and medications Medications are hard to code, as they are (usually) meaningless unless matched with doses. I have a dataframe with vectors (Drug1, Drug2. Drug 16) and individual patients are represente

Re: [R] Unbalanced Mixed Linear Models With Nested Stratum

Hi, I have found Pinheiro and Bates an absolute godsend for mixed modelling, http://www.amazon.com/Mixed-Effects-Models-S-S-Plus/dp/0387989579 alternatively, there is a REALLY good chapter in the 3rd edition of DAAG: http://www.amazon.com/Data-Analysis-Graphics-Using-Example-Based/dp/05217629

Re: [R] Iterate over a list of input files?

Quick and dirty answer: Try ?paste in conjunction with a "for" loop which iterates over i:length(Files), and you can paste "i" onto the end of your filename, then add the data to cd=list(), cd[i], then "unlist" the data or as.data.frame(cd). I'm sure there are more elegant alternatives. Ross --

[R] Help using "tm" text mining package - preprocessing

Thanks all for your help. I fear text mining is an abstract little corner of "R". I have imported 3228 text (.txt) files, each a news story, into R using [tm]: textd <- Corpus(DirSource("other/docs"), readerControl = list(reader =readPlain)) I can pre-process each individual document using tolo

[R] Reshape - novel error message and fail

Hello all, thanks in advance for the trouble. I can't reshape a particular large data.frame. I have used melt and reshape successfully before, alot. I am trying to go from wide to long format. I wish to separate The named variables below into 3 groups : SYST, DIAST and PULSE, and have them vary

[R] argument "x" is missing, with no default - Please help find argument x

R programming question, not machine learning, although that's the content. Apologies to all for whom the following code is eye-burning. I am using foreach() to run a simulation on a randomForest model (actually conditional randomForest ... "party" package). The simulation is in two dimensions. exa

[R] How do I extract coefficient standard errors /CI for a "coxme" model

Hello, and thanks for your time I'm trying to extract standard errors to produce confidence intervals from a multivariable coxme model object so I can write a function that will print a summary for some reproducible research. As far as I can glean, the SE is produced on-the-fly by the print meth

[R] Multilevel Survival Analysis - Cox PH Model

Hello all, thanks for your time and patience. I'm looking for a method in R to analyse the following data: Time to waking after anaesthetic for medical procedures repeated on the same individual. > str(mysurv) labelled [1:740, 1:2] 20 20 15 20 30+ 40+ 50 30 15 10 ... - attr(*, "dimname

Re: [R] Multilevel Survival Analysis - Cox PH Model

There is indeed right censoring, but I obviously didn't explain it very well. Patients are either fully oriented or not (1 or 2) after an hour. If they're not, then the data is right censored. However, I don't feel that "coxme" is overkill at all, as I may also have to account for repeated COURS

[R] number of items to replace is not a multiple of replacement length

Please help with this error message drugbook is an 885 x 32 dataframe >names(drugbook) [1] "DRUG1" "DRUG2" "DRUG3" "DRUG4" "DRUG5" [6] "DRUG6" "DRUG7" "DRUG8" "DRUG9" "DRUG10" [11] "DRUG11" "DRUG12" "DRUG13" "DRUG14" "DRUG15"

[R] glmmBUGS fails to accept higher level covariates

Hello I am using glmmBUGS to fit a multilevel model. Treatments are nested in Course are nested in Patients. The predicted variable in total EEG duration. The predictors are: at the observation level : Medication dose at the Course level: Weight in KG and Age at the Patient level: Weight in KG a