[R] R crash

[Duncan Murdoch]

On 13-04-09 2:41 PM, Tasnuva Tabassum wrote:
> 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?

Presumably geepack has a bug.  You should simplify your example as much 
as you can and send it to the maintainer of that package.

Duncan Murdoch

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