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

[Fernando Antonio de souza]

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