mysql的seq2_DESeq2处理TCGA数据库Seq-count数据

1、DESeq2需要导入两个数据集：mycounts, colData。先说mycounts，这就是处理完的TCGA数据RNAmatrix.txt，直接读入即可。

library(tidyverse)

library(DESeq2)

#导入数据

setwd("E:/2.Hitseq_counts/")

mycounts

head(mycounts)

#这里有个x，需要去除，先把第一列当作行名来处理

rownames(mycounts)

#把带X的列删除

mycounts

head(mycounts)

2、colData就是对每个样本的一个情况说明。这个可以生成，也可以自己写一个保存为csv格式。我一般自己写。

#载入colData文件

colData

head(colData)

3、构建矩阵

#构建数据矩阵

dds

#查看dds的内容

dds

#接下来，我们要查看treat VS control的总体结果，并根据p-value进行重新排序。利用summary命令统计显示一共多少个genes上调和下调(FDR0.1)

res = results(dds, contrast=c("condition", "control", "treat")) ##或者res= results(dds)

res = res[order(res$pvalue),]

head(res)

summary(res)

4、输出结果

#所有结果先进行输出

write.csv(res,file="All_results.csv")

table(res$padj<0.05)

#获取padj(p值经过多重校验校正后的值)小于0.05，表达倍数取以2为对数后大于1或者小于-1的差异表达基因。代码如下

diff_gene_deseq2 1.5) ##或者diff_gene_deseq2 1 | log2FoldChange < -1))

dim(diff_gene_deseq2)

head(diff_gene_deseq2)

write.csv(diff_gene_deseq2,file= "DEG_treat_vs_control.csv")

mysql的seq2_DESeq2处理TCGA数据库Seq-count数据相关推荐

linux下载TCGA数据,TCGA数据库讲解与数据下载
原标题:TCGA数据库讲解与数据下载前几期,生信草堂和大家分享了一篇TCGA数据下载的文章,有小白表示TCGA还是好深奥,找不到下载地址的,看不懂数据格式.今天小编就从头开始,娓娓道来,再给大家介绍 ...
mysql 字段操作_Mysql：数据库操作、数据表操作、字段操作整理
mysql 登录本机 mysql -u root -p password 数据库操作创建数据库命令:create database 例如:mysql> create database xhk ...
mysql表的组成_数据库中的数据表由什么组成？
在数据库中,数据表是由表名.表中的字段和表的记录三个部分组成的.在建立表之前都必须先设计它的结构,表结构描述了一个表的框架.设计表结构实际上就是定义组成一个表的字段个数,每个字段的名称.数据类型和长度 ...
jmeter mysql查询结果提取_Jmeter-从数据库中获取数据并作为变量使用
一.导入数据库驱动包: (1).在网上搜索下载数据库驱动包(MySQL的地址https://dev.mysql.com/downloads/file/?id=480090) (2).解压,获取到jar ...
mysql migrations_Code First Migrations更新数据库结构（数据迁移）
背景 code first起初当修改model后,要持久化至数据库中时,总要把原数据库给删除掉再创建(DropCreateDatabaseIfModelChanges),此时就会产生一个问题,当我们的 ...
MySQL官方提供的测试数据库脚本和数据文件下载
网盘下载链接: [sakila database]链接: https://pan.baidu.com/s/1qo0ybhid4DpyVVrhtxjwCQ 提取码: 5dqb [world.sql]链接 ...
mysql select详解_MySQL数据库 select 数据查询详解
数据查询数据查询是数据库系统应用的主要内容,也是用户对数据库最频繁.最常见的基本操作请求.数据查询可以根据用户提供的限定条件,从已存在的数据表中检索用户需要的数据.MySQL使用SELECT语句从数 ...
tcga数据下载_给你tcga数据库过万病人的原始测序数据你可以做什么
最近有两年前的学生过来寻求合作,让我想想给我tcga数据库过万病人的原始测序数据,我可以做什么方法学的创新.我想把这个问题抛给粉丝: 假设给你tcga数据库过万病人的原始测序数据你可以做什么??? 大 ...
java mysql 文本导入数据语句_Java利用MYSQL LOAD DATA LOCAL INFILE实现大批量导入数据到MySQL...
Mysql load data的使用数据库中,最常见的写入数据方式是通过SQL INSERT来写入,另外就是通过备份文件恢复数据库,这种备份文件在MySQL中是SQL脚本,实际上执行的还是在批量IN ...

mysql的seq2_DESeq2处理TCGA数据库Seq-count数据

mysql的seq2_DESeq2处理TCGA数据库Seq-count数据相关推荐

最新文章

热门文章