[Rd] formula(CO2)
Gavin Simpson
gavin.simpson at ucl.ac.uk
Mon Jul 16 16:21:44 CEST 2007
On Mon, 2007-07-16 at 14:57 +0100, ted.harding at nessie.mcc.ac.uk wrote:
> On 16-Jul-07 13:28:50, Gabor Grothendieck wrote:
> > The formula attribute of the builtin CO2 dataset seems a bit strange:
> >
> >> formula(CO2)
> > Plant ~ Type + Treatment + conc + uptake
> >
> > What is one supposed to do with that? Certainly its not suitable
> > for input to lm and none of the examples in ?CO2 use the above.
>
> I think one is supposed to ignore it! (Or maybe be inspired to
> write a mail to the list ... ).
>
> I couldn't find anything that looked like the above formula from
> str(CO2). But I did spot that the order of terms in the formula:
> Plant, Type, treatment, conc, uptake, is the same as the order
> of the "columns" in the dataframe.
CO2 is a groupedData object not a data.frame per se.
> class(CO2)
[1] "nfnGroupedData" "nfGroupedData" "groupedData" "data.frame"
What Gabor saw was the result of formula.data.frame being applied to
CO2, which, as you surmise Ted, is being produced from the columns of
CO2.
> formula(CO2)
Plant ~ Type + Treatment + conc + uptake
> stats:::formula.data.frame(CO2)
Plant ~ Type + Treatment + conc + uptake
But if we load nlme, we see that now (via a groupedData method for
formula) a more useful formula is displayed.
> require(nlme)
Loading required package: nlme
[1] TRUE
> formula(CO2)
uptake ~ conc | Plant
as this can then be directly used in lme
> lme(CO2)
Linear mixed-effects model fit by REML
Data: CO2
Log-restricted-likelihood: -283.1447
Fixed: uptake ~ conc
(Intercept) conc
19.50028981 0.01773059
<snip />
But like Gabor, I'm struggling to see where this [formula(data.frame)]
might be used (useful)?
G
--
%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%
Gavin Simpson [t] +44 (0)20 7679 0522
ECRC, UCL Geography, [f] +44 (0)20 7679 0565
Pearson Building, [e] gavin.simpsonATNOSPAMucl.ac.uk
Gower Street, London [w] http://www.ucl.ac.uk/~ucfagls/
UK. WC1E 6BT. [w] http://www.freshwaters.org.uk
%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%
More information about the R-devel
mailing list