Caros amigos Estou analisando dados em um modelo misto linear utilizando o pacote nlme contendo uma matrix de correlação do tipo compSymm e necessito utilizar como argumento "weights" o* inverso do quadrado do erro padrão* de cada média. Alguem sabe como faço isto. Sei que o argumento weigue utiliza algus objetos varfunction mas não compreendo muito bem seus objetos e não sei se algum deles fazem o que desejo. Alguém pode me ajudar? consumo<-lme(Fixed=Y~X, data=stpierre,random=X|Estudo) consumo<-update(consumo, random=compSymm(~X), Weight=?) Preciso usar o W2 como peso sendo: w2 =w1/ w sendo w1= o inverso do quadrado médio do erro w= média dos quadrados médios
Nas funções da nlme o argumento weights é usado para modelar a heterocedasticidade. É diferente do argumento weights da lm(), onde se passa um vetor de pesos que são fixos e de mesma dimensão das observações. Nas da nlme é passada uma função que corresponde ao modelo que descreve a variância, onde uma estimação de parâmetros é feita. Consulte a documentação para detalhes. À 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 twitter: @walmeszeviani homepage: http://www.leg.ufpr.br/~walmes linux user number: 531218 ==========================================================================
participantes (2)
-
Fernando Antonio de souza -
Walmes Zeviani