gwsem (version 2.1.4)

isSuspicious: Determine which results are suspicious

Description

maturing The GWAS function writes all results, both valid and invalid, to a log file. This function uses heuristics to try to classify rows as suspicious or unsuspicious.

Usage

isSuspicious(result, pars = attr(result, "focus"))

Arguments

result

output from loadResults

pars

names of the parameters available in result

Value

a vector of logicals for each row of result indicating suspicion (if TRUE)

Details

OpenMx reports exceptions in the ‘catch1’ column. Any error message in the ‘catch1’ column is suspicious. Any optimizer status code besides ‘OK’ is suspicious. It is suspicious if the focal parameter or its standard error is NA. If ‘signAdj’ was requested and it is NA then suspicion is also aroused.

See Also

Other reporting: loadResults(), plot.gwsemResult(), signifGxE(), signif()

Examples

# NOT RUN {
tdir <- tempdir()
dir <- system.file("extdata", package = "gwsem")
pheno <- data.frame(anxiety=rnorm(500))
m1 <- buildItem(pheno, 'anxiety')
GWAS(m1, file.path(dir,"example.bgen"),
    file.path(tdir,"out.log"))
r1 <- loadResults(file.path(tdir,"out.log"), "snp_to_anxiety")
r1[isSuspicious(r1, "snp_to_anxiety"),]
# }