## ----initialize, echo = FALSE, message = FALSE--------------------------------
library("Ternary")
library("CongreveLamsdell2016")
data("clBremPartitions", "clBremQuartets", "clMkvPartitions", "clMkvQuartets",
     "clBootFreqPartitions", "clBootFreqQuartets", 
     "clBootGcPartitions", "clBootGcQuartets",
     "clJackFreqPartitions", "clJackFreqQuartets",
     "clJackGcPartitions", "clJackGcQuartets", 
     "clCI")

## ----initialize-variables, echo = FALSE---------------------------------------
PCH <- c(
  brem  = 2,
  bootG = 0,
  bootF = 5,
  jackG = 3,
  jackF = 4,
  none  = NA,
  
  mk     = 1,
  markov = 1,
  equal  = 61, # "="
  eq     = 61,
  dot    = 183, # "."
  
  k1         = 4,
  implied1   = 4,
  k2         = 183, 
  implied2   = 183,
  k3         = 183, 
  implied3   = 183,
  k5         = 3,
  implied5   = 3,
  kX         = 183,
  implied10  = 183,
  implied20  = 183,
  implied200 = 183,
  kC = 17,
  impliedC = 17  # triupfilled
)

COL <- c(
  black      = paste0(cbPalette8[1], "99"),
  markov     = paste0(cbPalette8[4], "99"),
  mk         = paste0(cbPalette8[4], "99"),
  
  equal      = paste0(cbPalette8[8], "99"),
  eq        = paste0(cbPalette8[8], "99"),
  implied1   = paste0(cbPalette8[6], "42"),
  k1         = paste0(cbPalette8[6], "42"),
  implied2   = paste0(cbPalette8[6], "42"),
  k2         = paste0(cbPalette8[6], "42"),
  implied3   = paste0(cbPalette8[6], "42"),
  k3         = paste0(cbPalette8[6], "42"),
  implied5   = paste0(cbPalette8[6], "42"),
  k5         = paste0(cbPalette8[6], "42"),
  implied10  = paste0(cbPalette8[6], "99"),
  kX         = paste0(cbPalette8[6], "99"),
  implied20  = paste0(cbPalette8[6], "42"),
  implied200 = paste0(cbPalette8[6], "42"),
  kC         = paste0(cbPalette8[2], "99"),
  impliedC   = paste0(cbPalette8[2], "99")
)

MARGINS <- c(2.8, 0.3, 0.3, 0.3)
ROWS <- c(1, 3)

TernaryQuarts<-function(TREE = TREE, zoom = 1, padding = 0.1) {
  clInitializeTernaryQuarts(zoom, padding)
  clPlotQuartets(clBootGcQuartets, TREE, cex = 1.1, pch = PCH)
  clPlotQuartets(clMkvQuartets, TREE, cex = 1.1, pch = PCH["mk"])

  # Return:
  invisible()
}


TernaryParts<-function(TREE = TREE) {
  clInitializeTernarySplits()
  clPlotSplits(clBootGcPartitions, TREE, cex = 1.1, pch = PCH)
  clPlotSplits(clMkvPartitions, TREE, cex = 1.1, pch = PCH["mk"])
  
  # Return:
  invisible()
}

AddLegend <- function(pos = "bottomright") {
  legend(pos, cex = 0.8, bty = "n",
         lty = 1,
         pch = PCH[c("mk", "eq", "kX", "k5", "k3", "k1", "kC")], pt.cex = 1.1,
         col = COL[c("mk", "eq", "kX", "k5", "k3", "k1", "kC")],
         legend = c("Markov", "Equal weights", 
                  paste0("Implied, k = ", c(10, 5, "2, 3", 1, "2..10"))))
}

Plottem <- function (i) { 
  TernaryQuarts(TREE = i)
  
  title(main = paste0("\nQuartets"), cex.main = 0.8)
  arrows(sqrt(3/4) * 0.5, 0.5, sqrt(3/4) * 0.8, 0.5, length = 0.1)
  text  (sqrt(3/4) * 0.65, 0.5, pos = 3, "Decreasing resolution", cex = 0.8)
  
  arrows(sqrt(3/4) * 0.98, 0.40, sqrt(3/4) * 0.98, 0.20, length = 0.1)
  text  (sqrt(3/4) * 1.01, 0.30, pos = 3, "Increasing divergence",
         cex = 0.8, srt = 270)
  
  TernaryQuarts(TREE = i, zoom = 3, padding = 0.01)
  title(main = paste0("\nDataset ", i, ": CI = ",round(clCI[i], 2)), cex.main = 1.2)
  
  TernaryParts(TREE = i)
  
  arrows(sqrt(3/4) * 0.98, 0.40, sqrt(3/4) * 0.98, 0.20, length = 0.1)
  text  (sqrt(3/4) * 1.01, 0.30, pos = 3, "Increasing RF distance",
         cex = 0.8, srt = 270)
  
  AddLegend()
}

## ----Summary, echo = FALSE, fig.height = 3, fig.width = 9, warning = FALSE----
x <- Sys.setlocale("LC_ALL", "C") # allows PCH = 183
par(mfrow = c(1, ROWS[2]), mar = MARGINS)

clInitializeTernaryQuarts()
clPlotAverageQuartets(clBootGcQuartets, cex = 1.1, pch = PCH)
clPlotAverageQuartets(clMkvQuartets, cex = 1.1, pch = PCH["mk"])

title(main = paste0("\nQuartets"), cex.main = 0.8)
arrows(sqrt(3/4) * 0.5, 0.5, sqrt(3/4) * 0.8, 0.5, length = 0.1)
text  (sqrt(3/4) * 0.65, 0.5, pos = 3, "Decreasing resolution", cex = 0.8)

arrows(sqrt(3/4) * 0.98, 0.40, sqrt(3/4) * 0.98, 0.20, length = 0.1)
text  (sqrt(3/4) * 1.01, 0.30, pos = 3, "Increasing divergence",
       cex = 0.8, srt = 270)


clInitializeTernaryQuarts(zoom = 3, padding = 0.01)
clPlotAverageQuartets(clBootGcQuartets, cex = 1.1, pch = PCH)
clPlotAverageQuartets(clMkvQuartets, cex = 1.1, pch = PCH["mk"])
title("\nResults for all datasets (means)", cex.main = 1.2)

clInitializeTernarySplits()
clPlotAverageSplits(clBootGcPartitions, cex = 1.1, pch = PCH)
clPlotAverageSplits(clMkvPartitions, cex = 1.1, pch = PCH["mk"])

arrows(sqrt(3/4) * 0.98, 0.40, sqrt(3/4) * 0.98, 0.20, length = 0.1)
text  (sqrt(3/4) * 1.01, 0.30, pos = 3, "Increasing RF distance",
       cex = 0.8, srt = 270)

AddLegend()

## ----echo = FALSE, fig.width = 9, fig.height = 6, warning = FALSE-------------
par(mfrow = ROWS, mar = MARGINS); x <- Plottem(1)

## ----echo = FALSE, fig.width = 9, fig.height = 6, warning = FALSE-------------
par(mfrow = ROWS, mar = MARGINS); x <- Plottem(2)

## ----echo = FALSE, fig.width = 9, fig.height = 6, warning = FALSE-------------
par(mfrow = ROWS, mar = MARGINS); x <- Plottem(3)

## ----echo = FALSE, fig.width = 9, fig.height = 6, warning = FALSE-------------
par(mfrow = ROWS, mar = MARGINS); x <- Plottem(4)

## ----echo = FALSE, fig.width = 9, fig.height = 6, warning = FALSE-------------
par(mfrow = ROWS, mar = MARGINS); x <- Plottem(5)

## ----echo = FALSE, fig.width = 9, fig.height = 6, warning = FALSE-------------
par(mfrow = ROWS, mar = MARGINS); x <- Plottem(6)

## ----echo = FALSE, fig.width = 9, fig.height = 6, warning = FALSE-------------
par(mfrow = ROWS, mar = MARGINS); x <- Plottem(7)

## ----echo = FALSE, fig.width = 9, fig.height = 6, warning = FALSE-------------
par(mfrow = ROWS, mar = MARGINS); x <- Plottem(8)

## ----echo = FALSE, fig.width = 9, fig.height = 6, warning = FALSE-------------
par(mfrow = ROWS, mar = MARGINS); x <- Plottem(9)

## ----echo = FALSE, fig.width = 9, fig.height = 6, warning = FALSE-------------
par(mfrow = ROWS, mar = MARGINS); x <- Plottem(10)

## ----echo = FALSE, fig.width = 9, fig.height = 6, warning = FALSE-------------
par(mfrow = ROWS, mar = MARGINS); x <- Plottem(11)

## ----echo = FALSE, fig.width = 9, fig.height = 6, warning = FALSE-------------
par(mfrow = ROWS, mar = MARGINS); x <- Plottem(12)

## ----echo = FALSE, fig.width = 9, fig.height = 6, warning = FALSE-------------
par(mfrow = ROWS, mar = MARGINS); x <- Plottem(13)

## ----echo = FALSE, fig.width = 9, fig.height = 6, warning = FALSE-------------
par(mfrow = ROWS, mar = MARGINS); x <- Plottem(14)

## ----echo = FALSE, fig.width = 9, fig.height = 6, warning = FALSE-------------
par(mfrow = ROWS, mar = MARGINS); x <- Plottem(15)

## ----echo = FALSE, fig.width = 9, fig.height = 6, warning = FALSE-------------
par(mfrow = ROWS, mar = MARGINS); x <- Plottem(16)

## ----echo = FALSE, fig.width = 9, fig.height = 6, warning = FALSE-------------
par(mfrow = ROWS, mar = MARGINS); x <- Plottem(17)

## ----echo = FALSE, fig.width = 9, fig.height = 6, warning = FALSE-------------
par(mfrow = ROWS, mar = MARGINS); x <- Plottem(18)

## ----echo = FALSE, fig.width = 9, fig.height = 6, warning = FALSE-------------
par(mfrow = ROWS, mar = MARGINS); x <- Plottem(19)

## ----echo = FALSE, fig.width = 9, fig.height = 6, warning = FALSE-------------
par(mfrow = ROWS, mar = MARGINS); x <- Plottem(20)

## ----echo = FALSE, fig.width = 9, fig.height = 6, warning = FALSE-------------
par(mfrow = ROWS, mar = MARGINS); x <- Plottem(21)

## ----echo = FALSE, fig.width = 9, fig.height = 6, warning = FALSE-------------
par(mfrow = ROWS, mar = MARGINS); x <- Plottem(22)

## ----echo = FALSE, fig.width = 9, fig.height = 6, warning = FALSE-------------
par(mfrow = ROWS, mar = MARGINS); x <- Plottem(23)

## ----echo = FALSE, fig.width = 9, fig.height = 6, warning = FALSE-------------
par(mfrow = ROWS, mar = MARGINS); x <- Plottem(24)

## ----echo = FALSE, fig.width = 9, fig.height = 6, warning = FALSE-------------
par(mfrow = ROWS, mar = MARGINS); x <- Plottem(25)

## ----echo = FALSE, fig.width = 9, fig.height = 6, warning = FALSE-------------
par(mfrow = ROWS, mar = MARGINS); x <- Plottem(26)

## ----echo = FALSE, fig.width = 9, fig.height = 6, warning = FALSE-------------
par(mfrow = ROWS, mar = MARGINS); x <- Plottem(27)

## ----echo = FALSE, fig.width = 9, fig.height = 6, warning = FALSE-------------
par(mfrow = ROWS, mar = MARGINS); x <- Plottem(28)

## ----echo = FALSE, fig.width = 9, fig.height = 6, warning = FALSE-------------
par(mfrow = ROWS, mar = MARGINS); x <- Plottem(29)

## ----echo = FALSE, fig.width = 9, fig.height = 6, warning = FALSE-------------
par(mfrow = ROWS, mar = MARGINS); x <- Plottem(30)

## ----echo = FALSE, fig.width = 9, fig.height = 6, warning = FALSE-------------
par(mfrow = ROWS, mar = MARGINS); x <- Plottem(31)

## ----echo = FALSE, fig.width = 9, fig.height = 6, warning = FALSE-------------
par(mfrow = ROWS, mar = MARGINS); x <- Plottem(32)

## ----echo = FALSE, fig.width = 9, fig.height = 6, warning = FALSE-------------
par(mfrow = ROWS, mar = MARGINS); x <- Plottem(33)

## ----echo = FALSE, fig.width = 9, fig.height = 6, warning = FALSE-------------
par(mfrow = ROWS, mar = MARGINS); x <- Plottem(34)

## ----echo = FALSE, fig.width = 9, fig.height = 6, warning = FALSE-------------
par(mfrow = ROWS, mar = MARGINS); x <- Plottem(35)

## ----echo = FALSE, fig.width = 9, fig.height = 6, warning = FALSE-------------
par(mfrow = ROWS, mar = MARGINS); x <- Plottem(36)

## ----echo = FALSE, fig.width = 9, fig.height = 6, warning = FALSE-------------
par(mfrow = ROWS, mar = MARGINS); x <- Plottem(37)

## ----echo = FALSE, fig.width = 9, fig.height = 6, warning = FALSE-------------
par(mfrow = ROWS, mar = MARGINS); x <- Plottem(38)

## ----echo = FALSE, fig.width = 9, fig.height = 6, warning = FALSE-------------
par(mfrow = ROWS, mar = MARGINS); x <- Plottem(39)

## ----echo = FALSE, fig.width = 9, fig.height = 6, warning = FALSE-------------
par(mfrow = ROWS, mar = MARGINS); x <- Plottem(40)

## ----echo = FALSE, fig.width = 9, fig.height = 6, warning = FALSE-------------
par(mfrow = ROWS, mar = MARGINS); x <- Plottem(41)

## ----echo = FALSE, fig.width = 9, fig.height = 6, warning = FALSE-------------
par(mfrow = ROWS, mar = MARGINS); x <- Plottem(42)

## ----echo = FALSE, fig.width = 9, fig.height = 6, warning = FALSE-------------
par(mfrow = ROWS, mar = MARGINS); x <- Plottem(43)

## ----echo = FALSE, fig.width = 9, fig.height = 6, warning = FALSE-------------
par(mfrow = ROWS, mar = MARGINS); x <- Plottem(44)

## ----echo = FALSE, fig.width = 9, fig.height = 6, warning = FALSE-------------
par(mfrow = ROWS, mar = MARGINS); x <- Plottem(45)

## ----echo = FALSE, fig.width = 9, fig.height = 6, warning = FALSE-------------
par(mfrow = ROWS, mar = MARGINS); x <- Plottem(46)

## ----echo = FALSE, fig.width = 9, fig.height = 6, warning = FALSE-------------
par(mfrow = ROWS, mar = MARGINS); x <- Plottem(47)

## ----echo = FALSE, fig.width = 9, fig.height = 6, warning = FALSE-------------
par(mfrow = ROWS, mar = MARGINS); x <- Plottem(48)

## ----echo = FALSE, fig.width = 9, fig.height = 6, warning = FALSE-------------
par(mfrow = ROWS, mar = MARGINS); x <- Plottem(49)

## ----echo = FALSE, fig.width = 9, fig.height = 6, warning = FALSE-------------
par(mfrow = ROWS, mar = MARGINS); x <- Plottem(50)

## ----echo = FALSE, fig.width = 9, fig.height = 6, warning = FALSE-------------
par(mfrow = ROWS, mar = MARGINS); x <- Plottem(51)

## ----echo = FALSE, fig.width = 9, fig.height = 6, warning = FALSE-------------
par(mfrow = ROWS, mar = MARGINS); x <- Plottem(52)

## ----echo = FALSE, fig.width = 9, fig.height = 6, warning = FALSE-------------
par(mfrow = ROWS, mar = MARGINS); x <- Plottem(53)

## ----echo = FALSE, fig.width = 9, fig.height = 6, warning = FALSE-------------
par(mfrow = ROWS, mar = MARGINS); x <- Plottem(54)

## ----echo = FALSE, fig.width = 9, fig.height = 6, warning = FALSE-------------
par(mfrow = ROWS, mar = MARGINS); x <- Plottem(55)

## ----echo = FALSE, fig.width = 9, fig.height = 6, warning = FALSE-------------
par(mfrow = ROWS, mar = MARGINS); x <- Plottem(56)

## ----echo = FALSE, fig.width = 9, fig.height = 6, warning = FALSE-------------
par(mfrow = ROWS, mar = MARGINS); x <- Plottem(57)

## ----echo = FALSE, fig.width = 9, fig.height = 6, warning = FALSE-------------
par(mfrow = ROWS, mar = MARGINS); x <- Plottem(58)

## ----echo = FALSE, fig.width = 9, fig.height = 6, warning = FALSE-------------
par(mfrow = ROWS, mar = MARGINS); x <- Plottem(59)

## ----echo = FALSE, fig.width = 9, fig.height = 6, warning = FALSE-------------
par(mfrow = ROWS, mar = MARGINS); x <- Plottem(60)

## ----echo = FALSE, fig.width = 9, fig.height = 6, warning = FALSE-------------
par(mfrow = ROWS, mar = MARGINS); x <- Plottem(61)

## ----echo = FALSE, fig.width = 9, fig.height = 6, warning = FALSE-------------
par(mfrow = ROWS, mar = MARGINS); x <- Plottem(62)

## ----echo = FALSE, fig.width = 9, fig.height = 6, warning = FALSE-------------
par(mfrow = ROWS, mar = MARGINS); x <- Plottem(63)

## ----echo = FALSE, fig.width = 9, fig.height = 6, warning = FALSE-------------
par(mfrow = ROWS, mar = MARGINS); x <- Plottem(64)

## ----echo = FALSE, fig.width = 9, fig.height = 6, warning = FALSE-------------
par(mfrow = ROWS, mar = MARGINS); x <- Plottem(65)

## ----echo = FALSE, fig.width = 9, fig.height = 6, warning = FALSE-------------
par(mfrow = ROWS, mar = MARGINS); x <- Plottem(66)

## ----echo = FALSE, fig.width = 9, fig.height = 6, warning = FALSE-------------
par(mfrow = ROWS, mar = MARGINS); x <- Plottem(67)

## ----echo = FALSE, fig.width = 9, fig.height = 6, warning = FALSE-------------
par(mfrow = ROWS, mar = MARGINS); x <- Plottem(68)

## ----echo = FALSE, fig.width = 9, fig.height = 6, warning = FALSE-------------
par(mfrow = ROWS, mar = MARGINS); x <- Plottem(69)

## ----echo = FALSE, fig.width = 9, fig.height = 6, warning = FALSE-------------
par(mfrow = ROWS, mar = MARGINS); x <- Plottem(70)

## ----echo = FALSE, fig.width = 9, fig.height = 6, warning = FALSE-------------
par(mfrow = ROWS, mar = MARGINS); x <- Plottem(71)

## ----echo = FALSE, fig.width = 9, fig.height = 6, warning = FALSE-------------
par(mfrow = ROWS, mar = MARGINS); x <- Plottem(72)

## ----echo = FALSE, fig.width = 9, fig.height = 6, warning = FALSE-------------
par(mfrow = ROWS, mar = MARGINS); x <- Plottem(73)

## ----echo = FALSE, fig.width = 9, fig.height = 6, warning = FALSE-------------
par(mfrow = ROWS, mar = MARGINS); x <- Plottem(74)

## ----echo = FALSE, fig.width = 9, fig.height = 6, warning = FALSE-------------
par(mfrow = ROWS, mar = MARGINS); x <- Plottem(75)

## ----echo = FALSE, fig.width = 9, fig.height = 6, warning = FALSE-------------
par(mfrow = ROWS, mar = MARGINS); x <- Plottem(76)

## ----echo = FALSE, fig.width = 9, fig.height = 6, warning = FALSE-------------
par(mfrow = ROWS, mar = MARGINS); x <- Plottem(77)

## ----echo = FALSE, fig.width = 9, fig.height = 6, warning = FALSE-------------
par(mfrow = ROWS, mar = MARGINS); x <- Plottem(78)

## ----echo = FALSE, fig.width = 9, fig.height = 6, warning = FALSE-------------
par(mfrow = ROWS, mar = MARGINS); x <- Plottem(79)

## ----echo = FALSE, fig.width = 9, fig.height = 6, warning = FALSE-------------
par(mfrow = ROWS, mar = MARGINS); x <- Plottem(80)

## ----echo = FALSE, fig.width = 9, fig.height = 6, warning = FALSE-------------
par(mfrow = ROWS, mar = MARGINS); x <- Plottem(81)

## ----echo = FALSE, fig.width = 9, fig.height = 6, warning = FALSE-------------
par(mfrow = ROWS, mar = MARGINS); x <- Plottem(82)

## ----echo = FALSE, fig.width = 9, fig.height = 6, warning = FALSE-------------
par(mfrow = ROWS, mar = MARGINS); x <- Plottem(83)

## ----echo = FALSE, fig.width = 9, fig.height = 6, warning = FALSE-------------
par(mfrow = ROWS, mar = MARGINS); x <- Plottem(84)

## ----echo = FALSE, fig.width = 9, fig.height = 6, warning = FALSE-------------
par(mfrow = ROWS, mar = MARGINS); x <- Plottem(85)

## ----echo = FALSE, fig.width = 9, fig.height = 6, warning = FALSE-------------
par(mfrow = ROWS, mar = MARGINS); x <- Plottem(86)

## ----echo = FALSE, fig.width = 9, fig.height = 6, warning = FALSE-------------
par(mfrow = ROWS, mar = MARGINS); x <- Plottem(87)

## ----echo = FALSE, fig.width = 9, fig.height = 6, warning = FALSE-------------
par(mfrow = ROWS, mar = MARGINS); x <- Plottem(88)

## ----echo = FALSE, fig.width = 9, fig.height = 6, warning = FALSE-------------
par(mfrow = ROWS, mar = MARGINS); x <- Plottem(89)

## ----echo = FALSE, fig.width = 9, fig.height = 6, warning = FALSE-------------
par(mfrow = ROWS, mar = MARGINS); x <- Plottem(90)

## ----echo = FALSE, fig.width = 9, fig.height = 6, warning = FALSE-------------
par(mfrow = ROWS, mar = MARGINS); x <- Plottem(91)

## ----echo = FALSE, fig.width = 9, fig.height = 6, warning = FALSE-------------
par(mfrow = ROWS, mar = MARGINS); x <- Plottem(92)

## ----echo = FALSE, fig.width = 9, fig.height = 6, warning = FALSE-------------
par(mfrow = ROWS, mar = MARGINS); x <- Plottem(93)

## ----echo = FALSE, fig.width = 9, fig.height = 6, warning = FALSE-------------
par(mfrow = ROWS, mar = MARGINS); x <- Plottem(94)

## ----echo = FALSE, fig.width = 9, fig.height = 6, warning = FALSE-------------
par(mfrow = ROWS, mar = MARGINS); x <- Plottem(95)

## ----echo = FALSE, fig.width = 9, fig.height = 6, warning = FALSE-------------
par(mfrow = ROWS, mar = MARGINS); x <- Plottem(96)

## ----echo = FALSE, fig.width = 9, fig.height = 6, warning = FALSE-------------
par(mfrow = ROWS, mar = MARGINS); x <- Plottem(97)

## ----echo = FALSE, fig.width = 9, fig.height = 6, warning = FALSE-------------
par(mfrow = ROWS, mar = MARGINS); x <- Plottem(98)

## ----echo = FALSE, fig.width = 9, fig.height = 6, warning = FALSE-------------
par(mfrow = ROWS, mar = MARGINS); x <- Plottem(99)

## ----echo = FALSE, fig.width = 9, fig.height = 6, warning = FALSE-------------
par(mfrow = ROWS, mar = MARGINS); x <- Plottem(100)