[R-br] Função estimate para o R

[Walmes Zeviani]

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 em ufpr.br
twitter: @walmeszeviani
homepage: http://www.leg.ufpr.br/~walmes
linux user number: 531218
==========================================================================
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20110923/698f8839/attachment.html>