[R] problem with geepack

[wiersmap]

Hi,

I think in such cases you need to use the 'wave' parameter (wave = Date),
but I must admit I found contradictory comments on how missing values are
treated in geepack. Does someone know what geepack is doing with missing
values?

Popko



Giacomo Santini wrote:
> 
> Hi all
> 
> I am analyzing a data set containing information about the behaviour of 
> marine molluscs on a vertical wall. Since I have replicate observations 
> on the same individuals I was thinking to use the geepack library.
> 
> The data are organised in a dataframe with the following variables
> 
> Date = date of sampling,
> Size = dimensions (mm)
> Activity duration of activity (min)
> Water = duration of splashing by waves
> Hgt = resting eight of each specimen before activity begin
> Individual = a code indicating the id of the specimen.
> 
> I have up to 12 replicate observations for individual. Some observation 
> are missing and I organized the data frame to have exactly 12 rows for 
> each specimen, with NAs where there is a missing observation.
> 
> The following model worked fine:
> 
> gee1<-geese(Activity~Water, id=Individual, data=dataF, family=gaussian)
> 
> but when I use other variables e.g
> 
> gee2<-geese(Activity~Hgt+Size+Water, id=Individual, data=dataF, 
> family=gaussian)
> 
> I get the error message
> 
> "Error in geese.fit(x, y, id, offset, soffset, w, waves, zsca, zcor, 
> corp,  :
>   nrow(zsca) and length(y) not match"
> 
> which I am not able to understand.
> 
> The same problem has been reported in the list in 2006, but I have not 
> found any response to it.
> 
> Any suggestion?
> 
> 
> Giacomo
> 
> ______________________________________________
> 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.
> 
> 

-- 
View this message in context: http://www.nabble.com/problem-with--geepack-tp15851828p21601979.html
Sent from the R help mailing list archive at Nabble.com.