Displaying 1 result from an estimated 1 matches for "3bclp".
Did you mean:
hbclp
2004 May 27
1
Getting the same values of adjusted mean and standard errors as SAS
...e adjusted means for the males and females at meanAge using
predict
# set up data frame to get predicted values at meanAge
evalDF <- data.frame(Age = meanAge, Gender = c("F","M"),
Group = c("1NC","1NC", "2UCLP", "2UCLP", "3BCLP", "3BCLP",
"4ICP", "4ICP", "5CLPP", "5CLPP"))
mod <-lm(Measurement ~ Age + Gender + Group, data=dd)
pred <- predict(mod, evalDF, se.fit = TRUE)
adjDF <- data.frame(evalDF,fit = pred$fit, se = pred$se.fit)
adjMeanMale <- mean(adjDF$...