jackknife: A blocked jacknife estimator

In andremrsantos/htspop: High-Troughput Sequencing Population Genetics Statistics

Description

It performs a blocked jacknife estimator of mean and standard error of the input vector. Given a matrix as input, it considers the first column as numerator and second column as denominator. It can also perform a weighted estimation by counting the number of non-NA values per block. This implementation is largely inspired on (Treemix)https://bitbucket.org/nygcresearch/treemix/wiki/Home implementation of the algorithm.

Usage

jackknife(x, block = 10, weighted = TRUE)

Arguments

x	A numeric vector or matrix to estimate mean and standard error.
block	A integer scalar indicating the size of blocks (greater than zero).
weighted	A boolean scalar indicating if should perform a weighted estimation.

andremrsantos/htspop documentation built on May 14, 2020, 11:40 a.m.