maxlikboot: Bootstrap estimates for parameters and confidence intervals...

View source: R/maxlikboot.R

maxlikbootR Documentation

Bootstrap estimates for parameters and confidence intervals for heteroplasmy data

Description

Compute bootstrap estimates for parameters and confidence intervals for a given heteroplasmy set We can do this while imposing a specific h0 as an argument or by allowing a search over h0 values.

Usage

maxlikboot(h, nboot = 1000, conf.level = 0.95, h0 = F)

Arguments

h

A vector containig heteroplasmy measurements. Every observation should be in [0,1].

nboot

The number of bootstrap samples. Default value is 1000

conf.level

The preferred confidence interval calculation, Default value is 0,95 (95%).

h0

Logical parameter. A particular h0 value Default is to treat h0 as a fit parameter

Value

The maximum likelihood for the input data according to the Kimura distribution (using bootstrapping)

Examples

 X.1 = rnorm(50,0.5,0.1)
maxlik(X.1,nboot=10000,0.95)

kostasgian21/heteroplasmy documentation built on Jan. 30, 2024, 12:30 a.m.