R/suest.R

In Haplin: Analyzing Case-Parent Triad and/or Case-Control Data with SNP Haplotypes

suest <- function(reslist){

if(!is(reslist[[1]], "haplin")){
	stop('Each element of "reslist" must be the output from a "haplin" run. \n\nNote in particular that "haplinSlide" must be run with the argument "table.output" set to "FALSE".', call. = FALSE)
}

## CHECK FOR ERRORS IN INDIVIDUAL HAPLIN RUNS. REMOVE AND REPORT.
.fjern <- sapply(reslist, function(x) !is.null(attr(x, "error.message")))
.reslist <- reslist[!.fjern]
if(sum(.fjern) > 0){
	cat('The following estimation results were removed \nfrom the list before applying suest \n(due to incomplete estimation results):\n', '"', paste(names(reslist)[.fjern], collapse = '", "'), '"', '\n', sep = "")
}

.suest <- f.suest(.reslist, debug = F, diag.plot = F)
class(.suest) <- "suest"

return(.suest)

}