※ 引述《yhliu.bbs@bbs.wretch.cc (老怪物)》之銘言:
: ※ 引述《bcs.bbs@ptt.cc (= ="frailty..gggg XD)》之銘言:
: > y是0,1變數,誤差項是標準指數分配。
: > 想利用mle求算b係數,一階條件後b不見了。
: > 請問式中推導有哪出錯了?? 大謝^^"
: > y=bx+u, y~binomial,u~exp(λ=1)=> f(u)=exp(-u) F(u)=1-exp(-u)
: ^^^^^^^^^^^^^^^^^^^^^^這可行嗎?
對不起,前文吃掉。更正為:
p = P[u>-bx] = 1-F(-bx)=1-[1-exp(bx)]=exp(bx) {by F(u)=1-exp(-u)}
L = Π{y=0}pΠ{y=1}(1-p)
logL = Σy*log(p)+ Σ(1-y)*log(1-p)
= Σy*log(exp(bx))+Σ(1-y)log(1-exp(bx))
= Σy*bx+Σ(1-y)[log(1)-log(exp(bx))]
= bΣyx+Σ(1-y)[0-bx]
= bΣyx-bΣ(1-y)x
foc Σyx-Σ(1-y)x=0 <=???
--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 140.112.247.214