RStudio + sleuth! Error in if (sum(filter_bool) == 0) { : missing value where TRUE/FALSE needed
0
0
Entering edit mode
2.3 years ago
R • 0

hello my friends! i have a problem with Sleuth. RStudio output that error: Error in if (sum(filter_bool) == 0) { : missing value where TRUE/FALSE needed

sample_id <- dir(
                  file.path('~/Desktop/xxxxxx/kallisto_output')
)

kal_dirs <- file.path(
                  "~/Desktop/xxxxxx/kallisto_output", 
                  sample_id
)

s2c <- read.table(
    file.path("~/Desktop/xxxxxx/metadata/meta.tsv"), 
    header = TRUE, 
    stringsAsFactors = FALSE, 
    sep = " "
)

s2c <- dplyr::mutate(
                  s2c, 
                  path = kal_dirs
)

s2c$genotype_variation_s <- as.factor(
                  s2c$gene
)

s2c$genotype_variation_s <- relevel(
                  s2c$gene, 
                  ref = "c-like"
)

so <- sleuth_prep(s2c, full_model = ~gene,
                  read_bootstrap_tpm = TRUE, 
                  extra_bootstrap_summary = TRUE, 
                  num_cores = 4)


OUTPUT:
reading in kallisto results
dropping unused factor levels
...........................
normalizing est_counts
Error in if (sum(filter_bool) == 0) { : 
  missing value where TRUE/FALSE needed

s2c contain:
sample gene path
M_u_Kidney_clike c-like ~/path/to/kallisto/output
M_u_Kidney_x1 isoform ~/path/to/kallisto/output
M_u_Kidney_x2 isoform ~/path/to/kallisto/output
M_u_Kidney_x3 isoform ~/path/to/kallisto/output
M_u_Kidney_x4 isoform ~/path/to/kallisto/output
R_b_Kidney_clike c-like ~/path/to/kallisto/output
R_b_Kidney_x1 isoform ~/path/to/kallisto/output
R_b_Kidney_x2 isoform ~/path/to/kallisto/output
.................
R_b_Kidney_xX isoform ~/path/to/kallisto/output
_______________________________________________________________________________
27 objects
RNAseq gene • 1.5k views
ADD COMMENT
0
Entering edit mode

You probably aren't receiving any replies because your code is poorly formatted (it's very difficult to read it) and you aren't providing enough information (e.g. print out the first few lines of what each of your variables contain; what does t2g look like? what does s2c look like? etc.). Unless those things are fixed, it's hard for people to help you.

ADD REPLY
0
Entering edit mode

Hello my friend! thanks for answering! I redid the text to be more readable, I hope it got better.

ADD REPLY
0
Entering edit mode

There's only one sample in your "c-like" condition. That will be problematic for sleuth. You need replicates in all your conditions. You can't compare 1 sample with 26 other samples.

ADD REPLY
0
Entering edit mode

no, c-like 8 pieces.

for each tissue (kidney, skin, liver) 1 gene, somewhere even 2 genes

ADD REPLY
0
Entering edit mode

Thanks for clarifying. However, your entire path column has ~/path/to/kallisto/output for every single sample. Every sample should belong in a different directory.

ADD REPLY

Login before adding your answer.

Traffic: 2504 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6