[R] summary.rms help

Andras Farkas motyocska at yahoo.com
Wed Jan 3 20:57:10 CET 2018


Dear All,
using the example from the help of summary.rms

library(rms)
n <- 1000    # define sample size 
set.seed(17) # so can reproduce the results 
age            <- rnorm(n, 50, 10) 
blood.pressure <- rnorm(n, 120, 15) 
cholesterol    <- rnorm(n, 200, 25) 
sex            <- factor(sample(c('female','male'), n,TRUE)) 
label(age)            <- 'Age'      # label is in Hmisc 
label(cholesterol)    <- 'Total Cholesterol' 
label(blood.pressure) <- 'Systolic Blood Pressure' 
label(sex)            <- 'Sex' 
units(cholesterol)    <- 'mg/dl'   # uses units.default in Hmisc 
units(blood.pressure) <- 'mmHg' 
# Specify population model for log odds that Y=1 
L <- .4*(sex=='male') + .045*(age-50) + 
(log(cholesterol - 10)-5.2)*(-2*(sex=='female') + 2*(sex=='male')) 
# Simulate binary y to have Prob(y=1) = 1/[1+exp(-L)] 
y <- ifelse(runif(n) < plogis(L), 1, 0) 
ddist <- datadist(age, blood.pressure, cholesterol, sex) 
options(datadist='ddist') 
fit <- lrm(y ~ blood.pressure + sex * (age + rcs(cholesterol,4)))
s <- summary(fit) 
plot(s)
as you will see the plot will by default include the low and high values from the summary printed on the plot to the right of the variable name... Any thoughts on how printing these low and high values can be suppressed, ie: prevent them from being printed?

 
appreciate your help,
Andras



More information about the R-help mailing list