Description Usage Arguments Value Author(s) References Examples
The LPmerge algorithm is called for all linkage groups with an option for multicore processing. The mapfuser package further adds automatic selection of max.interval giving the lowest RMSE
1 2 | LPmerge_par(MF.obj, n.cores = 2, max.interval = 1:3, max.int_sel = "auto",
weights = NULL)
|
MF.obj |
A mapfuser object genetics maps loaded and optionally a reference map |
n.cores |
number of cores |
max.interval |
A whole number specifying the maximum interval size between bins to include in the objective function. An array of numbers can be passed to test different values (one consensus map is produced for each value in the array). |
max.int_sel |
Either automatically select the max.interval with the lowest RMSE for each linkage group or specify manually a vector of values of max.interval to select for the output |
weights |
Optional vector of length T containing the weights for each map in the objective function (see details). If not passed, the maps are given equal weight. |
Three items added to the mapfuser object under the results slot. 1) The consensus map at the selected max.interval, either manually specified or automatically select 2) A list of calculated Root Mean Square Error (RMSE) per chromosome with RMSE for each linkage group ID compared to the consensus map. THe mean and standard deviation over the linkage group IDs is also provided. 3) A list with length equal to the length of the max.interval parameter. Each entry in the list is a data frame containing the consensus map and the component linkage maps.
Dennis van Muijen
Endelman, JB, and C Plomion. 2014. LPmerge: An R package for merging genetic maps by linear programming. Bioinformatics 30:1623-1624.
1 2 3 4 5 6 7 8 9 | ## Not run:
MF.obj <- LPmerge_par(MF.obj = MF.obj, n.cores = 2,
max.interval = 1, max.int_sel = "auto", weights = NULL)
# Plot result
plot(MF.obj, which = "single_map", maps = "consensus")
# Access RMSE table
MF.obj$result$RMSE
## End(Not run)
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.