Comments (4)
Sorry to not have gotten to this sooner. The confidence intervals appear to be correct but the hazard ratios are exp(-beta). I'll take a look at the code to try to debug. If you have time to look at it earlier than I get to it that would be great (code for orm
and summary.rms
).
from rms.
Dear Prof. Harrell,
Thank you so much for your reply! I have noticed that the loglog link function in orm.fit.s
is given by -log(-log(x )) and not log(-log(x )) as stated in your excellent Regression Modeling Strategies book (Table 15.1, page 362). If the loglog link function is -log(-log(x )), exp(beta) will be given by logP1/logP2 while exp(-beta) will be given by logP2/logP1. If so, does exp(-beta) give the correct Hazard Ratio?
from rms.
The logic in summary.orm
was unnecessarily complex. Thanks for pointing this out. For the PH model (log-log link) I was negating the linear predictor before anti-logging, but was leaving the confidence limits in their original form. So the confidence limits were correct all along, and the point estimate was wrong. I'll have that logic removed for the next submission to CRAN. Changes will be committed to Github momentarily. Let me know if this change does not makes sense.
from rms.
Thank you for making those changes.
from rms.
Related Issues (20)
- After plot with nomogram(), the outputs are crowded together and the lines don't strehch HOT 1
- `poma` fails for response variables with >= 10 levels coded as a factor
- "Error in solve.default system is computationally singular" after restricted cubic splines with df > 3 (Cox Regression)
- the abscissa of the intersection
- datadist cannot find in function environment but can fine in global environment.
- latex works for cph, generates .tex but does not work for anova tinytex or anova problem
- .Rbuildignore contains broken symlink
- riskdist of val.prob couldn't be set to FALSE and remove risk distribution in plot
- dead link at http://hbiostat.org/R/rms/
- program logic error: variables vary in length but subset= was not given when event is not a logical HOT 1
- Issue installing on R version 4.0.0 HOT 2
- Fatal error in Gls(B > 0)
- allow to pass argument to solver (SparseM)
- Anova command does not run with lapply
- How to include time-vary covariates using tt=function from coxph in rcs code?
- print.Gls() using corSymm() returns error message.
- lrm and Rcmdr conflict HOT 2
- rexVar in survival model (cph) HOT 1
- How to estimate OR and 95%CI using Predict funtion for svyglm model?
- RCS curve in competing risk model
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.
from rms.