You are here:
Vanderbilt Biostatistics Wiki
>
Main Web
>
StatComp
>
RmS
>
Rms2Errata
(12 Feb 2017,
FrankHarrell
)
(raw view)
E
dit
A
ttach
---+ Errata for First Printing (Sept. 2015) of Regression Modeling Strategies, 2nd Edition For replacements for R code that no longer works due to changes in the R system or R packages, to go Rms2AltR |*Abbreviation*|*Meaning*| |Pg|Page| |L|Line| |Pr|Paragraph| |b|From Bottom (ignoring equations)| |Eq|Equation| |%Y%|Corrected in master file| |Pg _i_ | %RED%Feinberg%ENDCOLOR% should be %RED%Fienberg%ENDCOLOR% | | |Pg 150 | Changes in =ggplot2= require code changes (4) below | %Y% | |Pg 174 | Dotted blue line is missing in Figure 8.5. It should be at y=3958 | %Y% | |Pg 179 | Spline model in first large box - see (2) |%Y%*| |Pg 184 Fig. 9.2 | Caption should have log L not just L |%Y%| |Pg 240 | Table 10.4 %RED%114.51%ENDCOLOR% should be %RED%114.41%ENDCOLOR% |%Y%| |Pg 256 R code|%RED%x=TRUE, y=TRUE%ENDCOLOR% should be %RED%x=TRUE, y=TRUE, data=sex.age.response%ENDCOLOR%|%Y%| |Pg 256 R code|%RED%which.influence(f, .4)%ENDCOLOR% should be %RED%round(resid(f, 'dfbetas'), 1)%ENDCOLOR%|%Y%| |Pg 306|%RED%five completed datasets%ENDCOLOR% should be %RED%twenty completed datasets%ENDCOLOR% plus 2 other places|%Y%| |Pg 360 2 lines under Eq 15.1|%RED%X%ENDCOLOR% should be %RED%X subscript i%ENDCOLOR%| | |Pg 362 Eq. 15.7 | %RED%n%ENDCOLOR% should be %RED%k%ENDCOLOR%|%Y%| |Pg 373 Fig. 15.4 | Was run with %RED%default family%ENDCOLOR% (logistic) instead of intended %RED%family='loglog'%ENDCOLOR%; agreement increases |%Y%| |Pg 378 | %RED%family="loglog"%ENDCOLOR% can be %RED%family=loglog%ENDCOLOR% |%Y%| |Pg 474 | %RED%g 0.033 gr 1.959%ENDCOLOR% should be %RED%g=1.959 gr=7.095%ENDCOLOR% |%Y%| |Pg 508 Fig. 20.8|Omitted text (1) and index under figure|%Y%| |Pg 516 top|Code for Fig. 20.7: see Rms2AltR| | |Pg 516 bottom|Code for Fig. 20.10: see Rms2AltR| | |Pg 526| 3 paragraphs (3) omitted before %RED%rcspline.plot(lvef, d.time, event=cdeath, nk=3)%ENDCOLOR% due to omission of @ in previous chunk |%Y%| ---------------------------- * Denotes places in the master text where changes are indicated with TODO and lines need to be uncommented. <br> (1) You can see that the ratio is not stable from 160 to 190 days, but is mostly flat thereafter. \index{hazard ratio!interval-specific} Interval-specific estimates of the group 2 : group 1 log hazard ratios are shown in Table~\ref{tab:cox-ratios-group} with intervals specified to yield equal numbers of deaths. The Efron likelihood was used. [had omitted @ to close =knitr= chunk above] <br> (2) To make the model converge, and to actually show a dotted blue line, the model had to have some splines reduced from 5 knots to 4 knots: <highlight> f <- cph(S ~ rcs(sz,4) + rcs(sg,4) + rcs(log(ap),5) + rcs(sbp,4) + rcs(dbp,4) + rcs(age,3) + rcs(wt,4) + rcs(hg,4) + ekg + pf + bm + hx, tol=1e-14, data=imputed) </highlight> New text interpreting the results: The model allowing all continuous predictors to be nonlinear is better than the linear full model but not b etter than the 5 component model. <br> (3) The standard errors of log survival probabilities computed by \fu{cph} and the confidence intervals derived by the associated functions \fu{survest} and \fu{survplot} are only proper when the model contains only stratification factors or when predictions are made at the means of all covariables, unless you store in the fit object the data used in fitting the model using the options \Co{x=TRUE, y=TRUE}. The \fu{survfit.cph} and \fu{survest} functions are modifications of Therneau's \fu{survfit.coxph} function that can obtain predicted survival probabilities and confidence intervals for any setting of predictors. \fu{survfit.cph} and \fu{survest} are called as needed when you request plots. <br> The \fu{survest} function is used to obtain predicted survival curves or survival at a given vector of times, for a desired set of predictor values and a fit produced by \fu{cph}. Confidence limits may be obtained also, as can the standard error of the log--log survival. <br> The \fu{survplot} function is similar in some ways to the \fu{Predict} and \Co{plot}, but it plots follow-up time on the $x$-axis. It plots step-function survival curves with confidence bands or confidence bars for Cox models. <br> Figure~\ref{fig:cox-ef-spline} was produced using ... (4) <highlight> require(data.table) cdystonia <- data.table(cdystonia) cdys <- cdystonia[, j=as.list(quantile(twstrs, (1 : 3)/4)), by = list(treat, week)] cdys <- upData(cdys, rename=c('25%'='Q1', '50%'='Q2', '75%'='Q3')) ggplot(cdys, aes(x=week, y=Q2)) + xl + yl + ylim(0, 70) + geom_line() + facet_wrap(~ treat, nrow=2) + geom_ribbon(aes(ymin=Q1, ymax=Q3), alpha=0.2) </highlight>
E
dit
|
A
ttach
|
P
rint version
|
H
istory
: r14
<
r13
<
r12
<
r11
|
B
acklinks
|
V
iew topic
|
Edit
w
iki text
|
M
ore topic actions
Topic revision: r14 - 12 Feb 2017,
FrankHarrell
Main
Department Home Page
Biostatistics Graduate Program
Vanderbilt University Medical Center
Main Web
Main Web Home
Search
Recent Changes
Changes
Topic list
Biostatistics Webs
Archive
Main
Sandbox
System
Register
|
Log In
Copyright © 2013-2022 by the contributing authors. All material on this collaboration platform is the property of the contributing authors.
Ideas, requests, problems regarding Vanderbilt Biostatistics Wiki?
Send feedback