I get an error message when trying to specify a scale offset using
geese() in library(geepack). Here is a dummy example:
v <- rep( 1:2, each = 100 )
x <- rnorm( 200 ) * sqrt( v )
i <- 1:200
summary( geese( x~1, id=i ) ) # works ok
summary( geese( x~1, id=i, soffset=v ) ) # fails
Error in geese.fi
Dear Celine,
Please excuse the delay in responding over the holiday period. Here
are my answer to your questions. If you are having further
difficulties or these answers don't seem right do let me know.
(1) Yes, I recommend splitting your data by chromosome. If you perform
a genomewide analysis,
2 matches
Mail list logo