Displaying 4 results from an estimated 4 matches for "gee2".
Did you mean:
gee
2010 Feb 10
1
using step() with package geepack
...t;)).
The package does provide an anova method for its model objects, and update() seems to work:
library(geepack)
data(dietox) # the example from ?geeglm
# GEE model, using gaussian instead of poisson:
gee1 <- geeglm(Weight~1, data=dietox, id=Pig, family=gaussian, corstr="ar1")
gee2 <- update(gee1, ".~.+Time")
anova(gee1, gee2) # gives a p-value for adding the Time variable
step(gee2, direction="forward") # fails
cheers,
Steve
[[alternative HTML version deleted]]
2010 May 20
0
Installation problem with "geepack"
...or begin() const { return v_;}
^
../inst/include/tnt/vec.h(153): warning #858: type qualifier on return type is meaningless
const iterator end() const { return v_ + n_; }
^
icpc -I/usr/global/R/2.11.0/lib64/R/include -I../inst/include -I/usr/local/include -mp -fpic -g -O2 -c gee2.cc -o gee2.o
icpc: command line remark #10010: option '-mp' is deprecated and will be removed in a future release. See '-help deprecated'
../inst/include/tnt/vec.h(152): warning #858: type qualifier on return type is meaningless
const iterator begin() const { return v_;}...
2010 Apr 24
0
'geepack' and 'gee' package outputs
Hi, having used both the gee pacakge and the geepack package, i am unsure of
how to interpret the results.
Here are the results from the geeglm function from the geepack package
> gee2<-geeglm(data$erythema~data$product, data = data, id=subject,
> family=binomial, corstr="independence")
Warning message:
In model.response(mf, "numeric") :
using type="numeric" with a factor response will be ignored
> summary(gee2)
Call:
geeglm(formula = dat...
2008 Mar 05
1
problem with geepack
...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 rep...