[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

[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.

[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

[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 [16] DRUG16 DrugDose1

[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

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

[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(*,

[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

[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

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:

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] 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 13444185600

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] 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: MRN and

[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,

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