markovMSM (version 0.1.3)

plot.markovMSM: 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 markovMSM
plot(x = object, to=2, quantileOrder = NULL, 
axis.scale = c(-1, 1), difP = FALSE, ...)

Value

No value is returned.

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.

Author

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)
# \donttest{
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))
# }