PAGWAS (version 2.0)

SNAL: Sparse Normal/Adaptive lasso method for finding associated pathways

Description

Sparse Normal/Adaptive lasso method applied for finding the associated pathways. The iterative algorithm suggested by Wipf and Nagarajan (2008) is applied. A vector equal to the number of tested pathways is returned, the zero entries of the vector correspond to the pathways that are not associated. The posterior estimates of the beta coefficients are also returned as they are described by Wipf and Nagarajan (2008).

Usage

SNAL(y, G, P, a, s2)

Arguments

y

Response vector of length N

G

Genotype matrix, with N rows and L columns (number of tested SNPs)

P

Pathway matrix, with L columns and M columns (number of tested pathways)

a

Hyper-parameter of the variance assumed for the integrated out SNP effects

s2

Variance assumed for the response variable, the tuning parameter of adaptive lasso

Value

gamma.star

Estimates of gamma hyper-parameters

ARD

Posterior estimates of beta coefficients

References

Evangelou, M., Dudbridge, F., Wernisch, L. (2014). Two novel pathway analysis methods based on a hierarchical model. Bioinformatics, 30(5), 690 - 697.

Wipf, D. and Nagarajan, S. (2008). A new view of automatic relevance determination. Advances in Neural Information Processing Systems, 20

See Also

SNAL.calculation

Examples

## Not run: 
# 	data(genotypes)
# 	G=genotypes
# 	data(pathways)
# 	data(SNPs)
# 	data(genes)
# 	snps.genes=snps.to.genes(SNPs,genes,distance=0)
# 	snps.paths=snps.to.pathways(pathways,snps.genes)
# 	P=create.pathway.df(G,snps.paths)
# 	y=rnorm(nrow(G),mean=0,sd=10)
# 	SNAL(y,G,P,a,s2)## End(Not run)