当前位置: 首页 > news >正文

在线网站代理浏览驻马店市住房和城乡建设局网站

在线网站代理浏览,驻马店市住房和城乡建设局网站,江苏鑫圣建设工程有限公司网站,做外汇那个网站好上节#xff0c;我们使用结核病基因数据#xff0c;做了一个数据预处理的实操案例。例子中结核类型#xff0c;包括结核#xff0c;潜隐进展#xff0c;对照和潜隐#xff0c;四个类别。本节延续上个数据#xff0c;进行了差异分析。 差异分析 计算差异指标step12 加载… 上节我们使用结核病基因数据做了一个数据预处理的实操案例。例子中结核类型包括结核潜隐进展对照和潜隐四个类别。本节延续上个数据进行了差异分析。 差异分析 计算差异指标step12 加载数据 load(dataset_TB_LTBI_step8.Rdata) 构建差异比较矩阵 #样本列表 group_listgroup_data_TB_LTBI$group_more #构建分组 designmodel.matrix(~0factor(group_list)) colnames(design)levels(factor(group_list))#head(dataset_TB_LTBI)row.names(design)colnames(dataset_TB_LTBI) design #得到分组矩阵0代表不是1代表是#str(design)library(limma) ##差异比较矩阵 contrast_matrixmakeContrasts(paste0(c(LTBI,TB),collapse -),levels design) 计算差异基因指标 #step:lmFit fitlmFit(dataset_TB_LTBI,design) fit2contrasts.fit(fit,contrast_matrix) #step:eBayes fit3eBayes(fit2)#step3:topTable tempoutputtopTable(fit3,coef 1,nInf) DEG_Mna.omit(tempoutput) #得到差异分析矩阵重点看logFC和P值 head(DEG_M) #查看数据logFC AveExpr t P.Value adj.P.Val B ASPHD2 -1.452777 8.415563 -12.38370 5.885193e-22 5.868863e-18 39.30255 C1QC -3.978887 5.971935 -12.34993 6.954041e-22 5.868863e-18 39.14037 GBP1P1 -4.075057 5.607978 -12.24397 1.174622e-21 6.608814e-18 38.63087 GBP6 -3.225604 4.393248 -11.93968 5.320543e-21 1.692866e-17 37.16200 SDC3 -2.374911 7.388880 -11.92896 5.612049e-21 1.692866e-17 37.11012 LHFPL2 -1.705514 8.411180 -11.91494 6.017652e-21 1.692866e-17 37.04225#绘制前40个基因在不同样本之间的热图 library(pheatmap) #绘制前40个基因在不同样本之间的热图 f40_genehead(rownames(DEG_M),40) f40_subset_matrixdataset_TB_LTBI[f40_gene,] head(f40_subset_matrix) f40_subset_matrixxt(scale(t(f40_subset_matrix))) #数据标准化。。。数据标准化和归一化的区别平移和压缩 pheatmap(f40_subset_matrixx) #出图 差异分析 结果过滤筛选step13 res DEG[,c(logFC,P.Value,adj.P.Val)]colnames(res)-c(logFC,PValue,padj)colnames(res) library(dplyr) FC_filter 0.585 P_filter0.05 all_diff res %% filter(abs(logFC)FC_filter) %% filter(padjP_filter) res$id rownames(res) resselect(res,id,everything()) #write.table(res,all_diff.txt,sep\t,quoteF)up_diffres %% filter(logFCFC_filter) %% filter(padjP_filter) up_diff$id rownames(up_diff) up_diffselect(up_diff,id,everything()) #write.table(up_diff,up_diff.txt,sep\t,quoteF)down_diffres %% filter(logFC -FC_filter ) %% filter(padjP_filter) down_diff$id rownames(down_diff) down_diffselect(down_diff,id,everything()) #write.table(down_diff,down_diff.txt,sep\t,quoteF)group_data_clean -function(data){# colnames(data)[c(9,10,11)] c(logFC,PValue,padj)data[which(data$padj %in% NA),sig] - no diffdata[which(data$logFC FC_filter data$padj 0.05),sig] - updata[which(data$logFC -FC_filter data$padj 0.05),sig] - downdata[which(abs(data$logFC) FC_filter | data$padj 0.05),sig] - no diffcat( 上调,nrow(data[data$sig %in% up, ]))cat( 下调,nrow(data[data$sig %in% down, ]))cat( no fiff,nrow(data[data$sig %in% no diff, ]))# filter_data subset(data, data$sig up | data$sig down)# filter_data$Geneid - rownames(filter_data)return(data) } limma_clean_res group_data_clean(res)#上调 1381 下调 1432 no fiff 14066rownames(all_diff)dataset_TB_LTBI_DEG dataset_TB_LTBI[rownames(all_diff),] dim(dataset_TB_LTBI_DEG) #[1] 2813 102# #数据保存 # save(DEG,res,all_diff,limma_clean_res,dataset_TB_LTBI_DEG,file DEG_TB_LTBI_step13.Rdata) # #数据保存 #差异分析 绘制火山图step14 library(ggplot2)data - limma_clean_res################# # ggplot2绘制火山图 data$label - c(rownames(data)[1:10],rep(NA,nrow(data) - 10)) #sizeGrWindow(12, 9) pdf(file差异基因火山图step14.pdf, width 9, height 6) ggplot(data,aes(logFC,-log10(PValue),color sig)) xlab(log2FC) geom_point(size 0.6) scale_color_manual(valuesc(#00AFBB,#999999,#FC4E07)) geom_vline(xintercept c(-1,1), linetype dashed) geom_hline(yintercept -log10(0.05), linetype dashed) theme(title element_text(size 15), text element_text(size 15)) theme_classic() geom_text(aes(label label),size 3, vjust 1,hjust -0.1)dev.off()差异基因分析完毕下面我们可以观察一下这些基因富集在哪些通路之上。
http://www.dnsts.com.cn/news/36653.html

相关文章:

  • 网站推广广告 优帮云wordpress3.0手机版
  • 微信公众号申请网站商城网站制作方案
  • 三亚做网站服务中国对外贸易公司排名
  • 网站建设莱州哪家强?校园网站开发的目的
  • 手机网站一般做多大尺寸wordpress中文页面
  • 保亭整站优化青海省住房和城乡建设厅官方网站
  • 网站开发背景绪论软文推广案例500字
  • 网站建设的意义怎么写房地产微网站
  • 网络公司做的网站被告图片侵权网站开发 策划书
  • 做vi设计的网站企业标准网上备案网站
  • 温州seo建站wordpress找不到jquery
  • 门户网站推广方案点击器免费版
  • php导航网站flash 3d 网站源码
  • 手机端网站开发框架网站建设不用虚拟主机
  • 2017年网站推广怎么做网站开发需要几个域名
  • 做奢侈品回收网站特点新手学做网站下载
  • 万维网中文网站到期子域名大全
  • 做简报的网站软件开发兼职网站
  • 如何做网站跳转登入wordpress自建站哪里换logo
  • 太原做app网站建设查找网站后台的软件
  • 网站备案拍照是什么visual studio怎么新建网站
  • 网站建设代码结构做面膜的网站
  • 网站显示系统建设中wix和wordpress区别
  • 网站建设多少钱实惠湘潭磐石网络wordpress教程 吾爱破解
  • 南京网站建设设计新月直播大全免费下载手机版官网
  • 网页设计与网站建设+pdf网站添加微信
  • 网站开发需要什么人才网站备案收录下降
  • 哈尔滨站建好了吗小程序商城需要办什么证
  • 赢卡购网站建设一般通过什么确定险情的严重程度
  • 一半都有哪些做影视外包的网站恐怖音乐怎么做的视频网站