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]]

