-
Notifications
You must be signed in to change notification settings - Fork 0
/
OverlapPermutation.Rmd
39 lines (30 loc) · 1.1 KB
/
OverlapPermutation.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
---
title: "OverlapScorePermutationTest"
author: "Gregory Jefferis"
date: "21/06/2017"
output: html_document
---
```{r setup, include=FALSE}
knitr::opts_chunk$set(echo = TRUE)
```
## R Markdown
This is an R Markdown document. Markdown is a simple formatting syntax for authoring HTML, PDF, and MS Word documents. For more details on using R Markdown see <http://rmarkdown.rstudio.com>.
When you click the **Knit** button a document will be generated that includes both content as well as the output of any embedded R code chunks within the document. You can embed an R code chunk like this:
# Read in Data
```{r}
lhon_overlap <- readRDS("data_val.rds")
```
Let's pick 50% of total input as our threshold for interesting PNs
```{r}
# figure out 50% of distribution
threshold=sum(lhon_overlap$value)*0.50
lhon_overlap$cs=cumsum(lhon_overlap$value)
subset(lhon_overlap, cs<sum(lhon_overlap$value)*0.5)
```
So it turns out that all 6 inputs above that threshold are food. What's the likelihood
of that by random draw
```{r}
vv=replicate(10000,all(sample(lhon_overlap$BehaviouralSignificance, 6)=="Food"))
table(vv)
```
Low!