get.bot {mixstock}R Documentation

Masuda-Pella interchange functions

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


[Package mixstock version 0.9 Index]