bootstrap {sojourner}R Documentation

bootstrap

Description

Bootstrap confidience intervals with standard errors. bootstrap resamples dataset (e.g. diffusion coefficients) to calculate confidience intervals for a statistic measure of dataset.

Usage

bootstrap(fittedObj,n.reps=100)

Arguments

n.reps

number of replicates for bootstrapping

fittedObj

output from fitNormDistr

Details

A wrapper of boot::boot and mixtools::boot.se adapted for data format in sojourner package. Also returns stderr information by running boot.se from mixtools and an additional method for one-component distribution calculates the stderr separately. For multi-component distributions, the boot.se function from the mixtools package was used. For single-component distributions, a separate function was used to calculate the stderr and confidence interval.

Value

List of length 2 containing:

Examples

# read in track files
folder=system.file('extdata','SWR1',package='sojourner')
trackll=createTrackll(folder=folder, input = 3)

# filter track based length
trackll.flt=filterTrack(trackll,filter=c(min=5,max=Inf))
MSD=msd(trackll.flt,dt=6,summarize=FALSE,plot=TRUE)
dcoef=Dcoef(MSD=MSD,method='static',plot=FALSE)
# fit the dcoef result
fittedObj=fitNormDistr(dcoef)

# bootstrap new datasets
d.boot=bootstrap(fittedObj)
# manually set the number of bootstrap samples to 50
d.boot=bootstrap(fittedObj, n.reps=50)


[Package sojourner version 1.0.2 Index]