search for: rrms

A new version of rms is now available on CRAN for Linux/UNIX. I expect Mac and Windows versions to be available in a day or so. This version works with and requires the newest version of Therneau's survival package. More information is at http://biostat.mc.vanderbilt.edu/Rrms Changes in version 3.2-0 (2011-02-14) * Changed to be compatible with survival 2.36-3 which is now required * Added logLik.rms and AIC.rms functions to be compatible with standard R * Fixed oos.loglik.Glm * Fixed bootcov related to nfit='Glm' * Fixed (probably) old b...

...* lrm.fit, lrm: linear predictor stored in fit object now uses first intercept and not middle one (NOT DOWNWARD COMPATIBLE but makes predict work when using stored linear.predictors) * Fixed argument consistency with validate methods More information is at http://biostat.mc.vanderbilt.edu/Rrms -- Frank E Harrell Jr Professor and Chairman School of Medicine Department of Biostatistics Vanderbilt University _______________________________________________ R-packages mailing list R-packages at r-project.org https://stat.ethz.ch/mailman/listinfo/r-packages

...* lrm.fit, lrm: linear predictor stored in fit object now uses first intercept and not middle one (NOT DOWNWARD COMPATIBLE but makes predict work when using stored linear.predictors) * Fixed argument consistency with validate methods More information is at http://biostat.mc.vanderbilt.edu/Rrms -- Frank E Harrell Jr Professor and Chairman School of Medicine Department of Biostatistics Vanderbilt University _______________________________________________ R-packages mailing list R-packages at r-project.org https://stat.ethz.ch/mailman/listinfo/r-packages

A new version of rms is now available on CRAN for Linux/UNIX. I expect Mac and Windows versions to be available in a day or so. This version works with and requires the newest version of Therneau's survival package. More information is at http://biostat.mc.vanderbilt.edu/Rrms Changes in version 3.2-0 (2011-02-14) * Changed to be compatible with survival 2.36-3 which is now required * Added logLik.rms and AIC.rms functions to be compatible with standard R * Fixed oos.loglik.Glm * Fixed bootcov related to nfit='Glm' * Fixed (probably) old b...

The code library(Design) f <- lrm(y~x1+x2+x1*x2, data=data) plot(f) produces a plot of log odds vs x2 with 0.95 confidence intervals. How do I get a plot of odds ratios vs x2 instead? Thanks -- View this message in context: http://r.789695.n4.nabble.com/Odds-ratios-from-lrm-plot-tp4033340p4033340.html Sent from the R help mailing list archive at Nabble.com.

...as to analyze the data as a repeated-measures design where I consider the four rankings as four observations for the same person (one response per option). This would be an easy case of ordinal regression (e.g. proportional odds model) using the excellent package http://biostat.mc.vanderbilt.edu/Rrms of Frank Harrel (thank you Frank). However, these are not four independent observations per person, because if a respondent assigned the values to the first three objects, then the rank of the fourth object is pre-determined. Does anyone know if there is a package that can fit a more appropriate...

Hello, I would like to use the aregImpute and fit.mult.impute to impute missing values for my dataset and then conduct logistic regression analyses on the data, taking into account that we imputed values. I have no problems imputing the values using aregImpute, but I am getting an error at the fit.mult.impute stage. Here is some sample code (I actually have more observations and variables to

...ssion, accelerated failure time models, ordinary linear models, the Buckley-James model, generalized least squares for serially or spatially correlated observations, generalized linear models, and quantile regression. More details are at http://biostat.mc.vanderbilt.edu/Rrms ----- Frank Harrell Department of Biostatistics, Vanderbilt University -- View this message in context: http://r.789695.n4.nabble.com/New-version-of-rms-package-on-CRAN-tp3570655p3570655.html Sent from the R help mailing list archive at Nabble.com.

This is to announce a new package rms on CRAN. rms goes along with my book Regression Modeling Strategies. The home page for rms is http://biostat.mc.vanderbilt.edu/rms, or go directly to http://biostat.mc.vanderbilt.edu/Rrms for information just about the software. rms is a re-write of the Design package that has improved graphics and that duplicates very little code in the survival package. In particular, rms does not use low-level C language interfaces to other packages and will be easier to maintain. rms also...

This is to announce a new package rms on CRAN. rms goes along with my book Regression Modeling Strategies. The home page for rms is http://biostat.mc.vanderbilt.edu/rms, or go directly to http://biostat.mc.vanderbilt.edu/Rrms for information just about the software. rms is a re-write of the Design package that has improved graphics and that duplicates very little code in the survival package. In particular, rms does not use low-level C language interfaces to other packages and will be easier to maintain. rms also...

Hello, I am using Dr. Harrell's design package to make a nomogram. I was able to make a beautiful one. However, I want to change 5-year survival probability to 5-year failure probability. I couldn?t get hazard rate from Hazard(f1) because I used cph for the model. Here is my code: f1 <- cph(Surv(retime,dfs) ~ age+her2+t_stage+n_stage+er+grade+cytcyt+Cyt_PCDK2 , data=data11, surv=T,

...- html methods for model fit summaries especially when using Rmarkdown html notebooks - instead of running print(fit, latex=TRUE) use options(prType='latex') print(fit) Or: options(prType='html') A complete list of changes is below. See http://biostat.mc.vanderbilt.edu/Rrms for an html notebook showing examples of new features. Changes in version 5.0-0 (2016-10-31) * plot.summary.rms: implemented plotly interactive plots if options(grType='plotly') in effect * plot.anova.rms: implemented plotly interactive plots if options(grType='plotly') in effec...

...- html methods for model fit summaries especially when using Rmarkdown html notebooks - instead of running print(fit, latex=TRUE) use options(prType='latex') print(fit) Or: options(prType='html') A complete list of changes is below. See http://biostat.mc.vanderbilt.edu/Rrms for an html notebook showing examples of new features. Changes in version 5.0-0 (2016-10-31) * plot.summary.rms: implemented plotly interactive plots if options(grType='plotly') in effect * plot.anova.rms: implemented plotly interactive plots if options(grType='plotly') in effec...

Dear friends Hello I'm a PhD student in Edinburgh University and interested in survival analysis. By chance, I found design library in Splus software. I have some questions about it and it's highly appreciated if you can help me. In its help (Design library help) I found the topics(not the commands), which are available in the main software. for instance, survival analysis is

Dear all, i have re-installed Debian on my computer and for some reason I cannot install some packages anymore using the install.package command from the R console. Some packages I can, and some I cannot (Design for instance) this is what I get: > install.packages('Design') Installing package(s) into ‘/usr/local/lib/R/site-library’ (as ‘lib’ is unspecified) --- Please select a CRAN

Dear all, I am using Windows and R 2.9.2 for my analyses. I have a large dataset and I am particularly interested in looking at time to an event for a continuous variable. I would like to produce a plot of log(relative risk) or relative risk (also known as hazard ratio) against the continuous variable. I have spent a long time looking for advice on how to do this but my search has proved

Dear all, I try to show a subset of coefficients in my presentation. It seems that a "standard" table is not a good way to go. I found figure 9 (page 9) in this file ( http://www.destatis.de/jetspeed/portal/cms/Sites/destatis/Internet/DE/Content/Wissenschaftsforum/Kolloquien/VisualisierungModellierung__Beitrag,property=file.pdf ) looks pretty good. I wonder if there is any function for

Dear R-help, I am having a problem with the interpretation of result from validate.cph in the Design package. My purpose is to fit a cox model and validate the Somer's Dxy. I used the hypothetical data given in the help manual with modification to the cox model fit. My research problem is very similar to this example. This is the model without stratification: > library(Design) > f1

Hi all I have run into a case where I don't understand why predict.lrm and predict.glm don't yield the same results. My data look like this: set.seed(1) library(Design); ilogit <- function(x) { 1/(1+exp(-x)) } ORDER <- factor(sample(c("mc-sc", "sc-mc"), 403, TRUE)) CONJ <- factor(sample(c("als", "bevor", "nachdem",