modelo linear com `lm`: como obter variação de previsão da soma dos valores previstos
Estou somando os valores previstos de um modelo linear com vários preditores, como no exemplo abaixo, e quero calcular a variação combinada, o erro padrão e possivelmente os intervalos de confiança para essa soma.
lm.tree <- lm(Volume ~ poly(Girth,2), data = trees)
Suponha que eu tenha um conjunto deGirths
:
newdat <- list(Girth = c(10,12,14,16)
para o qual eu quero prever o totalVolume
:
pr <- predict(lm.tree, newdat, se.fit = TRUE)
total <- sum(pr$fit)
# [1] 111.512
Como posso obter a variação paratotal
?
Questões semelhantes sãoaqui (para GAMs), mas não sei como proceder com ovcov(lm.trees)
. Ficaria muito grato por uma referência para o método.