Description Usage Arguments Value Author(s) References See Also Examples
First, the recombination rates per segment are computed based on the regression model (generalized additive models) as well as the bias correction. Consequently, we apply SMUCE (simultaneous multiscale change-point estimator) of Frick (2014) and Futschik et al. (2014) to estimate locations and breakpoints in the recombination map. Under a specific type-I error probability alpha the number of distinct segments with respect to the recombination rate is not overestimated.
1 |
help |
a matrix containing a set of summary statistics is calculated in the function summary_statistics. These values are used in the regression model to calculate the (constant) recombination rates. |
segs |
A (non-negative) integer which reflects the number of segments considered. It is calculated in the program based on the user-defined |
alpha |
A value from the interval (0,1) for the type-I error probability used in the segmentation algorithm. We recommend to use 0.05. We enabled to estimate the recombination map efficiently (without recalculating all summary statistics) under several type-I errors when |
ll |
A (non-negative) integer which reflects the total sequence length of the sequences under study. |
quant |
A value between 0.1 and 0.5 with 0.05 distances in between which reflects the quantile used in the quantile regression. We recommend to use the 0.35 quantile. |
rescale |
an optional logical value: if |
constant |
an optional logical value: by default |
demography |
an optional character value: by default an empty string ("") indicates that the recombination rate estimation is estimated under neutrality. If |
regMod |
an optional character string: for the default empty string "" *LDJump* uses an existing regression model (constant population size or simple demography example, depending on |
seq.full.cor |
The final estimate of the recombination map. Depiction with plot-function of stepR package. |
pr.full.cor |
A vector of (constant) estimates of the recombination rate per segment. |
Philipp Hermann philipp.hermann@jku.at, Andreas Futschik
Frick, K., Munk, A., and Sieling, H. (2014). Multiscale change-point inference. Journal of the Royal Statistical Society: Series B, 76(3), 495–580.
Futschik, A., Hotz, T., Munk, A., and Sieling, H. (2014). Multiscale DNA partitioning: Statistical evidence for segments. Bioinformatics, 30(16), 2255–2262.
Hermann, P., Heissl, A., Tiemann-Boege, I., and Futschik, A. (2019), LDJump: Estimating Variable Recombination Rates from Population Genetic Data. Mol Ecol Resour. doi:10.1111/1755-0998.12994.
LDJump
, vcfR_to_fasta
, getPhi
, summary_statistics
, stepFit
, rq
, gam
1 2 3 | ##### Do not run these examples #####
##### In LDJump.R the function is called as follows #####
##### get_smuce(help, segs, alpha,ll,list.quantile.regs) #####
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.