Re: [問題] R program的指令
※ 引述《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函數)
--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 59.117.119.222
討論串 (同標題文章)
本文引述了以下文章的的內容：
完整討論串 (本文為第 2 之 3 篇)：