Nothing
## ----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)
Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.