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.2 Index]