TCGA_不同gleason评分的 gene 表达

library("stringr")
library("dplyr")
library("ggplot2")
library("readr")

加载数据(将第1列转为行名):

if(T){
  dataFilt=read_csv(file = 'dataFilt.csv') 
  dataFilt=as.data.frame(dataFilt) 
  rownames(dataFilt)=dataFilt[,1] 
  dataFilt=dataFilt[,-1]
  dataFilt=log(dataFilt+1) #取对数
} #测序数据

if(T){
  data_cl=read_csv(file = 'data_cl_index.csv') 
  data_cl=as.data.frame(data_cl) 
  rownames(data_cl)=data_cl[,1] 
  data_cl=data_cl[,-1]
} #临床数据

一、合并临床和测序数据

1 更改dataFilt使其能与data_cl匹配

dataFilt_t=t(dataFilt) #转置

dataFilt_t=as.data.frame(dataFilt_t)
dataFilt_t <- add_rownames(dataFilt_t, "Sample") #add_rownames()使row.names变为第一列(data.frame)
colnames(data_cl)[1]="Sample" 

dataFilt_t$Sample=str_sub(dataFilt_t$Sample,1,12) #取病人编号
</
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

郁柳_Fudan

你的鼓励将是我创作的最大动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付:¥1
获取中
扫码支付

您的余额不足,请更换扫码支付或充值

打赏作者

实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值