Displaying 7 results from an estimated 7 matches for "edema".
Did you mean:
edeca
2008 Nov 21
1
Discrepancy in the regression coefficients for Cox regression - PBC data set
...set (given in the "survival" package), the regression coefficients
do not agree with the results presented in Table 4.6.3 (p. 195) of Fleming &
Harrington's book.
library(survival)
data(pbc)
ans.cox <- coxph(Surv(time, status) ~ log(bili) + log(alb) + age +
log(protime) + edema)
ans.cox
> ans.cox <- coxph(Surv(time, status) ~ log(bili) + log(alb) + age +
log(protime) + edema)
> ans.cox
Call:
coxph(formula = Surv(time, status) ~ log(bili) + log(alb) + age +
log(protime) + edema)
coef exp(coef) se(coef) z p
log(bili) 0.8975...
2008 Nov 24
1
Discrepancy in the PBC data set
...at/upload/therneau_upload/pbc.d
at"
pbc <- read.table(pbcurl, header=F,
col.names=c('id', 'time', 'status', 'trt', 'age', 'sex',
'ascites', 'hepato', 'spiders', 'edema',
'bili', 'chol', 'albumin', 'copper',
'alk.phos', 'ast', 'trig', 'platelet',
'protime', 'stage'),...
2009 Jul 13
0
pbc data
...with the pbc data set. I just want to apply simple cox regression in the data set. I am a beginner in R but I don't think I am doing anything wrong.
I have the book of Fleming and Harrington 1990. I perform cox regression by typing:
out<- coxph(Surv(times/365,status)~log(bili)+log(proth)+edema+log(albumin)+age)
out
Call:
coxph(formula = Surv(times/365, status) ~ log(bili) + log(proth) +
edema + log(albumin) + age)
coef exp(coef) se(coef) z p
log(bili) 0.8636 2.3716 0.08294 10.41 0.0e+00
log(proth) 2.3868 10.8791 0.76851 3.11 1.9e-03
edema...
2003 Aug 04
1
coxph and frailty
...t;=250 & id<300] <- 6
group[id>=300] <- 7
I estimate the following model, using the Pbc data (with time-varying
covariates) from Therneau and Grambsch's book:
fitf <- coxph(Surv(start,stop,event==2) ~ age + log(bili) + log(protime) +
log(albumin) + edema + frailty(group),
na.action=na.exclude, data=Pbcseq)
Then I obtain:
> fitf[10]
$frail
[1] 0.06273372 0.16192093 0.10050877 0.37716999 -0.20853156 -0.71887977
[7] -0.10922275
And:
> fitf[11]
$fvar
[1] 0.03631634 0.03261972 0.03750465 0.04459689 0.05113802 0.08...
2006 Sep 03
2
Running cox models
....edu/~heagerty/Books/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...
2010 Dec 13
1
Multivariate binary response analysis
...The setting: Data from an assay, wherein 6-hours-post-fertilization zebrafish embryos (n=20, say) are exposed in a vial to a chemical (replicated 3 times, say), and 5 days later observed for the presence/absence (1/0) of defects in several organ systems (yolk sac, jaw, snout, body axis, pericardial edema, etc.) for each fish. The assay can be used as a screen (in which case a single response 'any' is first analyzed) as well as for comparing the responses of different classes of chemicals (a MANOVA-type analysis). Interest is focused on where response(s) are occurring, any associations among...
2008 Mar 02
0
new to latex to pdf
...='constant', f=0)$y
}
describe.survival <- function(y) {
km <- kmsurv(y, c(2,5))
c('2 Year'=km[1], '5 Year'=km[2], 'Mean, y'=sum(y[,1])/sum(y[,2]))
}
S <- with(pbc, Surv(fu.days/365.25, status))
s1 <- summary(S ~ age + albumin + ascites + bili + drug + edema + chol,
fun=describe.survival, data=pbc)
for(w in 1:2) {
if(w==1) setpdf(f1a,sublines=1,h=5.25) else
setpdf(f1b,sublines=1,h=5)
plot(s1, which=if(w==1)1:2 else 3,
cex.labels=.7, cex.group.labels=.7*1.15, subtitles=T, main='',
pch=if(w==2) 16 else c('2','5'), # 16=solid cir...