R/QMiRKAT.R

Defines functions QMiRKAT

QMiRKAT <-
function(r,ker,model){
	if(model =='gaussian'){
		Q = as.numeric((1/2)*(1/var(r))*t(r)%*%ker%*%r)
	}else{
		Q <- as.numeric((1/2)*t(r)%*%ker%*%r)
	}
	Q	
}
JiyuanHu/massMap documentation built on May 26, 2020, 2:20 a.m.