亚洲免费av电影一区二区三区,日韩爱爱视频,51精品视频一区二区三区,91视频爱爱,日韩欧美在线播放视频,中文字幕少妇AV,亚洲电影中文字幕,久久久久亚洲av成人网址,久久综合视频网站,国产在线不卡免费播放

        ?

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

        2013-12-18 06:58:34賈洪艷
        關(guān)鍵詞:對模型后驗參數(shù)估計

        賈洪艷

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

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

        1 函數(shù)模型的經(jīng)典統(tǒng)計學估計

        本文針對如下函數(shù)模型討論

        (1)

        1.1 最小二乘估計法

        由回歸模型(1),設

        (2)

        或者

        1.2 極大似然估計法

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

        利用

        可計算出參數(shù)a,b的估計值,其估計結(jié)果表達式為

        依據(jù)上述參數(shù)a,b的估計公式,對參數(shù)進行估計.

        2 貝葉斯估計法

        2.1 貝葉斯統(tǒng)計原理

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

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

        進一步有

        (3)

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

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

        (4)

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

        2.2 函數(shù)模型參數(shù)的后驗分布

        (5)

        (6)

        (7)

        (8)

        (9)

        (10)

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

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

        3)根據(jù)式(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)),計算各參數(shù)的估計值.

        4 實例分析

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

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

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

        5 結(jié)束語

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

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

        [2]茆詩松.貝葉斯統(tǒng)計[M].北京:中國統(tǒng)計出版社,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]莫惠棟.農(nóng)業(yè)試驗統(tǒng)計[M].上海:上海科學技術(shù)出版社,1984.

        猜你喜歡
        對模型后驗參數(shù)估計
        基于新型DFrFT的LFM信號參數(shù)估計算法
        基于對偶理論的橢圓變分不等式的后驗誤差分析(英)
        光源對模型貼圖的影響——3種人造光源在三維數(shù)字化采集中的應用
        廣州文博(2020年0期)2020-06-09 05:15:44
        貝葉斯統(tǒng)計中單參數(shù)后驗分布的精確計算方法
        一種基于最大后驗框架的聚類分析多基線干涉SAR高度重建算法
        雷達學報(2017年6期)2017-03-26 07:53:04
        Logistic回歸模型的幾乎無偏兩參數(shù)估計
        基于向前方程的平穩(wěn)分布參數(shù)估計
        蒙藥特潤舒都樂對模型小鼠脾臟NK細胞活性的影響
        蒙醫(yī)開竅補腎針刺法對模型大鼠胸腺、脾臟指數(shù)的影響
        蒙醫(yī)催熟療法對模型大鼠炎癥因子影響的實驗研究
        黄色录像成人播放免费99网| 亚洲成人色区| 国产美女遭强高潮网站| 亚洲日本无码一区二区在线观看 | av在线入口一区二区| 一本丁香综合久久久久不卡网站| 嫖妓丰满肥熟妇在线精品| 国产亚洲AV天天夜夜无码| 亚洲一区二区三区福利久久蜜桃| 亚洲色偷偷偷综合网| 国产亚洲情侣一区二区无| 亚洲国产成人精品激情资源9| 国产一区二区三区成人av| 92午夜少妇极品福利无码电影| 亚洲国产成人久久一区www| 2021精品综合久久久久| 中文字幕日韩高清乱码| 欧美黑寡妇特a级做爰| japanesehd中国产在线看| av毛片一区二区少妇颜射| 国产成年人毛片在线99| 中文字幕在线观看亚洲日韩| 亚洲AV无码乱码1区久久| 丝袜美腿亚洲综合一区| 久久天天躁夜夜躁狠狠| 成人免费ā片在线观看| 亚洲青涩在线不卡av| 中文字幕本久久精品一区| 中文日韩亚洲欧美制服| 一区二区三区不卡在线| 一区二区三区日本视频| 天堂а在线中文在线新版| 日本高清www午色夜高清视频| 蜜桃视频中文字幕一区二区三区| 中文字幕人妻饥渴浪妇| 国产97在线 | 中文| 欧美性爱一区二区三区无a| 久久午夜av一区二区| 免费毛片a线观看| 97中文字幕在线观看| 水蜜桃视频在线观看入口|