im3sanger / dndscv

dN/dS methods to quantify selection in cancer and somatic evolution
GNU General Public License v3.0
212 stars 48 forks source link

Problems encountered in using DNDSCV:[E::faidx_fetch_seq2] Failed to retrieve block. (Seeking in a compressed, .gzi unindexed, file?) Error in h(simpleError(msg, call)) : 在为'strsplit'函数选择方法时评估'x'参数出了错: record 1 (9:12285795-12285926) failed #88

Open manchengyanyu opened 1 year ago

manchengyanyu commented 1 year ago

path_cds_table = system.file("extdata", "5.txt", package = "dndscv", mustWork = TRUE) path_genome_fasta = system.file("extdata", "ref.fa", package = "dndscv", mustWork = TRUE) buildref(cdsfile=path_cds_table, genomefile=path_genome_fasta, outfile = "example_output_refcds.rda", excludechrs="MT") [1/3] Preparing the environment... [2/3] Building the RefCDS object... [E::faidx_fetch_seq2] Failed to retrieve block. (Seeking in a compressed, .gzi unindexed, file?) Error in h(simpleError(msg, call)) : 在为'strsplit'函数选择方法时评估'x'参数出了错: record 1 (9:12285795-12285926) failed file: E:/R/R-4.2.1/library/dndscv/extdata/ref.fa In addition: Warning message: In buildref(cdsfile = path_cds_table, genomefile = path_genome_fasta, : 35670 unique gene IDs (column 1) found. 16245 unique gene names (column 2) found. Consider combining gene names and gene IDs or replacing gene names by gene IDs to avoid losing genes (see useids argument in ? buildref)