R/SNAL.R

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

Any scripts or data that you put into this service are public.

PAGWAS documentation built on May 30, 2017, 6:40 a.m.