I have a generalized linear model to solve. I used package "geepack". When
I use the correlation structure "unstructured", I get a messeage that- R
GUI front-end has stopped working. Why this happens? What is the solution?
The r codes are as follows:
a<-read.table("d:/bmt.txt",header=T")
cutoffs<-c(1.75,3.34,5.09,6.93,9.11,11.8,15.3,20.1,23.7,27.3,30.9,34.3,37.6,41.5,46.5,51.1,57,63,73.1,93.68)
a$icr <- REL
library(pseudo)
pseudo <- pseudoci(time=a$TIME,event=a$icr,tmax=cutoffs)
#rearranging data
b <- NULL
for(it in 1:length(pseudo$time)){
b <- rbind(b,cbind(a,pseudo = pseudo$pseudo[[1]][,it],
tpseudo = pseudo$time[it],id=1:nrow(a)))
}
b <- b[order(b$id),]
library(geepack)
fit3 <- geese(pseudo ~ as.factor(AGE) + as.factor(WEIGHT) +
as.factor(HEIGHT) ,
data =b, id=id, jack = TRUE, scale.fix=TRUE,
mean.link = "logit", corstr="unstructured")
#The results using the AJ variance estimate
cbind(mean = round(fit3$beta,4), SD = round(sqrt(diag(fit3$vbeta.ajs)),4),
Z = round(fit3$beta/sqrt(diag(fit3$vbeta.ajs)),4),
PVal = round(2-2*pnorm(abs(fit3$beta/sqrt(diag(fit3$vbeta.ajs)))),4))
[[alternative HTML version deleted]]
______________________________________________
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.
