
Plot results of a permutation test performed with ptest
# S3 method for clusbootptest
plot(x, pcol = "red", pty = 1, mfrow = c(1, 1), ...)
object of class clusbootptest
color of vertical line indicating the observed Welch t test statistic
type of vertical line indicating the observed Welch t test statistic
vector of length 2 indicating the numbers of rows and columns in which the histograms will be drawn on the device.
other arguments to be passed into the hist
function.
# NOT RUN {
medication <- medication[medication$time %% 1 == 0,]
set.seed(1)
permtest.1 <- ptest(data = meds, outcome = pos, within = time, between = treat,
at.within = c(0,2,4,6), at.between = c(0,1), pn = 2000)
plot(permtest.1, pcol = "red", pty=2, mfrow = c(2,2), breaks="FD")
# }
Run the code above in your browser using DataLab