affymetrix芯片未标准化数据分布的示例分析

这篇文章给大家分享的是有关affymetrix芯片未标准化数据分布的示例分析的内容。小编觉得挺实用的，因此分享给大家做个参考，一起跟随小编过来看看吧。

获取芯片原始数据得到为标准化的数据，需要查看原数据的分布情况，可以通过多种方法获得分布图片。

以affy包ReadyAffy()读取CEL结果：AffyBatch对象为例（GSE66196）

1、箱线图：基于boxplot查看，在加载affy之后，利用boxplot可以直接对AffyBatchd

n.sample = length(AffyBatch$sample)
cols = rainbow(n.sample)
par(cex = 0.7)
if(n.sample>40) par(cex = 0.5)
boxplot(AffyBatch,col = cols,  las = 2)

2、hist曲线，绘图类似boxplot，加载affy之后，可以直接进行绘图

> SampleNames [1] "GSM1616746" "GSM1616747" "GSM1616748" "GSM1616749" "GSM1616750"
 [6] "GSM1616751" "GSM1616752" "GSM1616753" "GSM1616754" "GSM1616755"
[11] "GSM1616756" "GSM1616757"

par(cex=0.7)
if(n.sample>40) par(cex = 0.5)
hist(AffyBatch,col = cols, las = 1,lty=1)
legend("topright",legend =SampleNames,col=cols,lty=1,box.lty=0,bg=NA)

感谢各位的阅读！关于“affymetrix芯片未标准化数据分布的示例分析”这篇文章就分享到这里了，希望以上内容可以对大家有一定的帮助，让大家可以学到更多知识，如果觉得文章不错，可以把它分享出去让更多的人看到吧！