Hello,
I have been using the function lme() from package 'nlme' for several months
now without any problems. Suddenly, it cannot find a factor in my data.
Is this a new bug of some kind? My code and output are below.
Thanks for your help!
-Steve Pfeiffer
require(xlsx)
require(nlme)
Dat.middle<-read.xlsx(
"C:\\Users\\S\\Google Drive\\RESEARCH (flash
drive)\\Data\\sm reorganized 2.xlsx",
sheetName="Middle R Friendly",colIndex=1:5)
class(Dat.middle)
[1] "data.frame"
dim(Dat.middle)
[1] 380 5
names(Dat.middle)
[1] "Honeysuckle" "Trenching" "Plot" "DaysSinceEvent"
"soil.moisture"
Dat.middle[1:10,]
Honeysuckle Trenching Plot DaysSinceEvent soil.moisture
1 Ctrl tr TE.tr
0 15.21690
2 Ctrl tr TE.tr
1 16.68770
3 Ctrl tr TE.tr 2
16.11459
4 Ctrl tr TE.tr 3
14.57441
5 Ctrl tr TE.tr 4
15.16364
6 Ctrl tr TE.tr 5
15.32327
7 Ctrl tr TE.tr 8
13.26730
8 Ctrl tr TE.tr 9
11.81761
9 Ctrl tr TE.tr 10
13.87007
10 Ctrl tr TE.tr
11 12.26743
fit.1<-lme(soil.moisture ~ Trenching + DaysSinceEvent,
+ Trenching:DaysSinceEvent,
random = ~ DaysSinceEvent | Plot,
data=Dat.middle,
method="ML"
)
Error in lme.formula(sm ~ Trenching + DaysSinceEvent,
+Trenching:DaysSinceEvent, :
object 'Trenching' not found
class(Dat.middle$Trenching) #Oh, look, here is the factor it said it
couldn't find!
[1] "factor"
length(Dat.middle$Trenching)
[1] 380
Dat.middle$Trenching[1:10]
[1] tr tr tr tr tr tr tr tr tr tr
Levels: tr un
[[alternative HTML version deleted]]
______________________________________________
R-help@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.
