Displaying 20 results from an estimated 3000 matches similar to: "fitting a gamma frailty model (coxph)"
2005 Jul 21
1
output of variance estimate of random effect from a gamma frailty model using Coxph in R
Hi,
I have a question about the output for variance of random effect from a gamma
frailty model using coxph in R. Is it the vairance of frailties themselves or
variance of log frailties? Thanks.
Guanghui
2010 Apr 26
1
Interpreting output of coxph with frailty.gamma
Dear all,
this is probably a very silly question, but could anyone tell me what the
different parameters in a coxph model with a frailty.gamma term mean?
Specifically I have two questions:
(1) Compared to a "normal" coxph model, it seems that I obtain two standard
errors [se(coef) and se2].
What is the difference between those?
(2) Again compared to a "normal" coxph model,
2004 Nov 08
1
coxph models with frailty
Dear R users:
I'm generating the following survival data:
set.seed(123)
n=200 #sample size
x=rbinom(n,size=1,prob=.5) #binomial treatment
v=rgamma(n,shape=1,scale=1) #gamma frailty
w=rweibull(n,shape=1,scale=1) #Weibull deviates
b=-log(2) #treatment's slope
t=exp( -x*b -log(v) + log(w) ) #failure times
c=rep(1,n) #uncensored indicator
id=seq(1:n) #individual frailty indicator
2007 Apr 20
1
Approaches of Frailty estimation: coxme vs coxph(...frailty(id, dist='gauss'))
Dear List,
In documents (Therneau, 2003 : On mixed-effect cox
models, ...), as far as I came to know, coxme penalize
the partial likelihood (Ripatti, Palmgren, 2000) where
as frailtyPenal (in frailtypack package) uses the
penalized the full likelihood approach (Rondeau et al,
2003).
How, then, coxme and coxph(...frailty(id,
dist='gauss')) differs? Just the coding algorithm, or
in
2011 Jun 25
2
cluster() or frailty() in coxph
Dear List,
Can anyone please explain the difference between cluster() and
frailty() in a coxph? I am a bit puzzled about it. Would appreciate
any useful reference or direction.
cheers,
Ehsan
> marginal.model <- coxph(Surv(time, status) ~ rx + cluster(litter), rats)
> frailty.model <- coxph(Surv(time, status) ~ rx + frailty(litter), rats)
> marginal.model
Call:
coxph(formula =
2002 Oct 08
2
Frailty and coxph
Does someone know the rules by which 'coxph' returns 'frail', the
predicted frailty terms? In my test function:
-----------------------------------------------
fr <- function(){
#testing(frailty terms in 'survival'
require(survival)
dat <- data.frame(exit = 1:6,
event = rep(1, 6),
x = rep(c(0, 1), 3),
2003 Aug 04
1
coxph and frailty
Hi:
I have a few clarification questions about the elements returned by
the coxph function used in conjuction with a frailty term.
I create the following group variable:
group <- NULL
group[id<50] <- 1
group[id>=50 & id<100] <- 2
group[id>=100 & id<150] <- 3
group[id>=150 & id<200] <- 4
group[id>=200 & id<250] <- 5
group[id>=250
2009 Aug 31
2
How to extract the theta values from coxph frailty models
Hello,
I am working on the frailty model using coxph functions. I am running
some simulations and want to store the variance of frailty (theta)
values from each simulation result. Can anyone help me how to extract
the theta values from the results. I appreciate any help.
Thanks
Shankar Viswanathan
2007 Apr 17
3
Extracting approximate Wald test (Chisq) from coxph(..frailty)
Dear List,
How do I extract the approximate Wald test for the
frailty (in the following example 17.89 value)?
What about the P-values, other Chisq, DF, se(coef) and
se2? How can they be extracted?
######################################################>
kfitm1
Call:
coxph(formula = Surv(time, status) ~ age + sex +
disease + frailty(id,
dist = "gauss"), data = kidney)
2009 Jan 07
0
Frailty by strata interactions in coxph (or coxme)?
Hello,
I was hoping that someone could answer a few questions for me (the background is given below):
1) Can the coxph accept an interaction between a covariate and a frailty term
2) If so, is it possible to
a) test the model in which the covariate and the frailty appear as main terms using the penalized likelihood (for gaussian/t frailties)
b)augment model 1) by stratifying on the variable that
2006 Sep 22
0
$theta of frailty in coxph
Dear all,
Does the frailty.object$history[[1]]$theta returns the Variance of random
effect?
Why is the value different? Here is an example with kidney data:
> library(survival)
> data(kidney)
> frailty.object<-coxph(Surv(time, status)~ age + sex + disease +
frailty(id), kidney)
> frailty.object
Call:
coxph(formula = Surv(time, status) ~ age + sex + disease + frailty(id),
data
2006 Sep 21
0
frailty in coxph
Dear all,
I have been doing some frailty calculations and been facing some
difficulties.
I can extract coefficients, value of theta and the following things
library(survival)
fit<-coxph(Surv(time,status)~covariate+frailty(group), data=simulated.data)
fit$coef
fit$history[[1]]$theta
fit$history[[1]]$c.loglik
fit$var
fit$var2
from a frailty included coxph object:
but how can i know what other
2003 May 19
1
survit function and cox model with frailty
Hi:
I have a question about the use of the survfit function after the
estimation of a cox proportional hazard model with a frailty term. My goal
is to estimate expected survival probabilities while controlling for the
group-specific frailty term.
First, I estimate a model of the following form:
model1 <- coxph(Surv(t0, t, d) ~ x1 + x2 + frailty(id), na.action=na.exclude,
2009 Jun 24
1
Coxph frailty model counting process error X matrix deemed singular
Hello,
I am currently trying to simulate data and analyze it using the frailty option in the coxph function. I am working with recurrent event data, using counting process notation. Occasionally, (about 1 in every 100 simulations) I get the following warning:
Error in coxph(Surv(start, end, censorind) ~ binary + uniform + frailty(subject, :
X matrix deemed to be singular; variable 2
My
2007 Mar 14
0
Wald test and frailty models in coxph
Dear R members,
I am new in using frailty models in survival analyses and am getting
some contrasting results when I compare the Wald and likelihood ratio
tests provided by the r output.
I am testing the survivorship of different sunflower interspecific
crosses using cytoplasm (Cyt), Pollen and the interaction Cyt*Pollen
as fixed effects, and sub-block as a random effect. I stratified
2005 Sep 07
1
Survival analysis with COXPH
Dear all,
I would have some questions on the coxph function for survival analysis,
which I use with frailty terms.
My model is:
mdcox<-coxph(Surv(time,censor)~ gender + age + frailty(area, dist='gauss'),
data)
I have a very large proportion of censored observations.
- If I understand correctly, the function mdcox$frail will return the random
effect estimated for each group on the
2006 Nov 07
1
Extracting parameters for Gamma Distribution
I'm doing a cox regression with frailty:
model <- coxph(Surv(Start,Stop,Terminated)~ X + frailty(id),table)
I understand that model$frail returns the group level frailty
terms. Does this mean this is the average of the frailty
values for the respective groups? Also, if I'm fitting it to
a gamma frailty, how do I extract the rate and scale
parameters for the different gamma
2013 Nov 04
0
Fwd: Re: How to obtain nonparametric baseline hazard estimates in the gamma frailty model?
-------- Original Message --------
Subject: Re: How to obtain nonparametric baseline hazard estimates in the gamma frailty model?
Date: Mon, 04 Nov 2013 17:27:04 -0600
From: Terry Therneau <therneau.terry at mayo.edu>
To: Y <yuhanusa at gmail.com>
The cumulative hazard is just -log(sfit$surv).
The hazard is essentially a density estimate, and that is much harder. You'll notice
2004 Nov 17
1
frailty and time-dependent covariate
Hello,
I'm trying to estimate a cox model with a frailty variable and time-dependent covariate (below there is the statement I use and the error message). It's seems to be impossible, because every time I add the time-dependent covariate the model doesn't converge. Instead, if I estimate the same model without the time-dependent covariate it's converge. I'd like knowing if
2006 Sep 19
0
How to interpret these results from a simple gamma-frailty model
Dear R users,
I'm trying to fit a gamma-frailty model on a simulated dataset, with 6 covariates, and I'm running into some results I do not understand. I constructed an example from my simulation code, where I fit a coxph model without frailty (M1) and with frailty (M2) on a number of data samples with a varying degree of heterogeneity (I'm running R 2.3.1, running takes ~1 min).