pachterlab / sleuth

Differential analysis of RNA-Seq
http://pachterlab.github.io/sleuth
GNU General Public License v3.0
304 stars 95 forks source link

Issue with simpleLoess #278

Open asan-emirsaleh opened 1 year ago

asan-emirsaleh commented 1 year ago

Hello! I am enjoyed on the abilities the sleuth provides on RNA-Seq data analysis. So I try to use it in my current project. I have a large dataset, where from I am interested in looking the DE for 15 genes (which corresponds to 224 DE transcripts previously identified by lrt analysis. I ran an analysis again, with gene aggregation, target mapping and filtering by a vector containing the 224 transcript ids.

In sleuth_prep(s2c, num_cores = max(1L, parallel::detectCores() -  :
  325270 target_ids are missing annotations for the aggregation_column: gene.
These target_ids will be dropped from the gene-level analysis.
If you did not expect this, check your 'target_mapping' table for missing values.
fitting measurement error models
shrinkage estimation
Error in simpleLoess(y, x, w, span, degree = degree, parametric = parametric,  : 
  invalid 'x'
Calls: sleuth_fit -> shrink_df -> eval -> loess -> simpleLoess
Execution halted

The commands used:

t2g <- read.csv('DEG_vs_comb_prot_dedup.map.tsv', sep="\t", stringsAsFactors = FALSE)
head(t2g)
filter_id <- t2g$target_id[t2g$gene != '']

so <- sleuth_prep(s2c, 
   num_cores = max(1L, parallel::detectCores() - 10L), 
   extra_bootstrap_summary = TRUE, 
   read_bootstrap_tpm = TRUE, 
   target_mapping = t2g, 
   aggregation_column = "gene", 
   gene_mode = TRUE, 
   filter_target_id = filter_id, 
   max_bootstrap = 500)
apsteinberg commented 8 months ago

Were you able to resolve this in the end? I am getting a similar error.

kylagelev commented 2 months ago

I am also running into this error.

asan-emirsaleh commented 2 months ago

Were you able to resolve this in the end? I am getting a similar error.

Hello @apsteinberg ! I finally resolved the issue, but I don't remember the solution. It might be a dataset issue, and I just filtered out all transcripts that did not correspond the genes of interest.

I have not recieved suggestions from the developer team for any of my questions, but some community members guided me to a right answer. Developers did not response to email, on github or google group forums. The Sleuth seems to be deprecated by developers.

apsteinberg commented 2 months ago

Hi @asan-emirsaleh,

Got it, thanks this is all useful to know!

Best, Asher