Re: [R] Extract Variance Components

Well, I'm going to reply to my own thread with a solution here, turns out one attempt we made last week nearly had it, slight adjustment made it work. For anyone that is interested / in the future wants to achieve the same thing > *varcomp <- matrix(nrow=0, ncol=3) for (i in 1:nlevels(na

[R] Extract Variance Components

Hi, I'm still having problems putting the variance components of my model in to a data frame, it is a continuation of this discussion, http://r.789695.n4.nabble.com/ANOVA-problem-td4609062.html, but now focussed on the problem of extracting variance components. I have got my mixed effects model now

Re: [R] extract variance components

See the error message. Thanks. a.fit<-lme(distance~age, data=aaa, random=~day/subject) getVarCov(a.fit) Error in getVarCov.lme(a.fit) : Not implemented for multiple levels of nesting On Wed, Sep 10, 2008 at 9:53 PM, Lyman, Mark <[EMAIL PROTECTED]> wrote: > Take a look again at help(getVarCov

Re: [R] extract variance components

Take a look again at help(getVarCov). Mark Lyman, Statistician ATK Launch Systems [EMAIL PROTECTED] From: huang min [mailto:[EMAIL PROTECTED] Sent: Tuesday, September 09, 2008 10:02 PM To: r-help@r-project.org Cc: [EMAIL PROTECTED] Subject: Re: extract variance components Hi, Thanks

Re: [R] extract variance components

Hi, Thanks for all. I now know how to extract the \sigma's. For the unbalanced model y_{ijk}=x\beta+\alpha_i+\beta_{ij}+e_{ijk} i=1,2,\dots,a, j=1,2,\dots,b_i, k=1,2,\dots,n_{ij} How can I extract the variance matrix $V$? The variance for the ith group is also of help. Suppose the ith group has t

Re: [R] extract variance components

huang min gmail.com> writes: > > HI, > > I would like to extract the variance components estimation in lme function > like > > a.fit<-lme(distance~age, data=aaa, random=~day/subject) > > There should be three variances \sigma_day, \sigma_{day %in% subject } and > \sigma_e. > > I can extract

Re: [R] extract variance components

huang min gmail.com> writes: > I would like to extract the variance components estimation in lme function > like > > a.fit<-lme(distance~age, data=aaa, random=~day/subject) > Try VarCor(a.fit) Using lme for the type of problems you have is just fine; in many respects, lme4 currently is a wel

[R] extract variance components

HI, I would like to extract the variance components estimation in lme function like a.fit<-lme(distance~age, data=aaa, random=~day/subject) There should be three variances \sigma_day, \sigma_{day %in% subject } and \sigma_e. I can extract the \sigma_e using something like a.fit$var. However, I