Displaying 20 results from an estimated 600 matches similar to: "How to interpret these results from a simple gamma-frailty model"
2012 Jan 05
2
Bayesian estimate of prevalence with an imperfect test
Hi all!
I'm new to this forum so please excuse me if I don't conform perfectly to
the protocols on this board!
I'm trying to get an estimate of true prevalence based upon results from an
imperfect test. I have various estimates of se/sp which could inform my
priors (at least upper and lower limits even if with a uniform distribution)
and found the following code on this website..
2001 Jul 28
2
Re: [S] Labels wrong with lrm
Dear Jan,
Thank you very much for your excellent description of the
problem and the self-contained test code. This is a
problem that I've been meaning to either document better
or solve for some time. The root of the problem is with
the builtin S-Plus terms.inner function:
> attr(terms.inner(asthma ~ pol(age,kx) + smok),'variables')
expression(age, kx, smok)
You can see that
2005 Mar 16
1
Code to replace nested for loops
Dear list members,
How can I replace the nested for loops at then end of the script
below with more efficient code?
# Begin script__________________________________________________
# Dichotomous scores for 100 respondents on 3 items with
# probabilities of a correct response = .6, .4, and .7,
# respectively
x1 <- rbinom(100,1,.6)
x2 <- rbinom(100,1,.4)
x3 <- rbinom(100,1,.7)
#
2009 Jun 04
0
Dropping terms from regression w/ poly()
Hello r-help,
I'm fitting a model with lm() and using the orthogonal polynomials
from poly() as my basis:
dat <- read.csv("ConsolidatedData.csv", header=TRUE)
attach(dat)
nrows <- 1925
Rad <- poly(Radius, 2)
ntheta <- 14
Theta <- poly(T.Angle..deg., ntheta)
nbeta <- 4
Beta <- poly(B.Beta..deg., nbeta)
model.1 <- lm( Measurement ~ Block + Rad + Theta + Beta
2009 Feb 26
1
error message and convergence issues in fitting glmer in package lme4
I'm resending this message because I did not include a subject line in my first posting.
Apologies for the inconvenience!
Tanja
> Hello,
>
> I'm trying to fit a generalized linear mixed model to estimate diabetes prevalence at US county level. To do this I'm using the glmer() function in package lme4. I can fit relatively simple models (i.e. few covariates) but when
2006 Mar 28
3
Running text app without X
I'm sure this question comes up waaaay to often in this
list, and I apologize if I've missed the obvious answer.
I did spend the last two hours looking for a solution
and trying various things, but to no avail.
I'm trying to run a command-line app (text-only). It's a
cross-compiler tool for which we only have Windows binaries.
The app works fine, but when I run it (i.e. to
2003 May 07
0
frailty models in survreg() -- survival package (PR#2933)
I am confused on how the log-likelihood is calculated in a parametric
survival problem with frailty. I see a contradiction in the frailty() help
file vs. the source code of frailty.gamma(), frailty.gaussian() and
frailty.t().
The function frailty.gaussian() appears to calculate the penalty as the
negative log-density of independent Gaussian variables, as one would
expect:
>
2009 Feb 26
1
(no subject)
Hello,
I'm trying to fit a generalized linear mixed model to estimate diabetes prevalence at US county level. To do this I'm using the glmer() function in package lme4. I can fit relatively simple models (i.e. few covariates) but when expanding the number of covariates I usually encounter the following error message.
gm8 <-
2003 May 07
0
Re: frailty models in survreg() -- survival package (PR#2934)
On Tue, 6 May 2003, Jerome Asselin wrote:
>
> I am confused on how the log-likelihood is calculated in a parametric
> survival problem with frailty. I see a contradiction in the frailty() help
> file vs. the source code of frailty.gamma(), frailty.gaussian() and
> frailty.t().
>
> The function frailty.gaussian() appears to calculate the penalty as the
> negative
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
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
2011 Apr 05
0
frailty
Hi R-users
I spend a lot of time searching on the web but I didn?t found a clear
answer.
I have some doubts with 'frailty' function of 'survival' package.
The following model with the function R ?coxph? was fitted:
modx <- coxph(Surv(to_stroke, stroke) ~ age + sbp + dbp + sex +
frailty(center,distribution = "gamma", method='aic'), data=datax)
Then I get
2007 Apr 08
0
Simulation of the Frailty of the Cox PH model
Dear R-list users,
I am trying to do simulation of survival data to
enable it to run under frailty option. Below is the
function a that I am using. My questions are:
1. How do I modify it to get bigger (hopefully
significant) value of Variance of random effect?
2. What changes do I have to make in the function to
run it under correlated frailty model? (may be in
kinship package)
3. Is there
2006 Sep 21
0
Any examples of a frailty model actually used for prediction ?
Hi everyone,
I'm looking for any examples of useful frailty models, in particular any situation in which a cox proportional hazards model with frailty outperforms a regular cox proportional hazards model with respect to prediction of the time to event (or the X-year risk of an event). I have defined my own gamma-frailty cox PH model in R but on my simulated data sample it does not predict any
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
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 =
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
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),
2005 May 31
1
Shared Frailty in survival package (left truncation, time-dep. covariates)
Dear list,
I want o fit a shared gamma frailty model with the frailty specification in the survival package.
I have partly left-truncated data and time-dependent covariates. Is it possible to
combine these two things in the frailty function. Or are the results wrong if I use data in the start-stop-formulation which account for delayed entry?
Is the frailty distribution updated in the