Walmes,
Problema resolvido, valeu..
Ate mais
Em 4 de junho de 2011 19:46, Walmes Zeviani <[email protected]>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@Zt)
> t(as.matrix(fm1@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: [email protected]
> twitter: @walmeszeviani
> homepage: http://www.leg.ufpr.br/~walmes
> linux user number: 531218
> ==========================================================================
>
> _______________________________________________
> R-br mailing list
> [email protected]
> https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
>
>
_______________________________________________
R-br mailing list
[email protected]
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
