cum.residuals | R Documentation |

Computes cumulative residuals and approximative p-values based on resampling techniques.

cum.residuals( object, data = parent.frame(), modelmatrix = 0, cum.resid = 1, n.sim = 500, weighted.test = 0, max.point.func = 50, weights = NULL )

`object` |
an object of class 'aalen', 'timecox', 'cox.aalen' where the residuals are returned ('residuals=1') |

`data` |
data frame based on which residuals are computed. |

`modelmatrix` |
specifies a grouping of the data that is used for cumulating residuals. Must have same size as data and be ordered in the same way. |

`cum.resid` |
to compute residuals versus each of the continuous covariates in the model. |

`n.sim` |
number of simulations in resampling. |

`weighted.test` |
to compute a variance weighted version of the test-processes used for testing constant effects of covariates. |

`max.point.func` |
limits the amount of computations, only considers a max of 50 points on the covariate scales. |

`weights` |
weights for sum of martingale residuals, now for cum.resid=1. |

returns an object of type "cum.residuals" with the following arguments:

`cum` |
cumulative residuals versus time for the groups specified by modelmatrix. |

`var.cum` |
the martingale based pointwise variance estimates. |

`robvar.cum` |
robust pointwise variances estimates of cumulatives. |

`obs.testBeq0` |
observed absolute value of supremum of cumulative components scaled with the variance. |

`pval.testBeq0` |
p-value covariate effects based on supremum test. |

`sim.testBeq0` |
resampled supremum value. |

`conf.band` |
resampling based constant to construct robust 95% uniform confidence bands for cumulative residuals. |

`obs.test` |
absolute value of supremum of observed test-process. |

`pval.test` |
p-value for supremum test statistic. |

`sim.test` |
resampled absolute value of supremum cumulative residuals. |

`proc.cumz` |
observed cumulative residuals versus all continuous covariates of model. |

`sim.test.proccumz` |
list of 50 random realizations of test-processes under model for all continuous covariates. |

Thomas Scheike

Martinussen and Scheike, Dynamic Regression Models for Survival Data, Springer (2006).

data(sTRACE) # Fits Aalen model and returns residuals fit<-aalen(Surv(time,status==9)~age+sex+diabetes+chf+vf, data=sTRACE,max.time=7,n.sim=0,residuals=1) # constructs and simulates cumulative residuals versus age groups fit.mg<-cum.residuals(fit,data=sTRACE,n.sim=100, modelmatrix=model.matrix(~-1+factor(cut(age,4)),sTRACE)) par(mfrow=c(1,4)) # cumulative residuals with confidence intervals plot(fit.mg); # cumulative residuals versus processes under model plot(fit.mg,score=1); summary(fit.mg) # cumulative residuals vs. covariates Lin, Wei, Ying style fit.mg<-cum.residuals(fit,data=sTRACE,cum.resid=1,n.sim=100) par(mfrow=c(2,4)) plot(fit.mg,score=2) summary(fit.mg)

Embedding an R snippet on your website

Add the following code to your website.

For more information on customizing the embed code, read Embedding Snippets.