friendmiao
目的:采用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
yihui
你这两个问题我没时间仔细看了,最小值就是MA=min(y),后面的线性模型式子是对的(x后面增加-MA就好了)。
friendmiao
谢谢回复,问题没有这么简单MA=min(y),而是通过现有的数据你和一个对称的四次模型来求得模型y最小时的x.
friendmiao
确切来说,已知一组数据(R,DMA),来模拟得到PMA
模型建立为:
R=exp(A)(DMA-PMA)^4+exp(B)(DMA-PMA)^2+C
其中R,DMA是已知量,设参数A,B,C,且A,B>0,求PMA
请大侠指教!
rtist
MA is random, if it's a function of data.
yihui
不好意思,现在才大致明白你的意思。
你这个问题恐怕不太容易转换为线性模型来做,我觉得应该用非线性模型吧(PMA也作为一个参数一起估计)
friendmiao
[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"
friendmiao
[quote]引用第5楼谢益辉于2007-08-31 16:19发表的“”:
不好意思,现在才大致明白你的意思。
你这个问题恐怕不太容易转换为线性模型来做,我觉得应该用非线性模型吧(PMA也作为一个参数一起估计)[/quote]
谢谢,请问有constrained非线性的mode拟合的packagel吗?
yihui
friendmiao
Thanks very much~ Let me try!