set term pngcairo size 800,600 nocrop enhanced font 'Verdana,11'#define axis
set style line 11 lc rgb '#808080' lt 1
set border 3 back ls 11
set tics nomirror

#define key
#set key opaque
set key outside

# define grid
set style line 12 lc rgb '#808080' lt 0 lw 1
set grid back ls 12

# define linecolors

set style line 1 lc rgb '#0060ad' pt 7 ps 0.5 pi -1 lt 0 lw 2 # --- blue
set style line 2 lc rgb '#8b1a0e' pt 7 ps 0.5 pi -1 lt 0 lw 2 # --- red
set style line 3 lc rgb '#5e9c36' pt 7 ps 0.5 pi -1 lt 0 lw 2 # --- green
set style line 4 lc rgb '#ffa500' pt 7 ps 0.5 pi -1 lt 0 lw 2 # --- orange
set style line 5 lc rgb '#40e0d0' pt 7 ps 0.5 pi -1 lt 0 lw 2 # --- turquoise
set style line 6 lc rgb '#9400d3' pt 7 ps 0.5 pi -1 lt 0 lw 2 # --- darkviolet
set style line 7 lc rgb '#ff00ff' pt 7 ps 0.5 pi -1 lt 0 lw 2 # --- magenta
set style line 8 lc rgb '#c0c0c0' pt 7 ps 0.5 pi -1 lt 0 lw 2 # --- silver
set style line 9 lc rgb '#e6e6Fa' pt 7 ps 0.5 pi -1 lt 0 lw 2 # --- lavender

set pointintervalbox 1

set samples 300
#Start of user script
#---------------------

single_0       = 2
multi_0        = 3
single_default = 4
multi_default  = 5
ptp_default    = 6
real           = 7




# Kassian Score

set title "Average Kassian Score similar GMYC taxa"
set xlabel "Set number"
set ylabel "Average tree score"

ExtData1 = 'workfile_tree_scores'
set output 'plots/average_tree_scores.png'

plot ExtData1 using 1:2 title 'delimit single minbr 0' with linespoints ls 1, ExtData1 using 1:3 title 'delimit multi minbr 0' with linespoints ls 2, ExtData1 using 1:4 title 'delimit single minbr default' with linespoints ls 3, ExtData1 using 1:5 title 'delimit multi minbr default' with linespoints ls 4, ExtData1 using 1:6 title 'PTP minbr default' with linespoints ls 5



set xlabel "index"
set ylabel "Tree Score"

ExtData1_10000 = 'workfile_tree_scores_Ne10000'
ExtData1_100000 = 'workfile_tree_scores_Ne100000'
ExtData1_500000 = 'workfile_tree_scores_Ne500000'
ExtData1_1000000 = 'workfile_tree_scores_Ne1000000'

set title "Kassian Tree Score similar GMYC taxa delimit single minbr 0"
set output 'plots/tree_scores_delimit_single_minbr_0.png'
plot ExtData1_10000 using 1:single_0 title 'Ne10000' with linespoints ls 1, ExtData1_100000 using 1:single_0 title 'Ne10000' with linespoints ls 2, ExtData1_500000 using 1:single_0 title 'Ne500000' with linespoints ls 3, ExtData1_1000000 using 1:single_0 title 'Ne1000000' with linespoints ls 4

set title "Kassian Tree Score similar GMYC taxa delimit multi minbr 0"
set output 'plots/tree_scores_delimit_multi_minbr_0.png'
plot ExtData1_10000 using 1:multi_0 title 'Ne10000' with linespoints ls 1, ExtData1_100000 using 1:multi_0 title 'Ne10000' with linespoints ls 2, ExtData1_500000 using 1:multi_0 title 'Ne500000' with linespoints ls 3, ExtData1_1000000 using 1:multi_0 title 'Ne1000000' with linespoints ls 4

set title "Kassian Tree Score similar GMYC taxa delimit single minbr default"
set output 'plots/tree_scores_delimit_single_minbr_default.png'
plot ExtData1_10000 using 1:single_default title 'Ne10000' with linespoints ls 1, ExtData1_100000 using 1:single_default title 'Ne10000' with linespoints ls 2, ExtData1_500000 using 1:single_default title 'Ne500000' with linespoints ls 3, ExtData1_1000000 using 1:single_default title 'Ne1000000' with linespoints ls 4

set title "Kassian Tree Score similar GMYC taxa delimit multi minbr default"
set output 'plots/tree_scores_delimit_multi_minbr_default.png'
plot ExtData1_10000 using 1:multi_default title 'Ne10000' with linespoints ls 1, ExtData1_100000 using 1:multi_default title 'Ne10000' with linespoints ls 2, ExtData1_500000 using 1:multi_default title 'Ne500000' with linespoints ls 3, ExtData1_1000000 using 1:multi_default title 'Ne1000000' with linespoints ls 4

set title "Kassian Tree Score similar GMYC taxa PTP minbr default"
set output 'plots/tree_scores_PTP_minbr_default.png'
plot ExtData1_10000 using 1:ptp_default title 'Ne10000' with linespoints ls 1, ExtData1_100000 using 1:ptp_default title 'Ne10000' with linespoints ls 2, ExtData1_500000 using 1:ptp_default title 'Ne500000' with linespoints ls 3, ExtData1_1000000 using 1:ptp_default title 'Ne1000000' with linespoints ls 4



# NMI score

set title "Average NMI Score similar GMYC taxa"
set xlabel "Set number"
set ylabel "Average NMI score"

ExtData1 = 'workfile_nmi_scores'
set output 'plots/average_nmi_scores.png'

plot ExtData1 using 1:2 title 'delimit single minbr 0' with linespoints ls 1, ExtData1 using 1:3 title 'delimit multi minbr 0' with linespoints ls 2, ExtData1 using 1:4 title 'delimit single minbr default' with linespoints ls 3, ExtData1 using 1:5 title 'delimit multi minbr default' with linespoints ls 4, ExtData1 using 1:6 title 'PTP minbr default' with linespoints ls 5



set xlabel "index"
set ylabel "NMI Score"

ExtData1_10000 = 'workfile_nmi_scores_Ne10000'
ExtData1_100000 = 'workfile_nmi_scores_Ne100000'
ExtData1_500000 = 'workfile_nmi_scores_Ne500000'
ExtData1_1000000 = 'workfile_nmi_scores_Ne1000000'

set title "NMI Score similar GMYC taxa delimit single minbr 0"
set output 'plots/nmi_scores_delimit_single_minbr_0.png'
plot ExtData1_10000 using 1:single_0 title 'Ne10000' with linespoints ls 1, ExtData1_100000 using 1:single_0 title 'Ne10000' with linespoints ls 2, ExtData1_500000 using 1:single_0 title 'Ne500000' with linespoints ls 3, ExtData1_1000000 using 1:single_0 title 'Ne1000000' with linespoints ls 4

set title "NMI Score similar GMYC taxa delimit multi minbr 0"
set output 'plots/nmi_scores_delimit_multi_minbr_0.png'
plot ExtData1_10000 using 1:multi_0 title 'Ne10000' with linespoints ls 1, ExtData1_100000 using 1:multi_0 title 'Ne10000' with linespoints ls 2, ExtData1_500000 using 1:multi_0 title 'Ne500000' with linespoints ls 3, ExtData1_1000000 using 1:multi_0 title 'Ne1000000' with linespoints ls 4

set title "NMI Score similar GMYC taxa delimit single minbr default"
set output 'plots/nmi_scores_delimit_single_minbr_default.png'
plot ExtData1_10000 using 1:single_default title 'Ne10000' with linespoints ls 1, ExtData1_100000 using 1:single_default title 'Ne10000' with linespoints ls 2, ExtData1_500000 using 1:single_default title 'Ne500000' with linespoints ls 3, ExtData1_1000000 using 1:single_default title 'Ne1000000' with linespoints ls 4

set title "NMI Score similar GMYC taxa delimit multi minbr default"
set output 'plots/nmi_scores_delimit_multi_minbr_default.png'
plot ExtData1_10000 using 1:multi_default title 'Ne10000' with linespoints ls 1, ExtData1_100000 using 1:multi_default title 'Ne10000' with linespoints ls 2, ExtData1_500000 using 1:multi_default title 'Ne500000' with linespoints ls 3, ExtData1_1000000 using 1:multi_default title 'Ne1000000' with linespoints ls 4

set title "NMI Score similar GMYC taxa PTP minbr default"
set output 'plots/nmi_scores_PTP_minbr_default.png'
plot ExtData1_10000 using 1:ptp_default title 'Ne10000' with linespoints ls 1, ExtData1_100000 using 1:ptp_default title 'Ne10000' with linespoints ls 2, ExtData1_500000 using 1:ptp_default title 'Ne500000' with linespoints ls 3, ExtData1_1000000 using 1:ptp_default title 'Ne1000000' with linespoints ls 4



# number of species

set title "Average Number of Species similar GMYC taxa"
set xlabel "Set number"
set ylabel "Average number of species"

ExtData1 = 'workfile_num_species'
set output 'plots/average_num_species.png'

plot ExtData1 using 1:2 title 'delimit single minbr 0' with linespoints ls 1, ExtData1 using 1:3 title 'delimit multi minbr 0' with linespoints ls 2, ExtData1 using 1:4 title 'delimit single minbr default' with linespoints ls 3, ExtData1 using 1:5 title 'delimit multi minbr default' with linespoints ls 4, ExtData1 using 1:6 title 'PTP minbr default' with linespoints ls 5, ExtData1 using 1:7 title 'real' with linespoints ls 7



set xlabel "index"
set ylabel "Number of Species"

ExtData1_10000 = 'workfile_num_species_Ne10000'
ExtData1_100000 = 'workfile_num_species_Ne100000'
ExtData1_500000 = 'workfile_num_species_Ne500000'
ExtData1_1000000 = 'workfile_num_species_Ne1000000'

set title "Number of Species similar GMYC taxa delimit single minbr 0"
set output 'plots/num_species_delimit_single_minbr_0.png'
plot ExtData1_10000 using 1:single_0 title 'Ne10000' with linespoints ls 1, ExtData1_100000 using 1:single_0 title 'Ne10000' with linespoints ls 2, ExtData1_500000 using 1:single_0 title 'Ne500000' with linespoints ls 3, ExtData1_1000000 using 1:single_0 title 'Ne1000000' with linespoints ls 4

set title "Number of Species similar GMYC taxa delimit multi minbr 0"
set output 'plots/num_species_delimit_multi_minbr_0.png'
plot ExtData1_10000 using 1:multi_0 title 'Ne10000' with linespoints ls 1, ExtData1_100000 using 1:multi_0 title 'Ne10000' with linespoints ls 2, ExtData1_500000 using 1:multi_0 title 'Ne500000' with linespoints ls 3, ExtData1_1000000 using 1:multi_0 title 'Ne1000000' with linespoints ls 4

set title "Number of Species similar GMYC taxa delimit single minbr default"
set output 'plots/num_species_delimit_single_minbr_default.png'
plot ExtData1_10000 using 1:single_default title 'Ne10000' with linespoints ls 1, ExtData1_100000 using 1:single_default title 'Ne10000' with linespoints ls 2, ExtData1_500000 using 1:single_default title 'Ne500000' with linespoints ls 3, ExtData1_1000000 using 1:single_default title 'Ne1000000' with linespoints ls 4

set title "Number of Species similar GMYC taxa delimit multi minbr default"
set output 'plots/num_species_delimit_multi_minbr_default.png'
plot ExtData1_10000 using 1:multi_default title 'Ne10000' with linespoints ls 1, ExtData1_100000 using 1:multi_default title 'Ne10000' with linespoints ls 2, ExtData1_500000 using 1:multi_default title 'Ne500000' with linespoints ls 3, ExtData1_1000000 using 1:multi_default title 'Ne1000000' with linespoints ls 4

set title "Number of Species similar GMYC taxa PTP minbr default"
set output 'plots/num_species_PTP_minbr_default.png'
plot ExtData1_10000 using 1:ptp_default title 'Ne10000' with linespoints ls 1, ExtData1_100000 using 1:ptp_default title 'Ne10000' with linespoints ls 2, ExtData1_500000 using 1:ptp_default title 'Ne500000' with linespoints ls 3, ExtData1_1000000 using 1:ptp_default title 'Ne1000000' with linespoints ls 4

set title "Number of Species similar GMYC taxa real"
set output 'plots/num_species_real.png'
plot ExtData1_10000 using 1:real title 'Ne10000' with linespoints ls 1, ExtData1_100000 using 1:real title 'Ne10000' with linespoints ls 2, ExtData1_500000 using 1:real title 'Ne500000' with linespoints ls 3, ExtData1_1000000 using 1:real title 'Ne1000000' with linespoints ls 4



set xlabel "index"
set ylabel "Delta Number of Species"

ExtData1_10000 = 'workfile_delta_species_Ne10000'
ExtData1_100000 = 'workfile_delta_species_Ne100000'
ExtData1_500000 = 'workfile_delta_species_Ne500000'
ExtData1_1000000 = 'workfile_delta_species_Ne1000000'

set title "Delta Number of Species similar GMYC taxa delimit single minbr 0"
set output 'plots/delta_species_delimit_single_minbr_0.png'
plot ExtData1_10000 using 1:single_0 title 'Ne10000' with linespoints ls 1, ExtData1_100000 using 1:single_0 title 'Ne10000' with linespoints ls 2, ExtData1_500000 using 1:single_0 title 'Ne500000' with linespoints ls 3, ExtData1_1000000 using 1:single_0 title 'Ne1000000' with linespoints ls 4

set title "Delta Number of Species similar GMYC taxa delimit multi minbr 0"
set output 'plots/delta_species_delimit_multi_minbr_0.png'
plot ExtData1_10000 using 1:multi_0 title 'Ne10000' with linespoints ls 1, ExtData1_100000 using 1:multi_0 title 'Ne10000' with linespoints ls 2, ExtData1_500000 using 1:multi_0 title 'Ne500000' with linespoints ls 3, ExtData1_1000000 using 1:multi_0 title 'Ne1000000' with linespoints ls 4

set title "Delta Number of Species similar GMYC taxa delimit single minbr default"
set output 'plots/delta_species_delimit_single_minbr_default.png'
plot ExtData1_10000 using 1:single_default title 'Ne10000' with linespoints ls 1, ExtData1_100000 using 1:single_default title 'Ne10000' with linespoints ls 2, ExtData1_500000 using 1:single_default title 'Ne500000' with linespoints ls 3, ExtData1_1000000 using 1:single_default title 'Ne1000000' with linespoints ls 4

set title "Delta Number of Species similar GMYC taxa delimit multi minbr default"
set output 'plots/delta_species_delimit_multi_minbr_default.png'
plot ExtData1_10000 using 1:multi_default title 'Ne10000' with linespoints ls 1, ExtData1_100000 using 1:multi_default title 'Ne10000' with linespoints ls 2, ExtData1_500000 using 1:multi_default title 'Ne500000' with linespoints ls 3, ExtData1_1000000 using 1:multi_default title 'Ne1000000' with linespoints ls 4

set title "Delta Number of Species similar GMYC taxa PTP minbr default"
set output 'plots/delta_species_PTP_minbr_default.png'
plot ExtData1_10000 using 1:ptp_default title 'Ne10000' with linespoints ls 1, ExtData1_100000 using 1:ptp_default title 'Ne10000' with linespoints ls 2, ExtData1_500000 using 1:ptp_default title 'Ne500000' with linespoints ls 3, ExtData1_1000000 using 1:ptp_default title 'Ne1000000' with linespoints ls 4

set title "Delta Number of Species similar GMYC taxa real"
set output 'plots/delta_species_real.png'
plot ExtData1_10000 using 1:real title 'Ne10000' with linespoints ls 1, ExtData1_100000 using 1:real title 'Ne10000' with linespoints ls 2, ExtData1_500000 using 1:real title 'Ne500000' with linespoints ls 3, ExtData1_1000000 using 1:real title 'Ne1000000' with linespoints ls 4



# single lambda score

set title "Average Single Lambda Score similar GMYC taxa"
set xlabel "Set number"
set ylabel "Average single lambda score"

ExtData1 = 'workfile_single_scores'
set output 'plots/average_single_scores.png'

plot ExtData1 using 1:2 title 'delimit single minbr 0' with linespoints ls 1, ExtData1 using 1:3 title 'delimit multi minbr 0' with linespoints ls 2, ExtData1 using 1:4 title 'delimit single minbr default' with linespoints ls 3, ExtData1 using 1:5 title 'delimit multi minbr default' with linespoints ls 4, ExtData1 using 1:6 title 'PTP minbr default' with linespoints ls 5, ExtData1 using 1:7 title 'real minbr 0' with linespoints ls 7, ExtData1 using 1:8 title 'real minbr default' with linespoints ls 8




set xlabel "index"
set ylabel "Single Lambda Score"

ExtData1_10000 = 'workfile_single_scores_Ne10000'
ExtData1_100000 = 'workfile_single_scores_Ne100000'
ExtData1_500000 = 'workfile_single_scores_Ne500000'
ExtData1_1000000 = 'workfile_single_scores_Ne1000000'

set title "Single Lambda Score similar GMYC taxa delimit single minbr 0"
set output 'plots/single_scores_delimit_single_minbr_0.png'
plot ExtData1_10000 using 1:single_0 title 'Ne10000' with linespoints ls 1, ExtData1_100000 using 1:single_0 title 'Ne10000' with linespoints ls 2, ExtData1_500000 using 1:single_0 title 'Ne500000' with linespoints ls 3, ExtData1_1000000 using 1:single_0 title 'Ne1000000' with linespoints ls 4

set title "Single Lambda Score similar GMYC taxa delimit multi minbr 0"
set output 'plots/single_scores_delimit_multi_minbr_0.png'
plot ExtData1_10000 using 1:multi_0 title 'Ne10000' with linespoints ls 1, ExtData1_100000 using 1:multi_0 title 'Ne10000' with linespoints ls 2, ExtData1_500000 using 1:multi_0 title 'Ne500000' with linespoints ls 3, ExtData1_1000000 using 1:multi_0 title 'Ne1000000' with linespoints ls 4

set title "Single Lambda Score similar GMYC taxa delimit single minbr default"
set output 'plots/single_scores_delimit_single_minbr_default.png'
plot ExtData1_10000 using 1:single_default title 'Ne10000' with linespoints ls 1, ExtData1_100000 using 1:single_default title 'Ne10000' with linespoints ls 2, ExtData1_500000 using 1:single_default title 'Ne500000' with linespoints ls 3, ExtData1_1000000 using 1:single_default title 'Ne1000000' with linespoints ls 4

set title "Single Lambda Score similar GMYC taxa delimit multi minbr default"
set output 'plots/single_scores_delimit_multi_minbr_default.png'
plot ExtData1_10000 using 1:multi_default title 'Ne10000' with linespoints ls 1, ExtData1_100000 using 1:multi_default title 'Ne10000' with linespoints ls 2, ExtData1_500000 using 1:multi_default title 'Ne500000' with linespoints ls 3, ExtData1_1000000 using 1:multi_default title 'Ne1000000' with linespoints ls 4

set title "Single Lambda Score similar GMYC taxa PTP minbr default"
set output 'plots/single_scores_PTP_minbr_default.png'
plot ExtData1_10000 using 1:ptp_default title 'Ne10000' with linespoints ls 1, ExtData1_100000 using 1:ptp_default title 'Ne10000' with linespoints ls 2, ExtData1_500000 using 1:ptp_default title 'Ne500000' with linespoints ls 3, ExtData1_1000000 using 1:ptp_default title 'Ne1000000' with linespoints ls 4


set title "Average Multi Lambda Score similar GMYC taxa"
set xlabel "Set number"
set ylabel "Average multi lambda score"

ExtData1 = 'workfile_multi_scores'
set output 'plots/average_multi_scores.png'

plot ExtData1 using 1:2 title 'delimit single minbr 0' with linespoints ls 1, ExtData1 using 1:3 title 'delimit multi minbr 0' with linespoints ls 2, ExtData1 using 1:4 title 'delimit single minbr default' with linespoints ls 3, ExtData1 using 1:5 title 'delimit multi minbr default' with linespoints ls 4, ExtData1 using 1:6 title 'PTP minbr default' with linespoints ls 5, ExtData1 using 1:7 title 'real minbr 0' with linespoints ls 7, ExtData1 using 1:8 title 'real minbr default' with linespoints ls 8





set xlabel "index"
set ylabel "Multi Lambda Score"

ExtData1_10000 = 'workfile_multi_scores_Ne10000'
ExtData1_100000 = 'workfile_multi_scores_Ne100000'
ExtData1_500000 = 'workfile_multi_scores_Ne500000'
ExtData1_1000000 = 'workfile_multi_scores_Ne1000000'

set title "Multi Lambda Score similar GMYC taxa delimit single minbr 0"
set output 'plots/multi_scores_delimit_single_minbr_0.png'
plot ExtData1_10000 using 1:single_0 title 'Ne10000' with linespoints ls 1, ExtData1_100000 using 1:single_0 title 'Ne10000' with linespoints ls 2, ExtData1_500000 using 1:single_0 title 'Ne500000' with linespoints ls 3, ExtData1_1000000 using 1:single_0 title 'Ne1000000' with linespoints ls 4

set title "Multi Lambda Score similar GMYC taxa delimit multi minbr 0"
set output 'plots/multi_scores_delimit_multi_minbr_0.png'
plot ExtData1_10000 using 1:multi_0 title 'Ne10000' with linespoints ls 1, ExtData1_100000 using 1:multi_0 title 'Ne10000' with linespoints ls 2, ExtData1_500000 using 1:multi_0 title 'Ne500000' with linespoints ls 3, ExtData1_1000000 using 1:multi_0 title 'Ne1000000' with linespoints ls 4

set title "Multi Lambda Score similar GMYC taxa delimit single minbr default"
set output 'plots/multi_scores_delimit_single_minbr_default.png'
plot ExtData1_10000 using 1:single_default title 'Ne10000' with linespoints ls 1, ExtData1_100000 using 1:single_default title 'Ne10000' with linespoints ls 2, ExtData1_500000 using 1:single_default title 'Ne500000' with linespoints ls 3, ExtData1_1000000 using 1:single_default title 'Ne1000000' with linespoints ls 4

set title "Multi Lambda Score similar GMYC taxa delimit multi minbr default"
set output 'plots/multi_scores_delimit_multi_minbr_default.png'
plot ExtData1_10000 using 1:multi_default title 'Ne10000' with linespoints ls 1, ExtData1_100000 using 1:multi_default title 'Ne10000' with linespoints ls 2, ExtData1_500000 using 1:multi_default title 'Ne500000' with linespoints ls 3, ExtData1_1000000 using 1:multi_default title 'Ne1000000' with linespoints ls 4

set title "Multi Lambda Score similar GMYC taxa PTP minbr default"
set output 'plots/multi_scores_PTP_minbr_default.png'
plot ExtData1_10000 using 1:ptp_default title 'Ne10000' with linespoints ls 1, ExtData1_100000 using 1:ptp_default title 'Ne10000' with linespoints ls 2, ExtData1_500000 using 1:ptp_default title 'Ne500000' with linespoints ls 3, ExtData1_1000000 using 1:ptp_default title 'Ne1000000' with linespoints ls 4


reset;
