search for: rms

Hi I am trying to install "rms" package but while installing it shows following error package 'survival' 2.36-2 is loaded, but >= 2.36.3 is required by 'rms' what to do? i am using linux OS I have tried by updated r-base-core but it didnt work regards GRR [[alternative HTML version deleted]]

I am cleaning up the rms package to not export functions not to be called directly by users. rms uses generic functions defined in other packages. For example there is a latex method in the Hmisc package, and rms has a latex method for objects of class "anova.rms" so there are anova.rms and latex.anova.rms func...

I have updated the rms package to extensively use NAMESPACE. I cannot get certain S3 methods to dispatch. For example I have in NAMESPACE S3method(anova, rms) S3method(latex, anova.rms) anova.rms produces an object of class "anova.rms" and there is a latex.anova.rms function in rms. But when I do latex(ano...

rms version 3.3-1 has been installed on CRAN. New features/bug fixes are below. * Added new example for anova.rms for making dot plots of partial R^2 of predictors * Defined logLik.ols (calls logLik.lm) * Fixed and cleaned up logLik.rms, AIC.rms * Fixed residuals.psm to allow other t...

Dear R-users, I am using R 2.13.0 and rms 3.3-0 , but can not reproduce figure 7.8 of the handouts *Regression Modeling Strategies* ( http://biostat.mc.vanderbilt.edu/wiki/pub/Main/RmS/course2.pdf) by the following code. Could any one help me figure out how to solve this? setwd('C:/Rharrell') require(rms) load('data/counties....

Hi, I'm using the rms package to do regression analysis using the lrm function. Retrieving odds ratios is possible using summary.rms. However, I could not find any information on how exactly the odds ratios for continuous variables are calculated. It doesn't appear to be the odds ratio at 1 unit increase, because th...

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 graph...

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 graph...

The rms ("Regression Modeling Strategies") package has undergone a massive update. The entire list of updates is at the bottom of this note. CRAN has the update for linux and will soon have it for Windows and Mac - check http://cran.r-project.org/web/packages/rms/ for availability. This rm...

...' e <- ifelse(dt <= cens,1,0) dt <- pmin(dt, cens) units(dt) <- "Year" dd <- datadist(age, sex) options(datadist='dd') S <- Surv(dt,e) library(Design) f <- cph(S ~ age, surv=TRUE,x=T,y=T) plot(f,age=NA,time=5) But the same code won't work if I used rms package: detach(package:Design) library(rms) f <- cph(S ~ age, surv=TRUE,x=T,y=T) plot(f,age=NA,time=5) Error in xy.coords(x, y, xlabel, ylabel, log) : ? 'x' and 'y' lengths differ Is there a way to plot the same graph using rms package. I like to use Frank Harrell's ne...

How can I extract the adjusted R^2 value from an ols object (using rms package)? library(rms) x <- rnorm(10) y <- x + rnorm(10) ols1 <- ols(y ~ x) Typing "ols1" displays adjusted R^2 among other things, but how can I assign it to a variable? I tried str(ols1) but couldn't see where to go from there. Thanks, Mark Seeto

The integration of plot methods for various outputs from rms packages is a great appreciated aspect of the rms package. I particularly like to use: plot(summary(model)) for my own purposes, but... for publication/presentation I need to modify details like variable names, or the number of signficant digits used in the figure annotations. Is there a sim...

Dear All, I wonder if someone can point me in the right direction here. I'm working with the rms library, R 2.9.2 under Windows XP. I'm trying to arrange two plots side by side for a colleague. mfrow or mfcol do not seem to work, however, so I am obviously missing something important. I know that there have been changes in the graphics from Design to rms, but am just not sure where to f...

...n from: (NULL) (132.183.12.87) It seems that the the standard error of prediction of the linear regression, caclulated with predict.lm is incorrect. Consider the following example where the standard error is first calculated with predict.lm, then using delta method. and finally, using the formula rms*sqrt(1+1/n+(xp-x0)^2/Sxx). Marek Ancukiewicz > n <- 10 > x <- 1:n > y <- x > y[c(2,4,6)] <- y[c(2,4,6)] + 0.1 > y[c(3,5,7)] <- y[c(3,5,7)] - 0.1 > a <- lm(y~x) > rms <- sqrt(sum(a$residuals^2)/(n-2)) > s <- covmat(a)*rms^2 > xp <- 3 > xm &...

Dear R-help, I don't understand why Gls gives me an error when trying to fit a model with AR(2) errors, while gls (from nlme) does not. For example: library(nlme) library(rms) set.seed(1) d <- data.frame(x = rnorm(50), y = rnorm(50)) gls(y ~ x, data=d, correlation = corARMA(p=2)) #This works Gls(y ~ x, data=d, correlation = corARMA(p=2)) # Gives error # Error in `coef<-.corARMA`(`*tmp*`, value = value[parMap[, i]]) : # NA/NaN/Inf in foreign function call (arg...

Version 2.2-0 of the rms package is now available. This is a somewhat major update. One major change is not downward compatible: Instead of specifying predictor=. or predictor=NA to Predict, summary, nomogram, survplot, gendata, you just specify the name of the predictor. For example, to get predictions for the defa...

Version 2.2-0 of the rms package is now available. This is a somewhat major update. One major change is not downward compatible: Instead of specifying predictor=. or predictor=NA to Predict, summary, nomogram, survplot, gendata, you just specify the name of the predictor. For example, to get predictions for the defa...

The answer to this may be obvious, but I was wondering in the rms package and the survfit(), how you can plot the censored time points as ticks. Take for example, library(survival) library(rms) foo <- data.frame(Time=c(1,2,3,4,5,6,10), Status=c(1,1,0,0,1,1,1)) answer <- survfit(Surv(foo$Time, foo$Status==1) ~1) # this shows the censored time points as t...

Hello, I am using Dr. Harrell's rms 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: library(rms) f1 <- cph(Surv(retime,dfs) ~ age+her...

...83 obs. of 5 variables: ..$ DATE : int [1:83] 2001081500 2001081512 2001081600 2001081612 2001081700 2001081712 2001081800 2001081812 2001081900 2001081912 ... ..$ logrho: num [1:83] 1.16 -1.30 -1.30 -1.30 -1.30 ... ..$ w2 : num [1:83] 1.01 1.27 1.24 1.31 1.28 ... ..$ rms : num [1:83] 5.001 0.630 0.616 0.685 0.655 ... ..$ maxi : num [1:83] 8.66 3.39 3.83 3.35 3.23 ... $ log30:'data.frame': 71 obs. of 5 variables: ..$ DATE : int [1:71] 2001081500 2001081512 2001081600 2001081612 2001081700 2001081712 2001081800 2001081812 2001081900 2...