※ 引述《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