I am new to bioinformatics and I am not sure what the erorr could be.
I want to analyze single cell data with Seurat, I have filtered my genes.tsv file, to only contain protein coding genes, since then I can not create a seurat object when merging my files (barcode.tsv, genes.tsv, matrix.tsv)
here is how the code and output looks like
> for (sample in files){
+ tmp_data = Read10X(data.dir = paste0(working_dir,"/raw_only_protein_coding/",sample,"/longreads_genes_seurat"))
+ tmp_seu = CreateSeuratObject(counts = tmp_data,
+ project = sample)
+ tmp_seu$sample = sample
+ tmp_seu$disease = sample_meta[sample,"disease"]
+ tmp_seu$region = sample_meta[sample,"region"]
+ assign(sample,tmp_seu)
+ }
Error in fixupDN.if.valid(value, x@Dim) :
length of Dimnames[[1]] (168) is not equal to Dim[1] (22517)
thank you
I made sure that the genes file did not have a different structure, added headers, only difference between original and it, is that the rows with genes with non coding protein function are deleted.
Diana A. is a new contributor to this site. Take care in asking for clarification, commenting, and answering.
Check out our Code of Conduct.