get.bot: Masuda-Pella interchange functions

In mixstock: functions for mixed stock analysis

Description

Read and write control files for Masuda and Pella's program

Usage

get.bot(fn)
get.frq(fn)
get.ctl(fn)
get.bse(fn)
put.ctl(fn,which=1,tot,ranseed=c(899271,4480026,90092812),
                 thin=1,title="auto-gen R input",sourcenames,fprior,startval,H)
put.bse(fn,sourcesamp)
put.mix(fn,mixsamp)
put.mp(data,tot=25000,title="Input file auto-gen from R",
                fn="test",ranseed=c(899271,4480026,90092812))
get.mp.input(fn,which=1)

Arguments

fn	file name (without extension)
which	numeric value of control file
data	data in mixstock.data format (list of sourcesamp and mixsamp)
tot	total number of Gibbs steps
ranseed	numeric vector for random-number seed
thin	thinning factor
title	title for input files
sourcenames	vector of source names
fprior	prior on source contributions (default=c(1/R,R)
startval	starting values for source contributions
H	number of markers
sourcesamp	matrix (HxR) of marker samples in sources
mixsamp	vector of marker sample in mixed population

Value

produces output files or returns a list of

sourcesamp	source samples (matrix)
mixsamp	pooled samples (vector)
a	prior strength
startiter	starting iteration
maxiter	maximum iteration
startfval	starting fvalue
thin	thinning factor
fprior	source prior
rptiter	how often to print output to the screen

Note

put.mp and get.mp.input are portmanteau functions

Author(s)

Ben Bolker

mixstock documentation built on May 2, 2019, 6:48 p.m.