Obrigado Walmes, Era exatamente isso que queria. abr cleber Em 09/11/2015 14:40, Walmes Zeviani escreveu:
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 <http://www.leg.ufpr.br/%7Ewalmes/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.
--- Este email foi escaneado pelo Avast antivírus. https://www.avast.com/antivirus