[Rd] lme: how to extract the variance components?

Jean Wu darmy16@hotmail.com
Wed, 14 Nov 2001 12:43:53 -0500


This is a multi-part message in MIME format.

------=_NextPart_000_0005_01C16D0A.0459F770
Content-Type: text/plain;
	charset="iso-8859-1"
Content-Transfer-Encoding: quoted-printable

      Dear all,=20
      Here is the question:
      For example, using the "petrol" data offered with R.
      pet3.lme<-lme(Y~SG+VP+V10+EP,random=3D~1|No,data=3Dpetrol)
      pet3.lme$sigma gives the residual StdDev.
      But I can't figure out how to extract the "(intercept) StdDev",=20
      although it is in the print out if I do "summary(pet3.lme)".

      In S-plus3.4 , $var.ran is the estimate of the variance of the =
random=20
      effects between clusters. what is the equivalent thing to do in R?

      Thank you very much!

    =20

Jean

------=_NextPart_000_0005_01C16D0A.0459F770
Content-Type: text/html;
	charset="iso-8859-1"
Content-Transfer-Encoding: quoted-printable

<!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 4.0 Transitional//EN">
<HTML><HEAD>
<META content=3D"text/html; charset=3Diso-8859-1" =
http-equiv=3DContent-Type>
<META content=3D"MSHTML 5.00.3018.900" name=3DGENERATOR>
<STYLE></STYLE>
</HEAD>
<BODY bgColor=3D#ffffff>
<DIV><FONT size=3D2>
<TABLE width=3D"100%">
  <TBODY>
  <TR>
    <TD>
      <P><FONT face=3DArial>Dear all, <BR>Here is the question:<BR>For =
example,=20
      using the "petrol" data offered with=20
      =
R.<BR>pet3.lme&lt;-lme(Y~SG+VP+V10+EP,random=3D~1|No,data=3Dpetrol)<BR>pe=
t3.lme$sigma=20
      gives the residual StdDev.<BR>But I can't figure out how to =
extract the=20
      "(intercept) StdDev", <BR>although it is in the print out if I do=20
      "summary(pet3.lme)".<BR><BR>In S-plus3.4 , $var.ran is the =
estimate of the=20
      variance of the random <BR>effects between clusters. what is the=20
      equivalent thing to do in R?<BR><BR>Thank you very=20
      =
much!<BR></FONT></FONT></P></TD></TR></TBODY></TABLE></FONT></DIV><FONT=20
face=3DArial size=3D2>Jean</FONT></BODY></HTML>

------=_NextPart_000_0005_01C16D0A.0459F770--
-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-
r-devel mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html
Send "info", "help", or "[un]subscribe"
(in the "body", not the subject !)  To: r-devel-request@stat.math.ethz.ch
_._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._