scanonevar.to.p.values: Convert Scanonevar from LODs to Empirical p-values

Description Usage Arguments Details Value Author(s) See Also Examples

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)


Search within the vqtl package
Search all R packages, documentation and source code

Questions? Problems? Suggestions? or email at ian@mutexlabs.com.

Please suggest features or report bugs with the GitHub issue tracker.

All documentation is copyright its authors; we didn't write any of that.