Hi guys, how can i get the equivalent of the estimate code used in proc mixed in SAS linear mixed-method models in R.
The estimate statement is used to show the expected change in the dependent variable given a value of the independent variable. This is what it looks like in SAS.
PROC MIXED DATA=exerclong;
CLASS id group(ref='2')/ref=first;
MODEL y = group day group*day/solution;
random intercept/subject=ID;
- We use the estimate statement to get predictions based on fixed covariate values;
estimate 'Weight group: day 6 change in strength intercept 0 group 0 0 day 6 group*day 0 6;
My R code is something like this.
model1= lmer(y ~ factor(group) + Day + factor(group):Day + (1 |id), data = exerclong)
I have tried the predict and fit contrast function but looks like they are not appropriate to predict the value of y at a particular day, say day 6.
Can someone help with this?