FFPE中的rds和h5ad聚类标签互相添加

百创智造FFPE空间转录组上游分析,我们再聚类的时候选择了三种方法:Seurat聚类,Kmeans聚类,Graph-based聚类。其中seurat聚类得到的数据类型是RDS格式,存放在BSTViewer_project/cluster/level/seurat/L*/object.RDS,kmean和Graph-based聚类得到到数据类型是h5ad格式,存放在BSTViewer_project/cluster/level/kmeans_graph/L*/L*_clustered.h5ad

目前这两个数据(rds,h5ad)都是只存放了各自的聚类,那如果我们想RDS中能添加上kmean和Graph-based聚类标签,或者h5ad中添加上seurat聚类标签,该如何去操作呢?这里我们提供了两个脚本供大家使用。

将kmeans 和graph的聚类标签加到rds脚本 : add_cluster_csv_to_seurat.r

复制代码
nohup Rscript add_cluster_csv_to_seurat.r  \
--out out \ 
--h5ad_cluster_path ./BSTViewer_project/cluster/level/kmeans_graph/L7 \ 
--rds ./BSTViewer_project/cluster/level/seurat/L7/object.RDS \ 
--prefix L7

将seurat聚类加到h5ad :seurat_cluster_to_h5ad.py

复制代码
python seurat_cluster_to_h5ad.py \ 
--out outdir \ 
--h5ad ./BSTViewer_project/cluster/level/kmeans_graph/L18/L18_clustered.h5ad \ 
--spatial_snn_path ./BSTViewer_project/cluster/level/seurat/L18 \ 
--prefix L18

至此,rds和h5ad中都含有这三种方发的聚类标签了(Spatial_snn是seurat聚类标签):

graph_r0.5 graph_r0.8 graph_r1.0 graph_r1.2 graph_r1.5 kmeans_K4 kmeans_K6 ... kmeans_K12 kmeans_K14 kmeans_K16 kmeans_K18 kmeans_K20 Spatial_snn_res0.8 Spatial_snn_res1.5 Spatial_snn_res0.5 Spatial_snn_res1 Spatial_snn_res1.2

脚本下载方式:

通过网盘分享的文件:FFPE中的rds和h5ad聚类标签互相添加

链接: https://pan.baidu.com/s/1Mq_MtTa19XPYGX7nwZE0Og?pwd=eurr 提取码: eurr

相关推荐
奕鼎竜瑆5 天前
无监督学习:聚类与降维
学习·数据挖掘·聚类
袖清暮雨7 天前
机器学习之其他聚类(层次聚类+DBSCAN算法)
人工智能·机器学习·ai·聚类
袖清暮雨7 天前
机器学习之K-means聚类
机器学习·kmeans·聚类
卷毛迷你猪11 天前
快速实验篇(B16)用户级预测可行性审计(否定性意见)
大数据·hadoop·数据挖掘·聚类
ZhangJun9512 天前
Mobike 共享单车分析项目
人工智能·python·算法·kmeans·聚类·knn
卷毛迷你猪14 天前
快速实验篇(B12)异常流量与刷单识别(方法论演示)
大数据·hive·hadoop·数据挖掘·聚类
卷毛迷你猪15 天前
快速实验篇(B11)用户 RFM 与聚类分群实战
机器学习·数据挖掘·聚类
这张生成的图像能检测吗15 天前
(论文速读)AdaCluster:让视频 DiT 的稀疏注意力学会“区别对待”Q 和 K
聚类·注意力机制·视频生成
梦帮科技17 天前
多任务权重流形融合:SLERP 球形线性插值、DARE 稀疏剪枝与多专家模型融合落地
人工智能·深度学习·算法·机器学习·tensorflow·聚类·剪枝