R/fcc.new.R

In GenABEL: genome-wide SNP association analysis

"fcc.new" <-
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_new",as.raw(data@gtps),as.integer(pheno),as.integer(data@nids),as.integer(data@nsnps), chi2 = double(6*data@nsnps), PACKAGE="GenABEL")$chi2
	out
}