關(guān)于本站
人大經(jīng)濟(jì)論壇-經(jīng)管之家：分享大學(xué)、考研、論文、會計(jì)、留學(xué)、數(shù)據(jù)、經(jīng)濟(jì)學(xué)、金融學(xué)、管理學(xué)、統(tǒng)計(jì)學(xué)、博弈論、統(tǒng)計(jì)年鑒、行業(yè)分析包括等相關(guān)資源。
經(jīng)管之家是國內(nèi)活躍的在線教育咨詢平臺!
經(jīng)管之家新媒體交易平臺
提供"微信號、微博、抖音、快手、頭條、小紅書、百家號、企鵝號、UC號、一點(diǎn)資訊"等虛擬賬號交易，真正實(shí)現(xiàn)買賣雙方的共贏�！菊堻c(diǎn)擊這里訪問】
TOP熱門關(guān)鍵詞
實(shí)在頭痛的不行了，用R軟件做的極大似然估計(jì)，程序沒有問題，數(shù)據(jù)沒有問題，為什么就是運(yùn)行不出來？坑爹啊程序和數(shù)據(jù)都附上下面是我的程序，求大神指教我到底錯哪里了>d<-read.csv(file=file.choose(),header=T) ...
壇友互助群 |
掃碼加入各崗位、行業(yè)、專業(yè)交流群 |
下面是我的程序，求大神指教我到底錯哪里了
- > d<-read.csv(file=file.choose(),header=T)
- > x1<-d$ln目的地面積
- > x2<-d$目的地人口密度
- > x3<-d$目的地GDP
- > x4<-d$目的地距離
- > x5<-d$目的地是否為中心城區(qū)
- > y1<-d$出發(fā)地人口密度
- > y2<-d$出發(fā)地GDP
- > y3<-d$出發(fā)地是否為中心城區(qū)
- > s1<-d$出行成立
- > s2<-d$出行不成立
- > params<-c(0,0,0,0,0,0,0,0,0)
- > #定義log-likelihood函數(shù)
- > LL<-function(params,d)
- + {#參數(shù)"params"是一個向量.
- + #參數(shù)"d"，是整理之后的數(shù)據(jù).
- + f<-s1*log(1+exp((params[1]*x1+params[2]*x2+params[3]*x3+params[4]*x4+params[5]*x5-params[6]*y1-params[7]*y2-params[8]*y3+params[9])))+s2*log(1+exp((params[1]*x1+params[2]*x2+params[3]*x3+params[4]*x4+params[5]*x5-params[6]*y1-params[7]*y2-params[8]*y3+params[9])))
- + ll<-sum(f)
- + #log-likelihood函數(shù)
- + return(ll)
- + }
- > ##最大似然值的計(jì)算
- > res<-optim(params,LL,method="BFGS",hessian=TRUE,control=list(fnscale=-1))
- > ##參數(shù)估計(jì)
- > b<-res$par
- > hhh<-res$hessian
- > ##t值的計(jì)算
- >e<-sqrt(diag(solve(hhh)))
- >t<-b/e
- >##最終似然值
- >LL<-res$value
- >##初始似然值
- >A<-sum(s1)
- >B<-sum(s2)
- >L0=A*log(2)+B*log(2)
- >##擬合度計(jì)算
- >##結(jié)果輸出
- >##p^2的值
- >cat("roh = ",(L0-LL)/L0)
- >##修正的p^2值
- >cat("rohbar=",(L0-(LL-length(b)))/L0)
- >print(res)
- >print(t)
- >print(L0)
掃碼或添加微信號：壇友素質(zhì)互助
「經(jīng)管之家」APP：經(jīng)管人學(xué)習(xí)、答疑、交友，就上經(jīng)管之家！
免流量費(fèi)下載資料----在經(jīng)管之家app可以下載論壇上的所有資源，并且不額外收取下載高峰期的論壇幣。
涵蓋所有經(jīng)管領(lǐng)域的優(yōu)秀內(nèi)容----覆蓋經(jīng)濟(jì)、管理、金融投資、計(jì)量統(tǒng)計(jì)、數(shù)據(jù)分析、國貿(mào)、財(cái)會等專業(yè)的學(xué)習(xí)寶庫，各類資料應(yīng)有盡有。
來自五湖四海的經(jīng)管達(dá)人----已經(jīng)有上千萬的經(jīng)管人來到這里，你可以找到任何學(xué)科方向、有共同話題的朋友。
經(jīng)管之家（原人大經(jīng)濟(jì)論壇），跨越高校的圍墻，帶你走進(jìn)經(jīng)管知識的新世界。
掃描下方二維碼下載并注冊APP
免流量費(fèi)下載資料----在經(jīng)管之家app可以下載論壇上的所有資源，并且不額外收取下載高峰期的論壇幣。
涵蓋所有經(jīng)管領(lǐng)域的優(yōu)秀內(nèi)容----覆蓋經(jīng)濟(jì)、管理、金融投資、計(jì)量統(tǒng)計(jì)、數(shù)據(jù)分析、國貿(mào)、財(cái)會等專業(yè)的學(xué)習(xí)寶庫，各類資料應(yīng)有盡有。
來自五湖四海的經(jīng)管達(dá)人----已經(jīng)有上千萬的經(jīng)管人來到這里，你可以找到任何學(xué)科方向、有共同話題的朋友。
經(jīng)管之家（原人大經(jīng)濟(jì)論壇），跨越高校的圍墻，帶你走進(jìn)經(jīng)管知識的新世界。
掃描下方二維碼下載并注冊APP
本文論壇網(wǎng)址：http://xalimeijing.com/thread-2445051-1-1.html
您可能感興趣的文章
人氣文章
本文標(biāo)題：關(guān)于極大似然估計(jì)
本文鏈接網(wǎng)址：http://xalimeijing.com/jg/qikan_qikanku_2445051_1.html
2.轉(zhuǎn)載的文章僅代表原創(chuàng)作者觀點(diǎn),與本站無關(guān)。其原創(chuàng)性以及文中陳述文字和內(nèi)容未經(jīng)本站證實(shí),本站對該文以及其中全部或者部分內(nèi)容、文字的真實(shí)性、完整性、及時性，不作出任何保證或承若；
3.如本站轉(zhuǎn)載稿涉及版權(quán)等問題,請作者及時聯(lián)系本站,我們會及時處理。
