Fwd: Medidas repetidas com estruturas de co-variâncias
Olá. Estou com dificuldade de fazer atribuir corretamente os argumentos no meu modelo de medidas repetidas com estruturas de co-variâncias. . Estou tentando fazer no nlmer, talvez possam me ajudar. O código no sas seria: proc mixed data=cruza order=data; class raca data animal; model peso = raca data raca*data/outp=p ddfm=satterth; repeated data / subject = animal(raca) type = sp(pow)(data) r rcorr; run; no *R*: dados <- data.frame(read.csv("dados.txt", sep=" ", header=TRUE)) dados$Peso <- as.numeric(cruza$Peso) library(nlme) modelo <- lme(Peso~Data*Raca, .............., data=dados) Os dados (parcial) estão em: https://www.dropbox.com/s/eajz72kmgk1f897/dados.txt Att, Mauricio
Seu exemplo não é reproduzível. Facilite a vida de quem quer te ajudar fornecendo código para leitura e ajuste dos teus dados com um read.table() a partir do link do teu arquivo de texto, por exemplo. Para solucionar sua dúvida sugiro que leia http://plantecology.syr.edu/fridley/bio793/mixed2.html http://plantecology.syr.edu/fridley/bio793/mixed1.html À 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 skype: walmeszeviani twitter: @walmeszeviani homepage: http://www.leg.ufpr.br/~walmes linux user number: 531218 ==========================================================================
participantes (2)
-
Mauricio Morgado Oliveira -
Walmes Zeviani