A preprocessing function for objects of class AVS.

Share:

Description

This is a generic function.

Usage

1
avs.preprocess(object, nrand=1000, mergeColinked=TRUE, reldata="RTNdata.LDHapMap.rel27", snpop="all", verbose=TRUE)

Arguments

object

an object. When this function is implemented as the S4 method of class AVS-class, this argument is an object of class 'AVS'.

nrand

a single integer value specifying the size to the random variant set.

mergeColinked

a single logical value specifying to merge co-linked markers eventually present in more than one cluster in the AVS (when mergeColinked=TRUE) or not (when mergeColinked=FALSE).

reldata

a single character value specifying the RTNdata package used to build the associated variant set. Current options are "RTNdata.LDHapMap.rel27" and "RTNdata.LD1000g.rel20130502" (for additional details on the Linkage Disequilibrium data, see the RTNdata package documentation).

snpop

a single character value specifying the universe size used to build the null distributions for the VSE/EVSE analysis (i.e. random AVSs). Options: 'all' or 'ld'. The 1st option represents all the SNPs genotyped in the Linkage Disequilibrium (LD) dataset, while the 2nd represents all the SNPs observed after the LD threshold is applyed. Alternativelly, 'snpop' can be a customized 'BED file' data frame with rs# identifiers.

verbose

a single logical value specifying to display detailed messages (when verbose=TRUE) or not (when verbose=FALSE).

Author(s)

Mauro Castro

See Also

AVS-class

Examples

1
2
3
4
5
6
7
8
9
## Not run: 
#This example requires the RTNdata package! (currently available under request)
library(RTNdata.LDHapMap.rel27)
data(bcarisk)
avs <- new("AVS", markers=bcarisk)
avs <- avs.preprocess(avs, nrand=100)
avs.get(avs)

## End(Not run)

Want to suggest features or report bugs for rdrr.io? Use the GitHub issue tracker.