Open mbartl13 opened 2 years ago
Hi- I've been using pafr and can generate the plots for each replicate but want to combine all and add error bars in one plot.
SINVd3r1<- read_paf("sinvdpi3_rep1.paf") SINVd3r2<- read_paf("sinvdpi3_rep2.paf") ggplot(SINVd3r1) + geom_point(aes(x = alen, y = dv, color = strand, colour="black"))+geom_smooth(aes(color = strand, x = alen, y = dv), method = "loess", se = FALSE, span = 1/10)+scale_color_brewer(palette="PRGn")+scale_x_continuous("Alignment length (kb)", label = function(x) x/1e3)+scale_y_continuous("Per base divergence") +geom_line(data=SINVd2, aes(x = alen, y = dv, color = strand, colour="black"))+geom_smooth(aes(color = strand, x = alen, y = dv), method = "loess", se = FALSE, span = 1/10)+scale_color_brewer(palette="PRGn")+scale_x_continuous("Alignment length (kb)", label = function(x) x/1e3)+scale_y_continuous("Per base divergence") There were 12 warnings (use warnings() to see them)
SINVd3r1<- read_paf("sinvdpi3_rep1.paf") SINVd3r2<- read_paf("sinvdpi3_rep2.paf") ggplot(SINVd3r1) + geom_point(aes(x = alen, y = dv, color = strand, colour="black"))+geom_smooth(aes(color = strand, x = alen, y = dv), method = "loess", se = FALSE, span = 1/10)+scale_color_brewer(palette="PRGn")+scale_x_continuous("Alignment length (kb)", label = function(x) x/1e3)+scale_y_continuous("Per base divergence") +geom_line(data=SINVd2, aes(x = alen, y = dv, color = strand, colour="black"))+geom_smooth(aes(color = strand, x = alen, y = dv), method = "loess", se = FALSE, span = 1/10)+scale_color_brewer(palette="PRGn")+scale_x_continuous("Alignment length (kb)", label = function(x) x/1e3)+scale_y_continuous("Per base divergence")
Any suggestions for how to combine and plot on one? TY!
Hi- I've been using pafr and can generate the plots for each replicate but want to combine all and add error bars in one plot.
SINVd3r1<- read_paf("sinvdpi3_rep1.paf") SINVd3r2<- read_paf("sinvdpi3_rep2.paf") ggplot(SINVd3r1) + geom_point(aes(x = alen, y = dv, color = strand, colour="black"))+geom_smooth(aes(color = strand, x = alen, y = dv), method = "loess", se = FALSE, span = 1/10)+scale_color_brewer(palette="PRGn")+scale_x_continuous("Alignment length (kb)", label = function(x) x/1e3)+scale_y_continuous("Per base divergence") +geom_line(data=SINVd2, aes(x = alen, y = dv, color = strand, colour="black"))+geom_smooth(aes(color = strand, x = alen, y = dv), method = "loess", se = FALSE, span = 1/10)+scale_color_brewer(palette="PRGn")+scale_x_continuous("Alignment length (kb)", label = function(x) x/1e3)+scale_y_continuous("Per base divergence")There were 12 warnings (use warnings() to see them)Any suggestions for how to combine and plot on one? TY!