Re: [問題] R program的指令
※ 引述《hppavilion.bbs@ptt.cc (hppavilion)》之銘言：
> ※ 引述《hppavilion (hppavilion)》之銘言：
> : 想請問一下R的指令
> : (1) 解估計方程式的解 (f1,...fn,x,y已知)
> : f1(x,y,β)=0,f2(x,y,β)=0,...fn(x,y,β)=0
> : 求β=?
> 舉例來說,在logistic迴歸中
> 若迴歸方程式為P(Y=1)=exp(β1+β2x)/(1+exp(β1+β2x)
> 則我們可列出估計方程式為
> Σ(yi-exp(β1+β2xi)/(1+exp(β1+β2xi))=0 ----(1)
> Σxi(yi-exp(β1+β2xi)/(1+exp(β1+β2xi))=0 ----(2)
> 利用牛頓法求連立解β1,β2
> 請問若現在只有(1),(2)怎麼利用R求出連立解
也還是用牛頓法呀!
> : (2) 若f(x,y,β)已知, 要怎麼令函數 g(x,y,β)=df/dx
> : (若f(x,y,β)不好自己先微好在令g=df/dx)
> 舉例來說,若f(x,y,β)=exp(xy)+1/exp((xβ)^(1/2)+y)
> 我要怎麼用程式去寫g(x,y,β)=df(x,y,β)/dx
> (請別跟我說自己先微分微好在令g函數)
用 D() 或是 deriv(), 看需求.
--
夫兵者不祥之器物或惡之故有道者不處君子居則貴左用兵則貴右兵者不祥之器非君子
之器不得已而用之恬淡為上勝而不美而美之者是樂殺人夫樂殺人者則不可得志於天下
矣吉事尚左凶事尚右偏將軍居左上將軍居右言以喪禮處之殺人之眾以哀悲泣之戰勝以
喪禮處之道常無名樸雖小天下莫能臣侯王若能守之萬物將自賓天地相合以降甘露民莫
之令而自均始制有名名亦既有夫亦將知止知止可以不殆譬道之在 140.116.152.108海
討論串 (同標題文章)
完整討論串 (本文為第 3 之 3 篇)：