草庐IT

单细胞测序并不一定需要harmony去除批次效应

大家好,今天我们分享的是单细胞的学习教程https://www.singlecellworkshop.com/analysis-tutorial.html 教程的作者使用了四个样本,但是没有使用harmony或者其他方法去整合去除批次效应。主要内容:SCTransform流程代码及结果harmony流程代码及结果seurat单样本标准流程代码及结果三种方法结果比较是不是这四个样本就不需要去批次效应呢?接下来我们探索一下1首先是把教程的代码跑一遍#loadSeuratpackagelibrary(Seurat)dir.create("~/gzh/harmony_sct",recursive=TR

用活人脑细胞造AI系统!语音识别已成功,可无监督学习

由真实人脑细胞构建的“迷你大脑”和微电极组成的AI系统,已经能够进行语音识别——从数百个声音片段中准确认出某个特定人的声音的那种。最近,一项颇为前沿的类脑研究登上了Nature子刊。这个特别的AI系统甚至可以进行无监督学习:研究人员只是一遍遍播放音频片段,不提供任何形式的反馈来告诉系统答对还是错。最终,该系统在两天的训练之后,准确率直接从最初的51%升到了78%。这,究竟是怎么实现的?类器官神经网络来了发明该系统的主要目的,是解决硅芯片的高能耗等问题。一般来说,这个问题的解题思路都是靠类脑计算。但这种思想下设计的“传统”类脑芯片大多数都是直接基于数字电子原理,完全模仿大脑功能的能力着实有限。在

如何将可见细胞(过滤的表单元)复制和粘贴到CSV

我只想使用复制和粘贴可见单元格将一个Excel复制到另一个Excel,因为我在设定范围之间有一个过滤表。我想通过保存CSV来做到这一点,但显然这是不可能的。下面的代码有效,但似乎并没有像我想要的那样复制粘贴可见的单元/过滤单元。提前致谢。或者,如果还有另一种推荐方法可以将桌子过滤到CSV,我希望听到如何听到。干杯。SubMacro2()''Macro2Macro''DimlastRowAsLongDimwsAsWorksheet,tblAsListObjectSetws=Sheets("Sheet1")Settbl=ws.ListObjects("Table1")Withtbl.ListCol

细胞实例分割:DoNet: Deep De-overlapping Network for Cytology Instance Segmentation 论文阅读笔记

细胞实例分割:DoNet:DeepDe-overlappingNetworkforCytologyInstanceSegmentation论文阅读笔记一、Abstract二、引言三、相关工作细胞学实例分割遮挡实例分割四、方法4.1预览问题概述工作流程粗糙的Mask分割4.2解耦合和重组策略双路径区域分割模块Dual-pathRegionSegmentationModule(DRM)语义一致性引导的重组模块SemanticConsistency-guidedRecombinationModule(CRM)4.3Mask引导的区域提议Mask-guidedRegionProposal4.4端到端学

生信小白学单细胞转录组(sc-RNA)测序数据分析——R语言

一、数据准备10X单细胞转录组理论上有3个文件才能被读入R进行seurat分析,分别是barcodes.tsv、genes.tsv和matrix.mtx,文件barcodes.tsv和genes.tsv,就是表达矩阵的行名和列名pbmc.data文件解读genes.tsv文件(有时也叫features.tsv文件)文件内容:有两列,第一列为基因ID,第二列为基因SymbolID,区分各个基因。barcodes.tsv文件文件内容:有一列,内容为测序时为了区分各个细胞的标记信息,称为Barcodesmatrix.mtx文件内容:有三列,数字的第一行是测序的汇总信息。第一行的第一个为测序的总基因数

单细胞seurat入门—— 从原始数据到表达矩阵

根据所使用的建库方法,单细胞的RNA序列(也称为读取(reads)或标签(tags))将从转录本的3'端(或5'端)(10XGenomics,CEL-seq2,Drop-seq,inDrops)或全长转录本(Smart-seq)获得。图片来源:PapalexiEandSatijaR.Single-cellRNAsequencingtoexploreimmunecellheterogeneity,NatureReviewsImmunology2018(https://doi.org/10.1038/nri.2017.76)我们可以根据自己感兴趣的生物学问题而选择不同的方法。这些方法具有以下优点:

Python Unet ++ :医学图像分割,医学细胞分割,Unet医学图像处理,语义分割

一,语义分割:分割领域前几年的发展图像分割是机器视觉任务的一个重要基础任务,在图像分析、自动驾驶、视频监控等方面都有很重要的作用。图像分割可以被看成一个分类任务,需要给每个像素进行分类,所以就比图像分类任务更加复杂。此处主要介绍DeepLearning-based相关方法。    主要介绍unet和unet++ 二,数据介绍---医学细胞分割任务原数据:标签数据:   三,代码部分模型包含以下文件:archs.py为模型的主体部分:importtorchfromtorchimportnn__all__=['UNet','NestedUNet']classVGGBlock(nn.Module):

【生物信息学】单细胞RNA测序数据分析:计算亲和力矩阵(基于距离、皮尔逊相关系数)及绘制热图(Heatmap)

文章目录一、实验介绍二、实验环境1.配置虚拟环境2.库版本介绍三、实验内容0.导入必要的库1.读取数据集2.质量控制(可选)3.基于距离的亲和力矩阵4.绘制基因表达的Heatmap5.基于皮尔逊相关系数的亲和力矩阵6.代码整合一、实验介绍  计算亲和力矩阵,一般按照以下步骤进行:导入数据:加载单细胞RNA测序数据集。数据预处理:根据需要对数据进行预处理,例如基因过滤、归一化等。计算亲和力:使用合适的算法(例如,欧几里德距离、Pearson相关系数或其他距离/相似度度量)计算样本之间的亲和力(可以使用现有的生物信息学工具包(如Scanpy)来执行此计算。构建亲和力矩阵:将计算得到的亲和力值组织成

单细胞差异基因可视化|Seurat和scanpy提取细胞簇主要差异基因(DEG)并绘制气泡图

1.背景单细胞数据分析在进行完细胞自聚类或者细胞类型注释后,一般需要对查到的差异基因可视化,用来显示基因和细胞群的相关性,进行后续分析。当然Seurat和scanpy本身可视化的方式有非常多,例如featureplot,violinplot,dotplot等,但是问题在于差异基因分析后,如何快速将每个细胞簇所对应的topdeg汇总,然后再对接函数绘制成图像。Seurat的操作比较简单,因为FindMarker()后自身生成的就是一个数据框,但scanpy的sc.tl.rank_genes_groups()就没有那么用户友好了。2.Seurat的实现library(Seurat)library(

【单细胞】monocle2拟时序分析图美化(2)

前面那个帖子我们讲了如何提取monocle2的结果,然后利用pheatmap自己可以进行多方位的美化。今天我们测试一下如何利用complexheatmap进行更多的美化,因为相对来说complexheatmap能做更多的控制和美化。annotation_col=data.frame( pseudotime=rescale(newdata$Pseudotime,to=c(-1,1)))row.names(annotation_col)annotation_rowrow.names(annotation_row)rowcolornames(rowcolor)ann_colors         C