回复 第16楼 的 wanguan2000：原来如此
惭愧，我就是用上面那个例子写的。
</p>
<p>par(cex.axis=0.8)</p>
<p>boxplot(Ct~gene,data=a,boxwex=0.15,at=1:12-0.6,<br />
subset=mygroup==1,col=grey(0.4),<br />
main="gene",xlab="gene",ylab="Ct",xlim=c(0.5,12),ylim=c(18,38),frame.plot=T,axes=F)</p>
<p>boxplot(Ct~gene,data=a,add=T,boxwex=0.15,at=1:12-0.4,<br />
subset=mygroup==2,col=grey(0.6))</p>
<p>boxplot(Ct~gene,data=a,add=T,boxwex=0.15,at=1:12-0.2,<br />
subset=mygroup==3,col=grey(0.8),axes=F)</p>
<p>boxplot(Ct~gene,data=a,add=T,boxwex=0.15,at=1:12,<br />
subset=mygroup==4,col=grey(1),axes=F)</p>
<p>legend("topleft",c("1","2","3","4"),fill=c(grey(0.4),grey(0.6),grey(0.8),grey(1)))</p>
<p>
</p>
[attachment=210774,645]
