Displaying 1 result from an estimated 1 matches for "lmmean".
Did you mean:
mmean
2018 Apr 09
1
linear regression with repeated measures
...3*dat[,4:7]+runif(48)
# some individual measures may be missing :
dat[2:3,4] <- NA
colnames(dat) <- paste(rep(c("a","b"),3:4),c(1:3,1:4),sep="")
# In analogy to the ample documentation of lm() :
datMean <- cbind(aM=rowMeans(dat[,1:3]),bM=rowMeans(dat[,4:7]))
(lmMean <- lm(bM ~ aM,data=as.data.frame(datMean)))
# I suppose the estimated parameters (intercet & slope) may be correct but
sice the degrees of freedom are not made of means I am convinced they are
incorrect and thus any statistics using them will be so, too ...
df.residual(lmMean)
summary(lmMean...