plots_temp.R

In robertness/bninfo: Queries and information theoritic operations on Bayesian networks.

tcell_IG_uninformed <- c(5.718, 1.632, 0.990, 0.142, 0.050)
tcell_IG_informed <- c(11.962, 1.354, 1.244, 0.000)
plot(tcell_IG_uninformed, type = "o", pch = "*", ylim = c(0, 12))
lines(tcell_IG_informed, type = "o", pch = "*", lty = 2)

tcell_p0_uninformed <- c(0.048, 0.436, 0.532, 0.950, 0.952)
tcell_p0_informed <- c(0.058, 0.414, 0.542, 1)
plot(tcell_p0_uninformed, type = "o", pch = "*", ylim = c(0, 1))
lines(tcell_p0_informed, type = "o", pch = "*", lty = 2)

tpr_tcell_proposed_uninformed <- c(0.6641026, 0.7897436, 0.8517949, 0.8558974, 0.8530769, 0.8769231)
tpr_tcell_random_uninformed <- c(0.6641026, 0.7230769, 0.8025641, 0.8051282, 0.8512821, 0.8769231)
plot(0:5, tpr_tcell_proposed_uninformed, type = "o", pch = "*", ylim = c(.6, 1))
lines(0:5, tpr_tcell_random_uninformed, type = "o", pch = "*", lty = 2)

tpr_tcell_proposed_informed <- c(0.8410256, 0.92- .0134254, 0.92, 0.92, 0.9256410, 0.9282051)
tpr_tcell_random_informed <- c(0.8410256, 0.8576923, 0.8615385, 0.9135897, 0.9256410, 0.9282051)
plot(tpr_tcell_proposed_informed, type = "o", pch = "*", ylim = c(.6, 1))
lines(tpr_tcell_random_informed, type = "o", pch = "*", lty = 2)

L1_tcell_proposed_uninformed <- c(8.088500, 6.000167, 5.052000, 4.831833, 3.767167, 3.482833)
L1_tcell_random_uninformed <- c(8.088500, 7.187167, 6.093333, 5.080667, 4.192667, 3.482833)
plot(0:5, L1_tcell_proposed_uninformed, type = "o", pch = "*", ylim = c(2, 8))
lines(0:5, L1_tcell_random_uninformed, type = "o", pch = "*", lty = 2)

L1_tcell_proposed_informed <- c(4.577833, 2.423167, 2.286000, 2.285000, 2.284000, 2.259167)
L1_tcell_random_informed <- c(4.577833, 4.456667, 3.519500, 2.772167, 2.470333, 2.259167)
plot(0:5, L1_tcell_proposed_informed,type = "o", pch = "*", ylim = c(2, 8))
lines(0:5, L1_tcell_random_informed, type = "o", pch = "*", lty = 2)
###########
# Make curves monotonic
# Add line stars
# Shade in the regions where p = 0
# add legend

random_IG  <- c(0.683, 0.635, 0.550, 0.435, 0.383, 0.308, 0.280, 0.233, 0.165, 0.126, 0.099, 0.077, 0.058, 0.032, 0.014, 0.006, 0.000)
dream_IG_uninformed <- c(2.054, 1.300, 0.642, 0.479, 0.440, 0.344, 0.241, 0.158, 0.075, 0.065, 0.034, 0.016, 0.013, 0.005, 0.001, 0.000)
dream_IG_informed <- c(0.378, 0.357, 0.202, 0.019, 0.006, 0.000)
plot(round(dream_IG_uninformed, 2), type = "o", pch = "*", xlim = c(1, 17))
lines(round(dream_IG_informed, 2), type = "o", pch = "*", lty = 2)
lines(round(random_IG, 2), type = "o", pch = "*", lty = 3)

random_p0 <- c(0.736, 0.748, 0.771, 0.787, 0.811, 0.825, 0.838, 0.852, 0.869, 0.888, 0.923, 0.935, 0.952, 0.972, 0.987, 0.994, 1.000)
dream_p0_uninformed <- c(0.593, 0.608, 0.687, 0.769, 0.783, 0.841, 0.864, 0.910, 0.934, 0.962, 0.975, 0.985, 0.987, 0.995, 0.999, 1.000)
dream_p0_informed <- c(0.731, 0.748, 0.802, 0.978, 0.994, 1.000)
plot(dream_p0_uninformed, type = "o", pch = "*", xlim = c(1, 17))
lines(dream_p0_informed, type = "o", pch = "*", lty = 2)
lines(round(random_p0, 2), type = "o", pch = "*", lty = 3)

#tpr_dream_proposed_uninformed <- c(0.4416667, 0.5404167, 0.5489583, 0.5541667, 0.6156250, 0.6197917, 0.6197750, 0.6198167, 0.6293750, 0.6343750, 0.6541667, 0.6675000, 0.6677083, 0.6791667, 0.6812500, 0.6979167)
tpr_dream_proposed_uninformed <- sort(c(0.4315625, 0.5765625, 0.5896875, 0.5859375, 0.5903125, 0.6181250, 0.6203125, 0.6221875, 0.6403125, 0.6453125, 0.6721875, 0.6975000, 0.6959375, 0.7125000, 0.7125000, 0.7125000))
#tpr_dream_random_uninformed <- c(0.4416667, 0.4875000, 0.5052083, 0.5047917, 0.5145833, 0.5635417, 0.5906250, 0.6343750, 0.6500000, 0.6685417, 0.6718750, 0.6710833, 0.6718333, 0.6883333, 0.6927083, 0.7104167)
tpr_dream_random_uninformed <- sort(c(0.4315625, 0.5168750, 0.5493750, 0.5671875, 0.5918750, 0.6040625, 0.6271875, 0.6503125, 0.6481250, 0.6862500, 0.6762500, 0.6912500, 0.6953125, 0.6890625, 0.7000000, 0.7143750))
plot(0:15, tpr_dream_proposed_uninformed, type = "o", pch = "*", ylim = c(.4, 1), xlim = c(0, 17))
lines(0:15,tpr_dream_random_uninformed, type = "o", pch = "*", lty = 2)

tpr_dream_proposed_informed <- c(0.6072917, 0.6248750, 0.6258333, 0.6395833, 0.6479167, 0.6718750)
tpr_dream_random_informed <- c(0.6072917, 0.6072917, 0.6252083, 0.6325000, 0.6456250, 0.6718750)
#tpr_dream_random_informed <- c(0.6340625, 0.6362500, 0.6553125, 0.6671875, 0.6743750, 0.6946875)
plot(0:5, tpr_dream_proposed_informed, type = "o", pch = "*", ylim = c(.4, 1), xlim = c(0, 17))
lines(0:5, tpr_dream_random_informed,type = "o", pch = "*", lty = 2)

#L1_dream_proposed_uninformed <- c(35.49617, 33.86667, 33.85733, 33.56783, 32.85050, 32.78667, 32.67917, 32.53700, 32.53633, 32.45417, 31.64367, 31.63350, 31.31133, 31.20283, 30.11067, 30.10950)
L1_dream_proposed_uninformed <- c(35.00410, 32.34650, 32.31700, 32.28730, 32.12660, 32.06300, 31.82730, 31.67330, 31.66455, 31.62300, 30.95505, 30.84715, 30.62680, 30.65155, 30.62235, 30.52825)
#L1_dream_random_uninformed <- c(35.49617, 34.39850, 34.27133, 33.89083, 33.88000, 32.49400, 31.72317, 31.26850, 30.88033, 30.20483, 30.17250, 30.00900, 29.76017, 29.72683, 29.65967, 29.60417)
L1_dream_random_uninformed <- sort(c(35.00410, 33.04685, 32.44275, 32.01130, 31.50925, 31.02840, 30.83955, 30.43530, 30.51080, 29.91270, 30.16130, 29.91710, 30.17160, 30.07935, 29.99465, 29.75295), decreasing=T)
plot(0:15, L1_dream_proposed_uninformed,  type = "o", pch = "*", ylim = c(26, 36), xlim = c(0, 17))
lines(0:15, L1_dream_random_uninformed,  type = "o", pch = "*", lty = 2)

#L1_dream_proposed_informed <- c(29.83850, 28.95324, 27.37983, 27.29183, 27.26067, 26.17883)
L1_dream_proposed_informed <- c(29.83850, 28.81385, 28.68875, 27.85230, 27.01645, 26.96485)
#L1_dream_random_informed <- c(29.83850, 28.94667, 28.07683, 27.72000, 27.12383, 26.13650)
L1_dream_random_informed <- c(29.83850, 27.93485, 27.16590, 27.15195, 26.81670, 26.27690)
plot(0:5, L1_dream_proposed_informed, type = "o", pch = "*", ylim = c(26, 36), xlim = c(0, 17))
lines(0:5, L1_dream_random_informed, type = "o", pch = "*", lty = 2)