# NOT RUN {
data(abide_9002_23)
res <- graph.flm(nsim = 19, # Increase nsim for serious analysis!
formula.full = Y ~ Group + Sex + Age,
formula.reduced = Y ~ Sex + Age,
curve_sets = list(Y = subset(abide_9002_23[['curve_set']], 1:50)),
factors = abide_9002_23[['factors']][1:50,],
contrasts = FALSE,
GET.args = list(type = "area"))
plot(res)
plot(res, what=c("obs", "hi"))
plot(res, what=c("hi", "lo"), fixedscales=1)
plot(res, what=c("obs", "lo", "hi"), fixedscales=FALSE)
if(requireNamespace("gridExtra", quietly=TRUE) && require("ggplot2", quietly=TRUE)) {
# Edit style of "fixedscales = 2" plots
plot(res, what=c("obs", "hi")) + theme_minimal()
# Edit style (e.g. theme) of "fixedscales = 1 or 0" plots
gs <- lapply(res, plot, what=c("obs", "hi"), main="")
gridExtra::grid.arrange(grobs=gs, ncol=1, top="My main")
gs <- outer(res, c("obs", "hi"), FUN=Vectorize(function(res, what)
list(plot(res, what=what, main="") + theme(axis.ticks=element_blank(),
axis.text=element_blank(), axis.title=element_blank()))))
gridExtra::grid.arrange(grobs=t(gs))
}
# }
Run the code above in your browser using DataLab