search for: age75

...ndouts *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.sav') older <- counties$age6574 + counties$age75 label(older) <- '% age >= 65, 1990' pdensity <- logb(counties$pop.density+1, 10) label(pdensity) <- 'log 10 of 1992 pop per 1990 miles^2' counties <- cbind(counties, older, pdensity) # add 2 vars. not in data frame dd <- datadist(counties) options(datadist='dd...

...te a case study described by Prof. Harrell in Chapter 7 of his Regression Modeling Strategies book, but failed on using anova.Design to reproduce his table 7.1, Following is the code: rm(list=ls()) library(Hmisc) library(Design) getHdata(counties) counties$older <- counties$age6574 + counties$age75 label(counties$older) <- '% age >= 65, 1990' counties$pdensity <- log10(counties$pop.density+1) label(counties$pdensity) <- 'log 10 of 1992 pop per 1990 miles^2' dd <- datadist(counties) options(datadist='dd') f <- ols(democrat ~ rcs(pdensity,4) + rcs(pop...

Dear R-users: I am trying to use the following code to reproduce the figures on page 140 of Prof. Frank Harrell's book 'Regression Modeling Strategies': rm(list=ls()) options(width=128) library(Hmisc) library(Design) getHdata(counties) counties$older <- counties$age6574 + counties$age75 label(counties$older) <- '% age >= 65, 1990' counties$pdensity <- log10(counties$pop.density+1) label(counties$pdensity) <- 'log 10 of 1992 pop per 1990 miles^2' dd <- datadist(counties) options(datadist='dd') f <- ols(democrat ~ rcs(pdensity,4) + rcs(pop...

You report an error message: > plot(zph.revasFit[1]) Error in plot.window(xlim, ylim, log, asp, ...) : need finite 'ylim' values I have never seen this error before, and I cannot guess what causes it. You need to provide more information, and likely a small data set that produces the problem. Perhaps you have an x variable that is a constant? Terry Therneau