Re: [R] Problems with lme

From: ONKELINX, Thierry <Thierry.ONKELINX_at_inbo.be>
Date: Sat, 24 May 2008 11:31:06 +0200

That's because you have not specified the random effect. I think you want

lme(Error ~ Sequence + Temperature + Tumour, random = ~1|ID, data = YourDataset)

HTH, Thierry



ir. Thierry Onkelinx
Instituut voor natuur- en bosonderzoek / Research Institute for Nature and Forest Cel biometrie, methodologie en kwaliteitszorg / Section biometrics, methodology and quality assurance Gaverstraat 4
9500 Geraardsbergen
Belgium
tel. + 32 54/436 185
Thierry.Onkelinx_at_inbo.be
www.inbo.be

To call in the statistician after the experiment is done may be no more than asking him to perform a post-mortem examination: he may be able to say what the experiment died of.
~ Sir Ronald Aylmer Fisher

The plural of anecdote is not data.
~ Roger Brinner

The combination of some data and an aching desire for an answer does not ensure that a reasonable answer can be extracted from a given body of data.
~ John Tukey

-----Oorspronkelijk bericht-----
Van: r-help-bounces_at_r-project.org [mailto:r-help-bounces_at_r-project.org] Namens Dani Valverde Verzonden: zaterdag 24 mei 2008 11:01
Aan: R Help
Onderwerp: [R] Problems with lme

Hello,
I want to perform an lme on a database with this structure:

 ID Sequence Temperature Tumour                Error
 1        5           0                  1                     8.721872e-08
 1        5           0                  2                     8.695348e-08
 1        5           0                  3                     2.019604e-13
 1        5          37                  1                     2.746476e-01
 1        5          37                  2                     1.580415e-02
 1        5          37                  3                     1.586857e-01
 1        2           0                  1                     8.721872e-08
 1        2           0                  2                     8.695348e-08
 1        2           0                  3                     2.019604e-13
 1        2          37                  1                    2.746476e-01

I try to perform an lme with the lme() function of the lme package (first I attach the data frame with attach()). This is the code: lme(Error ~ ID + Sequence + Temperature + Tumour)

If I use this code, I get the following error:

Error in getGroups.data.frame(dataMix, groups) :   Invalid formula for groups

Can anyone help me?
Best,

Dani

-- 

Daniel Valverde Saubí

Grup de Biologia Molecular de Llevats
Facultat de Veterinària de la Universitat Autònoma de Barcelona
Edifici V, Campus UAB
08193 Cerdanyola del Vallès- SPAIN

Centro de Investigación Biomédica en Red
en Bioingeniería, Biomateriales y
Nanomedicina (CIBER-BBN)

Grup d'Aplicacions Biomèdiques de la RMN
Facultat de Biociències
Universitat Autònoma de Barcelona
Edifici Cs, Campus UAB
08193 Cerdanyola del Vallès- SPAIN
+34 93 5814126

______________________________________________
R-help_at_r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.

______________________________________________
R-help_at_r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
Received on Sat 24 May 2008 - 09:35:06 GMT

Archive maintained by Robert King, hosted by the discipline of statistics at the University of Newcastle, Australia.
Archive generated by hypermail 2.2.0, at Sat 24 May 2008 - 12:30:43 GMT.

Mailing list information is available at https://stat.ethz.ch/mailman/listinfo/r-help. Please read the posting guide before posting to the list.

list of date sections of archive