plot.markovMSM {markovMSM} | R Documentation |
Plot for an object of class "markovMSM"
Description
Plots for the visual inspection for checking the Markov condition through the difference between the AJ and LMAJ from a starting time of the transition probabilities for each transition.
Usage
## S3 method for class 'markovMSM'
plot(x = object, to=2, quantileOrder = NULL,
axis.scale = c(-1, 1), difP = FALSE, ...)
Arguments
x |
A dataframe in long format in case of the plot with the differences between AJ and LMAJ estimators or a object of class "markovMSM" with the results of the global.test or local.test function. |
to |
The last receiving state considered for the estimation of the transition probabilities. Plot of the transition probabilities between the starting to this state is shown. |
quantileOrder |
Order of the quantil used in the formula of the AUC global test. |
axis.scale |
Limit of the y axis of the plots. |
difP |
Type of plot representing in case of x means the results of the Localtest. If difP=TRUE plot depicts the discrepancies between AJ and LMAJ estimators. If difP=FALSE plot show the AJ and the LMAJ estimates. By default difP=FALSE. |
... |
For future methods. |
Value
No value is returned.
Author(s)
Gustavo Soutinho and Luis Meira-Machado.
References
Soutinho G, Meira-Machado L (2021). Methods for checking the Markov condition in multi-state survival data. Computational Statistics.
Examples
library(markovMSM)
data(prothr)
res<-AUC.test(data=prothr, times=30, from=1, to=3, type='local',
replicas=10)
plot(res, to=3, axis.scale=c(-0.25,.25), difP=TRUE)
plot(res, to=2, axis.scale=c(0,.25), difP=FALSE)
res2<-AUC.test(data=prothr, from=1, to=3,
type='global', replicas=10, limit=0.90, quantiles=c(.05, .10, .20, .30, 0.40))
plot(res2, quantileOrder=3, 2, axis.scale=c(-0.05,.15))