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

        ?

        基于MCEM算法的多元正態(tài)分布均值向量估計(jì)

        2022-08-17 09:17:40殷雨晨陳兆榮
        關(guān)鍵詞:補(bǔ)法樣本量正態(tài)分布

        殷雨晨,陳兆榮

        (銅陵學(xué)院 經(jīng)濟(jì)學(xué)院,安徽 銅陵 244000)

        在數(shù)據(jù)收集和處理的過程中,往往會(huì)出現(xiàn)數(shù)據(jù)缺失的現(xiàn)象,不恰當(dāng)?shù)奶幚矸椒▽?huì)影響分析結(jié)果的可靠性.因此,如何在數(shù)據(jù)存在缺失的情況下減小估計(jì)的誤差成了亟待解決的問題.對(duì)于缺失的數(shù)據(jù),利用數(shù)據(jù)特征來進(jìn)行填補(bǔ)是目前統(tǒng)計(jì)學(xué)中最常用的方法,包括均值插補(bǔ)法、回歸插補(bǔ)法和最近距離插補(bǔ)法等.其中,均值插補(bǔ)法就是用所研究樣本的已觀測(cè)數(shù)據(jù)的均值作為缺失值的替代值,這種方法操作簡單,實(shí)用性強(qiáng).不過也存在著明顯的弊端,即均值填補(bǔ)忽略了數(shù)據(jù)的離散程度,可能會(huì)造成較大的推斷誤差.

        除了進(jìn)行傳統(tǒng)的方法進(jìn)行數(shù)據(jù)填補(bǔ)以外,EM(Expectation Maximization)算法也逐漸被應(yīng)用到經(jīng)濟(jì)、醫(yī)療和工程等不完全數(shù)據(jù)的處理中.孫大飛[1]將EM算法應(yīng)用到混合密度極大似然參數(shù)估計(jì)當(dāng)中,驗(yàn)證了此算法的收斂性和有效性.張夢(mèng)琇[2]將EM算法應(yīng)用到左截?cái)嘤覄h失數(shù)據(jù)的幾何分布參數(shù)估計(jì)中.此外,學(xué)者們也對(duì)EM算法不斷進(jìn)行研究和改進(jìn),先后提出了ECM(Expectation Conditional Maximization)、ECME(Expectation Conditional Maximization Either)、PX-EM(Parameter-Expanded EM)和MCEM(Monte Carlo EM)等算法.溫艷清[3]在區(qū)間型數(shù)據(jù)條件下利用ECM算法對(duì)Weibull分布進(jìn)行了極大似然估計(jì).嚴(yán)海芳[4]應(yīng)用MCEM算法給出了對(duì)數(shù)正態(tài)分布的參數(shù)的迭代公式,給出了比EM算法更有效、收斂速度更快的模擬結(jié)果.本文將MCEM算法引入到含缺失數(shù)據(jù)的多元正態(tài)數(shù)據(jù)的均值向量估計(jì)當(dāng)中,得到迭代公式,并據(jù)此進(jìn)行數(shù)值模擬,驗(yàn)證算法的合理性和準(zhǔn)確性.

        1 EM算法和MCEM算法簡介

        EM算法,也稱最大期望算法,是一種能夠在不完全數(shù)據(jù)下通過迭代得到參數(shù)最大似然估計(jì)的算法,最早由Dempster[5]提出,該算法迭代速度較快并且估計(jì)結(jié)果具有很好的穩(wěn)定性.它是在概率模型中尋找參數(shù)最大似然估計(jì)或者最大后驗(yàn)估計(jì)的算法,其中概率模型依賴于無法觀測(cè)的隱性變量.EM算法的每次迭代包括一個(gè)E步和一個(gè)M步,兩個(gè)步驟交替進(jìn)行計(jì)算.其基本思想是:首先根據(jù)已經(jīng)給出的觀測(cè)數(shù)據(jù),估計(jì)出模型參數(shù)的值;然后再依據(jù)上一步估計(jì)出的參數(shù)值估計(jì)缺失數(shù)據(jù)的值,再根據(jù)估計(jì)出的缺失數(shù)據(jù)加上之前已經(jīng)觀測(cè)到的數(shù)據(jù)重新再對(duì)參數(shù)值進(jìn)行估計(jì),然后反復(fù)迭代,直至最后收斂,迭代結(jié)束.

        E步:計(jì)算目標(biāo)Q函數(shù)

        Q(θ|θ(t))=Ef[l(θ|Yobs,Ymis)|Yobs,θ(t)]

        (1)

        M步:求目標(biāo)Q函數(shù)的最優(yōu)解

        θ(t+1)=argmaxQ(θ|θ(t));

        (2)

        E2步:計(jì)算目標(biāo)Q函數(shù),其中

        (3)

        M步:求目標(biāo)Q函數(shù)的最優(yōu)解

        θ(t+1)=argmaxQ(θ|θ(t));

        (4)

        2 含缺失數(shù)據(jù)的多元正態(tài)數(shù)據(jù)的MCEM算法

        2.1 多元正態(tài)分布的定義

        多元正態(tài)分布是一元正態(tài)分布的推廣[7].已知一元正態(tài)分布的密度函數(shù)為

        (5)

        將式(5)進(jìn)行推廣,可以給出多元正態(tài)分布的定義

        定義1p元正態(tài)變量X=(X1,X2,…,Xp)′的概率密度函數(shù)為

        (6)

        其中,|Σ|為協(xié)方差陣Σ的行列式.

        2.2 多元正態(tài)分布的條件分布

        記X~Np(μ,Σ),將X,μ和Σ剖分如下

        其中,X(1),μ(1)為q×1維,Σ11為q×q維,(X(1)|X(2))為給定X(2)時(shí)X(1)的條件分布[8].

        定理1 設(shè)X~Np(μ,Σ),Σ>0,則

        (X(1)|X(2))~Nq(μ1·2,Σ11·2),

        其中

        (7)

        (8)

        2.3 二元正態(tài)分布均值的MCEM估計(jì)

        (9)

        (10)

        最后實(shí)現(xiàn)M步,對(duì)目標(biāo)函數(shù)求極值

        (11)

        (12)

        3 MCEM算法的數(shù)值模擬

        由表1可知,相比于傳統(tǒng)的均值插補(bǔ)法,MCEM算法估計(jì)的均值離真實(shí)值的偏差更小,效果更好.并且隨著樣本量的增大,MCEM算法的估計(jì)精度也有所增加,即當(dāng)樣本量從100增加到1000時(shí),μ1估計(jì)偏差的絕對(duì)值從0.3938降低到0.0351,μ2估計(jì)偏差的絕對(duì)值從0.1198降低到0.0721.因此,本文提出的MCEM算法比均值插補(bǔ)法更有優(yōu)勢(shì),在多元正態(tài)分布的均值向量估計(jì)中是切實(shí)可行的.

        表1 不同樣本量下均值估計(jì)結(jié)果

        4 結(jié)論

        本文將MCEM算法應(yīng)用到含缺失數(shù)據(jù)的多元正態(tài)分布的均值向量估計(jì)中,既解決了傳統(tǒng)均值填補(bǔ)方法忽略數(shù)據(jù)離散程度的問題,同時(shí)也化解了 EM算法中計(jì)算E步期望的難題,將原本復(fù)雜的問題簡化為蒙特卡洛模擬求解積分.最后,數(shù)值模擬的結(jié)果也驗(yàn)證了MCEM算法在多元正態(tài)分布均值向量估計(jì)中是可行且高效的.

        猜你喜歡
        補(bǔ)法樣本量正態(tài)分布
        醫(yī)學(xué)研究中樣本量的選擇
        基于絡(luò)病理論探討絡(luò)虛通補(bǔ)法在氣虛血瘀型椎動(dòng)脈型頸椎病中的應(yīng)用
        基于少數(shù)類過采樣的傾向得分匹配插補(bǔ)法
        航空裝備測(cè)試性試驗(yàn)樣本量確定方法
        響應(yīng)傾向得分匹配插補(bǔ)法
        Sample Size Calculations for Comparing Groups with Binary Outcomes
        基于對(duì)數(shù)正態(tài)分布的出行時(shí)長可靠性計(jì)算
        正態(tài)分布及其應(yīng)用
        正態(tài)分布題型剖析
        χ2分布、t 分布、F 分布與正態(tài)分布間的關(guān)系
        日韩精品熟女中文字幕| 色婷婷七月| 精品人无码一区二区三区 | 成在人线av无码免费| 蜜臀av一区二区| 久久狠狠高潮亚洲精品暴力打 | 亚洲一本到无码av中文字幕| 嗯啊哦快使劲呻吟高潮视频| 高清在线亚洲中文精品视频| 国产精品久久久一本精品| 国产精品亚洲A∨无码遮挡| 女同性恋亚洲一区二区| 婷婷久久亚洲中文字幕| 自拍偷自拍亚洲精品第按摩| 亚洲午夜成人精品无码色欲| 男男啪啪激烈高潮cc漫画免费 | 久精品国产欧美亚洲色aⅴ大片| 综合激情网站| 蜜臀av国内精品久久久人妻| 99精品人妻少妇一区二区三区| 国产精品人伦一区二区三| 未发育成型小奶头毛片av| av综合网男人的天堂| 中文字幕人妻丝袜乱一区三区| 亚洲午夜精品a区| 国产中文字幕乱码在线| 一本色道久久综合中文字幕| 成人在线观看视频免费播放| 亚洲一区二区国产一区| 日韩夜夜高潮夜夜爽无码| 永久免费av无码入口国语片| 精品香蕉久久久午夜福利| 久热香蕉精品视频在线播放| 人妻精品人妻一区二区三区四五| 风流熟女一区二区三区| 成人欧美一区二区三区在线| 久久国产劲暴∨内射| 美丽人妻被按摩中出中文字幕 | 人人澡人人澡人人看添av| 曰本女人与公拘交酡免费视频| 亚洲综合自拍|