-
Notifications
You must be signed in to change notification settings - Fork 5
/
Copy pathREADME.Rmd
49 lines (33 loc) · 1.85 KB
/
README.Rmd
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
---
output: github_document
---
[](https://travis-ci.org/leekgroup/phenopredict)
[](https://ci.appveyor.com/project/leekgroup/phenopredict)
[](https://coveralls.io/r/leekgroup/phenopredict?branch=master)
<!-- README.md is generated from README.Rmd. Please edit that file -->
```{r, echo = FALSE}
knitr::opts_chunk$set(
collapse = TRUE,
comment = "#>",
fig.path = "README-"
)
```
# phenopredict
The goal of phenopredict is to build predictors from expression data (RNA-Seq) that will predict necessary phenotype information across samples.The phenotypes sex, tissue, sequencing strategy, and sample source have been predicted across the ~70,000 samples currently available in [recount](https://jhubiostatistics.shinyapps.io/recount/).
## Installation
You can install phenopredict from GitHub with:
```{r gh-installation, eval = FALSE}
# install.packages("devtools")
devtools::install_github("leekgroup/phenopredict")
```
## Examples using phenopredict
Details of how the predicted phenotypes previously generated for the [recount data](https://jhubiostatistics.shinyapps.io/recount/)) are [available on GitHub](https://github.com/ShanEllis/phenopredict_phenotypes).
## Adding predicted phenotypes in recount
If you're looking to find the phenotype information previously generated across the recount samples, an example of how to use `add_predictions()` within recount is provided below.
```{r example}
library('recount')
download_study('ERP001942', type='rse-gene')
load(file.path('ERP001942', 'rse_gene.Rdata'))
rse <- scale_counts(rse_gene)
rse_with_pred <- add_predictions(rse)
```