forked from elijahedmondson/NASA-GRSD-HS
-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathcloud.file.R
26 lines (21 loc) · 846 Bytes
/
cloud.file.R
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
library(BSgenome.Mmusculus.UCSC.mm10)
library(doParallel)
library(foreach)
library(Rsamtools)
library(VariantAnnotation)
library(DOQTL)
library(GenomicRanges)
library(regress)
library(MASS)
library(lmtest)
library(HZE)
options(stringsAsFactors = F)
load(file = "/home/ubuntu/CloudLOG.Rdata")
setwd("~/Dropbox/Rstudio.cloud/WD")
outdir = "~/Dropbox/Rstudio.cloud/WD"
perms <- GRSDassoc.perms(perms = 2, chr = 1:2, pheno = HZE, Xchr = F, addcovar = HZE.add,
pheno.col = "HCC", probs = probs, K = K, markers = markers,
snp.file = snp.file, outdir = "~/Dropbox/Rstudio.cloud/WD", tx = "Test",
sanger.dir = "/home/ubuntu/HS.sanger.files/")
source("https://bioconductor.org/biocLite.R")
biocLite("doParallel", "foreach", "Rsamtools", "VariantAnnotation", "GenomicRanges")