search for: alldata3

...to plot my actual averaged data points within 95 % confidence intervals from the model. Plotting the model is easy, but determining confidence bands is not. Here is my model: ratiomodel<-glmmPQL(ratio~as.factor(joint)*time, random = ~ 1 | subject, family = Gamma(link = "identity"),alldata3) I am used to seeing confidence intervals from models that increase, “flair out” in the y direction, at the beginning and ending time points (x values) of the simulated data. If I use 'lm' and pass the command 'int = "c" ' 'to create this model I can easily find and...