Skip to content

Commit

Permalink
remove Sleuth.
Browse files Browse the repository at this point in the history
  • Loading branch information
jianhong committed Jul 25, 2024
1 parent 45d0edf commit c3e67b6
Show file tree
Hide file tree
Showing 2 changed files with 1 addition and 23 deletions.
4 changes: 1 addition & 3 deletions README.md
Original file line number Diff line number Diff line change
Expand Up @@ -39,9 +39,7 @@ variable instead of plain text to pass along passwords and other secrets in dock

## pipeline for RNA-seq

The sample RNA-seq data will be analyzed with [kallisto](https://pachterlab.github.io/kallisto/about)
+ [sleuth](https://pachterlab.github.io/sleuth/about),
[kallisto](https://pachterlab.github.io/kallisto/about) or [Salmon](https://combine-lab.github.io/salmon/)
The sample RNA-seq data will be analyzed with [kallisto](https://pachterlab.github.io/kallisto/about) or [Salmon](https://combine-lab.github.io/salmon/)
+ [tximport](https://bioconductor.org/packages/tximport/) + [DESeq2](https://bioconductor.org/packages/DESeq2).

## pipeline for ChIP-seq
Expand Down
20 changes: 0 additions & 20 deletions vignettes/scripts.Rmd
Original file line number Diff line number Diff line change
Expand Up @@ -227,26 +227,6 @@ res.kallisto[!is.na(res.kallisto$padj), ]
```


### run Sleuth

The following code is showing how to use sleuth to do DE analysis.
However, the sleuth package is not updated. It may be broken.
```{r, eval=FALSE}
#BiocManager::install("pachterlab/sleuth")
library(sleuth)
samples <- dir("kallisto_quant")
s2c <- data.frame(path=dir("kallisto_quant", full.names = TRUE),
sample=samples,
condition=sub(".rep.", "", samples))
s2c
so <- sleuth_prep(s2c, ~-1+condition, target_mapping = t2g)
so <- sleuth_fit(so)
## which_beta must be in the colname of design table.
so <- sleuth_wt(so, which_beta = "conditionAblated")
sleuth_live(so)
```


## scripts for ChIPseq

### prepare index file for bwa
Expand Down

0 comments on commit c3e67b6

Please sign in to comment.