MASS

Running examples in ‘MASS-Ex.R’ failed
The error most likely occurred in:

> base::assign(".ptime", proc.time(), pos = "CheckExEnv")
> ### Name: motors
> ### Title: Accelerated Life Testing of Motorettes
> ### Aliases: motors
> ### Keywords: datasets
> 
> ### ** Examples
> 
> library(survival)
> plot(survfit(Surv(time, cens) ~ factor(temp), motors), conf.int = FALSE)
> # fit Weibull model
> motor.wei <- survreg(Surv(time, cens) ~ temp, motors)
> ## IGNORE_RDIFF_BEGIN
> summary(motor.wei)

Call:
survreg(formula = Surv(time, cens) ~ temp, data = motors)
               Value Std. Error     z       p
(Intercept) 16.31852    0.62296  26.2 < 2e-16
temp        -0.04531    0.00319 -14.2 < 2e-16
Log(scale)  -1.09564    0.21480  -5.1 3.4e-07

Scale= 0.334 

Weibull distribution
Loglik(model)= -147.4   Loglik(intercept only)= -169.5
	Chisq= 44.32 on 1 degrees of freedom, p= 2.8e-11 
Number of Newton-Raphson Iterations: 8 
n= 40 

> ## IGNORE_RDIFF_END
> # and predict at 130C
> unlist(predict(motor.wei, data.frame(temp=130), se.fit = TRUE))
   fit.1 se.fit.1 
33813.06  7506.36 
> 
> motor.cox <- coxph(Surv(time, cens) ~ temp, motors)
> summary(motor.cox)
Call:
coxph(formula = Surv(time, cens) ~ temp, data = motors)

  n= 40, number of events= 17 

        coef exp(coef) se(coef)     z Pr(>|z|)    
temp 0.09185   1.09620  0.02736 3.358 0.000786 ***
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

     exp(coef) exp(-coef) lower .95 upper .95
temp     1.096     0.9122     1.039     1.157

Concordance= 0.84  (se = 0.035 )
Likelihood ratio test= 25.56  on 1 df,   p=4e-07
Wald test            = 11.27  on 1 df,   p=8e-04
Score (logrank) test = 22.73  on 1 df,   p=2e-06

> # predict at temperature 200
> plot(survfit(motor.cox, newdata = data.frame(temp=200),
+      conf.type = "log-log"))
Error in a$specials[[i]] : subscript out of bounds
Calls: plot ... survfit.coxph -> model.matrix.coxph -> delete.response -> off
Execution halted