R/fcc.R

"fcc" <-
function(data,pheno) {
	if (is(data,"gwaa.data")) data <- data@gtdata else if (!is(data,"snp.data")) stop ("Wrong argument (should be \"gwaa.data\" or \"snp.data\")")
	out <- .C("fastcc",as.raw(data@gtps),as.integer(pheno),as.integer(data@nids),as.integer(data@nsnps), chi2 = double(7*data@nsnps), PACKAGE="GenABEL")$chi2
	out
}

Try the GenABEL package in your browser

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

GenABEL documentation built on May 30, 2017, 3:36 a.m.