selectLocalOptimalSNPs-methods: Select locally-optimized SNPs on one or more chromosomes

In selectSNPs: Optimal Selection of SNPs Using a Unified Local Function

Description

This function is used to select locally-optimal SNPs using an unified, local function uder varied scenarios. The idea is to divide each chromosome into equal segments, and within each segment, an local optimal SNP is selected as the one having the largest weighted score between the uniformmness (U score) and a function (average Shannon entropy, or E score) of minor allele frequency. If the weight for the E score is zero, it selects (approximately) uniformly-distributed SNPs. If the weight for the U score is zero, it selects SNPs with the largest minor allele frequencies within pre-defined local chromosome regions or bins.

Usage

selectLocalOptimalSNPs(object, ...)

## S4 method for signature 'numeric'
selectLocalOptimalSNPs(object, f, n, t1 = 1, t2 = 1, w1 = 0.5, w2 = 0.5)

## S4 method for signature 'Chrom'
selectLocalOptimalSNPs(object, n, t1 = 1, t2 = 1, w1 = 0.5, w2 = 0.5)

## S4 method for signature 'Map'
selectLocalOptimalSNPs(object, n, t1 = 1, t2 = 1, w1 = 0.5, w2 = 0.5)

Arguments

object	An input object, which can be a vector of map positions, a Chrom object, a Map object, or a data frame containing map information for all the candidate SNPs.
...	Extra input parameters, as needed.
f	A vector of minor allele frequencies for all the candidate SNPs.
n	NUmber of SNP positions to be filled.
t1	A penalty parameter (0~100) for the E score, which by default is 1.
t2	A shrinkage parameter (0-100) for the U score, which by default is 1.
w1	A weight for the impact of MAF on selecting markers, which by default is 0.5.
w2	A weight for the impact of map position on selecting markers, which by default is 0.5.

Value

A Map object with locally-optimal SNPs when the input is a Map object, or a Chrom object with locally-optimal SNPs when the input is a Chrom object, or a vector of map positions of locally-optimal SNPs when the input are a vector of map positions of SNPs.

Author(s)

Nick X-L Wu

References

1. Wu X-L, Li H, Ferretti R, Simpson B, Walker J, Parham J, Mastro L, Qiu J, Schultz T, Tait RG. Jr., and Bauck S. (2020) A unified local objective function for optimally selecting SNPs on arrays for agricultural genomics applications. Anim. Genet. (in press)

2. Wu XL, Xu J, Feng G, Wiggans GR, Taylor JF, He J, Qian C, Qiu J, Simpson B, Walker J, Bauck S. Optimal Design of Low-Density SNP Arrays for Genomic Prediction: Algorithm and Applications. PLoS One. 2016, 11(9):0161719. doi: 10.1371/journal.pone.0161719. eCollection 2016.

Examples

data(bov80K)
map1<-selectLocalOptimalSNPs(bov80K,bin=1,n=1000,w1=0,w2=1)
plot(map1)
scoreU(map1)
scoreE(map1)

map2<-selectLocalOptimalSNPs(bov80K,bin=2,n=1000,w1=0,w2=1)
plot(map2)
scoreU(map2)
scoreE(map2)

map3<-selectLocalOptimalSNPs(bov80K,bin=4,n=1000,w1=0.5,w2=0.5)
plot(map3)
scoreU(map3)
scoreE(map3)
 

selectSNPs documentation built on Feb. 28, 2020, 5:08 p.m.