廣義正交表相遇平衡性的等價條件

吳亞楨，楊　林，張應(yīng)山，廖靖宇，田　萍

(1.許昌學(xué)院 數(shù)學(xué)與統(tǒng)計學(xué)院，河南 許昌 461000；2.華東師范大學(xué) 統(tǒng)計學(xué)院，上海 200241)

?

廣義正交表相遇平衡性的等價條件

吳亞楨1，楊林1，張應(yīng)山2，廖靖宇1，田萍1

(1.許昌學(xué)院 數(shù)學(xué)與統(tǒng)計學(xué)院，河南 許昌 461000；2.華東師范大學(xué) 統(tǒng)計學(xué)院，上海 200241)

廣義正交表是一種類似于正交表的新設(shè)計,其數(shù)據(jù)分析保持了正交表的優(yōu)良性,即,各個因子的估計是無偏估計,并且方差達到最小,但試驗次數(shù)大幅減少.相遇平衡性是廣義正交表的兩個平衡條件之一, 也是正交相遇平衡區(qū)組設(shè)計必須滿足的一個基本要求.利用矩陣象技術(shù), 給出了相遇平衡性的一個等價條件, 借助于SAS 軟件可以方便快速的進行驗證.

正交相遇平衡區(qū)組設(shè)計；廣義正交表；相遇平衡；矩陣象

試驗設(shè)計的基本問題是要求試驗的分析結(jié)論具有再現(xiàn)性，即：對于同樣的試驗問題，不同的試驗工作者，無論采取什么樣的合理設(shè)計來收集數(shù)據(jù)，根據(jù)數(shù)據(jù)進行數(shù)據(jù)分析的分析結(jié)論應(yīng)該是唯一的，這種分析結(jié)論不但和試驗工作者的假設(shè)無關(guān)，而且和試驗工作者的操作行為無關(guān).可以證明：正交表試驗和相應(yīng)的數(shù)據(jù)分析結(jié)論具有這種再現(xiàn)性要求[1,2].廣義正交表是在保證試驗具有類似于正交表數(shù)據(jù)分析的基本要求下提出的一種新的設(shè)計，這種設(shè)計可以和正交表一樣用于試驗設(shè)計和相應(yīng)的數(shù)據(jù)分析，也可以證明其數(shù)據(jù)分析結(jié)論是具有再現(xiàn)性的[3-6].

廣義正交表要求設(shè)計表具有兩個最基本的平衡性質(zhì)：正交平衡和相遇平衡.

正交平衡是指多個試驗因子中的各個試驗因子的組合配搭之間的一種平衡.這種平衡性是保證對多個試驗因子中的各個試驗因子之間的處理是公平的，試驗工作者無論對于試驗因子賦予任何內(nèi)容，雖然試驗的數(shù)據(jù)分析結(jié)論將會隨著試驗工作者的賦值予以變動，但不會影響各個試驗因子的分析結(jié)論.如果試驗沒有正交平衡性的要求，那么對多個試驗因子中的各個因子的試驗效果比較將會失去公平，也將不能保證試驗的再現(xiàn)性要求.因此，正交平衡是進行多個試驗因子的試驗設(shè)計必須考慮的一種平衡要求.

相遇平衡是指每個試驗因子的各個試驗水平的組合配搭之間的一種平衡.這種平衡性保證試驗因子對各個試驗水平的處理是公平的，試驗工作者無論對于試驗水平賦予任何試驗內(nèi)容，雖然試驗的數(shù)據(jù)分析結(jié)論將會隨著試驗工作者的賦值予以變動，但不會影響各個試驗水平的分析結(jié)論.如果試驗沒有相遇平衡性的要求，那么對試驗因子的各個水平的試驗比較將會失去公平，也將不能保證試驗的再現(xiàn)性要求.因此，相遇平衡是試驗設(shè)計的最基本的一種平衡要求.但是，這種平衡性的組合判定方法較為復(fù)雜，本文利用矩陣象技術(shù),給出了相遇平衡性的一個等價條件, 借助于SAS軟件可以方便快速的進行驗證.

1　正交相遇平衡區(qū)組設(shè)計的定義

在試驗設(shè)計理論當(dāng)中，具有兩種形式的設(shè)計，其一為行列設(shè)計，如正交表、均勻設(shè)計等，其二為區(qū)組設(shè)計，如平衡不完全區(qū)組設(shè)計、正交拉丁方等.

先從任意形式的區(qū)組設(shè)計考慮起

定義1[3，4]任意具有v個水平b個區(qū)組

B1=(b11,…,b1k1)T,…,Bb=(bb1,…,bbkb)T

的設(shè)計都稱為區(qū)組設(shè)計，這里bij為整數(shù)滿足：1≤bij≤v,j=1,…,ki，而ki為區(qū)組Bi的大小，即向量Bi的維數(shù).記這樣的區(qū)組設(shè)計為

這里k=(k1,…,kb)T為向量.為了以后敘述方便,我們也把區(qū)組設(shè)計

這種形式的區(qū)組設(shè)計僅僅考慮設(shè)計的形狀，并不具有設(shè)計的優(yōu)良性.一個設(shè)計的優(yōu)良性，主要體現(xiàn)在各種各樣的平衡性.

定義3[3，4]將v個水平安排到b個區(qū)組B1,B2,…,Bb的一個區(qū)組設(shè)計Db×k(v)稱為相遇平衡的, 假若它滿足如下條件:

其中λ稱為相遇度，相遇度必須大于0才理解為平衡，這個條件稱為相遇平衡條件.

那么n×(m+1)矩陣

2　區(qū)組設(shè)計相遇平衡性的等價條件

相遇平衡條件的更一般的形式是

記k=(k1,…,kb)T,r=(r1,…,rb)T,那么一個相遇平衡的區(qū)組設(shè)計Db×k(v)都可以用參數(shù)(v,b,r,k,λ)刻畫.常?？紤]k=k11b的情況, 其中1b為元素全是1的b維列向量.例如：設(shè)B1=(1)、B2=(1,2)T,那么B1、B2就組成了一個區(qū)組設(shè)計, 記為

一般地，對于任意的水平數(shù)是2的區(qū)組設(shè)計，只要不是各個區(qū)組內(nèi)僅僅只有一個水平出現(xiàn)，那么相應(yīng)的區(qū)組設(shè)計都是相遇平衡的.這說明相遇平衡的條件是比較弱的一個條件.

引理1對于任意區(qū)組設(shè)計Db×k(v)，我們有

②由于區(qū)組設(shè)計Db×k(v)是相遇平衡的，即

所以

③在ki=k1的條件下，相遇次數(shù)為

相遇度λ(x,y)是常數(shù)當(dāng)且僅當(dāng)相遇次數(shù)λ*(x,y)是常數(shù).

WTW=CTτkC=λvτv，

此時區(qū)組設(shè)計D=Db×k(v)的矩陣象為

證明首先證明：區(qū)組設(shè)計D=Db×k(v)是相遇平衡的當(dāng)且僅當(dāng)CTτkC=λvτv.

我們利用關(guān)聯(lián)矩陣的定義,并結(jié)合相遇平衡性的定義及引理1之②中相遇平衡的等價條件分兩步來證明.

CTτkC的(x,y)號元素為

注意到關(guān)聯(lián)矩陣C=Design(Vec(D))的行和是1,所以C1v=1n,而τk1n=0,這樣得到τkCPvCTτk=0,從而定理1得證.

計算可知

其中Vec是矩陣CTτkC的上三角拉長，RC是Vec的極差，RC=0是區(qū)組設(shè)計為相遇平衡的判定條件.

用定理1可知這個區(qū)組設(shè)計是相遇平衡的.將這個相遇平衡的區(qū)組設(shè)計寫成廣義正交表的形式為

這是一個不飽和的廣義正交表，在此廣義正交表的基礎(chǔ)上，還可以增加一個2水平列，形成飽和的廣義正交表GL6(3121;31).

3　 結(jié)論

一般的廣義正交表GLn(b1v1…vl;vl+1…vm)，其數(shù)據(jù)分析的優(yōu)良性質(zhì)和正交表基本類似，但組合定義性質(zhì)和正交表的定義有較大差別.而用廣義正交表的矩陣象來分析，其矩陣象性質(zhì)和正交表的矩陣象性質(zhì)[7,8]又基本類似.實際上，定理1說明：廣義正交表的組合相遇平衡性質(zhì)是它的矩陣象為投影矩陣的一個充分條件.同時廣義正交表的組合正交平衡性質(zhì)是它的各個列矩陣象為相互正交的矩陣的一個充分條件.因此，廣義正交表是保證它的各個列的矩陣象是相互正交的投影矩陣的一個充分條件.利用矩陣象作為工具是判定一個設(shè)計是否為一個廣義正交表的方便和有效的方法.根據(jù)多元分析的二次型的克赫倫定理可知，這也是保證廣義正交表設(shè)計的數(shù)據(jù)分析結(jié)論具有再現(xiàn)性，相應(yīng)的估計具有無偏和方差最小性質(zhì)的基本保證.因此，我們建議用廣義正交表代替正交表來進行試驗設(shè)計.

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

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

[3]吳亞楨,廖靖宇,張應(yīng)山,等.正交平衡區(qū)組設(shè)計統(tǒng)計分析模型的參數(shù)估計[J].數(shù)學(xué)的實踐與認識,2012,42(2):200-208.

[4]吳亞楨,廖靖宇,張應(yīng)山,等.正交平衡區(qū)組設(shè)計統(tǒng)計分析模型參數(shù)估計的分布特征研究[J].數(shù)學(xué)的實踐與認識,2012,42(3):212-221.

[5]廖靖宇,張建軍,田萍,等.正交平衡區(qū)組設(shè)計矩陣象的基本概念及其基本定理[J].數(shù)學(xué)的實踐與認識,2012,42(17):170-177.

[6]羅純,潘長緣.窮舉法尋找正交平衡區(qū)組設(shè)計[J].應(yīng)用概率統(tǒng)計,2011,27(1):1-13.

[7]Zhang Yingshan, Lu Yiqiang and Pang Shanqi. Orthogonal arrays obtained by orthogonal decomposition of projection matrices[J].Statistica Sinica, 1999, 9: 595-604.

[8]Pang Shanqi, Zhang Yingshan and Liu Sanyang. Normal mixed difference matrix and the construction of orthogonal arrays[J].Stat. & Prob.Lett, 2004, 69: 431-437.

責(zé)任編輯：周倫

The Equivalent Condition of Meeting Balance of Generalized Orthogonal Arrays

WU Ya-zhen1, YANG Lin1, ZHANG Ying-shan2, LIAO Jing-yu1, TIAN Ping1

(1.SchoolofMathsandStatistics，XuchangUniversity，Xuchang461000，China;2.SchoolofStatistics,EastChinaNormalUniversity,Shanghai200241,China)

Generalized orthogonal arrays are similar to the orthogonal arrays,whose data analysis maintains the superiority as orthogonal arrays do, i.e., thatis to say estimation of each factor is unbiased estimation in the sense of variance minimizing, but test tim are remarkably decreased. Meeting balance is one of the two balanced conditions of generalized orthogonal arrays,it is also a basic requirement which orthogonal meeting balanced block design must meet. In this paper,using the matrix images technology, it is to prove an equivalent condition of meeting balance. With the help of SAS software verified quickly and easily.

orthogonal meeting balanced block designs；generalized orthogonal arrays；meeting balance; matrix image

2016-02-12

教育部高校博士點專項基金(44K55050);許昌學(xué)院科研基金(2015083)

吳亞楨(1980—)，男，河南許昌人，講師，碩士，研究方向:應(yīng)用統(tǒng)計，試驗設(shè)計.

1671-9824(2016)05-0021-05

O212.6

A