Cell发表的单细胞整合方法:LIGER,很好用!

之前,我们已经介绍过其他的整合方法:

  1. harmony整合单细胞数据-去除批次效应

  2. Seurat一键完成五种数据整合:harmony,CCA,RPCA,FastMNN,scVI,代码分享

今天来看看Cell发表的LIGER

介绍

LIGER(Linked Inference of Genomic Experimental Relationships,基因实验关系的链接推断)是一个用于整合和分析多个单细胞数据集的软件包,由Macosko实验室开发,并由Welch实验室维护和扩展。它依赖于整合的非负矩阵分解技术来识别共享的和数据集特定的因子。

LIGER可以用来在多种背景下比较和对比实验数据集,例如:

  • 实验批次间

  • 不同个体间

  • 不同性别间

  • 不同组织间

  • 不同物种间(例如,小鼠和人类)

  • 不同模态间(例如,单细胞RNA测序和空间转录组数据,单细胞甲基化或单细胞ATAC测序)

一旦多个数据集被整合,该软件包提供了进一步的数据探索、分析和可视化的功能。用户可以:

  1. 识别cluster

  2. 找到显著的共享(和数据集特定的)基因标记

  3. 与之前识别的细胞类型比较簇

  4. 使用t-SNE和UMAP可视化簇和基因表达

让人兴奋的是,作者设计了LIGER r包,使其能够与现有的单细胞分析软件包(包括Seurat)接入。


代码实战

这里的数据仅作演示使用,读者使用自己的数据可以试一试

如果你想使用我的pbmc:

  1. 20231126单细胞直播一理解pbmc处理流程与常用R技巧

  2. 单细胞数据pbmc数据下载,标准处理流程获取seurat对象

#remotes::install_github('satijalab/seurat-wrappers').libPaths(  c(    '/home/rootyll/seurat_v5/',    "/usr/local/lib/R/site-library",    "/usr/lib/R/site-library",    "/usr/lib/R/library"  ))#install.packages('rliger')library(rliger)library(Seurat)library(SeuratData)library(Seurat)library(SeuratWrappers)#下面正式开始-----------------------------------------------------# options(timeout = 9000)# #InstallData("pbmcsca")# data("pbmcsca")load("~/gzh/pbmc3k_final_v4.rds")pbmc$Method=pbmc$grouppbmcsca=pbmc# Please update your `liger` version to 0.5.0 or above before following this tutorialpbmcsca <- NormalizeData(pbmcsca)pbmcsca <- FindVariableFeatures(pbmcsca)pbmcsca <- ScaleData(pbmcsca, split.by = "Method", do.center = FALSE)#整合的过程就至需要下面这两句代码----pbmcsca <- RunOptimizeALS(pbmcsca, k = 20, lambda = 5, split.by = "Method")pbmcsca <- RunQuantileNorm(pbmcsca, split.by = "Method")# You can optionally perform Louvain clustering (`FindNeighbors` and `FindClusters`) after# `RunQuantileNorm` according to your needspbmcsca <- FindNeighbors(pbmcsca, reduction = "iNMF", dims = 1:20)pbmcsca <- FindClusters(pbmcsca, resolution = 0.3)# Dimensional reduction and plottingpbmcsca <- RunUMAP(pbmcsca, dims = 1:ncol(pbmcsca[["iNMF"]]), reduction = "iNMF")DimPlot(pbmcsca, group.by = c("Method", "ident", "cell.type"), ncol = 3)head(pbmc@meta.data)

LIGER还具有其他功能,后续我们再分享~

  • Iterative Single-Cell Multi-Omic Integration Using Online iNMF

  • Integrating unshared features with UINMF

  • scATAC and scRNA Integration using unshared features (UINMF)

  • Cross-species Analysis with UINMF

  • Performing Parameter Selection

  • Integrating spatial transcriptomic and transcriptomic datasets using UINMF (Click to Download)

  • Integrating Multiple Single-Cell RNA-seq Datasets

  • Jointly Defining Cell Types from scRNA-seq and scATAC-seq

  • Jointly Defining Cell Types from Single-Cell RNA-seq and DNA Methylation

  • Running Liger directly on Seurat objects using Seurat wrappers

生信小博士

【生物信息学】R语言开始,学习生信。Seurat,单细胞测序,空间转录组。 Python,scanpy,cell2location。资料分享

公众号

参考;https://github.com/welch-lab/ligerhttps://github.com/satijalab/seurat-wrappershttps://htmlpreview.github.io/?https://github.com/satijalab/seurat.wrappers/blob/master/docs/liger.html

看完记得顺手点个“在看”哦!

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.mzph.cn/news/761928.shtml

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

C++语言学习(二)—— C++语言的基本知识

目录 一、面向对象的三个核心概念 二、C语言中的I/O口 三、C语言中的数据类型​​​​​​​ 3.1 逻辑类型 3.2 引用类型 3.2.1 引用作为函数参数 3.2.2 引用作为函数返回值 3.2.3 引用作为类成员 3.3 类类型 四、 C语言中的内联函数 五、 函数重载 六、 带默认形参…

docker 容器挂掉,无法exec 进入bash 怎么修改容器里的文件

在使用tdengine 数据库时出现了 TDengine.Driver.TDengineError:“code:[0x334],error:Out of dnodes” 查找文档发现需要修改一个配置文件 。 /etc/taos/taos.cfg 中的 supportVnodes 参数 于是修改 保存。然后&#xff0c;运行出错。 03/21 06:56:27.986498 00000064 …

IDEA通过内网穿透实现固定公网地址远程SSH连接本地Linux服务器

文章目录 1. 检查Linux SSH服务2. 本地连接测试3. Linux 安装Cpolar4. 创建远程连接公网地址5. 公网远程连接测试6. 固定连接公网地址7. 固定地址连接测试 本文主要介绍如何在IDEA中设置远程连接服务器开发环境&#xff0c;并结合Cpolar内网穿透工具实现无公网远程连接&#xf…

[音视频学习笔记]六、自制音视频播放器Part1 -新版本ffmpeg,Qt +VS2022,都什么年代了还在写传统播放器?

前言 参考了雷神的自制播放器项目&#xff0c;100行代码实现最简单的基于FFMPEGSDL的视频播放器&#xff08;SDL1.x&#xff09; 不过老版本的代码参考意义不大了&#xff0c;我现在准备使用Qt VS2022 FFmpeg59重写这部分代码&#xff0c;具体的代码仓库如下&#xff1a; …

基于python+vue渔船出海及海货统计系统的设计与实现flask-django-php-nodejs

当今社会已经步入了科学技术进步和经济社会快速发展的新时期&#xff0c;国际信息和学术交流也不断加强&#xff0c;计算机技术对经济社会发展和人民生活改善的影响也日益突出&#xff0c;人类的生存和思考方式也产生了变化。传统渔船出海及海货统计采取了人工的管理方法&#…

抖音视频爬虫下载软件|可导出视频分享链接|视频批量采集工具

无水印视频关键词批量下载工具操作指南 我们自主开发了一款便捷的音视频批量下载工具&#xff0c;不仅支持单个视频链接提取&#xff0c;还可通过关键词搜索实现批量采集和选择性下载&#xff0c;让您轻松获取所需视频内容。 操作说明&#xff1a; 关键词批量采集视频 进入软件…

AI绘画可以稳定生成中文了:白嫖阿里云部署AnyText

长久以来&#xff0c;在AI绘画中书写文字一直是个难题。即使到了SDXL时代&#xff0c;我们也只能输出英文&#xff0c;而且还经常出现漏掉字母的情况。现在阿里达摩院搞出了一个解决方案&#xff0c;可以在Stable Diffusion生成的作品中稳定输出中、英、日、韩等多种文字&#…

python共享单车信息系统的设计与实现flask-django-php-nodejs

课题主要分为二大模块&#xff1a;即管理员模块和用户模块&#xff0c;主要功能包括&#xff1a;用户、区域、共享单车、单车租赁、租赁归还、报修信息、检修信息等&#xff1b; 语言&#xff1a;Python 框架&#xff1a;django/flask 软件版本&#xff1a;python3.7.7 数据库…

探索Code Llama 70B:Meta让AI辅助编程更易获得的举措

探索Code Llama 70B&#xff1a;Meta让AI辅助编程更易获取 在AI技术不断改变软件开发的今天&#xff0c;Meta推出了其最先进的开源基础模型&#xff0c;简化了软件开发流程。这个模型名为Code Llama 70B&#xff0c;旨在让AI辅助代码生成及其相关任务更容易被更广泛的受众获取…

【Unity】从0到1的横版2d制作笔记-DAY3

确定碰撞体积 选择rigidbody2d&#xff0c;创建player重力 创建player碰撞体积 创建瓦片地图碰撞体积 使平台变成一个整体 ​​​​​ 设置Body Type为Static&#xff08;避免平台也因为重力影响下落&#xff09; 回到Player&#xff0c;在Rigidbody2D中设置为冻结旋转 Player设…

python家政服务系统flask-django-php-nodejs

相比于以前的传统手工管理方式&#xff0c;智能化的管理方式可以大幅降低家政公司的运营人员成本&#xff0c;实现了家政服务的标准化、制度化、程序化的管理&#xff0c;有效地防止了家政服务的随意管理&#xff0c;提高了信息的处理速度和精确度&#xff0c;能够及时、准确地…

解决mini2440 LCD(型号:P43)驱动的背光失效问题

目录 概述 1 LCD(P43)背光问题 1 移植LCD驱动程序 1.1 编写驱动代码 1.2 编写LCD的配置参数 1.3 配置LCD驱动至内核 2 移植触摸功能程序 2.1 移植ADC驱动程序 2.1.1 编写ADC驱动代码 2.1.2 编写配置文件 2.2 移植P43的一线触摸程序 2.2.1 编写代码 2.2.2 编写配置链…

-bash: ./1.sh: /bin/bash^M: bad interpreter: No such file or directory解决方法

1、执行脚本 ./1.sh时报如下错误 -bash: ./1.sh: /bin/bash^M: bad interpreter: No such file or directory 2、在Windows编辑的脚本导入Linux系统中&#xff0c;执行报错问题 yum install -y dos2unix 3、或者本地安装 rpm -ivh /mnt/Packages/dos...... 4、然…

Flutter探索之旅:控制键盘可见性的神奇工具(flutter_keyboard_visibility)

随着移动应用的不断发展&#xff0c;用户体验的重要性愈发突显。而键盘的弹出和隐藏对于用户体验来说是至关重要的一环。在Flutter中&#xff0c;我们有幸拥有一个强大的工具——flutter_keyboard_visibility&#xff0c;它让我们能够轻松地监测键盘的可见性并做出相应的响应。…

设计模式(行为型设计模式——备忘录模式)

设计模式&#xff08;行为型设计模式——备忘录模式&#xff09; 备忘录模式 基本定义 在不破坏封装的前提下&#xff0c;捕获一个对象的内部状态&#xff0c;并在该对象之外保存这个状态&#xff0c;这样可以在以后将对象恢复到原先保存的状态。 模式结构 Originator&#x…

3527. 旋转矩阵 考研上机复试真题 暴力模拟

任意输入两个 9 阶以下矩阵&#xff0c;要求判断第二个是否是第一个的旋转矩阵&#xff08;顺时针&#xff09;&#xff0c;如果是&#xff0c;输出旋转角度&#xff08;0、90、180、270&#xff09;&#xff0c;如果不是&#xff0c;输出 −1。 输入格式 第一行包含整数 n&a…

spring boot @Autowired 为null 处理

spring boot 开发中&#xff0c;经常用到 Autowired 注入某些服务&#xff0c;但有时发现注入的服务为null。原因很简单&#xff0c;因为服务生成是有先后顺序的。当要注入的服务还没有生成时&#xff0c;自然为null。 解决的方法&#xff1a;通过getBean 手工获取服务。 1、定…

基于python+vue中医学习服务管理系统flask-django-php-nodejs

随着世界经济信息化、全球化的到来和互联网的飞速发展&#xff0c;推动了各行业的改革。若想达到安全&#xff0c;快捷的目的&#xff0c;就需要拥有信息化的组织和管理模式&#xff0c;建立一套合理、动态的、交互友好的、高效的中医学习服务管理系统。当前的信息管理存在工作…

利用Jmeter工具对服务器,数据库进行性能监控,压测,导出性能测试报告

Jmeter是Apache基金会旗下的一款免费,开源,轻量级的性能测试工具,主要针对web应用程序客户端/服务器进行性能测试.它可以分别测试静态、动态资源(Java Servlet,CGI Scripts,Java Object,数据库和FTP服务器等),它可以通过线程组来模拟数个用户,在一段时间内同时登录服务器,数个用…

Matlab使用教程(持续更新)

1. Matlab Matlab被广泛的应用在数据分析&#xff0c;汽车仿真&#xff0c;机器人以及医学研究等众多方面。 它可以帮助我们理解研究复杂的系统。 在60年代和70年代&#xff0c;计算机使得科学家和工程师完成了以前不可能进行的计算&#xff1b;但是需要懂得计算机编程。 C…