Displaying 1 result from an estimated 1 matches for "fit_vol_model1".
2012 Apr 01
1
Degrees of Freedom for lme.
...e ID (time they were seen).
I believe my data is set up the right way, but when I run it, I get DF for
Intercept is 49, and DF for slope is 13? Why?
lme.gort=lme(Volume ~ GORT_FLUENCY+AGE, random = ~ 1 | TIMEID, data =
subset(vol_data, INDEX1=='LH_FUSIFORM'), na.action=na.omit)
fit_vol_model1 <- function(df){
tryCatch(lme.gort <- lme(Volume ~ GORT_FLUENCY+AGE, random = ~ 1 | UID,
data=df, na.action=na.omit), error=function(err) tag <<-1)
data.frame(Term = rownames(anova(lme.gort)), anova(lme.gort))
+ }
models = list()
models$anova = ddply(vol_data, c("IND...