-
-
Notifications
You must be signed in to change notification settings - Fork 556
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
weibull_aft.plot_survival_function: object has no attribute #1562
Comments
What is the survival function of an AFT model? The AFT model is conditional (i.e requires covariates). Maybe you want |
I see an example here on page 11-12 |
That's the |
Ah! thank you.
How can I get the baseline curve alone?
At this time, my only covariate is a continuous one.
I will use the partial feature after augmenting the dataset.
…On Thu, Sep 28, 2023 at 5:26 PM Cameron Davidson-Pilon < ***@***.***> wrote:
That's the WeibullFitter, not WeibullAFTFitter
—
Reply to this email directly, view it on GitHub
<#1562 (comment)>,
or unsubscribe
<https://github.com/notifications/unsubscribe-auth/ANQJMZPRCFFESSINWQOOPODX4YIT7ANCNFSM6AAAAAA5LPF4UI>
.
You are receiving this because you authored the thread.Message ID:
***@***.***>
|
There's not really a baseline curve for AFT models (terminology isn't used), but I think you just want to set the covariate to 0 in |
Never mind, I added a relevant constant column T.
|
The first three questions can be answered by checking out the docs: https://lifelines.readthedocs.io/en/latest/fitters/regression/WeibullAFTFitter.html?highlight=plot_partial_effects_on_outcome#lifelines.fitters.weibull_aft_fitter.WeibullAFTFitter.plot_partial_effects_on_outcome I don't quite understand your 4th question, however |
That helps. Thanks Cam! I guess I can get the cumulative distribution function (CDF) as 1-survival from the plot collection of partials. Is there a confidence interval band fill between for the survival plot? If not, how can I get it? |
Unfortunately, not |
Under what conditions is the cumulative_hazard < hazard?
My plots have t = 0 to 6.
The Y for hazard rate goes up to 750 at t=6.
THe Y for cumulative hazard goes only up to 450 at t=6
The survival plot looks as expected..
…On Sat, Sep 30, 2023 at 3:15 PM Cameron Davidson-Pilon < ***@***.***> wrote:
Is there a confidence interval band fill between for the survival plot? If
not, how can I get it?
Unfortunately, not
—
Reply to this email directly, view it on GitHub
<#1562 (comment)>,
or unsubscribe
<https://github.com/notifications/unsubscribe-auth/ANQJMZNT4D363SP4OXZKTPDX5CKWVANCNFSM6AAAAAA5LPF4UI>
.
You are receiving this because you authored the thread.Message ID:
***@***.***>
|
cumulative hazard (t) = int_0^t hazard(s) ds so it's possible for cumulative_hazard < hazard. Think about how a short spike in a function might affect its integral. |
Thanks Cam.
If my input data itself is in log10 and the shape is:
coef exp(coef)
rho_ Intercept 2 9
Is it safe to say, my antilog(rho_) for interpretation is 10^(2) and not
10^(9)?
…On Tue, Oct 3, 2023 at 2:19 PM Cameron Davidson-Pilon < ***@***.***> wrote:
cumulative hazard (t) = int_0^t hazard(s) ds
so it's possible for cumulative_hazard < hazard. Think about how a short
spike in a function might affect its integral.
—
Reply to this email directly, view it on GitHub
<#1562 (comment)>,
or unsubscribe
<https://github.com/notifications/unsubscribe-auth/ANQJMZOHNEKIIZ4BQKY3MHDX5R6NDAVCNFSM6AAAAAA5LPF4UKVHI2DSMVQWIX3LMV43OSLTON2WKQ3PNVWWK3TUHMYTONBVG42DEMBXGQ>
.
You are receiving this because you authored the thread.Message ID:
***@***.***>
|
Sure, yea, but you should expect very small variation in your log10 variable as a consequence. |
Is the plot_partial_effects_on_outcome same as a regular weibull(1.75,2.24) plot conditional on T? (By the way, there is some literature on confidence intervals for Cox PH survival with covariates. Not sure if those are conditional on a covariate. |
For WeibullAFT, why does not predict_survival_function have partial
outcomes?
…On Thu, Oct 5, 2023 at 3:51 PM Cameron Davidson-Pilon < ***@***.***> wrote:
Sure, yea, but you should expect very small variation in your log10
variable as a consequence.
—
Reply to this email directly, view it on GitHub
<#1562 (comment)>,
or unsubscribe
<https://github.com/notifications/unsubscribe-auth/ANQJMZOGANFWD5USMNRK67DX542WXAVCNFSM6AAAAAA5LPF4UKVHI2DSMVQWIX3LMV43OSLTON2WKQ3PNVWWK3TUHMYTONBZG43DCOBRGY>
.
You are receiving this because you authored the thread.Message ID:
***@***.***>
|
Is this a known issue?
Is there a workaround?
Traceback (most recent call last):
File "lifelines_example.py", line 31, in
weibull_aft.plot_survival_function(ax=axes[0][0])
AttributeError: 'WeibullAFTFitter' object has no attribute 'plot_survival_function'
The text was updated successfully, but these errors were encountered: