[R-br] Matrix Modelo

[Eder David Borges da Silva]

Walmes,
Problema resolvido, valeu..
Ate mais

Em 4 de junho de 2011 19:46, Walmes Zeviani <walmeszeviani em gmail.com>escreveu:

> Éder,
>
> Se você quer obter essa matriz saturada, então é só remover os interceptos
> de cada porção e parear as matrizes (como o Benilton sugeriu), ó
>
> terms <- c(~-1+BLOCO, ~-1+CLONE, ~-1+LOCAL:CLONE)
> str(terms)
> X <- cbind("(Intercept)"=1, do.call(cbind, lapply(terms, model.matrix,
> data=dados)))
>
> Para extrair a matriz de ajustes da lmer tem use abaixo
>
> fm1 <- lmer(Reaction ~ Days + (Days|Subject), sleepstudy)
> summary(fm1)
> model.matrix(fm1) # efeitos fixos, X
> str(fm1)
> str(fm1 em Zt)
> t(as.matrix(fm1 em Zt)) # efeitos aleatórios, Z
>
> À 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
>
>
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20110604/38d29f11/attachment.html>