Convert Scanonevar from LODs to Empirical p-values

Description

scanonevar.to.p.values takes a scanonevar with LODs as units and maxes from permutation scans, estimates an extreme value distribution for the maxes, and returns the probability of observing the LOD scores in those EVDs.

Usage

1

Arguments

scanonevar

the scanonevar in LODs to be converted to emprirical p values

perm.scan.maxes

the tbl_df object returned by scanonevar.perm, the maximum LOD score observed on a per-scan, per-chromosome-type basis in permutation scans.

Details

none

Value

Returns a scanonevar object in terms of p-values, with attr(x, 'units') = 'emp.ps'.

Author(s)

Robert Corty rcorty@gmail.com

See Also

scanonevar, scanonevar.perm

Examples

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
## Not run: 
   my.scanonevar <- scanonevar.perm(cross = my.cross,

   my.perms <- scanonevar.perm(cross = my.cross,
                               n.perms = 10))

   scanonevar.to.p.values(scanonevar = my.scanonevar,
                          perm.scan.maxes = my.perms)
 
## End(Not run)