| 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)