[R] call works with gee and yags, but not geepack

2008-10-29 Thread Juliet Hannah
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,

Re: [R] call works with gee and yags, but not geepack

2008-10-29 Thread Juliet Hannah
Sorry. I did not output the NAs correctly. dat - read.table(clipboard, header = TRUE) 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

Re: [R] call works with gee and yags, but not geepack

2008-10-29 Thread Rolf Turner
On 30/10/2008, at 9:08 AM, Juliet Hannah wrote: 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)