Displaying 2 results from an estimated 2 matches for "obsno".
2012 Jan 19
2
Reading in tab (and space) delimited data within a script XXXX
...data <- read.table( textConnection(mystring),
header=TRUE, sep=",",
row.names="id", na.strings=" ")
closeAllConnections()
mydata
Can anyone suggest a similar approach for reading in tab-delimited or
single space delimited data? Example data:
data3<-
"OBSNO AGE SEX ALKPHOS LAB CAMMOL PHOSMMOL AGEGROUP
21 76 M 84 5 3.2 0.9 3
22 76 M 5 2.18 0.84 3
23 68 M 82 5 2.15 0.52 1
24 69 M 84 5 2.3 1.36 1
25 76 F 100 3 25.3 1.07 3
26 70 F 90 3 20 0.97 2
27 71 F 109 3 22.3 0.94 2
28 70 -99 65 3 24.3 1.42 2
29 74 F 61 3 25 0.87 2
30 74 F 62 3 23.3 0.94 2"
Th...
2012 Aug 07
1
Which R function for GLMM with binary response, nested random factors with temporal correlation?
...sponse is binary.
So I need a GLMM with a correlation structure. I've tried using GEE, but the R packages can't cope with multilevel nested data. The only R function I've found that can do this is glmmPQL.
m <- glmmPQL(y ~ f1 * f2 * f3 + (1|subj/day), correlation=corAR1(form =~obsno|subj/day))
f1 - f3 are fixed factors
However, PQL estimation is not recommended for binary response data. With no AIC and unreliable p values, model selection seems impossible! So my question is:
1) are there any other functions which are suitable for a GLMM with multilevel nested random effects...