wenjie2wang / intsurv Goto Github PK
View Code? Open in Web Editor NEWIntegrative Survival Models
Home Page: https://wwenjie.org/intsurv/
License: GNU General Public License v3.0
intsurv's Introduction
intsurv's People
Contributors
Stargazers
Watchers
intsurv's Issues
Different results from smcure for the lung data?
Received the following issue via emails.
Issue
Fitting a Cox cure rate model to the survival::lung dataset by using intsurv::cox_cure() and smcure::smcure() produces different coefficient estimates.
library(survival)
library(intsurv)
library(smcure)
data(lung)cox_cure(~ age + sex, ~ age + sex, time, status - 1, data = lung)## Call:
## cox_cure(surv_formula = ~age + sex, cure_formula = ~age + sex,
## time = time, event = status - 1, data = lung)
##
## Coefficient estimates from the survival part:
##
## coef exp(coef) se(coef) z Pr(>|z|)
## age 0.018277 1.018445 NA NA NA
## sex -0.519556 0.594785 NA NA NA
##
## Coefficient estimates from the cure rate part:
##
## coef exp(coef) se(coef) z Pr(>|z|)
## (Intercept) 2.686813 14.684801 NA NA NA
## age 0.016017 1.016146 NA NA NA
## sex -0.397289 0.672140 NA NA NA
smcure(Surv(time, status - 1) ~ age + sex, cureform = ~ age + sex,
data = lung, model = "ph", Var = FALSE)## Program is running..be patient... done.
## Call:
## smcure(formula = Surv(time, status - 1) ~ age + sex, cureform = ~age +
## sex, data = lung, model = "ph", Var = FALSE)
##
## Cure probability model:
## Estimate
## (Intercept) 8.20142875
## age -0.05227707
## sex -1.29707091
##
##
## Failure time distribution model:
## Estimate
## age 0.02754377
## sex -0.31575014
While the coefficient estimates are close for the bmt dataset in smcure.
data(bmt)
cox_cure(~ TRT, ~ TRT, Time, Status, data = bmt)## Call:
## cox_cure(surv_formula = ~TRT, cure_formula = ~TRT, time = Time,
## event = Status, data = bmt)
##
## Coefficient estimates from the survival part:
##
## coef exp(coef) se(coef) z Pr(>|z|)
## TRT 0.63639 1.88965 NA NA NA
##
## Coefficient estimates from the cure rate part:
##
## coef exp(coef) se(coef) z Pr(>|z|)
## (Intercept) 1.0567 2.8770 NA NA NA
## TRT 0.3578 1.4302 NA NA NA
smcure(Surv(Time, Status) ~ TRT, cureform = ~ TRT,
data = bmt, model = "ph", Var = FALSE)## Program is running..be patient... done.
## Call:
## smcure(formula = Surv(Time, Status) ~ TRT, cureform = ~TRT, data = bmt,
## model = "ph", Var = FALSE)
##
## Cure probability model:
## Estimate
## (Intercept) 1.0565750
## TRT 0.3579095
##
##
## Failure time distribution model:
## Estimate
## TRT 0.6363645
Recommend Projects
-
React
A declarative, efficient, and flexible JavaScript library for building user interfaces.
-
Vue.js
🖖 Vue.js is a progressive, incrementally-adoptable JavaScript framework for building UI on the web.
-
Typescript
TypeScript is a superset of JavaScript that compiles to clean JavaScript output.
-
TensorFlow
An Open Source Machine Learning Framework for Everyone
-
Django
The Web framework for perfectionists with deadlines.
-
Laravel
A PHP framework for web artisans
-
D3
Bring data to life with SVG, Canvas and HTML. 📊📈🎉
-
Recommend Topics
-
javascript
JavaScript (JS) is a lightweight interpreted programming language with first-class functions.
-
web
Some thing interesting about web. New door for the world.
-
server
A server is a program made to process requests and deliver data to clients.
-
Machine learning
Machine learning is a way of modeling and interpreting data that allows a piece of software to respond intelligently.
-
Visualization
Some thing interesting about visualization, use data art
-
Game
Some thing interesting about game, make everyone happy.
Recommend Org
-
Facebook
We are working to build community through open source technology. NB: members must have two-factor auth.
-
Microsoft
Open source projects and samples from Microsoft.
-
Google
Google ❤️ Open Source for everyone.
-
Alibaba
Alibaba Open Source for everyone
-
D3
Data-Driven Documents codes.
-
Tencent
China tencent open source team.