fit_athena function

Porting of 'ATHENA' linear combination fitting

The function can be used to check which combinations of standards produce a good fit and if output from 'ATHENA' is similar.

fit_athena(all.samples, all.standards, LC.vals, amoSTD, ex.smaller = NULL,
  file.output = NULL, best.fits = NULL)

Arguments

• all.samples: List of all samples
• all.standards: List of all standards
• LC.vals: The fitting range values for the linear combination fitting
• amoSTD: Use at most X standards
• ex.smaller: Exclude portions smaller than a given value (decimal form), default to NULL
• file.output: Possibility to have a file output, default to NULL
• best.fits: Possibility to output more than the best fit (e.g. the first 10 best fits), default to 1

Examples

data(stdmix)
corr.spec.standards  <- initial_load(specdat[1:4], 
  corr.norm = c(-36, -15, 37, 58))
corr.spec.samples    <- initial_load(specdat[5:8], 
  corr.norm = c(-36, -15, 37, 58))
athena.fit <- fit_athena(all.samples = corr.spec.samples, 
  all.standards = corr.spec.standards, LC.vals = c(-14, 46), amoSTD = 4)
## exclude portions smaller 5% = 0.05
athena.fit.exlcude <- fit_athena(all.samples = corr.spec.samples, 
  all.standards = corr.spec.standards, 
  LC.vals = c(-14, 46), amoSTD = 4, ex.smaller = 0.05)