confidence intervals or SE of gls {nlme} predicted values
-
02-10-2022 - |
質問
I'm running a multivariate gls model:
m <- gls(y ~ x + factor1 + factor2, cor = corPagel(1,phylogeny), weight= ~1/log(n))
I want to plot the results and I could get predicted values like this:
newdata <- data.frame(expand.grid( x = mean(x),
factor1= unique(factor1),
factor2 = unique(factor2)))
predvals <- predict(m,newdata)
Although two of the variables I want to plot are factors, and I want to get confidence intervals or standard errors for these predicted values.
Does anyone know how to do it?
Thanks in advance for your help! Cheers, Jasmine
正しい解決策はありません
所属していません StackOverflow