[R-es] extraer un valor de un objeto lmer

gmonteru gmonteru en balcarce.inta.gov.ar
Mie Nov 25 00:23:28 CET 2009


Hola Teresa,

Voy a llamar "modelo" a tu modelo. 
Explorá en modelo en deviance, en particular con
modelo en deviance["sigmaREML"] obtenés el desvío
estándar (estimación REML) de la componente que
buscás.
Saludos

Gloria

----- Original Message ----- From: "Teresa Boca"
<tboca en cnia.inta.gov.ar> To: <R-help-es en r-project.org>
Sent: Tuesday, November 24, 2009 3:25 PM
Subject: [R-es] extraer un valor de un objeto lmer



> Hola, a todos!
>


Estoy ajustando un modelo componentes de varianza utilizando
el paquete lme4, quiero extraer luego los valores estimados
de varianza de cada factor aleatorio de este objeto lmer y
guardarlo en un vector. Existe la función VarCorr,  Esto
es
lo que hago:

Mi ejemplo: tengo 3 factores NDVIaño;  Máx.IH; TV

# Ajusto el modelos
lmer(peligro100 ~ 1 + (1|NDVIaño) + (1|Máx.IH) + (1|TV),
data=Verano) #Extraigo las componentes de varianza
Com_Var<-VarCorr(var_verano_NDVIaño )

Com_Var$NDVIaño[1] # Este es el valor de varianza
deNDVIaño
Com_Var$Máx.IH[1] # Este es el valor de varianza
deMáx.IH
Com_Var$TV[1] # Este es el valor de varianza de TV


*Pero no se como extraer la componente de varianza
residual....* Desde ya agradezco alguna sugerencia

saludos
teresa









-- 
Ing. Agr. M.Sc.Teresa Boca
Instituto de Clima y Agua INTA Castelar.
Las Cabañas y Los Reseros s/n°
1712 Castelar Pcia de Buenos Aires, Argentina
Cel: 01115-3572-9584
Tel: (54-11)-4621-0125/1684
Fax:(54-11)-4621-5663
tboca en cnia.inta.gov.ar

[[alternative HTML version deleted]]




------------------------------------------------------------
--------------------


_______________________________________________
R-help-es mailing list
R-help-es en r-project.org
https://stat.ethz.ch/mailman/listinfo/r-help-es

____________________________________________________________
_______________

Aviso:
====
El contenido del presente e-mail y sus posibles adjuntos
pertenecen al INTA y pueden contener información
confidencial. Si usted no es el destinatario original de
este mensaje y por este medio pudo acceder a dicha
información, por favor solicitamos contactar al remitente
y
eliminar el mensaje de inmediato. Se encuentra prohibida la
divulgación, copia, distribución o cualquier otro uso de
la
información contenida en el presente e-mail por parte de
personas distintas al destinatario.

This e-mail contents and its possible attachments belong to
INTA and may contain confidential information. If this
message was not originally addressed to you, but you have
accessed to such information by this means, please contact
the sender and eliminate this message immediately.
Circulation, copy, distribution, or any other use of the
information contained in this e-mail is not allowed on part
of those different from the addressee.


Antes de imprimir este mensaje, asegúrese de que sea
necesario. Proteger el medio ambiente está también en su
mano.


___________________________________________________________________________

Aviso:
==
El contenido del presente e-mail y sus posibles adjuntos pertenecen al INTA y pueden contener información confidencial. Si usted no es el destinatario original de este mensaje y por este medio pudo acceder a dicha información, por favor solicitamos contactar al remitente y eliminar el mensaje de inmediato. Se encuentra prohibida la divulgación, copia, distribución o cualquier otro uso de la información contenida en el presente e-mail por parte de personas distintas al destinatario. 


This e-mail contents and its possible attachments belong to INTA and may contain confidential information. If this message was not originally addressed to you, but you have accessed to such information by this means, please contact the sender and eliminate this message immediately. Circulation, copy, distribution, or any other use of the information contained in this e-mail is not allowed on part of those different from the addressee.


Antes de imprimir este mensaje, asegúrese de que sea necesario. Proteger el medio ambiente está también en su mano.



Más información sobre la lista de distribución R-help-es