Nothing
### R code from vignette source 'GLIDE.Rnw'
###################################################
### code chunk number 1: loadLibrary
###################################################
library(GLIDE)
###################################################
### code chunk number 2: data
###################################################
data(simdata)
#The example dataset is a list composed of two dataframes.
#simat stores 20,000 observations of 81 variables,
#including outcome, 5 ajusting covariates, and 75 SNPs.
simdat=simdata$simdat
dim(simdat)
head(colnames(simdat),n=10)
#coeff stores the 75 external regression cofficients.
coeff=simdata$coeff
head(coeff)
###################################################
### code chunk number 3: genotype_columns
###################################################
genotype_columns=which(grepl("^SNP",colnames(simdat)))
###################################################
### code chunk number 4: formula
###################################################
formula=as.formula("outcome~age+sex+pc1+pc2+pc3")
formula
###################################################
### code chunk number 5: out
###################################################
out=glide(formula=formula,exposure_coeff=coeff,genotype_columns,data=simdat,
np=100000,qcutoff=0.2,parallel=TRUE,corenumber=1,verbose=TRUE)
head(out)
###################################################
### code chunk number 6: GLIDE.Rnw:120-124
###################################################
cat('\\begin{figure}[h]\n')
file = "./plot.png"
cat('\\includegraphics{', file, '}\n', sep = '')
cat('\\end{figure}\n')
###################################################
### code chunk number 7: sessionInfo
###################################################
sessionInfo()
Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.