看板Statistics
标 题Re: [问题] R program的指令
发信站无名小站 (Thu Nov 9 17:26:35 2006)
转信站ptt!Group.NCTU!grouppost!Group.NCTU!wretch
※ 引述《[email protected] (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海