10x单细胞数据处理各种报错记录

bio_yys

已于 2023-02-02 15:20:14 修改

阅读量764

点赞数

分类专栏：生物信息文章标签： linux

于 2023-02-02 14:26:50 首次发布

本文链接：https://blog.csdn.net/yangyongsan/article/details/128844385

版权

生物信息专栏收录该内容

4 篇文章 1 订阅

订阅专栏

1.SRR数据解压

SRR数据有时候解压出来就一个fastq文件，但是在geo上看原始数据是提供了I1,R1,R2文件的

解压命令如下：

fasterq-dump --split-files  --threads 128 --outdir /home/yangys/data/healthy_pbmc/fastq/SRR8712357 /home/yangys/data/healthy_pbmc/fastq/SRR8712357

解决办法：换解压命令

/home/yangys/pfastq-dump/bin/pfastq-dump --threads 96 --split-files --outdir /home/yangys/data/healthy_pbmc/fastq/SRR8712357 /home/yangys/data/healthy_pbmc/SRR/SRR8712357.sra

有时候--split-files参数也会有影响，可以和--split-3换着测试一下

2.cellranger count报错

解决办法：加上--chemistry=***参数，具体得看数据介绍去确定参数

3.GEO上提供的Original format是bam文件

解决办法：下载bam文件然后用cellranger bamtofastq命令解压

cellranger bamtofastq --nthreads 96 --traceback /home/yangys/data/healthy_pbmc/SRR/SRR17296837.bam /home/yangys/data/healthy_pbmc/fastq/test

优惠劵

bio_yys

关注关注

0
点赞
踩
2

收藏

觉得还不错? 一键收藏
1
评论
10x单细胞数据处理各种报错记录

10x单细胞数据处理各种报错记录
复制链接

扫一扫

专栏目录

cellranger配置环境变量失败怎么办？

Zmonster_ion的博客

09-30

2022

我在下载并解压完cellranger-6.1.1后配置环境变量 export PATH=/home/XXX/XXXX/cellranger-6.1.1:$PATH 根据官网的教程Installing Cell Ranger -Software -Single Cell Gene Expression -Official 10x Genomics Support，下一步应该直接输入 cellranger 以查看是否配置成功，但是报错提示 [User@localhost bin]$ cellra

单细胞数据资源之pbmc3k

最新发布

03-05

该资源为Seurat的单细胞基础分析的基础数据集：pbmc3k_filtered_gene_bc_matrices data/filtered_gene_bc_matrices/hg19

1 条评论您还未登录，请先登录后发表或查看评论

2024单细胞数据分析之拟时序分析monocle2

02-28

monocle数据

BD流式细胞仪脱机数据处理系统的建立.pdf

08-15

#资源达人分享计划#

sctools:用于大规模分布式单细胞数据处理的实用程序

05-22

单细胞工具 Single Cell Tools提供用于处理序列数据格式的实用程序，适用于分析大型生物数据集的分布式系统。下载及安装凳子包 sctools软件包提供了命令行实用程序和设计用于python程序的类。命令行实用程序 Attach10XBarcodes：附加的条形码存储在fastq文件中，以读取未对齐的bam文件 SplitBam：将bam文件拆分为多个块，确保将单元格包含在1个块中 CalculateGeneMetrics：计算有关实验或组块中基因的信息 CalculateCellMetrics：计算有关实验或块中细胞的信息 MergeGeneMetrics：合并从实验的不同部分计算得出的基因指标 MergeCellMetrics合并根据实验的不同部分计算得出的单元格指标主要包装类别平台：一个抽象类，为不同的3'测序格式定义了通用的数据结构。此程序包中所有旨在处理

Linux系统中已知SRR号如何从NCBI上下载SRA数据到服务器中

weixin_52660253的博客

08-01

1579

高通量组学数据目前已经成为生物研究的重要板块，对于一些文章中出现的数据的挖掘尤其是人体数据的再利用也成为探究科学问题的重要前沿组成。通常情况下文章的高通量数据需要上传到NCBI的SRA（Sequence Read Archive）供大家下载学习，而我们也可以通过多种方法对数据进行下载再挖掘。如果仍然出现问题，请再次检查PATH设置是否正确，并确保在PATH中包含了sra-toolkit的bin目录。首先，我们在下面的网址中进入SRA Toolkit下载的官网。我们在NCBI上找到我们需要的数据的。

count

weixin_63884580的博客

11-25

558

cellranger count

Cellranger故障排除

clayluo的专栏

12-15

903

cellranger

cellranger处理单细胞数据的一些错误记录

zengwanqin的博客

12-15

554

2、确认化学版本：确保您使用的是与您的测序实验相匹配的正确的化学版本。错误消息指出，FASTQ 文件 “T_C_5_S32_L4_R1_001.fastq.gz” 在第 0 行出现问题，可能是由于文件格式不正确或未被正确压缩。您可以使用 zcat T_C_5_S32_L4_R1_001.fastq.gz | head 命令来检查文件的前几行，确认它们是否以 ‘@’ 开头。4、重新下载或生成 FASTQ 文件：如果文件确实损坏或格式错误，您可能需要重新从您的测序平台获取或重新生成 FASTQ 文件。

基于10X平台的scRNA-seq数据准备

yangyongsan的博客

09-21

1031

基于10X平台的scRNA-seq数据准备

单细胞常规巴氏涂片图像数据集

12-09

单细胞常规巴氏涂片图像数据集，数据库由4049个分离细胞图像组成，这些细胞图像是从巴氏涂片的966个群集细胞图像中手工提取的。（包含：im_Dyskeratotic im_Koilocytotic im_Metaplastic im_Parabasal im_...

shell脚本：unexpected end of file错误解决

Brave_heart4pzj的博客

03-25

1963

shell

samtools sort -@ 4 -t CB -O BAM -o $new $idsamtools排序原理过程

生信小博士的博客

03-26

3476

读取每个小文件的第一个read，将它们按照位置大小排序，选出最小的一个read，将该read写入输出文件中，并从对应的小文件中读取下一个read，并继续进行排序，重复此过程直到所有read都被排序。因此，在进行单细胞的scvelo分析时，通常需要对BAM文件进行排序，以确保后续的数据分析结果的准确性和可靠性。读取每个小文件的第一个read，将它们按照位置大小排序，选出最小的一个read，将该read写入输出文件中，并从对应的小文件中读取下一个read，并继续进行排序，重复此过程直到所有read都被排序。

seurat 与geo geo单细胞导入seurat实战数据下载数据分析 scrna-seq

生信小博士的博客

01-21

6208

前面jimmy老师分享了两个祖传的单细胞转录组数据分析代码，非常给力，是标准流程：祖传的单个10x样本的seurat标准代码祖传的单个10x样本的seurat标准代码（人和鼠需要区别对待）其中有一个环节是需要比较seurat分群以及singleR的分群，这样就可以合理的命名啦。在jimmy老师的督促下，我使用老师的代码处理了GSE135927数据集，直接套用了jimmy老师的标准代码，希望对所有的初学者有帮助！首先进入GEO可以看到是两个10X的样本：教程目录大纲如下： 1、准备原始分析数据 2.

单细胞数据读取（二）之Read10X读不出来dgCMatrix报错

weixin_43949246的博客

11-18

1万+

前面我们也遇到过10x的数据读取不进去，如果大家遇到下面的报错，可以通过修改10x的原始重新读取，详细可以见链接https://blog.csdn.net/weixin_43949246/article/details/121225791 当然，这次跟大家说的是另一种10x的数据，同样的，我们先说报错，如果出现dgTMatrix的报错，我们该怎么办呢？这里给大家推荐一种解决方法，首先先定义个函数，修改读取的方式 read10x_ymc <- function( data.dir = NULL

单细胞10x的数据读取不进去怎么办？

weixin_43949246的博客

11-09

9758

面对越来越多的单细胞数据被上传至NCBI上，单细胞数据挖掘分析也逐渐走入大家的眼中。如何寻找一个合适的单细胞数据用于后续非常重要，这时候大家可能会经常遇到这么一个问题，一个标准的10x的数据，为什么我怎么也读取不进去呢？，这里我们以GSE163974数据为例。首先，进入NCBI，下载该数据集https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi 下载之后解压一个一个样本一个文件夹，如GSM4994385样本修改文件名这时候，打开Rstudio #修改工作

SRA文件的下载（prefetch）和解压SRA文件（fastq-dump）

热门推荐

yearstime的博客

03-24

2万+

sra文件下载方式 NCBI-SRA和EBI-ENA数据库 SRA数据库: Sequence Read Archive：隶属NCBI (National Center for Biotechnology Information)，它是一个保存大规模平行测序原始数据以及比对信息和元数据 (metadata) 的数据库，所有已发表的文献中高通量测序数据基本都上传至此，方便其他研究者下载及再研究。其中的数据则是通过压缩后以.sra文件格式来保存的，SRA数据库可以用于搜索和展示SRA项目数据，包括SRA主页和 E

10X cellranger count 后的数据认识和使用

高锦的博客

12-27

2626

对SRR7722937的三个fastq文件运行cellranger count ，完成后得到结果文件存于SRR7722937，打开其下的outs文件夹下游分析所需的是红框的文件，利用R读入以上文件，形成expression_matrix #!/usr/bin/R library(Matrix) cellbarcodes ...

10X单细胞merge

08-29

在使用10X单细胞测序数据时，merge是将多个样本的数据合并在一起的一种常见操作。具体来说，通过将不同样本的数据读取并创建为Seurat对象，然后使用merge函数将这些对象合并在一起。这可以使用以下步骤完成： 1. 读取数据：首先，使用Read10X函数读取每个样本的数据，并创建相应的Seurat对象。这可以使用lapply函数将每个样本的数据读取为一个列表。例如，使用system.time函数来记录运行时间，然后使用lapply函数循环读取数据并创建Seurat对象。最后，返回一个包含Seurat对象的列表。 2. 处理不规则的数据：有时候，可能会遇到一些不规则的数据，例如某个样本的基因和条码信息以txt格式而不是tsv格式给出。在这种情况下，需要将这些样本的数据分开读取。可以使用类似的方法读取不规则的数据，并将其创建为Seurat对象。 3. 合并数据：一旦所有样本的数据都已读取并创建为Seurat对象，可以使用merge函数将它们合并在一起。首先，使用names函数将每个Seurat对象的名称设置为对应的GSE编号。然后，使用merge函数将这些Seurat对象合并在一起，使用add.cell.ids参数指定样本的标识符。这将创建一个包含合并后数据的新的Seurat对象。综上所述，对于10X单细胞数据的merge操作，您可以按照上述步骤进行处理。123 #### 引用[.reference_title] - *1* *2* *3* [GEO10X单细胞scRNK-seq中数据的下载与整理读取](https://blog.csdn.net/m0_58549466/article/details/128302022)[target="_blank" data-report-click={"spm":"1018.2226.3001.9630","extra":{"utm_source":"vip_chatgpt_common_search_pc_result","utm_medium":"distribute.pc_search_result.none-task-cask-2~all~insert_cask~default-1-null.142^v93^chatsearchT3_1"}}] [.reference_item style="max-width: 100%"] [ .reference_list ]

“相关推荐”对你有帮助么？

非常没帮助
没帮助
一般
有帮助
非常有帮助

提交