I have included data at the bottom of this email. It can be read in by
highlighting the data and then using this command: dat <-
read.table("clipboard", header = TRUE,sep="\t")
I can obtain solutions with both of these:library(gee) fit.gee<-gee(score ~ chem + time, id=id, family=gaussian,corstr="exchangeable",data=dat) and library(yags) fit.yags <- yags(score ~ chem + time, id=id, family=gaussian,corstr="exchangeable",data=dat,alphainit=0.05) However, I am making a mistake with: library(geepack) fit.geese <- geese(score ~ chem + time, id=id, family=gaussian,corstr="exch",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 22 24 18 22 1000 1000 1055 950 3 18 1055 2 1 22 24 18 22 1000 1000 1055 950 6 22 950 3 1 14 10 24 10 1000 1999 800 1700 0 14 1000 3 1 14 10 24 10 1000 1999 800 1700 2 10 1999 3 1 14 10 24 10 1000 1999 800 1700 3 24 800 3 1 14 10 24 10 1000 1999 800 1700 6 10 1700 4 1 38 34 32 24 1000 1050 1400 0 38 1000 4 1 38 34 32 24 1000 1050 1400 2 34 1050 4 1 38 34 32 24 1000 1050 1400 3 32 4 1 38 34 32 24 1000 1050 1400 6 24 1400 5 1 25 29 25 29 1000 1020 1040 1045 0 25 1000 5 1 25 29 25 29 1000 1020 1040 1045 2 29 1020 5 1 25 29 25 29 1000 1020 1040 1045 3 25 1040 5 1 25 29 25 29 1000 1020 1040 1045 6 29 1045 6 1 30 28 26 14 1000 1100 1180 1500 0 30 1000 6 1 30 28 26 14 1000 1100 1180 1500 2 28 1100 6 1 30 28 26 14 1000 1100 1180 1500 3 26 1180 6 1 30 28 26 14 1000 1100 1180 1500 6 14 1500 7 0 20 15 21 20 1000 1200 1200 1000 0 20 1000 7 0 20 15 21 20 1000 1200 1200 1000 2 15 1200 7 0 20 15 21 20 1000 1200 1200 1000 3 21 1200 7 0 20 15 21 20 1000 1200 1200 1000 6 20 1000 8 0 21 27 19 27 1000 900 1075 900 0 21 1000 8 0 21 27 19 27 1000 900 1075 900 2 27 900 8 0 21 27 19 27 1000 900 1075 900 3 19 1075 8 0 21 27 19 27 1000 900 1075 900 6 27 900 9 0 15 22 22 20 1000 1000 700 0 15 1000 9 0 15 22 22 20 1000 1000 700 2 22 9 0 15 22 22 20 1000 1000 700 3 22 1000 9 0 15 22 22 20 1000 1000 700 6 20 700 10 0 39 39 34 1000 950 1033 1025 0 39 1000 10 0 39 39 34 1000 950 1033 1025 2 950 10 0 39 39 34 1000 950 1033 1025 3 39 1033 10 0 39 39 34 1000 950 1033 1025 6 34 1025 11 0 27 27 31 22 1000 950 910 1050 0 27 1000 11 0 27 27 31 22 1000 950 910 1050 2 27 950 11 0 27 27 31 22 1000 950 910 1050 3 31 910 11 0 27 27 31 22 1000 950 910 1050 6 22 1050 12 0 28 24 33 1000 1015 985 0 28 1000 12 0 28 24 33 1000 1015 985 2 24 1015 12 0 28 24 33 1000 1015 985 3 33 985 12 0 28 24 33 1000 1015 985 6 ______________________________________________ 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.
