Displaying 3 results from an estimated 3 matches for "hepmeg".
2007 May 31
0
Using MIcombine for coxph fits
...u might be able to help with. I will work through a
reproducible example to demonstrate the issue.
First, make a dataset from the pbc dataset in the survival package
---------------
# Make a dataset
library(survival)
d <- pbc[,c('time','status','age','sex','hepmeg','platelet',
'trt', 'trig')]
d[d==-9] <- NA
d[,c(4,5,7)] <- lapply(d[,c(4,5,7)], FUN=as.factor)
str(d)
summary(d)
---------------
Second, since there is missing data for several (but not all) of the
variables, investigate the patterns.
---------------
library(...
2006 Sep 03
2
Running cox models
...oks/Biostatistics/chapter16.html
(Primary Biliary Cirrhosis data link at top of the page),
I'm using the following code:
--------------- start of code
library(survival)
liver <- scan("liver2.txt",list(age=0,albumin=0,alkphos=0,ascites=0,bili=0,
cholest=0,edema=0,edmadj=0,hepmeg=0,obstime=0,platelet=0,protime=0,
sex=0,sgot=0,spiders=0,stage=0,status=0,treatmnt=0,
triglyc=0,urinecu=0))
fit<-coxph(Surv(obstime,status)~bili+edmadj+albumin+protime+age,data=liver)
summary(fit)
----------------- End of code
but the answer is rather different from that in the...
2007 May 17
1
MICE for Cox model
...--
The problem that I encounter occurs when I try to use the "pool"
function to pool the fitted models into one general model. Here is some
code that reproduces the error using the pbc dataset.
------------
d <- pbc[,c('time','status','age','sex','hepmeg','platelet', 'trt',
'trig')]
d[d==-9] <- NA
d[,c(4,5,7)] <- lapply(d[,c(4,5,7)], FUN=as.factor)
str(d)
imp <- mice(d, m=10, maxit=10, diagnostics=T, seed=500,
defaultImputationMethod=c('norm', 'logreg', 'polyreg'))
fit <- cox.mids...