999精品在线视频,手机成人午夜在线视频,久久不卡国产精品无码,中日无码在线观看,成人av手机在线观看,日韩精品亚洲一区中文字幕,亚洲av无码人妻,四虎国产在线观看 ?

基于SAS程序模擬運算的函數模型Y=(a+bX)/X+ε參數估計及其應用

2013-12-18 06:58:34賈洪艷
關鍵詞:程序分析方法

賈洪艷

(連云港職業技術學院 信息工程學院, 江蘇 連云港222006)

由一個或者一組非隨機變量來估計或預測某一個隨機變量的觀測值時,所建立的數學模型和進行的統計分析稱為回歸分析[1].長期以來,回歸分析在宏觀經濟預測、醫療診斷、農業產量分析等方面得到了廣泛應用.回歸分析本質上是借助于函數模型,分析并確定函數模型中的參數.一直以來,由于計算手段的限制,對模型參數的估計常常依賴單因素方差分析、最小二乘法等經典統計方法,這些經典統計方法都有一個共同的特點,那就是參數估計必須有明確的數學公式.近十幾年來,計算機技術的飛速發展,以及吉布斯抽樣方法(Gibbs方法)的引入,為回歸分析拓寬思路和方法提供了可能.本文依據函數模型Y=(a+bX)/X+ε,根據貝葉斯統計原理,對模型進行參數估計,然后利用SAS程序對參數進行概率抽取,從而實現對模型的參數估計.

1 函數模型的經典統計學估計

本文針對如下函數模型討論

(1)

1.1 最小二乘估計法

由回歸模型(1),設

(2)

或者

1.2 極大似然估計法

由于n組樣本觀測值可看作n個相互獨立的隨機變量,因此其似然函數為

利用

可計算出參數a,b的估計值,其估計結果表達式為

依據上述參數a,b的估計公式,對參數進行估計.

2 貝葉斯估計法

2.1 貝葉斯統計原理

設θ為包含所有未知參數的向量,p(θ|Y)是θ關于事件Y的條件概率,由統計學原理[2]知

p(θ,Y)=p(Y)p(θ|Y)=p(θ)p(Y|θ)

進一步有

(3)

其中:p(θ)和p(θ|Y)分別為先驗概率和后驗概率,p(Y|θ)為Y的概率密度函數.由于p(Y)無法計算且不依賴參數θ,故可略去,(3)式可進一步改寫為

p(θ|Y)∝p(θ)p(Y|θ)

(4)

依據(4)式可以推導出θ中各個參數的后驗分布.由于這些后驗分布并不是明確的數學公式,因此,需要編寫SAS程序,依據SAS程序,運用Gibbs抽樣產生后驗樣本,通過一定的方法,分析該后驗樣本數據,進而獲得該參數的估計.

2.2 函數模型參數的后驗分布

(5)

(6)

(7)

(8)

(9)

(10)

3 基于SAS程序的函數模型參數的Gibbs模擬算法

2)根據式(6)~式(7),運用SAS程序更新a;

3)根據式(8)~式(9),運用SAS程序更新b;

5)重復步驟2)~步驟4)t輪;

6)分析樣本(a(1),b(1),σ2(1)),(a(2),b(2),σ2(2)),…,(a(t),b(t),σ2(t)),計算各參數的估計值.

4 實例分析

選用文獻[5]中例12.1演示上述算法.

例研究“岱字棉”自播種至齊苗(以80%出苗為準)的天數(X)與日平均土溫(Y,℃)的關系,得表1結果,試作回歸分析.

表1 土溫(Y,℃)和岱字棉播種至齊苗天數(X)的關系

5 結束語

長期以來,對數學模型參數進行估計,一直采用最小二乘法、極大似然法等經典統計方法.這些方法的主要優點是利用偏導函數指導搜索尋優方向,搜索效率比較高,在適合的初始值條件下,可以較快實現目標函數的優化.但是,這些方法也明顯存在著諸多缺陷,主要體現在:(1)對于一個需要進行參數估計的非線性方程必須提供其偏導函數,但事實上某些方程本身較為復雜,不易獲得其偏導函數.(2)編制通用軟件難度較大.(3)對于不同數據進行估計時,必須根據數據的基本情況對模型提供合適的初始值,否則很容易陷入局部最優陷阱.盡管貝葉斯方法早就被提出,但是,由于計算上的困難,其統計方法很少提及.得益于計算機技術的快速發展,貝葉斯方法重新被提起.貝葉斯統計對所要進行參數估計的非線性方程,不需要提供其偏導函數,只是對方程模型的參數進行概率抽取,因此,無論方程復雜與否,均不受影響.同時,與經典統計不同的是,貝葉斯學派很重視先驗信息的收集、挖掘和加工,使它數量化,形成先驗分布,參加到統計推斷中來,從而極大地提高了統計推斷的質量.本研究嘗試將SAS程序和貝葉斯統計相結合對數學模型參數進行估計,為統計方法提供了新的手段和途徑.

[1]何蘊理,賀亞平,陳中和,等.概率論與數理統計[M].北京:高等教育出版社,2010.

[2]茆詩松.貝葉斯統計[M].北京:中國統計出版社,2008.

[3]Geman S,Geman D.Stochastic relaxation,Gibbs distribution and the Bayesian restoration of images[J].IEEE Tansactions on Pattern Analysis and Mechine Intelligence,1984,6:721~741.

[4]Gelfand A E,Smith A F M.Sampling based on approaches to calculating marginal densities[J].Journal of American Statistical Association,1990,85(2):339~355.

[5]莫惠棟.農業試驗統計[M].上海:上??茖W技術出版社,1984.

猜你喜歡
程序分析方法
隱蔽失效適航要求符合性驗證分析
試論我國未決羈押程序的立法完善
人大建設(2019年12期)2019-05-21 02:55:44
電力系統不平衡分析
電子制作(2018年18期)2018-11-14 01:48:24
“程序猿”的生活什么樣
英國與歐盟正式啟動“離婚”程序程序
環球時報(2017-03-30)2017-03-30 06:44:45
電力系統及其自動化發展趨勢分析
用對方法才能瘦
Coco薇(2016年2期)2016-03-22 02:42:52
創衛暗訪程序有待改進
中國衛生(2015年3期)2015-11-19 02:53:32
四大方法 教你不再“坐以待病”!
Coco薇(2015年1期)2015-08-13 02:47:34
捕魚
主站蜘蛛池模板: 久久无码av三级| 国产高清不卡| 精品国产一区二区三区在线观看| 国产视频一区二区在线观看| 亚洲精品卡2卡3卡4卡5卡区| 国产精品99在线观看| 国产免费福利网站| 欧美国产综合色视频| 九九精品在线观看| 日本五区在线不卡精品| 国产精品欧美在线观看| 亚洲三级色| 午夜不卡福利| 欧美亚洲中文精品三区| 欧美h在线观看| 91精品国产麻豆国产自产在线| 亚洲性日韩精品一区二区| 欧美精品影院| 亚洲AⅤ无码日韩AV无码网站| 在线网站18禁| 久久无码免费束人妻| a毛片在线播放| 日韩一级毛一欧美一国产| 国产导航在线| 久久a级片| 国产精品高清国产三级囯产AV| 国产精品一线天| 99精品国产高清一区二区| 亚洲精品第一页不卡| 国产精品999在线| 日韩在线欧美在线| 国产熟女一级毛片| 毛片在线播放网址| 久久久久人妻一区精品色奶水 | 极品私人尤物在线精品首页| 亚洲性日韩精品一区二区| 视频国产精品丝袜第一页| 亚洲天堂网2014| 亚洲永久色| 色综合天天视频在线观看| 亚洲最猛黑人xxxx黑人猛交| 亚洲国产中文精品va在线播放 | 色老头综合网| 国产噜噜在线视频观看| 内射人妻无码色AV天堂| 亚洲第一视频区| 亚洲精品在线91| 永久在线播放| 中文字幕人妻无码系列第三区| 久久久久久高潮白浆| www精品久久| 欧美成人午夜视频免看| 免费午夜无码18禁无码影院| 天堂成人av| 欧美国产日韩一区二区三区精品影视| 91小视频版在线观看www| 欧美一道本| 日韩av资源在线| 国产在线精品美女观看| 欧美午夜久久| 精品国产成人高清在线| 色综合中文综合网| 九九热在线视频| 国产一级精品毛片基地| 午夜三级在线| 国产精品一区二区国产主播| 亚洲福利一区二区三区| 国产精品9| 亚洲国产天堂在线观看| 国产真实乱子伦视频播放| 国产欧美日韩综合一区在线播放| 男人天堂伊人网| 成年女人18毛片毛片免费| 色亚洲成人| 国产呦精品一区二区三区下载 | 久久午夜夜伦鲁鲁片无码免费| 亚洲网综合| 日韩精品高清自在线| 伊人91视频| 久久国产高潮流白浆免费观看| 狠狠色综合网| 丁香婷婷久久|