R-tips
- rpkm->tpm
CCLE_TPM_ERV <- apply(CCLE_RPKM_ERV1, 2, function(x) x/sum(x))*1000000
- 读gct文件格式
rt<-read.table("GTEx_Analysis_2016-01-15_v7_RNASeQCv1.1.8_gene_tpm.gct", skip = 2, header = TRUE, sep = "\t",nrow=5)
- 找df中X1 列的值位于{1,2}范围内,X2 列的值位于{2,5,6}范围
data[which(((dataX1==1)|(dataX1==2))&
((dataX5==2)|(dataX2==5)| (dataX2==6))),]
a <- c(1,2)
b <- c(2,5,6)
data[data$X1 %in% a & data$X2 %in% b,]
- R语言如何去除矩阵中相同的列
df[,!duplicated(t(df))]
t(unique(t(df)) #转成行~
df[!duplicated(df),] #去除相同的行
- 用表达式调用变量