| proxyA {ClamR} | R Documentation |
Optimum Wilkinson Curve Fitting
Description
Runs one cycle of optimum Wilkinson curve fitting for a single sinusoid fitting.
Usage
proxyA(ax, ay, xin)
Arguments
ax |
x-axis values |
ay |
y-axis values |
xin |
starting model: c( Phs,Pos,Amp,Prd) |
Details
This program implements the Wilkinson and ivany approach to climate time series modeling. This is used in the more sophisticated proxyJK code.
Value
Optimum model, vector of 4 values
Note
Uses stats package routine optim for optimization
Author(s)
Jonathan M. Lees<jonathan.lees@unc.edu>
References
Wang, T., Surge, D., and Lees, J. M., (2015) ClamR: A Statistical Evaluation of Isotopic and Temperature Records in Sclerochronologic Studies. Palaeogeography, Palaeoclimatology, Palaeoecology, doi:10.1016/j.palaeo.2015.07.008.\
Wilkinson, B. H. and Ivany, L. C., Paleoclimatic inference from stable isotope profiles of accretionary biogenic hardparts; a quantitative approach to the evaluation of incomplete data, Palaeogeography, Palaeoclimatology, Palaeoecology, vol. 185, no. 1-2, pp.95-114, 01 Sep 2002
See Also
proxyJK
Examples
data(CLAM1)
x = CLAM1$x
y = CLAM1$y
dx = 3.392
Aout = proxyA(x, y, dx)