2011年第十三届全国大学生数学建模邀请赛C题---期待高人指点~~

赤霄

急需高手指点~~~~
) j/ ~! C& q$ q% N" n

黄窗帘

好东西大家分享

qinniyang

暴力求解有木有

一般不笑

多元约束Welsch-Kuh距离与
3 ~/ F$ U4 o# R3 o/ }多元约束Cook距离
, w7 ^) X) \5 P
唐年胜　韦博成
5 y8 @8 n; m* l8 }# k( n  I
摘　要
4 Z) G. C+ T7 o8 j3 \- |
! T# h8 r  B$ {对多元约束线性回归模型,本文讨论了约束最小二乘估计的影响问题,定义了度量影响的多元约束Cook距离和多元约束W-K统计量,给出了其分布,建立了约束W-K统计量与约束广义相关系数之间的联系。
9 C5 E8 _- z( |关键词　多元W-K统计量,多元Cook距离,影响,约束,广义相关系数
; Z# D+ S& X/ d+ j# v分类号　AMS(1991) 62F;　　　　CCL O212

! x$ b8 n3 s& F1 LMultivariate Restricted Welsch-Kuh Distance and
Multivariate Restricted Cook Distance
" F* Q" V* A2 G' Q1 G, z2 g% e
Tang Niansheng
# k! Y" m4 t9 c* E8 X6 Q7 A' {(Adult Education College,Yunnan University,　Kunming 650091)
Wei Bocheng
(Department of Applied Mathematics,Southeast University,Nanjing 210096)

Abstract
4 _" \$ i2 B( h
2 ]2 z$ q4 r+ h2 Z! `In this paper,the influence of the least squares estimates for the multivariate linear regression models under the restricted conditions is discussed,and multivariate restricted W-K distance and Cook distance is defined as measures of influence.Their distributions are given,and the relationships between the restricted W-K distances and the generalized correlation coefficients are also established.

　　考虑多元线性回归模型

; a# `4 @; D" sYn×p=Xn×mBm×p+En×p,　　　　　　　(1)

, i  M2 s2 [) W0 B2 n# N7 A其中Y=(y1,…,yn)T为因变量的n次观测数据阵,X=(x1,…,xn)T为第一列元素全是1的已知自变量设计阵且rank(X)=m,B=(β1,…,βm)T为回归系数矩阵,E=(ε1,…,εn)T为随机误差阵,ε1,…,εn相互独立,且E(εi)=0,Cov(εi)=V(i=1,…,n),V为正定阵,则由［1］p348知模型(1)中B的最小二乘估计(简记为LSE)为,协方差阵V的无偏估计为,其中P=X(XTX)-1XT,残差阵为。然而,在实际问题中,由于主观和客观的种种原因,人们或多或少对B有些先验信息。即是说B总有一定的约束。
) ?  l4 M# e0 f　　记模型(1)在约束条件AB=0(A为q×m阶行满秩固定阵)下的模型为模型(2)。
5 P4 x, M" B8 Y, u2 Z由［1］知模型(2)中B的LSE为,其中F=AT［A(XTX)-1AT］-1A,协方差阵V的无偏估计为,其中W=YTLY,L=I-P+Q,Q=X(XTX)-1F(XTX)-1XT,残差阵为。本文约定,若A=0,则F=0。
2 r. v2 N6 ~1 x0 }6 Z3 Q　　考虑在模型(1)中剔除k个数据点后所得的一新的多元线性回归模型。
2 D/ W- X: u- j6 A6 `
　　　　Y(I)=X(I)B+E(I),　　　　　　(3)

' _6 G* o* G! r: k  {* f这里Y(I),X(I)和E(I)的意义同［2］p350。
　　记模型(3)在条件AB=0下的模型为模型(4)。
/ \5 A* M! S# C记模型(3)和(4)的相应估计量分别为，，和，，。
　　对模型(1)剔除一组或多组数据对回归系数最小二乘估计的影响问题现已有许多学者做了大量研究工作,并得到了许多很好的结果,见文［2－5］。文［6］对模型(2)在p=1的情形下讨论了约束W-K统计量与广义相关系数之间的联系。本文对多元约束线性回归模型(2)考虑了类似的影响问题,§2将文［2－7］的结果推广到模型(2)的情形,得到了更一般的结果。

1　预备引理

　　引理1　在模型(4)的条件下,若E(I)的各行向量均为正态随机向量,且n＞m+p-q+k,则
5 {) [! t  |# D8 ]0 h; |　　①
, e- C+ g! P& O$ \# n3 n
' o" I# d+ I! @- V; P8 h/ l　　② 相互独立,
　　③

DTX,YI=DTY,D=(di1,…,dik)为n×k阵,dij为第ij个元素是1其余全为0的n维向量。
( g- }( x& u& v  k1 f　　证明　用同［6］完全类似的方法易得①成立。现证②,记R=L(I)Y(I),则



由文［1］知和R独立和Vec(R)独立
7 _7 ]3 B* b* J& H$ I1 q1 ^。而后者是显然成立的,故和R独立,从而与独立,即得②成立。最后,由W1=YT(I)L(I)Y(I),L2(I)=L(I),rank(L(I))=n-m+q-k,L(I)X(I)B=0,及文［1］知③成立。
2 j8 S4 U' E; [1 a& @
证毕

　　引理2　对模型(2)和(4),若记且假设E的各行向量均为正态随机向量,且n＞m+p-q+k,则
　　① W1与W2相互独立,且W=W1+W2,W2～Wp(k,V),W～Wp(n-m+q,V),
　　②
7 r: Y2 w  d' p4 b$ N特别地,当k=1时,
　　证明　记L1=L［I-D(DTLD)-1DT］L,　L2=LD(DTLD)-1DTL,W1=YTL1Y,W2=YTL2Y。则易证W=W1+W2,L=L1+L2。由此即知L1和L2都为投影阵,且rank(L1)=n-m+q-k,rank(L2)=k。因rank(L)=rank(L1)+rank(L2)且LB=L1B=L2B=0,故由多元Cochran定理知①成立。
　　由Wilks分布的定义及①和引理1③知②成立。
, \$ H, m1 I* P6 H, _# @- j
证毕
2 U7 Q' i, p- I% F( o; `
/ [$ r" E! b  y: G4 j& Q2　主要结论
1 w  g+ U4 i3 d# @- u
6 A1 t' O: r- L+ ?/ W& q　　定义1　对模型(2),称

　　　　　　(5)

2 H3 y2 F  q% l8 F, S为k个数据点的多元约束广义Cook距离。对不同的M和C,可得到DHI(M,C)的各种距离和统计意义。
　　定理1　设PI－QI的谱分解式为PI-QI=ΓΛΓT,Γ为正交阵,Λ为对角阵,其对角元素λ1,…,λk为PI-QI的特征根。则在引理2的条件下有

　　　　　　(6)

　　　　　　(7)
8 X! l3 F/ X( i! @$ g
　　　　　　(8)
) c0 U& ?) S2 Z  J$ F5 k; h$ h( I' j
其中h定义为
& l3 Y7 j( i6 ]3 p) b. g2 j
; F- \5 U0 {3 e　　　

5 N- y' L2 z: ~+ gh1,…,hk相互独立且服从p元正态分布。
　　证明　由引理1①和(5)式,并经直接运算即可得定理的结论。
0 F! t! ~, u( c$ z" G
2 g* v7 I, k  W7 ?6 t4 N; Z证毕
& Z: w+ Y5 _+ J9 F1 r; s
　　特别地在定理1中令A=0且k=1,则得
　　推论1　对模型(1),有
+ Q1 }, \# b1 A4 _5 [,其中。
. p' e' E% p6 A+ Z2 M4 N9 s　　这与文［2］、［4－7］的相应结果一致。因此定理1是文［2］、［4－7］相应结果的推广。
* p1 j3 G$ s) p5 y　　因为当k＞1时的分布较复杂,因此,这里只讨论k=1的情形。
　　定理2　在模型(2)的条件下,若ε1,…,εn均为正态随机向量,且n＞m+p-q+1,则,
- Q' ~0 z% K1 r' O
# O( ]2 R# ]: H+ D. T. {都服从Beta分布
7 L) j# G; g+ E/ o+ a5 f. r　　证明　记,则由定理1知上述统计量全等于τ2i／(n-m+q)。因此,我们只须证明。而由引理2知
0 H7 G+ y6 H( ~1 p


' f$ ~3 o! ^$ M5 l3 M* L由此即得。
6 V  x1 ~* [4 u- x
0 E" G: ~: ?  I- w/ w$ \. b2 P证毕

* U) J) A) T' ?& i　　特别地,若A=0,则得
　　推论2　对模型(1)在定理2的条件下,有
, M% A! A% @# E3 \+ T

9 O  ~3 x# L& s' c( g" y
都服从Beta分布。
　　定理3　在定理1的假设下,有
; x  L* s4 H3 E) E( G8 k  c
　　　　　　(9)
7 ]: b" |; F6 C0 L+ N
　　　　　　　(10)

　　　　　　(11)
& _; I+ ]9 w; t& {- e* k0 b9 o
' A! J6 a5 @) ^其中相互独立且都服从Hotelling T2(p,n-m+q-1)。
9 N8 ?! W: q9 O$ \% s$ c6 a7 h& R: m　　证明　由引理1①和(5)式,易得(9)、(10)和(11)式成立。为证定理的后一部分,须证h和相互独立。由h=(I-Λ)-1/2ΓDLY和q-k)YTL1Y,并用同引理1完全类似的方法即得h和相互独立。因此,由引理1③和定理1及文［1］知t*2j～HotellingT2(p,n-m+q-1),j=1,…,k且t*21,…,t*2k相互独立。

证毕
) O& q  y6 K( M& d: R' F+ }# m
+ y" ]+ e+ ~$ m6 x0 N  g　　特别地,当k=1时,有
  P/ T) G" J% ?4 P　　定理4　对模型(2),若ε1,…,εn均为正态随机向量,且n＞m+p-q+1,则



都服从F(p,n-m+q-p)。
　　证明　因为在定理3中令k=1即得上述统计量都等于,其中,故只需证T～F(p,n-m+q-p)。用同文［2］p362完全类似的方法易得与相互独立。因此,由和定理1及文［1］知τ*2i～HotellingT2(p,n-m+q-1)。由T2分布与F分布之间的关系即得

& |, Q+ l, C3 w) @6 a" E8 ^　　　　T～F(p,n-m+q-p)

证毕
8 v" k0 `" M# e3 r. U
　　现考虑k个数据点的剔除对XI处拟合值的影响。由一元线性回归分析方法得
9 X4 Y/ d& t! |; p: O" l1 s# t　　定义2　给定模型(2),k个数据点为XI处拟合值的影响可定义为
: l, e& k5 o0 s. D


) W( W3 P3 p' a+ w; x6 |称WKHI为多元约束Welsch-Kuh统计量。注意,此时WKHI为k×p阶矩阵,为使用方便,考虑W-K统计量的如下形式:
, z8 ?# z% f: }3 |  a! K, ?4 D# J

+ S  E" u" N" b% Y! j0 z- Y
由引理1①和定理3即得
; z$ @2 ?- q* [# f' d
! e" H( ^( u& z9 q$ s7 G
: |0 C" ~# [! W. D. Z' D0 z1 ^
# K5 q" v5 m& J, I, ^8 Z由此可知,多元约束W-K统计量与多元约束Cook距离具有相同的统计意义。这与一元线性回归的情形是一致的。
　　定义3　给定模型(2),k个数据点对XI处拟合值的影响亦可定义为
# u8 `0 T  {2 d/ i) t+ ^+ Q


. @% ]# ^# W# {2 U　　定理5　在定理1的假设下,有



其中GT2为Hotelling广义T2分布。
) Z; H( t" d5 m, f　　证明　记则
# w; \! ~# Y* c. Q' n
　　　　(n-m+q-k)-1WK2HI=tr［(h*Th*)U-1］。

由定理1知,故h*Th*～Wp(m,Ip)。又


/ Y$ I7 ^2 e0 K8 I
) g" U) _# q& W7 A' u且与h*独立,故U～Wp(n-m+q-k,Im)且U与h*Th*相互独立。因此,

　　　　tr［(h*Th*)U-1］～GT2(p,m,n-m+q-k)。
; [/ g& Q4 g, ~
  k& c+ O8 B0 V4 f/ |证毕
8 _7 @* ]: D9 }/ d
# f! S0 g0 G, J4 b3 d0 Q. R# \　　现考虑第i个数据点(yTi,xTi)的剔除对xj(j≠i)处拟合值的影响。
1 b6 |) Q& [  F" F0 E8 b* ^　　定义4　给定模型(2),则第i个数据点对xj(j≠i)处拟合值的影响定义为:
. A7 d# b8 `6 o. y- ~2 k

4 x# \% W3 f5 {6 b
此时,WKTHj为p维向量。因此,考虑其范数
! ?$ h6 P1 C) |0 H+ C
　　　　　(12)
' v7 ?$ S5 f. x% ]3 a& L
7 b! t' V) F8 R7 `　　定理6
7 F, e$ I* K! O7 b! _$ T1 i- g  v7 ]
' S- V; c* ?- W$ Y2 |" y6 B% \  a　　　　　　　(13)

, a$ M/ k& @( t* S$ A3 Y8 B其中ρ(1)H为Hotelling广义相关系数［8］。
　　证明　记pji=xTj(XTX)-1xi,qji=xTj(XTX)-1F(XTX)-1xi,ωji=pji-qji,
则由引理1得
! D+ ?; o' O* G
) B* o4 f& d! a' d
8 i) V0 n. [5 c  X% K) V( ?
将上式代入(12),并令ωjj=pj-qj,ωii=pi-qi,得
7 [' V6 r3 H0 |! _6 @" Q  N& O/ u( `& I
) l0 N3 E. s! O　　　　　　(14)

4 l( {7 L7 `, s& q1 v! K由于与在约束条件AB=0下的联合协方差矩阵为
9 P9 _" p: P0 b3 T$ G4 X6 r
- F8 H' ?8 |( n9 Z
$ s/ i+ K% @8 u3 m
6 F7 c& H5 T; T' Q' I6 l0 R则由［8］P74知与的Hotelling广义相关系数为

7 _% o1 I5 c4 Y( |: o  Q& J　　　　　(15)

8 R: k" w) ]  h/ v5 b# G8 ]6 ^1 ^因此,由(14)和(15)式即得(13)式。
2 J, I! Z/ x  s4 y  D- J5 M; ?6 L
7 r, @5 Q. O! b' Z3 V证毕

- e1 I3 `! ^3 B& X: J* r0 @　　由定理6可知:第一,第i个数据点(yTi,xTi)的剔除对xj(j≠i)处拟合值影响的大小仅依赖于(当i固定时)(j≠i)与的约束Hotelling广义相关系数ρ(1)H,｜ρ(1)H｜越大,其影响WK2Hj也越大,即与在AB=0下的相关程度越好,第i个数据点的剔除对xj(j≠i)处拟合值的影响就越厉害。第二,WK2Hj≤WK*Hi,当j=i时,WK2Hj达到最大值。这与一元约束线性回归［6］的情形是一致的。
* C( h/ F2 S7 E- m/ _; [6 Z9 Q　　致谢:作者感谢审稿同志提出的宝贵意见。
2 ~; g( l3 c. w) k( X5 \" l% i" Q# ?. b
基金项目：国家自然科学基金和云南省自然科学基金资助项目。
作者单位：唐年胜（云南大学成人教育学院,昆明 650091）
　　　　　韦博成（东南大学应用数学系,南京 210096）

! f- g" V: J, h8 V9 A参考文献

: Y8 _6 d2 J2 r9 i) Y, _  q［1］方开泰.实用多元统计分析.上海:华东师范大学出版社,1989;145～159
［2］韦博成,鲁国斌,史建清.统计诊断引论.南京:东南大学出版社,1991;358～362
2 F7 L) P! B8 H" `# {［3］鲁国斌.多元线性回归中度量影响的置信域体积比统计量.高校应用数学学报,1993;8(2):149～156
［4］Hossain A,Naik D N.Detection of influential observations in multivariate regression.Appl Statist,1989;16(1):25～37
0 _) e+ m; e; Y9 w［5］岳珠.多元线性回归中强影响点的判别方法.高校应用数学学报,1987;2(4):343～351
［6］唐年胜,王学仁.约束Welsch-Kuh统计量与广义相关系数.数理统计与应用概率,1998;13(3):239～246
2 w) w5 F# j( W5 p, U& h8 X/ E［7］于义良,吴诗怺.约束Welsch-Kuh统计量与约束Cook距离.应用概率统计,1981;7(2):136～142
  `2 \% b9 r* o" l' B2 ^' U［8］王松桂.线性模型的理论及其应用.安徽,安徽教育出版社,1987;71～74

7 m6 A- D: k* m/ n9 Q

霍水

正在努力求解！！！

qhbsyys

第13届华东数模C题解法推荐
http://www.madio.cn/forum.php?mo ... &fromuid=296190

biqian2012

上面的楼把结果说说呗

lihsky

回复 pbj1203 的帖子
3 ?  p$ m5 v* o4 f& [
求问详细思路

pbj1203

还是没有什么想法啊……