#plot multipanel figure
library(gridBase)
library(grid)
pdf("../mappcabarplot.pdf",width=10,height=10)
layout(matrix(c(1,2,3,3), 2, 3, byrow = TRUE))
## map
plot.new()
vps <- baseViewports()
pushViewport(vps$figure)
vp1 <-plotViewport(c(0.1,0,0,0.1))
print(p1, vp=vp1)
### pca ###
palette(mycols.list)
par(mar=c(4,5,0.5,0.5))
mycols.list<-c("#999999","#999999","#999999","#E69F00","#56B4E9","#999999","#56B4E9","#56B4E9","yellowgreen","#999999","#E69F00","#E69F00","#56B4E9","#56B4E9","#999999","#999999","#999999","#999999","#E69F00","#56B4E9","#E69F00","#56B4E9","#56B4E9")
palette(mycols.list)
plot(g.pca$x[,1], g.pca$x[,2],cex.lab=3, cex.axis=1, col="grey20", pch=21,bg=family, xlab= "", ylab = "", cex=1.5)
legend(-5, 62, legend=c(expression(italic("L. idas")),expression(italic("L. melissa")),"Old hybrids", "New hybrids-DBS"), fill=c("#E69F00", "#56B4E9", "#999999","yellowgreen"), cex=1.3)
title(xlab = "PC1 = 2.4%", cex.lab = 2,line = 3)
title(ylab = "PC2 = 0.98%", cex.lab = 2,line = 3)
### barplot #####
par(mar=c(5,7,3,1))
cols=c(rep("#999999",10),"yellowgreen")
boxplot(hyb$pcscores~hyb$factor, col=cols, ylab="", cex.lab=2, cex.axis=1.5)
title(ylab="PC1 Scores", cex.lab=2, line=3.5)
dev.off()