Use o recurso de procurar palavras do seu navegador e procure por "se.fit" dentro da página http://www.leg.ufpr.br/~walmes/cursoR/fiocruz/fiocruz05hipot.html. Basicamente, para obter o erro padrão de um valor predito, você faz
sqrt(diag(L%*%V%*%t(L))) em que L é a matriz de coeficientes da função linear estimada, no caso X%*%beta, e V a matriz de covariância das estimativas dos parâmetros. No código eu uso decomposição de Cholesky porque é computacionalmente mais interessante. À disposição. Walmes.
_______________________________________________ R-br mailing list [email protected] https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forne�a c�digo m�nimo reproduz�vel.
