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
_______________________________________________ R-br mailing list [email protected] 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.
