bootstrap: Obtain a bootstrap sample of Genomic Architecture s

In danjlawson/robustarchitecture: Estimate Genomic Architecture parameters using robust statistics

Description

Resample the data dat with replacement a given number nbs times, for each running fullsearch and reporting the inferred parameter.

Usage

bootstrap(
  nbs,
  dat,
  bins = defaultbins(),
  thresholds = c(2.5, 3),
  sigma = NULL,
  verbose = TRUE,
  ...
)

Arguments

nbs	The number of bootstrap samples to obtain. We recommend 20 for exploration, 100 for estimating standard deviations and more for estimating the shape of the distribution.
dat	A data frame of (f,beta) pairs
bins	(default: defaultbins) the histogram bins to use
thresholds	(Default: defaultthresholds) the tail areas to compare (in s.d units).
sigma	(Default: NULL, meaning use the best value) Genomic architecture sigma
verbose	(Default: TRUE) Whether to print iteration progress.
...	Additional parameters to fullsearch.

Value

A vector of length nbs of estimates of s.

See Also

linesearch to explore a range of s, or fullsearch to search a single dataset.

danjlawson/robustarchitecture documentation built on Dec. 19, 2021, 8:09 p.m.