inst/doc/paper.R

In micompr: Multivariate Independent Comparison of Observations

## ----pphpc1, eval = FALSE-----------------------------------------------------
# # Load library
# library(micompr)
# 
# # Output names
# outputs <- c("$P^s$", "$P^w$", "$P^c$", "$\\mean{E}^s$",
#              "$\\overline{E}^w$", "$\\overline{C}$",
#              "$\\widetilde{A}$")
# 
# # Outputs from the NetLogo implementation
# dir_nl_ok <- paste0(dir_data, "nl_ok")
# # Outputs from the Java implementation, first configuration
# dir_jex_ok <- paste0(dir_data, "j_ex_ok")
# # Outputs from the Java implementation, second configuration
# dir_jex_noshuff <- paste0(dir_data, "j_ex_noshuff")
# # Outputs from the Java implementation, third configuration
# dir_jex_diff <- paste0(dir_data, "j_ex_diff")
# 
# # Files for model size 400, parameter set 1
# filez <- glob2rx("stats400v1*.txt")
# 
# # Perform the three comparison cases
# mic <- micomp(outputs,
#               ve_npcs = 0.75,
#               list(list(name = "I",
#                         folders = c(dir_nl_ok, dir_jex_ok),
#                         files = c(filez, filez),
#                         lvls = c("NLOK", "JEXOK")),
#                    list(name = "II",
#                         folders = c(dir_nl_ok, dir_jex_noshuff),
#                         files = c(filez, filez),
#                         lvls = c("NLOK", "JEXNS")),
#                    list(name = "III",
#                         folders = c(dir_nl_ok, dir_jex_diff),
#                         files = c(filez, filez),
#                         lvls = c("NLOK","JEXDIF"))),
#               concat = T)

## ----pphpc2, eval = FALSE-----------------------------------------------------
# toLatex(mic,
#         booktabs = T,
#         data_show = c("npcs-1", "mnvp-1", "parp-1", "scoreplot"),
#         data_labels = c("$\\#$PCs", "MNV", "$t$-test", "PCS"),
#         col_width = T,
#         pvalf_params = list(minval = 1e-8, na_str = "*"),
#         label = "tab:pphpc",
#         caption = paste("Comparison of a NetLogo implementation of",
#                         "the PPHPC model against three configurations",
#                         "of a parallel Java implementation."))

## ----sunspot1, results = 'hide', warning = FALSE------------------------------
# Load library
library(micompr)

# Months in the 1749-1859 interval (110 years)
# Months in the 1902-2012 interval (110 years)
m <- sunspot.month[c(1:1320, 1837:3156)]
m <- matrix(m, nrow = 20)

# Factor vector, two levels:
# a) ten 11-year cycles from 1749 to 1859
# b) ten 11-year cycles from 1902 to 2012
groups <- factor(c(rep("A", 10), rep("B", 10)))

# Compare the two groups, use 9 PCs for MANOVA
cmp <- cmpoutput("SunSpots", 9, m, groups)

## ----sunspot2, results = 'markup', warning = FALSE----------------------------
cmp

## ----sunspot3, results = 'markup', warning = FALSE----------------------------
assumptions(cmp)

## ----sunspot4, fig.show = 'asis', fig.env = 'figure', fig.cap = 'Plots produced by sunspots example.'----
plot(cmp)

## ----derma1, eval = FALSE-----------------------------------------------------
# # Load libraries
# library(bmp)
# library(micompr)
# 
# # Image definitions
# imgs <- dir(imgfolder)
# nimgs <- length(imgs)
# npixels <- 760 * 570
# 
# # Specify image groups (Common nevi, atypical nevi,
# # melanomas).
# f <- read.table(grpsfile, row.names = 1)
# grps <- f[order(row.names(f)), ]
# 
# # Read images from disk
# # Use different color channels as outputs, and also
# # use a concatenated output
# rimgs <- matrix(nrow = nimgs, ncol = npixels)
# gimgs <- matrix(nrow = nimgs, ncol = npixels)
# bimgs <- matrix(nrow = nimgs, ncol = npixels)
# rgbimgs <- matrix(nrow = nimgs, ncol = npixels * 3)
# 
# for (i in 1:nimgs) {
# 
#   cimg <- read.bmp(paste0(imgfolder, imgs[i]))
#   rimgs[i, ] <- c(cimg[ , , 1])
#   gimgs[i, ] <- c(cimg[ , , 2])
#   bimgs[i, ] <- c(cimg[ , , 3])
#   rgbimgs[i, ] <- c(cimg[ , , 1], cimg[ , , 2], cimg[ , , 3])
# 
# }
# 
# # Perform multivariate independent comparison of images
# mic <-
#   micomp(outputs = c("R", "G", "B", "RGB"),
#          ve_npcs = 0.9,
#          comps = list(
#            list(name = "1v2",
#                 grpout = list(
#                   data = list(R = rimgs[grps != 3, ],
#                               G = gimgs[grps != 3, ],
#                               B = bimgs[grps != 3, ],
#                               RGB = rgbimgs[grps != 3, ]),
#                   obs_lvls = factor(grps[grps != 3]))),
#            list(name = "1v3",
#                 grpout = list(
#                   data = list(R = rimgs[grps != 2, ],
#                               G = gimgs[grps != 2, ],
#                               B = bimgs[grps != 2, ],
#                               RGB = rgbimgs[grps != 2, ]),
#                   obs_lvls = factor(grps[grps != 2]))),
#            list(name = "2v3",
#                 grpout = list(
#                   data = list(R = rimgs[grps != 1, ],
#                               G = gimgs[grps != 1, ],
#                               B = bimgs[grps != 1, ],
#                               RGB = rgbimgs[grps != 1, ]),
#                   obs_lvls = factor(grps[grps != 1])))))

## ----derma2, eval = FALSE-----------------------------------------------------
# toLatex(mic,
#         booktabs = T,
#         data_show = c("parp-1", "nparp-1", "scoreplot"),
#         data_labels = c("$t$-test", "$U$ test", "PCS"),
#         pvalf_params = list(minval = 1e-8, na_str = "*"),
#         label = "tab:ph2",
#         caption = paste("Comparison of PH$^2$ dataset images",
#                         "grouped by lesion type."))