王劍紅,楊素芳

(山西藥科職業(yè)學(xué)院 基礎(chǔ)部,山西 太原 030031)

線性模型中總平方和分解公式的證明*

王劍紅,楊素芳

(山西藥科職業(yè)學(xué)院 基礎(chǔ)部,山西 太原 030031)

文中結(jié)合矩陣的知識(shí),通過(guò)正交變換,給出了正交飽和設(shè)計(jì)對(duì)應(yīng)的線性統(tǒng)計(jì)模型中總平方和分解公式的另一種證明方法，優(yōu)化了文獻(xiàn)[1]中對(duì)其的證明.

矩陣;冪等陣;正交矩陣;投影矩陣

對(duì)于正交飽和設(shè)計(jì)問(wèn)題，通?？捎萌缦碌木€性統(tǒng)計(jì)模型來(lái)描述

Y=β01n+β1x1+…+βmxm+ε=μ+ε

(1)

其中Y=(y1,y2,…,yn)T是觀察值向量；xj=(x1j,x2j,…,xnj)/，j=1,2,…,m由試驗(yàn)設(shè)計(jì)來(lái)確定；矩陣X=(x1,x2,…,xm)為上述正交表H所對(duì)應(yīng)的設(shè)計(jì)陣；β=(β0,`β1,…,βm)T；ε=(ε1,ε2,…,εn)/是誤差向量，且εi(i=1,2,…,n)是相互獨(dú)立同分布的隨機(jī)變量，且有ε～N(0,σ2In).

由于模型(1)是基于正交飽和設(shè)計(jì)，故此時(shí)的誤差平方和等于零，從而使總平方和SSj與各列的效應(yīng)平方和之間有如下總平方和分解公式

SST=SS1+SS2+…+SSm,fT=f1+f2+…+fm成立.

這篇文章我們給出了總平方和SST與各列的效應(yīng)平方和SSj之間的矩陣證明方法,優(yōu)化了文獻(xiàn)[1]中對(duì)其的證明.

1 預(yù)備知識(shí)

在一個(gè)試驗(yàn)設(shè)計(jì)中，當(dāng)被考慮因子(包括交互作用)個(gè)數(shù)多到使得需估計(jì)參數(shù)的個(gè)數(shù)達(dá)到可估計(jì)參數(shù)的最大個(gè)數(shù)時(shí)，這樣的試驗(yàn)設(shè)計(jì)稱為飽和設(shè)計(jì).當(dāng)一個(gè)飽和設(shè)計(jì)又為一個(gè)正交設(shè)計(jì)時(shí)，稱為正交飽和設(shè)計(jì).

2 主要研究結(jié)果及證明

又注意到

Tj為對(duì)應(yīng)設(shè)計(jì)陣第j列的置換矩陣，即Tj中只有0和1兩個(gè)元素，而且每一行、每一列有且只有一個(gè)1，其余的元素全是0.由張應(yīng)山的博士論文知,Tj具有存在性.

取

故

這樣,對(duì)于模型(1)，各列的效應(yīng)平方和

YTAjY,j=1,2,…,m.

且Pn,A1,…,Am為相互正交的投影陣.

因?yàn)?/p>

diag(Pr1,Pr2,…,Prm)Tdiag(Pr1,Pr2,…,Prm)=

diag(Pr1,Pr2,…,Prm)且Pndiag(Pr1,Pr2,…,Prm)=

Pn,PnAj=0,AiAj=0(i≠j).

所以Pn,A1,…,Am為相互正交的投影陣.

容易驗(yàn)證

τn=In-Pn=

所以SST=YTτnY=YT(A1+A2+…+Am)Y=

SS1+SS2+…+SSm.

總平方和自由度為fT=n-1,fj=pj-1為第j列平方和的自由度(j=1,2,…,m),而完全正交設(shè)計(jì)是指上述正交設(shè)計(jì)滿足如下的等式n-1=(p1-1)+(p2-1)+…+(pm-1).

故SST=SS1+SS2+…+SSm,fT=f1+f2+…+fm.對(duì)比文獻(xiàn)[1]和本文中證明的方法，顯然用矩陣知識(shí)證明更為簡(jiǎn)潔.

[1]張曉琴.正交飽和設(shè)計(jì)的統(tǒng)計(jì)分析[J].應(yīng)用概率統(tǒng)計(jì),2007,23(1):91-102.

[2]張應(yīng)山.正交表的數(shù)據(jù)分析及其構(gòu)造[D].上海:華東師范大學(xué)博士論文,2006.

[3]茆詩(shī)松,周紀(jì)薌,陳穎.試驗(yàn)設(shè)計(jì)[M].北京:中國(guó)統(tǒng)計(jì)出版社,2004.

[4]張金槐.線性模型參數(shù)估計(jì)及其改進(jìn)[M].長(zhǎng)沙:國(guó)防科技大學(xué)出版社,1992.

[5]張應(yīng)山.多邊矩陣?yán)碚揫M].北京:中國(guó)統(tǒng)計(jì)出版社,1993.

[6]程云鵬.矩陣論[M].西安:西北工業(yè)大學(xué)出版社,2002.

[7]余錦華,楊維權(quán).多元統(tǒng)計(jì)分析與應(yīng)用[M].中山大學(xué)出版社,2005..

(責(zé)任編輯:王宏志)

2014-01-05

王劍紅(1971-),山西汾陽(yáng)人,碩士,講師.

O

A

1008-7974(2014)01-00032-03