2008/06/15

在 R 底下畫 boxplot ,並把他們合併



rm(list=ls())
#read the data
sg01auc <- read.table("data/sg01_auc", header=T)
sg02auc <- read.table("data/sg02_auc", header=T)
sg03auc <- read.table("data/sg03_auc", header=T)
sg04auc <- read.table("data/sg04_auc", header=T)
sg05auc <- read.table("data/sg05_auc", header=T)
sg01kappa <- read.table("data/sg01_kappa", header=T)
sg02kappa <- read.table("data/sg02_kappa", header=T)
sg03kappa <- read.table("data/sg03_kappa", header=T)
sg04kappa <- read.table("data/sg04_kappa", header=T)
sg05kappa <- read.table("data/sg05_kappa", header=T)
postscript(file="boxplotRDo-GAM.ps")
op <- par(mfrow= c(1,2))
boxplot(cbind(sg01auc, sg02auc, sg03auc, sg04auc, sg05auc),
col="lightblue", axes=F, main="AUC")
box()
axis(side=2)
axis(side=1, 1:5, cbind("sg01", "sg02", "sg03", "sg04", "sg05"))
boxplot(cbind(sg01kappa, sg02kappa, sg03kappa, sg04kappa, sg05kappa),
col="lightgreen", axes=F, main="Kappa")
box()
axis(side=4)
axis(side=1, 1:5, cbind("sg01", "sg02", "sg03", "sg04", "sg05"))
dev.off()




boxplotRDo-GAM