基因家族分析套路(四).docVIP

基因家族分析套路（四） 广告：大家有基因组家族分析的需求，麻烦留言哦，小编会第一时间回复的。 今天是基因家族分析类文章最后的一部分，也是一个文章亮点所在的部分，小伙伴们仔细阅读学习吧！ 一、转录组及芯片原始数据下载网站 1、 GEO datesets/profile(/gds ).。用法见下图。GEO数据ID命名规则：GPL-gt;GSE-gt;GSM.GPL: platformGSE: multiple series.GSM: multiple samples.GDS ≈ GSE. Thedifference concentrated on the data labeled GDS can be analyzed for one geneonline. It is simple and easily.The data in the sameGPL can be used to compare inexperiment. 下面是在线分析转录组数据的用法： 2、EBI ArrayExpress(http://www.ebi.ac.uk/arrayexpress/) 该数据库下载数据用法如下： 3、PLEXdb(/).该数据库下载数据用法如下，注意用户名和密码！ 4、SRA db(/sra/) 5、DRA db（http://trace.ddbj.nig.ac.jp/DRASearch/） 二、数据处理 拿到原始数据，要进行处理，才能进行后续数据分析。1、芯片数据。原始数据格式“.cel”格式。以AffyMicroarray数据处理为例讲述主要的命令如下：gt; library(affy); gt;library(makecdfenv); gt;library……gt; barleyGenome = make.cdf.env(“barleyGenome.cdf)gt;mydata lt;- ReadAffy() ##choose “.cel “ file analyzed.gt;eset lt;- rma(mydata);gt;write.exprs(eset,file=mydata.txt)gt;design lt;- model.matrix(~-1 factor(c(1,1,2,2,3,3))) # Createsappropriate design matrix. gt;colnames(design) lt;-c(group1, group2, group3) # Assigns column names.gt;fit lt;- lmFit(eset, design) # Fits a linear model for each gene based onthe given series of arrays.gt;contrast.matrix lt;- makeContrasts(group2-group1,group3-group2, group3-group1, levels=design) # Creates appropriate contrast matrix toperform all pairwise comparisons.gt;fit2 lt;- contrasts.fit(fit, contrast.matrix)# Computes estimatedcoefficients and standard errors for a given set of contrasts.gt;fit2 lt;- eBayes(fit2) # Computes moderated t-statistics and log-oddsof differential expression by empirical Bayes gt;topTable(fit2, coef=1,adjust=fdr, sort.by=B, number=10) # Generates list of top 10 (number=10)differentially expressed genes sorted by B-values (sort.by=B) for firstcomparison group.gt;write.table(topTable(fit2, coef=1,adjust=fdr, sort.by=B, number=500),file=limma_complete.xls, s=F, sep=\t) # Exports complete limma statistics table forfirst comparison group.gt;results lt;- decideTests(fit2,