typoq

lib/scripts/plot_hyp_simu_validation.R

@@ -79,7 +79,7 @@ plotD = plotD + facet_grid(hyp ~ tree)
 alpha = 0.7
 x_labs = " # of AA in NonConvLeaves"
 y_labs = "# of sites"
-plotE = ggplot(df, aes(x=NbnonConvAA, fill=hyp)) + theme_bw() + labs(x=x_labs, y=y_labs)
+plotE = ggplot(df, aes(x=NbNonConvAA, fill=hyp)) + theme_bw() + labs(x=x_labs, y=y_labs)
 plotE = plotE + geom_histogram(binwidth = 0.05)
 plotE = plotE + facet_grid(hyp ~ tree)
 

lib/scripts/plot_sens_spe_all_trees.R

@@ -140,11 +140,11 @@ df = do.call(rbind, apply(files_df_ok, 1, read_dir))
 rownames(df) <- NULL
 
 hyps_from = unique(df$hyp)
-hyps_to = rep("Specificity",length(hyps))
+hyps_to = rep("Specificity",length(hyps_from))
 hyps_to[grep("Ha",hyps_from)] = "Sensitivity"
 
-df$variable = mapvalues(df$hyp, from=hyps_from,
-                                  to = hyps_to)
+df$variable = mapvalues(df$hyp, from=hyps_from, to = hyps_to)
+
 df$P[df$variable == "Specificity"] = 1-df$P[df$variable == "Specificity"]