补充材料 supplementary
补充 材料
2023-09-14 09:09:48 时间
load("G:/silicosis/geo/GSE104154_scRNA-seq_fibrotic MC_bleomycin/normalized/only_AMs.rds")
markers=FindAllMarkers(subset_data,only.pos = T)
#openxlsx::write.xlsx(markers,file = "G:/silicosis/geo/GSE104154_scRNA-seq_fibrotic MC_bleomycin/normalized/makers_for_AM1_AM2_AM3.xlsx")
markers=openxlsx::read.xlsx("G:/silicosis/geo/GSE104154_scRNA-seq_fibrotic MC_bleomycin/normalized/makers_for_AM1_AM2_AM3.xlsx")
head(markers)
mymarkers=list()
for (eachcluster_markers in unique(markers$cluster)) {
#eachcluster_markers="AM1"
mymarkers[[paste0(eachcluster_markers)]]=markers[markers$cluster==eachcluster_markers,"gene"]
}
head(mymarkers)
fibrosis_gene=openxlsx::read.xlsx("G:/silicosis/fibrosis_geness/harmonizome_gfibrosis.xlsx")
head(fibrosis_gene)
fibrosis_gene=fibrosis_gene$Symbol
#intersect_fibrosis_gene_with_AM
intersected_fibrosis_genes=list()
for (eachcluster_AM in names(mymarkers)) {
intersected_fibrosis_genes[[eachcluster_AM]]=intersect(mymarkers[[eachcluster_AM]],fibrosis_gene)
}
intersected_fibrosis_genes
lapply(intersected_fibrosis_genes,length)
intersected_fibrosis_genes_as_a_vector=as.vector(unlist(intersected_fibrosis_genes))
marker_exp=AverageExpression(subset_data, features = intersected_fibrosis_genes_as_a_vector, return.seurat = TRUE)
DoHeatmap(marker_exp, features = intersected_fibrosis_genes_as_a_vector,
label=TRUE, group.bar = TRUE, draw.lines = FALSE)
getwd()
pdf('heatmap.pdf',width = 7,height = 14)
p=DoHeatmap(marker_exp, features = intersected_fibrosis_genes_as_a_vector,
label=TRUE, group.bar = TRUE, draw.lines = FALSE)
print(p)
dev.off()
#supplementary table S7
PF_silicosis_bleo=openxlsx::read.xlsx("G:\\silicosis\\sicosis\\Supplementary Table.xlsx",
sheet = "S7_hPF_Silica_Bleomycin")
head(PF_silicosis_bleo)
library(Hmisc)
library(dplyr)
disease_markers_for_AM3=list()
for (eachdisease in colnames(PF_silicosis_bleo)) {
#eachdisease="Silicosis,.Cluster.AM3"
disease_markers_for_AM3[[paste0(eachdisease)]]=PF_silicosis_bleo[,eachdisease] %>%
capitalize()%>% tolower() %>% na.omit() %>% unique() %>%capitalize()
}
head(disease_markers_for_AM3)
#https://www.jianshu.com/p/58f5429b5402
genes <- Reduce(intersect, list(disease_markers_for_AM3[[1]],
disease_markers_for_AM3[[2]],
disease_markers_for_AM3[[3]])) %>% capitalize()
length(genes)
getwd()
library(VennDiagram)#https://www.jianshu.com/p/b5a4c40c3a33
venn_list <- list(Human_PF = disease_markers_for_AM3[[1]],
Silicosis = disease_markers_for_AM3[[2]],
Bleomycin=disease_markers_for_AM3[[3]])
venn.diagram(venn_list, filename = 'venn2_for_humanPF-Silicosi-Bleomycin.png', imagetype = 'png',
fill = c('red', 'blue','orange'), alpha = 0.50, cat.col = rep('black', 3),
col = 'black', cex = 1.3, fontfamily = 'serif',
cat.cex = 1.3, cat.fontfamily = 'serif')
getwd()
#继续以上述4个分组为例,组间交集元素获得
inter <- get.venn.partitions(venn_list)
for (i in 1:nrow(inter)) inter[i,'values'] <- paste(inter[[i,'..values..']], collapse = ', ')
write.table(inter[-c(5, 6)], 'venn4_inter.txt', row.names = FALSE, sep = '\t', quote = FALSE)
for (i in 1:nrow(inter)) inter[i,'values'] <- paste(inter[[i,'..values..']], collapse = ', ')
openxlsx::write.xlsx(inter[-c(5, 6)], 'venn4_inter.xlsx', row.names = FALSE, sep = '\t', quote = FALSE)
相关文章
- 对sssp项目搭建的补充,总错误处理。
- html5知识点补充—mark元素的使用
- flask-本地线程-请求上下文补充
- Django-form组件补充
- dorado spring知识补充
- jQuery补充,基于jQuery的bxslider轮播器插件
- 【Python 补充01】Python运算符
- Atitit 学习方法 补充 艾龙 著 attilax著 1. Atitit 学习的方法 attilax总结1 1.1. 2. 基于学习策略的分类2机械 示教 演绎 类比 解释 归纳2 1.
- paip.自适应网页设计 跟 响应式 设计方法与工具补充(2)o54
- 习题 11.3 将例11.3 的程序修改、补充,写成一个完整、正确的程序,用保护继承方式。在程序中应包括输入数据的函数。
- 集团信息化之路—补充在线销售和支付平台解决方案
- Appium自动化测试基础--补充:C/S架构和B/S架构说明
- DNS隧道工具汇总——补充,还有IP over DNS的工具NSTX、Iodine、DNSCat
- FPGA | 深入理解复位(补充)
- 00-ES6补充
- 〖Python语法进阶篇㉑〗- 综合实战 - 抽奖系统之优化补充