R/SNAL.R

In PAGWAS: Pathway Analysis Methods for Genomewide Association Data

Defines functions SNAL

Documented in SNAL

SNAL <-
function(y,G,P,a,s2){
	P=data.matrix(P)
	Z=data.matrix(G)
	X=Z%*%P
	ZZ=Z%*%t(Z)
	V=(a*ZZ+diag(length(y)))
	W=matrix.invroot.calculation(V)
	Y=W%*%y
	Phi=W%*%X
	results=SNAL.calculation(Y,Phi,s2)
	results
}

Try the PAGWAS package in your browser