[R] lmer() and "$ operator is invalid for atomic vectors"

ONKELINX, Thierry Thierry.ONKELINX at inbo.be
Fri Jul 31 10:28:05 CEST 2009


Dear Angela,

I think you have to omit the last lamba in the lmer code. Is there a
reason for putting it there?

lmer.lamda <- lmer(lamda ~ Age*Item + ( 1 | Subject), data=new4) 

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 Angela Radulescu
Verzonden: donderdag 30 juli 2009 23:38
Aan: r-help at r-project.org
Onderwerp: [R] lmer() and "$ operator is invalid for atomic vectors"


Hi all,

I am a bit mystified by this error message that I get when I try to
apply
lmer() to a simple dataset with one between factor (age) and one within
factor (item): "$ operator is invalid for atomic vectors"

I'll just provide the code, because I don't see where the problem is: 


library(lme4)
options(contrasts=c("contr.helmert","contr.poly"))


data =
read.csv("C:/Users/Angela/Desktop/Summer_at_CDS/aging_JDM_Martine/datase
ts/lamda.data_07.27.09.csv",
header=TRUE)
attach(data)

new<-data.frame(lamda, age.o_1, Item.l2_2,Subject) new

Subject <- factor(Subject)
Age <- factor(age.o_1,levels=c("0","1"),labels=c("young","old"))
Item <- factor(Item.l2_2,levels=c("1","2"),labels=c("lamda1","lamda2"))

new2<-data.frame(lamda, Subject, Age, Item)
new2

new3 <- na.omit(new2)
new3
new4=as.data.frame(new3)


> (lmer.lamda <- lmer(lamda ~ Age*Item + ( 1 | Subject), lamda, 
> data=new4))
Error in family$family : $ operator is invalid for atomic vectors


Here are the first 25 values of the dataframe lmer() should operate on: 

 lamda Subject   Age   Item
1     1.6000000       1   old lamda1
2     1.0434783       2   old lamda1
3     0.9600000       3   old lamda1
4     3.4285714       4   old lamda1
5     1.1428571       5   old lamda1
6     2.1818182       6   old lamda1
7     0.9600000       7   old lamda1
8     1.8461538       8   old lamda1
9     2.1818182       9   old lamda1
10    1.8461538      10   old lamda1
11    1.8461538      11   old lamda1
12    1.8461538      12   old lamda1
13    2.6666667      13   old lamda1
14    2.6666667      14   old lamda1
15    1.8461538      15   old lamda1
16    1.0434783      16   old lamda1
17    1.8461538      17   old lamda1
18    2.6666667      18   old lamda1
19    2.6666667      19   old lamda1
20    8.0000000      20   old lamda1
21    1.8461538      21   old lamda1
22    8.0000000      22   old lamda1
23    1.8461538      23   old lamda1
24    2.6666667      24   old lamda1
25    2.1818182      25   old lamda1


Any help would be much appreciated this one. I am in the dark and tried
pretty much everything... 

Many thanks,
Angela 

--
View this message in context:
http://www.nabble.com/lmer%28%29-and-%22%24-operator-is-invalid-for-atom
ic-vectors%22-tp24747770p24747770.html
Sent from the R help mailing list archive at Nabble.com.

______________________________________________
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.

Dit bericht en eventuele bijlagen geven enkel de visie van de schrijver weer 
en binden het INBO onder geen enkel beding, zolang dit bericht niet bevestigd is
door een geldig ondertekend document. The views expressed in  this message 
and any annex are purely those of the writer and may not be regarded as stating 
an official position of INBO, as long as the message is not confirmed by a duly 
signed document.




More information about the R-help mailing list