## Not run:
# #The examples uses the OBrienKaiser dataset from car and needs reshape.
# # This extends the examples from raw.means.plot
# require(reshape)
# require(car)
# data(OBrienKaiser)
# OBKnew <- cbind(factor(1:nrow(OBrienKaiser)), OBrienKaiser)
# colnames(OBKnew)[1] <- "id"
# OBK.long <- melt(OBKnew)
# OBK.long[, c("measurement", "time")] <-
# t(vapply(strsplit(as.character(OBK.long$variable), "\\."), "[", c("", "")))
#
# # For this example the position at each x-axis are within-subject comparisons!
# raw.means.plot2(OBK.long, "id", "measurement", "gender", "value")
# add.ps(OBK.long, "id", "measurement", "gender", "value", paired = TRUE)
# #reference is "fup"
#
# raw.means.plot2(OBK.long, "id", "measurement", "gender", "value")
# add.ps(OBK.long, "id", "measurement", "gender", "value", ref.offset = 2,
# paired = TRUE) #reference is "post"
#
# # Use R's standard (i.e., Welch test)
# raw.means.plot2(OBK.long, "id", "treatment", "gender", "value")
# add.ps(OBK.long, "id", "treatment", "gender", "value",
# prefixes = c("p(control vs. A)", "p(control vs. B)"))
#
# # Use standard t-test
# raw.means.plot2(OBK.long, "id", "treatment", "gender", "value")
# add.ps(OBK.long, "id", "treatment", "gender", "value", var.equal = TRUE,
# prefixes = c("p(control vs. A)", "p(control vs. B)"))
#
# ## End(Not run)
Run the code above in your browser using DataLab