Caro walmes, Mais uma vez meu muito obrigado realmente era o que eu desejava fazer. até a próxima Em 23 de setembro de 2011 11:03, Walmes Zeviani <walmeszeviani@gmail.com>escreveu:
Fernando,
Para isso você precisa saber qual a definição dele de intercepto-inclinação geral. Penso que seja aquela obtida usando a restrição soma zero nas estimativas, mas pode ser outra. Para conseguir essas estimativas você nem precisa de pacote extra, é sou definir corretamente o tipo de contraste (ou restrição) que você quer. Evite enviar dados em anexo ou colados na mensagem. Isso não deixa o seu código reproduzível com um copia e cola. Eu teria que fazer download, blá blá, isso enrola e eu prefiro gerar dados do que ter que todo esse trabalhão. Para facilitar a vida de quem tá disposto a te ajudar, passe o dput() do seus dados, que assim com um copia e cola reproduzímos o seu código. Mostro isso no final.
da <- expand.grid(estudo=gl(5,1), x=1:10) da$y <- round(rnorm(nrow(da),0,1),2)
m0 <- lm(y~estudo*x, data=da) coef(m0) m1 <- lm(y~-1+estudo/x, data=da) # mesmo espaço coluna do anterior coef(m1) # são os interceptos e inclinações para cada estudo
m2 <- lm(y~estudo*x, data=da, contrast=list(estudo=contr.sum)) coef(m2) coef(m2)[c(1,6)] # intercepto e inclinação geral coef(m2)[2:5] # interceptos desvios do geral para estudos de 1 a 4 -sum(coef(m2)[2:5]) # estudo 5 coef(m2)[7:10] # inclinações desvios da geral para estudos de 1 a 4 -sum(coef(m2)[7:10]) # estudo 5
summary(m2)
# o resultado dessa função desse ver enviado para um # código reproduzível com copia e cola dput(da)
À disposição. Walmes.
========================================================================== Walmes Marques Zeviani LEG (Laboratório de Estatística e Geoinformação, 25.450418 S, 49.231759 W) Departamento de Estatística - Universidade Federal do Paraná fone: (+55) 41 3361 3573 VoIP: (3361 3600) 1053 1173 e-mail: walmes@ufpr.br twitter: @walmeszeviani homepage: http://www.leg.ufpr.br/~walmes linux user number: 531218 ==========================================================================
_______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br 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.