[R] Rare behaviour for nlme::reStruct example and question about ?lmeObject

Martin Maechler m@ech|er @end|ng |rom @t@t@m@th@ethz@ch
Fri Nov 18 10:45:50 CET 2022


>>>>> Iago  
>>>>>     on Thu, 17 Nov 2022 11:53:31 +0100 writes:

    > Thank you Martin,
    > Regarding my question about `terms`, I meant the `terms` component of 
    > the `lme` output. For example, for

    > fm1 <- lme(distance ~ age, data = Orthodont)

    > I can see through str(fm1), the next:

    > $ terms       :Classes 'terms', 'formula'  language distance ~ age
    >   .. ..- attr(*, "variables")= language list(distance, age)
    >   .. ..- attr(*, "factors")= int [1:2, 1] 0 1
    >   .. .. ..- attr(*, "dimnames")=List of 2
    >   .. .. .. ..$ : chr [1:2] "distance" "age"
    >   .. .. .. ..$ : chr "age"
    >   .. ..- attr(*, "term.labels")= chr "age"
    >   .. ..- attr(*, "order")= int 1
    >   .. ..- attr(*, "intercept")= int 1
    >   .. ..- attr(*, "response")= int 1
    >   .. ..- attr(*, ".Environment")=<environment: R_GlobalEnv>
    >   .. ..- attr(*, "predvars")= language list(distance, age)
    >   .. ..- attr(*, "dataClasses")= Named chr [1:2] "numeric" "numeric"
    >   .. .. ..- attr(*, "names")= chr [1:2] "distance" "age"

    > However, there is no paragraph for this component among those in 
    > ?lmeObject (in the Value section) as there is for fixDF, apVar, call, 
    > etc.. If not needed, I assume the explanation would be that `terms` does 
    > not need to be included in a legitimate ´"lme"´ object'.

    > Kind regards,
    > Iago

I see; now I understand what you meant. I've taken time to check and
found that the 'terms' part *is* needed, even though not for the
simplest methods such as summary(), residuals(), ..  and so have
added it to the  ?lmeObject  help page for the next version of
nlme.

Thank you, Iago!

Best,
Martin



More information about the R-help mailing list