如何使用R估算數據點的對數形式？

Question

我有代表對數函數的數據點。

有沒有一種方法可以估算描述這些數據的函數？

謝謝。

Answer 1

如果我理解正確，你想估計一個給定的（x，y）值的函數。如果是，請檢查以下鏈接。

閱讀一下：

http://en.wikipedia.org/wiki/Spline_%28mathematics%29 
http://en.wikipedia.org/wiki/Polynomial_interpolation 
http://en.wikipedia.org/wiki/Newton_polynomial 
http://en.wikipedia.org/wiki/Lagrange_polynomial 

GOOGLE了它：

http://www.stat.wisc.edu/~xie/smooth_spline_tutorial.html 
http://stat.ethz.ch/R-manual/R-patched/library/stats/html/smooth.spline.html 
http://www.image.ucar.edu/GSP/Software/Fields/Help/splint.html 

我從來沒有使用[R所以我不知道是否可行與否，但如果你有Matlab的我可以更說明你。

Answer 2

我假定你的意思是你有矢量y和x，你試着做一個功能y(x)=Alog(x)。
首先，擬合日誌是一個壞主意，因爲它表現不好。幸運的是我們有x(y)=exp(y/A)，所以我們可以擬合一個更方便的指數函數。我們可以用非線性最小二乘做到這一點：

nls(x~exp(y/A),start=list(A=1.),algorithm="port") 

其中start是A初始估計值。這種方法是數值優化，因此可能會失敗。
更穩定的方式是將其變換成線性函數，log(x(y))=y/A和適合使用lm直線：需要

lm(log(x)~y)