Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Issue with simpleLoess #278

Open
asan-emirsaleh opened this issue Jun 7, 2023 · 1 comment
Open

Issue with simpleLoess #278

asan-emirsaleh opened this issue Jun 7, 2023 · 1 comment

Comments

@asan-emirsaleh
Copy link

asan-emirsaleh commented Jun 7, 2023

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
Copy link

apsteinberg commented Feb 27, 2024

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

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

2 participants