保存及加載數(shù)據(jù):

save(a,file="GSE17215.Rdata")
load("GSE17215.Rdata")

清除內(nèi)存變量况鸣、釋放內(nèi)存：

rm(list = ls())  #清除內(nèi)存變量
gc()  #釋放內(nèi)存

畫一個(gè)熱圖：

#參考http://www.reibang.com/p/4e3a94bd155a 第五部分
a<-read.table("GSE17215_series_matrix.txt",comment.char = "!",header = T,sep = "\t")
write.csv(a,"a1.csv")

rownames(a)<-a[,1]
a<-a[-1]

a<-log2(a)
library(pheatmap)
pheatmap(a)

OncoLnc網(wǎng)站 獲取某基因在TCGA數(shù)據(jù)庫中的表達(dá)情況

繪制生存曲線：

#從網(wǎng)頁(例如OncoLnc)獲取數(shù)據(jù)民镜，繪制生存曲線
a<-read.csv("LIHC_93663_50_50.csv",header = T,sep = ",",fill = T)
colnames(a)
dat<-a
library(ggstatsplot)
ggbetweenstats(data=dat,x=Group,y=Expression)
ggbetweenstats(data=dat,x=Status,y=Expression)
library(ggplot2)
library(survival)
library(survminer)
table(dat$Status)
dat$Status<-ifelse(dat$Status=="Dead",1,0)
sfit<-survfit(Surv(Days,Status)~Group,data=dat)
sfit
summary(sfit)
ggsurvplot(sfit,conf.int=F,pval=T)
ggsave("survival_ARHGAP18_in_LGG.png")

ggsurvplot(sfit,palette = c("#E7B800", "#2E9FDF"),
           risk.table=T,pval=T,conf.int=T,
           xlab="Time in months",ggtheme = theme_light(),ncensor.plot=T)
ggsave("survival_ARHGAP18_in_LGG_2.png")

cBioPortal 網(wǎng)站獲取臨床信息

表達(dá)量與臨床分期畫圖：

a<-read.table("plot.txt",header = T,sep="\t",fill = T)
colnames(a)<-c("id","stage","gene","mut")
dat<-a
library(ggstatsplot)
ggbetweenstats(data = dat,x = stage,y = gene)
ggsave("ov_ARHGAP18.png")

#改變統(tǒng)計(jì)方法
res.aov <- aov(gene~stage,data=dat)
summary(res.aov)
TukeyHSD(res.aov)

相關(guān)性cor(a,b) --->a,b之間的相關(guān)性,cor(exprSet)可以查看exprSet矩陣所有列之間的相關(guān)性；

dim()`可以看某數(shù)據(jù)的維度帆锋。

表達(dá)矩陣樣本的相關(guān)性：

rm(list = ls())
options(stringsAsFactors = F)
library(airway)
data("airway")
exprSet<-assay(airway)  #h獲得表達(dá)矩陣

colnames(exprSet)
cor(exprSet)
pheatmap::pheatmap(col(exprSet))

group_list<-colData(airway)[,3]
tmp<-data.frame(g<-group_list)
rownames(tmp)<-colnames(exprSet)
pheatmap::pheatmap(col(exprSet),annotation_col = tmp)
dim(exprSet)
exprSet<-exprSet[apply(exprSet,1,function(x) sum(x>1) >5),]
dim(exprSet)

exprSet<-log(edgeR::cpm(exprSet)+1)   #去除文庫大小差異，并得到新矩陣
exprSet<-exprSet[names(sort(apply(exprSet,1, mad),decreasing = T)[1:500]),]
dim(exprSet)
M<-cor(log2(exprSet+1))
tmp<-data.frame(g<-group_list)
rownames(tmp)<-colnames(M)
pheatmap::pheatmap(M,annotation_col = tmp,filename = "cor.png")