这篇文章给大家分享的是有关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芯片未标准化数据分布的示例分析”这篇文章就分享到这里了,希望以上内容可以对大家有一定的帮助,让大家可以学到更多知识,如果觉得文章不错,可以把它分享出去让更多的人看到吧!
免责声明:本站发布的内容(图片、视频和文字)以原创、转载和分享为主,文章观点不代表本网站立场,如果涉及侵权请联系站长邮箱:is@yisu.com进行举报,并提供相关证据,一经查实,将立刻删除涉嫌侵权内容。