R: Plot method for a MarkovTest object

plot.MarkovTest {mstate} | R Documentation

## Plot method for a MarkovTest object

### Description

Plot method for an object of class 'MarkovTest'. It plots the trace of the log-rank statistics provided by `MarkovTest`.

### Usage

    ## S3 method for class 'MarkovTest'
    plot(
      x,
      y,
      what = c("states", "overall"),
      idx = NULL,
      quantiles = TRUE,
      qsup,
      states,
      xlab,
      ylab,
      main,
      ...
    )

### Arguments

`x` | Object of class 'MarkovTest'
`y` | The grid at which `MarkovTest` was calculated
`what` | Choose "states" for plotting state-specific traces, and "overall" for the overall chi-squared trace
`idx` | Vector of indices of wild bootstrap traces to plot
`quantiles` | Boolean whether or not to plot the 2.5 and 97.5 percent quantiles, default is `TRUE`
`qsup` | The index of the function in either `fn` (when plotting state-specific) or `fn2` (when plotting overall) to plot along with the traces; when missing this line is not included
`states` | Number of the qualifying state(s) to plot trace for
`xlab` | Text for x-axis label
`ylab` | Text for y-axis label
`main` | Text for title (main)
`...` | Further arguments to plot

### Value

No return value

### Author(s)

Hein Putter H.Putter@lumc.nl

### See Also

`MarkovTest`

### Examples

    ## Not run:
    # Example provided by the prothrombin data
    data("prothr")
    # Apply Markov test to grid of monthly time points over the first 7.5 years
    year <- 365.25
    month <- year / 12
    grid <- month * (1:90)
    # Markov test for transition 1 (wild bootstrap based on 100 replications)
    MT <- MarkovTest(prothr, id = "id", transition = 1,
                     grid = grid, B = 100)

    plot(MT, grid, what="states", idx=1:50, states=rownames(attr(prothr, "trans")),
         xlab="Days since randomisation", ylab="Log-rank test statistic",
         main="Transition Normal -> Low")

    plot(MT, grid,what="overall", idx=1:50,
         xlab="Days since randomisation", ylab="Chi-square test statistic",
         main="Transition Normal -> Low")

    plot(MT, grid, what="states", quantiles=FALSE) # only trace
    plot(MT, grid, what="states") # trace plus quantiles (default)
    plot(MT, grid, what="states", idx=1:10) # trace plus quantiles, plus first 10 bootstrap traces

    plot(MT, grid, what="overall", quantiles=FALSE) # only trace
    plot(MT, grid, what="overall") # trace plus quantiles (default)
    plot(MT, grid, what="overall", idx=1:10) # trace plus quantiles, plus first 10 bootstrap traces

    ## End(Not run)

[Package mstate version 0.3.2 Index]