|
Maria@70
|
1 source("MetadataPlots.R")
|
|
Maria@70
|
2
|
|
Maria@70
|
3 PlotCountryOutliers(df=read.csv("data/global_outliers.csv",header=TRUE), output="data/global_outliers.pdf")
|
|
Maria@70
|
4 PlotCountryOutliers(df=read.csv("data/global_outliers_rhy.csv",header=TRUE), output="data/global_outliers_rhy.pdf")
|
|
Maria@70
|
5 PlotCountryOutliers(df=read.csv("data/global_outliers_mel.csv",header=TRUE), output="data/global_outliers_mel.pdf")
|
|
Maria@70
|
6 PlotCountryOutliers(df=read.csv("data/global_outliers_mfc.csv",header=TRUE), output="data/global_outliers_mfc.pdf")
|
|
Maria@70
|
7 PlotCountryOutliers(df=read.csv("data/global_outliers_chr.csv",header=TRUE), output="data/global_outliers_chr.pdf")
|
|
Maria@70
|
8 PlotCountryOutliers(df=read.csv("data/spatial_outliers.csv",header=TRUE), output="data/spatial_outliers.pdf")
|
|
Maria@70
|
9 #PlotCountryOutliers(df=read.csv("data/global_outliers_rhy_1band.csv",header=TRUE))
|
|
Maria@70
|
10
|
|
Maria@70
|
11 require(graphics)
|
|
Maria@70
|
12 par(mfrow=c(2,2))
|
|
Maria@70
|
13 g1<-PlotCountryOutliers(df=read.csv("data/global_outliers_rhy.csv",header=TRUE))
|
|
Maria@70
|
14 g2<-PlotCountryOutliers(df=read.csv("data/global_outliers_mel.csv",header=TRUE))
|
|
Maria@70
|
15 g3<-PlotCountryOutliers(df=read.csv("data/global_outliers_mfc.csv",header=TRUE))
|
|
Maria@70
|
16 g4<-PlotCountryOutliers(df=read.csv("data/global_outliers_chr.csv",header=TRUE))
|
|
Maria@70
|
17 #do.call(addMapLegend, c(g3,labelFontSize=0.7, legendWidth=0.5, tcl=0.3, legendMar = 7, legendLabels="all",horizontal=T, legendIntervals="page"))
|
|
Maria@70
|
18 #legend("bottomleft", legend = c(paste(seq(100,1,-10),'%'), 'missing countries'), fill = c(heat.colors(10, alpha = 1), 'grey'), cex = 0.56, bty = "n")
|
|
Maria@70
|
19 legend("right", legend = c(paste(seq(90,0,-10),'-',seq(100,10,-10),'%'), 'NA'), fill = c(heat.colors(10, alpha = 1), 'grey'), cex = 0.56, bty = "o",bg="white",box.lwd=0,box.col="white")
|
