【R作图】lattice包，画多个分布柱形图，hist图纵轴转换为百分比

一开始用lattice包，感觉在多元数据的可视化方面，确实做得非常好。各种函数，可以实现任何想要实现的展示。

barchart(y ~ x)	y对x的直方图
bwplot(y ~ x)	盒形图
densityplot(~ x)	密度函数图
dotplot(y ~ x)	Cleveland点图(逐行逐列累加图)
histogram(~ x)	x的频率直方图
qqmath(~ x)	x的关于某理论分布的分位数-分位数图
stripplot(y ~ x)	一维图，x必须是数值型，y可以是因子
qq(y ~ x)	比较两个分布的分位数，x必须是数值型，y可以是数值型，字符型，或者因子，但是必须有两个"水平"
xyplot(y ~ x)	二元图（有许多功能）
levelplot(z ~ x*y) contourplot(z ~ x*y)	在x，y坐标点的z值的彩色等值线图（x，y和z等长）
cloud(z ~ x*y)	3-D透视图（点）
wireframe(z ~ x*y)	同上（面）
splom(~ x)	二维图矩阵
parallel(~ x)	平行坐标图

下面详细介绍下xyplot，xyplot相当于plot，不过可以画多个子图。

一.不过如果想要在每个子图里展示几种点或线，就要想想办法了。

1)这是一个利用公式关系的方法：Sepal.Length + Sepal.Width ~ Petal.Length + Petal.Width能在一个子图中画出

Sepal.Length ~ Petal.Length，Sepal.Width ~ Petal.Length，Sepal.Length ~ Petal.Width，Sepal.Width ~ Petal.Width四种关系

*代码：

library(lattice)

require(stats)

xyplot(Sepal.Length + Sepal.Width ~ Petal.Length + Petal.Width | Species,
data = iris, scales = "free", layout = c(2, 2),
auto.key = list(x = .6, y = .7, corner = c(0, 0)))

2)下面是另一种方法，利用groups参数：

*代码：

xyplot(freq~abun |KO ,groups=country,data =zz, type="o", layout = c(4, 6),xlab="Abundance (log10)",ylab="Frequency",scales="free")

 

不过我想要画hist分布图，所以尝试了几种不同方法：

1）barchart：可惜柱子的透明度无法调节

*代码：

barchart(freq~abun |KO ,groups=country,data =zz,stack =T,strip=T, layout = c(2, 2),xlab="Abundance (log10)",ylab="Frequency",scales="free",col=rainbow(2,alpha=0.5),box.ratio =100,horizontal=F)

2）想用xyplot 函数，把type改成“h”柱形，但是 lwd 这个参数没法每个子图不一样，lend参数无效 。

*代码：

xyplot(freq~abun |KO ,groups=country,data =zz, type="h", lwd = 10, lend =2, layout = c(4, 6),xlab="Abundance (log10)",ylab="Frequency",scales="free",col=rainbow(2,alpha=0.5))

3）经过一番研究，用histogram函以及函数中的panel（子图中加函数）可以实现，但是发觉xlim无法各个子图自己调整，最后无奈选择最原始的方法。

*代码：

histogram(~abun |KO ,data=zz , layout = c(4, 6),xlab="Abundance (log10)",ylab="Frequency",scales="free",col=rainbow(2,alpha=0.5)[1],
panel = function(x, ...) {
panel.histogram(x[1:100],breaks=NULL,type="count",col = rainbow(2,alpha=0.5)[1])
panel.histogram(x[101:160],breaks=NULL,type="count",col = rainbow(2,alpha=0.5)[2])
}
)

4）最原始的方法：

*代码：

par(mfrow=c(6,4),mar=c(2.1,2.1,0.1,0.1))
for (i in 1:length(id))
{
tmp=log10(as.numeric(data[id[i],2:161]))
a=hist(tmp,breaks=10,plot=F)
b=hist(tmp[1:100],breaks=a$breaks,plot=F)
c=hist(tmp[101:160],breaks=a$breaks,plot=F)
if((max(b$count)/100)>(max(c$count)/60)){
hist(tmp[1:100],breaks=a$breaks,col=rainbow(2,alpha=0.5)[1],xlab="",ylab="",main="");
e=hist(tmp[101:160],breaks=a$breaks,plot=F)
e$counts = e$counts*100/sum(e$counts)
plot(e,add=T,col=rainbow(2,alpha=0.5)[2],xlab="",ylab="",main="")
}else{
e=hist(tmp[101:160],breaks=a$breaks,plot=F)
e$counts = e$counts*100/sum(e$counts)
plot(e,col=rainbow(2,alpha=0.5)[2],xlab="",ylab="",main="",border=0)
hist(tmp[1:100],breaks=a$breaks,add=T,col=rainbow(2,alpha=0.5)[1],xlab="",ylab="",main="");
plot(e,add=T,col=rainbow(2,alpha=0.5)[2],xlab="",ylab="",main="")
}
}

最后终于大功告成，里面还有很多地方值得仔细研究。

例如：

　　i.我两组样品的数量不一样，但是我想在一个hist图中比较，怎么办？hist纵轴百分比显示！！！

e=hist(tmp[101:160],breaks=a$breaks,plot=F)#先hist存到变量中
e$counts = e$counts*100/sum(e$counts) #将counts转换成百分比
plot(e,add=T,col=rainbow(2,alpha=0.5)[2],xlab="",ylab="",main="")#再画图

　　ii.我画两个具有一定透明度的hist，发觉画图先后顺序对重叠部分的颜色有影响，例如我先画红色再画蓝色vs.先画蓝色再画红色，结果就不一样了。

但是我又需要让x轴和y轴依据最大的图来确定。

所以我先根据所有数据来取breaks，保持两个hist的break一致，这样x轴问题就解决了；

然后判断哪个hist的y值最大，那么我先根据最大的那个画一个空白的hist，y轴问题解决；

然后再按照先红后蓝的顺序画上两个hist。