Displaying 2 results from an estimated 2 matches for "pseudoci".
Did you mean:
pseudoc
2013 Apr 09
2
R crash
...s 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) +...
2013 Mar 07
1
Comparing Cox model with Competing Risk model
...model the cumulative
incidence function is used directly. I used the jackknife (pseudovalue) of
the cumulative incidence function for each cause (AIDS or Cancer) in a
generalized estimating equation. I used the following code:
>library(pseudo)
>cutoffs<-c(0.5,5.5,10.5,15.5)
>pseudo<-pseudoci(time,event,tmax=cutoffs)
>library(geepack)
>fit<-geese(pseudo~as.factor(Age)+as.factor(Sex),data=b,id=id, jack = TRUE,
scale.fix=TRUE,
mean.link = "logit", corstr="independence").
I want to know whether I am doing the right thing?
[[alternative HTML version deleted...