similar to: lognormal frailty in frailtypack

Dear All, We are happy to announce, after a long gestation, the release of the new version of FRAILTYPACK (version 2.2-9.5) which is now available from CRAN. The package fit general frailty models using penalized likelihood estimation, for clustered or recurrent events. For instance : -- ADDITIVE FRAILTY MODELS for proportional hazards models with two correlated random effects (intercept

Dear All, We are happy to announce, after a long gestation, the release of the new version of FRAILTYPACK (version 2.2-9.5) which is now available from CRAN. The package fit general frailty models using penalized likelihood estimation, for clustered or recurrent events. For instance : -- ADDITIVE FRAILTY MODELS for proportional hazards models with two correlated random effects (intercept

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

Dear R users, We are pleased to tell you that "FRAILTYPACK" has been updated. "FRAILTYPACK" stands now for general frailty models estimated with a semi-parametrical penalized likelihood, but also with a parametrical approach. In case of comments/corrections/remarks/suggestions -- which are very welcome --please contact the maintainer directly. Kind regards, The

Dear R users, We are pleased to tell you that "FRAILTYPACK" has been updated. "FRAILTYPACK" stands now for general frailty models estimated with a semi-parametrical penalized likelihood, but also with a parametrical approach. In case of comments/corrections/remarks/suggestions -- which are very welcome --please contact the maintainer directly. Kind regards, The

Dear, I want to analyze two-level survival data using a shared frailty model, for which I want to use the R package 'Frailtypack", proposed by Rondeau et al. The dataset was built using SAS software. I also tried to change the format using SPSS and Excell. My (reduced) dataset has following column names: ID entry time status family var1 I used following command: >

Hi, when i run the following code i get this massege: "The instruction at 0x######## reference memory at 0x#######, the memory cannot be "read". and then i have to close R. what is the problem and how can i solve it? thanks in advance Avi my code # frailtypack library(frailtypack) cgd.ag <- read.csv("C:/rfiles/RE/cgd.csv") cgd.nfm <-frailtyPenal(Surv(TStart,

Dear R-listers, I have a very strange problem. I made a package (under Windows and Linux). The package passed the R CMD Check without problem. Then, I installed the package and executed a function which calls to a 'dll' mod<-frailtyPenal(Surv(time,status)~sex+age+cluster(id), + n.knots=8,kappa1=10000,data=kidney) mod Call: frailtyPenal(formula = Surv(time, status)

A new version of the package FRAILTYPACK is now available on CRAN. -- possibility to fit now a Shared and a Joint Frailty model with a log-normal distribution for the random effects. -- possibility to deal with interval-censored data (for a shared frailty model) -- possibility to fit a joint frailty model for clustered data For more details see the corresponding NEWS files in the pkgs. We are

A new version of the package FRAILTYPACK is now available on CRAN. -- possibility to fit now a Shared and a Joint Frailty model with a log-normal distribution for the random effects. -- possibility to deal with interval-censored data (for a shared frailty model) -- possibility to fit a joint frailty model for clustered data For more details see the corresponding NEWS files in the pkgs. We are

Dear R-help, I am trying to estimate a Cox model with nested effects, or better h(t,v,w)=v*w*h0(t)*exp(B'x) where h(t,v,w) is the individual hazard function w and v are both frailty terms (gamma or normal distributed) I have 12 clusters and for each one of them I would like to associate a realization of v, while w is a random effect for the whole population. At the population level

I can't comment on frailtypack issues, but would like to mention that coxme will handle nested models, contrary to the statement below that "frailtypack is perhaps the only .... for nested survival data". To reprise the original post's model cgd.nfm <- coxme(Surv(Tstart, Tstop, Status) ~ Treatment + (1 | Center/ID), data=cgd.ag) And a note to the poster-- you should

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: >

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

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

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

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

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

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).

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