数模中国里头的人，是英雄？还是虚伪？

陈华良

! X. F! C/ `: Q数据组：
data awt5_2;
( p- X# b: a  C% Ainput x1 x2 @@;
% p1 c0 L2 X$ `. K1 q0 C4 |cards;
3   2.30
5   1.90
/ p# i+ J8 O% N, O7 w5   1.00
& c4 |. }2 j* U8 y9 }  m$ L4 x7   0.70
7   0.30
/ l" j2 q( D# h. C$ O7   1.00
8   1.05
( z! n3 |8 c3 E2 o9   0.45
. g2 W' V. A8 a10  0.70
11  0.30
# [4 p+ m5 Z" g+ D% u) x' M) L% X;
% J3 o  ^+ J4 @4 \" a! C9 m
% A& u3 y; |9 j

( _' Z% M- M, E  y) W5 ^6 i

/ L1 C( w5 S! k) J1）计算每一对观测值到样本均值向量的统计距离平方？
2）计算样本点落在二元正态50％置信区域内的比率。（χ22（0.5）=1.39）
3）作（X1，X2）的Chi-square图。
4）判断（X1，X2）的二元正态性。

* ^  ]8 H6 d# p* N9 _, H& ^6 Q, w3 x, V高手请指教哈？？
# J  J: P7 D2 ~5 e+ U, R在此谢过啦
( \5 k3 h6 [6 S5 J针对第一问！
proc sql;
& O' R# s1 w" K- a$ Y! h  select x1, x2, mean(x1) as x1bar, mean(x2) as x2bar,
+ a3 e1 V- ^. k$ ]" E# n( L( K         sqrt( (x1-mean(x1))**2 + (x2-mean(x2))**2 ) as dist
  from awt5_2;
  quit;
* P- D8 w, l* T2 ?9 d+ G! m
**********************
2 a; z, s6 C6 A$ ?        x1        x2     x1bar     x2bar      dist
# f/ B5 s" S8 q  ƒƒƒƒƒƒƒƒƒƒƒƒƒƒƒƒƒƒƒƒƒƒƒƒƒƒƒƒƒƒƒƒƒƒƒƒƒƒƒƒƒƒƒƒƒƒƒƒ
         3       2.3       7.2      0.97  4.405553
         5       1.9       7.2      0.97  2.388493
" h4 k8 r/ M$ A3 f0 q; \5 r         5         1       7.2      0.97  2.200205
. k7 m- W5 V/ Y1 M, r- o$ V( b8 i# L         7       0.7       7.2      0.97  0.336006
* x1 ^5 [8 y/ {) w         7       0.3       7.2      0.97  0.699214
' m/ T' U: H- w0 H. t         7         1       7.2      0.97  0.202237
7 x) Z- t6 S5 {" R! T9 e         8      1.05       7.2      0.97   0.80399
         9      0.45       7.2      0.97  1.873606
        10       0.7       7.2      0.97  2.812988
5 s* R! D4 T6 W; D        11       0.3       7.2      0.97  3.858614
* h# ]* A; q' d0 }- B* q) u
: Y* H) M' x) Y0 n针对第二问：
7 O) n, a; C0 ^# I0 U5 w  tproc lifetest method=pl width=2;
time  x2*x2(0);
, I$ K: u4 D. |8 Qfreq x2;
9 P9 F: E0 i- t; Z* `6 Orun;
proc lifetest method=pl width=2;
5 `3 h, v/ Z9 ]2 z) dtime  x1*x2(0);
! z. p2 P( K0 E! o9 yfreq x2;
run;


数模中国里头的人，看样子都只会索取资源？不会解决问题啊，太让人失望了........

" u- _+ _/ Z7 D0 q# L

---

  f: E4 k  Z8 N. I# o3 g5 e  u% N解答详情见三楼（数学中国总策划）

---

2 g& V$ P  ~8 k+ k0 t1 P

魔力枫叶

厚积薄发

chaoccqiang

淡定，做事要从别人角度出发考虑！

陈华良

; F0 c$ \: w% @8 z

厚积薄发 发表于 2012-4-12 11:55

8 c- ~- X6 u  U; a2 X
上面答案不是很正确，下面给出正确答案，以提供交流探讨
: G! w) i% p% l( e
所有程序见下文：
) j' |' K6 |& U. {: b3 C) n
5 d$ D: B4 c. A  ?( c9 N# O; M+ \1）初始数据部
6 Y: G; I& O4 a3 R- w1 ]options nodate nonumber ls=80;
! @7 l$ H$ {& f3 G' P+ Pdata awt5_2;
input x1 x2 @@;
) l( ]$ P! ~( H  A% H: E& [5 Ycards;
4 R) l# |7 H; I, Z  w7 c' l1 V4 ?3   2.30  
5   1.90  
5   1.00  
7   0.70  
7   0.30  
7   1.00  
8   1.05  
9   0.45  
& U* a3 _( {9 L8 O( |10  0.70  
, j) F& J: |0 J' i11  0.30  
;
ods html;
4 _  \$ j( _5 b* L! L1 n
2）具体操作：
. I$ d) }* F: Y2 _( u统计距离是计算马氏距离！
' b7 {% A3 ?) |% D  e! x: u) ~proc iml;
n=10;p=2;
& O5 Z2 K! \! p* I' u- w" @* cxx={x1 x2};
use awt5_2;
read all var xx into x;
e={[10] 1};
x0=(e*x)/n;
mm=i(10)-j(10,10,1)/n;
a=x`*mm*x;
s=a/(n-1);
) H7 \2 W# q' {4 B% T+ P" u( o( }5 esi=inv(s);print x0 s si; /*si为s的逆矩阵*/
, K8 z) `9 F& S2 @  Xuse awt5_2(obs=1);
read all var xx into xx1;
' f5 h6 }0 h, @% ed1=(xx1-x0)*si*(xx1-x0)`; /*d 为马氏距离*/
use awt5_2(firstobs=2 obs=2);
4 N' K& F# O$ E9 qread all var xx into xx2;
d2=(xx2-x0)*si*(xx2-x0)`;
% ?) Q( @: T4 _8 H& C: huse awt5_2(firstobs=3 obs=3);
0 S: }3 y# ?! J  e. E  m; Jread all var xx into xx3;
" T& y/ Y: C3 B$ Q$ Q5 B# |5 _d3=(xx3-x0)*si*(xx3-x0)`;
* n& p4 @8 F, v9 ?use awt5_2(firstobs=4 obs=4);
read all var xx into xx4;
d4=(xx4-x0)*si*(xx4-x0)`;
use awt5_2(firstobs=5 obs=5);
2 v0 l7 p, I, l$ A& hread all var xx into xx5;
d5=(xx5-x0)*si*(xx5-x0)`;
use awt5_2(firstobs=6 obs=6);
read all var xx into xx6;
, P0 w5 A/ o. W9 _2 Y2 |0 r0 K$ ^d6=(xx6-x0)*si*(xx6-x0)`;
use awt5_2(firstobs=7 obs=7);
read all var xx into xx7;
& {9 J4 j5 e. K4 c0 ~d7=(xx7-x0)*si*(xx7-x0)`;
" U4 P% ~2 n7 ]' [/ p5 fuse awt5_2(firstobs=8 obs=8);
read all var xx into xx8;
/ x4 m8 l3 i. r4 E: qd8=(xx8-x0)*si*(xx8-x0)`;
use awt5_2(firstobs=9 obs=9);
0 ?2 ]3 }5 \+ s1 J! r2 q3 @) qread all var xx into xx9;
d9=(xx9-x0)*si*(xx9-x0)`;
use awt5_2(firstobs=10 obs=10);
read all var xx into xx10;
& d  D' X' |, r- r! \d10=(xx10-x0)*si*(xx10-x0)`;
print d1 d2 d3 d4 d5 d6 d7 d8 d9 d10;
run;
9 A  y. @3 _. K; D* D( X5 V# ddata md;
8 u! _; _7 b+ Ninput n d @@;
0 z; |( N- f$ i! T' J! ocards;
1  4.0586824
1 y. z5 t! L+ C! p2  2.1095808
+ H- j2 B0 r' `9 g3  2.1074318
9 {- L1 w! |9 s! s4  0.6361144
. }+ V  Q0 n* K7 y* ^5  3.2654794
6  0.0079034
7 L! b4 |4 v  \- ~, U7  0.5218616
8  0.6479336
$ W7 T9 [+ v" u, D5 K7 E# R9  2.0590803
10 2.5859323
;
5 ?* z% O5 w& W$ G+ Qrun;
proc sort data=md;
; p$ e: f! m: V$ C; v+ G5 m5 |by d;
run;
" k$ ^( \, f- t" n6 W8 }% uproc print data=md;
& E6 x( o0 U7 ]/ ]" erun;
* h6 D- Z7 C2 vproc means  data=md noprint;
var d;
. _! [' B; I1 i9 N, i4 K5 qoutput out=chiqn n=totn;
run;
# o/ ~% t  ^, K- fdata chiqq;
; B6 q( q3 t2 O5 pif (_n_=1) then set chiqn;
) o& C6 x" r: t/ A" rset md;
- C, J+ H$ @' h* J; f9 A2 lnovar=2;
$ ^$ l6 P& ]/ K/ T* N, I: Ychisq=cinv(((_n_-0.5)/totn),novar);
# s& m$ n9 B: g  [prop=0;
- _" [  c6 c3 T2 b8 u1 y" V4 Id0=cinv(0.5,novar);
if d <=d0 then prop=1;
proc univariate data=chiqq;
var prop;
/ h( B9 t" g# v1 Crun;
+ U, d" p5 ]$ Y1 B3 T( xproc gplot;
plot d*chisq;
3 i( I+ K1 S8 plabel d='Mahalanobis Distance'
1 `) z. P; C# K9 B, J! V1 `; c/ W      chisq='Chi-Square Quantile';
          symbol1 v=star;
8 ^  h) _, Z( s, G" z. i          *symbol2 i=join v=+;
+ q1 e! i8 X* k  Crun;
4 X7 R$ J1 Y9 O3 y* I0 a
程序结果：
4 r8 Y& w# N" A2 X6 e' }                                    SAS 系统
) F* G8 H* z# j/ ~6 [) k5 e
                 x0                   s                  si
% F: Y" j6 k0 ?" N4 F9 O
& T& L2 p! Q( x: C                7.2      0.97 5.9555556   -1.298889   0.4722924  1.3956265
                                   -1.298889    0.4395556  1.3956265  6.3991086

1 |" T5 s9 Z* W
$ Y- r9 {0 M) ?        d1               d2            d3                 d4               d5             d6                d7        
! A5 m2 x! k! `% ^# ?( T( s# g4.0586824  2.1095808  2.1074318  0.6361144  3.2654794  0.0079034  0.5218616
) e9 \! I$ _7 E2 a+ ^       d8                d9             d10
0.6479336  2.0590803  2.5859323
$ u! }; `$ q* E9 s                                    SAS 系统

' O  h7 p9 I9 \4 ]8 d                              Obs     n       d

$ q; P; q1 ~5 @: o0 |/ D                                1     6    0.00790
                                2     7    0.52186
                                3     4    0.63611
                                4     8    0.64793
! I: n1 y( G. h& L! A, c8 ~2 D2 q                                5     9    2.05908
                                6     3    2.10743
9 S, G* F( ]# g! p& `# }                                7     2    2.10958
                                8    10    2.58593
/ z2 k# a0 A" F! \& g3 @, Y; N                                9     5    3.26548
                               10     1    4.05868
                                    SAS 系统

7 Q; G$ Q% q2 P                              UNIVARIATE PROCEDURE
: }* M# Z4 M( ^! w7 U* o                                  变量:   prop
7 w: d+ I+ G) k% J5 K
                                       矩

2 p! Z  N' d2 ^! e- @: ~* n             N                     10    权重总和                10
7 w9 o+ z( h" }6 j% J. x  [             均值                 0.4    观测总和                 4
. F" O4 h; _2 D             标准差        0.51639778    方差            0.26666667
4 q3 D9 X! {* R             偏度          0.48412292    峰度            -2.2767857
$ c: G  }3 n8 R& G2 e             未校平方和             4    校正平方和             2.4
             变异系数      129.099445    标准误差均值    0.16329932
$ F! c- S7 v- z  B: q

                                 基本统计测度

                       位置                    变异性
% p$ X( D* }. k# M( _. g0 l
. f0 |# Y, ?* N1 O) y( W/ `9 Y                 均值     0.400000     标准差          0.51640
                 中位数   0.000000     方差            0.26667
* r6 |/ g  o' _3 @) H8 M% J, Q                 众数     0.000000     极差            1.00000
: ?# `3 t* D% G% j                                       四分位极差      1.00000

9 X5 h' e" F  Y. ^" b% T
                                位置检验: Mu0=0

                  检验      --统计量---    -------P 值-------

7 H& s% l8 E6 ]* e0 c7 O                  学生 t    t   2.44949    Pr > |t|    0.0368
                  符号      M         2    Pr >= |M|   0.1250
+ j9 o' w/ R7 T/ u4 A* Q                  符号秩    S         5    Pr >= |S|   0.1250

/ u+ i( D' J$ l. A2 J+ r
& e# N2 j8 x6 a& q1 [                                分位数（定义 5）

4 r" G7 |! d" p) z9 A                             分位数          估计值
# u  [) O2 n# d9 j( O, H+ k
' H/ E  o" J  d2 h' D                             100%  最大值    1
4 n  ^# n) G7 m# y* t3 d                             99%                   1
                             95%                   1
! l. w# x5 T2 P  v$ Q( j) I                             90%                   1
                             75% Q3             1
                             50%  中位数       0
: @" G/ K0 @. L  B$ n+ V; Z5 }6 z                             25% Q1              0
                             10%                   0
                             5%                     0
- \3 A3 h8 m: P3 n* D9 l" t                             1%                     0
                             0%  最小值         0
9 L& H  h  w$ k
) d7 h0 S. D5 Y  w4 `
                                    极值观测

$ @1 N8 a3 y! r9 e- f  B) H. O                         ---最小值--        ---最大值--

; s) H- e' q$ {                         值     观测        值     观测

/ `, c7 [9 n& a9 L6 x0 n0 D2 r4 D                          0       10         0       10
                          0        9         1        1
, c6 y% n; J* B. N; w                          0        8         1        2
                          0        7         1        3
0 Z6 F6 ?, s9 P% w7 o% y                          0        6         1        4
) [$ Z7 O% F3 c4 A
. |9 j  I+ B* F+ I& C

厚积薄发

陈华良 发表于 2012-4-17 13:34
& E6 @# L- Q" E9 }3 w上面答案不是很正确，下面给出正确答案，以提供交流探讨
: H0 b$ V* w, n6 S& v; Z0 n! L
所有程序见下文：

  F, @9 |: e; B/ O错误的话，指出错误所在

陈华良

厚积薄发 发表于 2012-4-17 13:35
错误的话，指出错误所在

' w2 R: ?5 h- V# a8 _7 v没说错误，只是有点答非所问的感觉，正确答案详见5楼

（一路向北）

看不懂，请解析一下

sdccumcm

pcyaoqiang

要换位思考。