Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Support for estimation of expected survival time and similar measures #136

Open
adibender opened this issue Feb 6, 2020 · 4 comments
Open
Assignees
Labels
feature feature requests priority: low

Comments

@adibender
Copy link
Owner

see also: https://journal.r-project.org/archive/2019/RJ-2019-042/index.html

@adibender adibender self-assigned this Feb 6, 2020
@fabian-s
Copy link
Collaborator

fabian-s commented Feb 6, 2020

see also: https://journal.r-project.org/archive/2019/RJ-2019-042/index.html

o wow.... these trivial tiny weird things survival people get published.... ;)

ok, but, (numerically) integrating S(t) as mentioned in the comments to PR #131 would seem the way to go since this is fairly easy to do based on existing functionality in pammtools? that other weird stuff (put a GAM on top of our rankings and times etc) only works for PH models anyway, right?
these coxed people simply integrate to max(event times, censoring times). That seems simplistic/dangerous to me.

@adibender
Copy link
Owner Author

adibender commented Feb 6, 2020

Yeah, i was thinking that their stuff is not necessary for us, but wanted to link for reference.

these coxed people simply integrate to max(event times, censoring times). That seems simplistic/dangerous to me

true, but we also only have S(t) estimates until the last interval (which using default cuts will end at max(event times, censoring times).

Also, once we have the estimates for E(t|x) we also need an estimate for SE (potentially for all 3 SE types) ...

@fabian-s
Copy link
Collaborator

fabian-s commented Feb 6, 2020

true, but we also only have S(t) estimates until the last interval (which using default cuts will end at max(event times, censoring times).

true, too, all i'm saying it really does not make sense to set the upper limit of the integration to a Tmax for which S(Tmax) is not already almost zero, or you obviously get a large downward bias.

Also, once we have the estimates for E(t|x) we also need an estimate for SE (potentially for all 3 SE types) ...

seems exceedingly easy for sim...? for direct, we'd just integrate the upper/lower pointwise CI of S(t)? for delta, we need to run Maple I think 🤕

adibender added a commit that referenced this issue Feb 8, 2020
This was referenced Feb 8, 2020
@adibender adibender removed the addition label Feb 8, 2020
@adibender
Copy link
Owner Author

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
feature feature requests priority: low
Projects
None yet
Development

No branches or pull requests

2 participants