作者qqwerty (love ocean)
看板Statistics
标题Re: [问题] 使用R来解估计方程式
时间Thu Apr 10 22:54:51 2008
※ 引述《[email protected] ()》之铭言:
: ※ 引述《[email protected] (love ocean)》之铭言:
: > 最近刚刚接触R,之前都是使用Gauss来写程式
: > 之前的程式老师希望我能用R重新在写一次
: > 但我碰到一个问题是,我写出了一个估计方程式来解参数(a,b)
: > 例如
: > g1(a,b,x,y)=0
: > g2(a,b,x,y)=0
: > (x,y) is data
: > 之前在Gauss可以使用eqSolve来解这个连立方程式
: > 但在R里却不知道该用什麽指令来解
: > 不知道有没有了解的人可以帮忙一下?...感谢了!
: 先令你的 objective function 为
: g1(a,b,x,y)^2 + g2(a,b,x,y)^2
: 再利用 optim() 来解未知参数 a, b
嗯~首先感谢大家热心的帮忙,提供了相关的指令能参考
不过我有一些问题想请教一下
大部份推文所提到的指令都是对於某个f(a,b,x,y)找出使
这个函数最小的a,b的值
例如在R中的Example所提的f(x1,x2)=100(x2-x1^2)^2+(1-x1)^2
找出x1,x2使f(x1,x2)最小
用optim()可以求出为(1,1)
转成微积分的问题可视为在解
f1(x1,x2)=x2-x1^2=0
f2(x1,x2)=200x1(x2-x1^2)+(1-x1)=0
若是以上面回覆所说的
我令g(x1,x2)=f1^2+f2^2去求x1,x2使g有最小值
是否跟原来要找f有最小值的问题不同?
除此之外,我也能够令g*=|f1|+|f2|来解
这样看来似乎与我原先的问题又不大相同了
我的问题比较像是解f1,f2的连立方程解
--
※ 发信站: 批踢踢实业坊(ptt.cc)
◆ From: 220.133.117.103