startvec0 {mixstock} | R Documentation |
Provides raw (untransformed) starting vector of source contribution parameters for mixed stock analysis.
startvec0(sourcesamp, mixsamp=NULL, type="equal", sd=1, lmin=0.001) startvec(sourcesamp, mixsamp=NULL, type="equal", marktype="sample",a=1,cond=FALSE,transf="full", fuzz=0,sd=1) lsolve(n, s, tol = 1e-05, warn = FALSE)
sourcesamp |
Marker frequencies in sources, or a list with
source and mixed samples, or a turtle.data object. |
mixsamp |
Marker frequencies in mixed stock |
type |
Various options for setting starting contributions.
equal : equal contributions from all sources. random :
random multinomial sample with equal multinomial
probabilities. rand2 : random sample from transformed normal
variates with mean 0 and standard deviation sd .
A number n between 1 and the number of sources
inclusive gives a starting condition with 95% of the contribution
from source n and the other 5% evenly split between the
other sources. The default is to attempt a solution of the linear
equation (sourcesamp*f=mixsamp ) and use these values as the
starting contributions. |
sd |
standard deviations for starting type rand2 |
lmin |
When doing linear solutions, the boundary values are
{lmin , 1-lmin . |
marktype |
method for starting marker frequencies:
sample uses the observed sample probabilities; random
used observed sample probabilities; weighted does a
Bayes-weighted start (a la Masuda and Pella) |
cond |
Conditional likelihood? (i.e. cond=TRUE gives just
the parameters for source contributions, not parameters for source
marker frequencies) |
a |
prior strength parameter |
transf |
transform ("full","part", or "none") |
fuzz |
fuzz parameter for moving parameters away from the boundary |
n |
source samples (HxR matrix) |
s |
mixed stock samples (vector) |
tol |
tolerance for linear fit |
warn |
warn if numeric problems with solution |
A parameter vector of the contributions from each of the sources:
just the raw source contributions in the case of startvec0
, or
source contributions and possibly marker frequencies (transformed
or untransformed) in the case of startvec
.
lsolve
attempts to get a starting value by solving
the linear equation (solve(n,s,tol=tol)
)
Ben Bolker
data(simex) startvec0(simex) startvec(simex,transf="part") startvec(simex)