Displaying 2 results from an estimated 2 matches for "chem3".
Did you mean:
chem
2010 Dec 13
0
Help with lme for split plot
...icate ponds for each treatment
data table looks like this
pond side sampperiod response
1 cont 1 x1
1 cont 2 x2
1 cont 3 x3
1 chem1 1 x1
1 chem1 2 x2
1 chem1 3 x3
2 cont 1 x1
2 cont 2 x2
2 cont 3 x3
2 chem2 1 x1
2 chem2 2 x2
2 chem2 3 x3
3 cont 1 x1
3 cont 2 x2
3 cont 3 x3
3 chem3 1 x1
3 chem3 2 x2
3 chem3 3 x3
I can fit a model using aov
m1<-aov(response~side*sampperiod + Error(pond/side/sampperiod,
data=datafile)
Since there is unequal replication, 4 replicates for one chemical treatment
and 3 for the other two, and there is one missing value, lme from the nlme
pa...
2008 Oct 29
2
call works with gee and yags, but not geepack
...t;,data=dat)
I obtain the following error:
Error in geese.fit(x, y, id, offset, soffset, w, waves, zsca, zcor, corp, :
nrow(zsca) and length(y) not match
Could someone tell me what I have done incorrectly. Thanks for your
time, Juliet.
Data Below:
id treat time1 time2 time3 time4 chem1 chem2 chem3 chem4 time score chem
1 1 20 18 15 15 1000 1100 1200 1300 0 20 1000
1 1 20 18 15 15 1000 1100 1200 1300 2 18 1100
1 1 20 18 15 15 1000 1100 1200 1300 3 15 1200
1 1 20 18 15 15 1000 1100 1200 1300 6 15 1300
2 1 22 24 18 22 1000 1000 1055 950 0 22 1000
2 1 22 24 18 22 1000 1000 1055 950 2 24 1000
2 1...