• R语言

  • 请教：模型拟合问题（四次对称函数拟合）

目的：采用y=a(x-MA)^4+b(x-MA)^2+c模型对temp数据进行拟合，求a,b,c,MA（MA是当y最小时的值）

问题：如何code?



本人尝试着用lm(y~I(x^2)+I(x^4))但又不会求MA，请达人指教？

万分感谢？





temp

y    x

20.26008    -0.08

19.72613    -0.065

19.501619 -0.05

18.72662    -0.03

18.58769    -0.015

18.592199 0.015

18.88372    0.03

19.5453    0.05

19.88743    0.065

20.9914    0.08

18.12336    0
你这两个问题我没时间仔细看了，最小值就是MA=min(y)，后面的线性模型式子是对的（x后面增加-MA就好了）。
谢谢回复，问题没有这么简单MA=min(y)，而是通过现有的数据你和一个对称的四次模型来求得模型y最小时的x.
确切来说，已知一组数据（R,DMA），来模拟得到PMA

模型建立为：



R=exp(A)(DMA-PMA)^4+exp(B)(DMA-PMA)^2+C



其中R,DMA是已知量，设参数A,B,C,且A,B>0，求PMA



请大侠指教！
MA is random, if it's a function of data.
不好意思，现在才大致明白你的意思。



你这个问题恐怕不太容易转换为线性模型来做，我觉得应该用非线性模型吧（PMA也作为一个参数一起估计）
[quote]引用第4楼rtist于2007-08-31 05:03发表的“”:

MA is random, if it's a function of data. [/quote]



Thanks for reply



pls refer to my quesion at "3-Lou"
[quote]引用第5楼谢益辉于2007-08-31 16:19发表的“”:

不好意思，现在才大致明白你的意思。



你这个问题恐怕不太容易转换为线性模型来做，我觉得应该用非线性模型吧（PMA也作为一个参数一起估计）[/quote]



谢谢，请问有constrained非线性的mode拟合的packagel吗？