HI, I tried to processed the fragments data in the ATAC-seq results, but I cannot find the idx information. If I generate the idx file for each fragment file, I found that seurat cannot correctly read it: Error in ecdf(x = object$TSS.enrichment): 'x' must have 1 or more non-missing values Traceback: 1. TSSEnrichment(pbmc) 2. TSSFast(object = object, assay = assay, tss.positions = tss.positions, . process_n = process_n, verbose = verbose, region_extension = region_extension) 3. ecdf(x = object$TSS.enrichment) 4. stop("'x' must have 1 or more non-missing values") Therefore, what is the correct approach to process the fragment data? Can I find the processed data in ROSMAP? Thanks.

Created by Tianyu Liu tl688
Hi Tianyu, Thank you for your post. Please provide the synIDs to the files or folders you are experiencing these issues with so that I may investigate further. Kind regards, Victor

How to correctly processed the atac-seq data in MIT_MULTIOME page is loading…