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

陈华良

数据组：
data awt5_2;
! j: Q1 B* [- tinput x1 x2 @@;
cards;
3   2.30
5   1.90
5   1.00
7   0.70
7   0.30
7   1.00
# w4 }3 ]: g  ~" x8 X& n6 H7 S+ g8   1.05
  w7 N1 n9 g* e+ L' _- V6 B9   0.45
10  0.70
11  0.30
;


* K4 R, O) D6 J3 {


) v2 M& _+ m( R2 [* L- f! a" a  S1）计算每一对观测值到样本均值向量的统计距离平方？
2 Y; e$ t( a8 z3 N( k0 z' X2）计算样本点落在二元正态50％置信区域内的比率。（χ22（0.5）=1.39）
/ F7 a' o) C6 {3）作（X1，X2）的Chi-square图。
4）判断（X1，X2）的二元正态性。

高手请指教哈？？
, S/ C) ]7 h. W7 X7 D' }在此谢过啦
针对第一问！
5 @+ W' Q3 b, D% W9 n( gproc sql;
  select x1, x2, mean(x1) as x1bar, mean(x2) as x2bar,
* O; q; X4 W! G, r8 W* }. v5 }         sqrt( (x1-mean(x1))**2 + (x2-mean(x2))**2 ) as dist
  from awt5_2;
  quit;
& w( Z0 I, C  K
+ C) u- w7 A; ]2 c0 z**********************
        x1        x2     x1bar     x2bar      dist
  ƒƒƒƒƒƒƒƒƒƒƒƒƒƒƒƒƒƒƒƒƒƒƒƒƒƒƒƒƒƒƒƒƒƒƒƒƒƒƒƒƒƒƒƒƒƒƒƒ
         3       2.3       7.2      0.97  4.405553
         5       1.9       7.2      0.97  2.388493
         5         1       7.2      0.97  2.200205
         7       0.7       7.2      0.97  0.336006
         7       0.3       7.2      0.97  0.699214
0 S) I; u! X8 H4 O2 B" |, ^         7         1       7.2      0.97  0.202237
3 A+ [5 ^3 i# N. j1 {. q         8      1.05       7.2      0.97   0.80399
3 V% q' M# k) z# ?8 y; l         9      0.45       7.2      0.97  1.873606
        10       0.7       7.2      0.97  2.812988
        11       0.3       7.2      0.97  3.858614

针对第二问：
5 W- `6 `9 _  c2 T4 }0 Uproc lifetest method=pl width=2;
% o0 a0 A  K; P) y: N* p) }time  x2*x2(0);
4 a* D$ Z; ~+ T1 U0 ~: ?freq x2;
run;
proc lifetest method=pl width=2;
% [( G: _3 X/ b  I2 ltime  x1*x2(0);
& P2 f2 x/ X' b. |freq x2;
run;

# N# P# Y; \* b5 a* H4 K
5 q1 l& q* ]9 n5 E数模中国里头的人，看样子都只会索取资源？不会解决问题啊，太让人失望了........

% X0 P4 i% U6 [6 v6 }" c

---

; L$ T; ]8 _& H9 H# l0 ?4 U解答详情见三楼（数学中国总策划）
9 Q3 f+ _* |2 S' V6 I& s

---

魔力枫叶

厚积薄发

chaoccqiang

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

陈华良

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

0 B! N; A: |5 v  h* m
% |3 r! h6 @/ Y- h3 i上面答案不是很正确，下面给出正确答案，以提供交流探讨
  s1 Z+ a* ~% M/ \- F0 g' q
% y3 l4 n+ k- w  i7 \/ m% @7 u所有程序见下文：
8 C4 ~  B- e3 {. V6 I  m# s
1）初始数据部
options nodate nonumber ls=80;
data awt5_2;
8 I, B7 a: D1 s8 k  ?input x1 x2 @@;
cards;
3   2.30  
- ^7 I1 \, W8 s9 a5   1.90  
5   1.00  
! k5 l7 {8 X; Q$ F+ [7   0.70  
$ n( W" E+ F5 M  {% k- p9 E7   0.30  
7   1.00  
) E, E( E) n- L/ Q* F* L8   1.05  
/ }( v# a! x1 W1 B4 P+ a9   0.45  
. n8 K6 ^, d  u# i, m10  0.70  
11  0.30  
;
: O1 _/ C! C8 r5 `; ~ods html;
! p7 h. z+ u- z( q6 K( Y1 E" e
2）具体操作：
统计距离是计算马氏距离！
proc iml;
' W& Z6 q( ?6 T, h! W- Gn=10;p=2;
xx={x1 x2};
6 |8 Y% d0 L" I; c# Uuse awt5_2;
8 _$ ~: Q4 y3 ]+ Y! Eread all var xx into x;
) ?4 j9 x. ?2 ie={[10] 1};
x0=(e*x)/n;
+ N- E0 Q4 F1 mmm=i(10)-j(10,10,1)/n;
, F" @( z. Z8 x& b& {a=x`*mm*x;
* h( g- N9 |) T& [s=a/(n-1);
3 [# j0 a4 H- _- ssi=inv(s);print x0 s si; /*si为s的逆矩阵*/
/ B: C! n$ S) B0 Buse awt5_2(obs=1);
read all var xx into xx1;
d1=(xx1-x0)*si*(xx1-x0)`; /*d 为马氏距离*/
% \2 M4 }5 Z* c. D; Puse awt5_2(firstobs=2 obs=2);
read all var xx into xx2;
d2=(xx2-x0)*si*(xx2-x0)`;
8 N" c' q' E" M, v, A" P) S* Tuse awt5_2(firstobs=3 obs=3);
read all var xx into xx3;
d3=(xx3-x0)*si*(xx3-x0)`;
: A; O. \8 C. l, H, y7 Ouse awt5_2(firstobs=4 obs=4);
read all var xx into xx4;
d4=(xx4-x0)*si*(xx4-x0)`;
- e) L* a: ?/ l/ P8 U6 H  kuse awt5_2(firstobs=5 obs=5);
read all var xx into xx5;
) l4 s1 t3 G6 h/ c+ J7 x: h: ?d5=(xx5-x0)*si*(xx5-x0)`;
use awt5_2(firstobs=6 obs=6);
8 p  I5 W2 C& ?: \read all var xx into xx6;
d6=(xx6-x0)*si*(xx6-x0)`;
# h" v( H) |+ y9 V# f0 tuse awt5_2(firstobs=7 obs=7);
; M- e0 g& P5 q2 x+ oread all var xx into xx7;
d7=(xx7-x0)*si*(xx7-x0)`;
& p' u) k  s& k* juse awt5_2(firstobs=8 obs=8);
2 T. i( d3 j0 J8 rread all var xx into xx8;
d8=(xx8-x0)*si*(xx8-x0)`;
use awt5_2(firstobs=9 obs=9);
, q7 d& s2 W; B  }/ G& T1 Vread all var xx into xx9;
. R% G7 p5 v) [! T& q3 Bd9=(xx9-x0)*si*(xx9-x0)`;
use awt5_2(firstobs=10 obs=10);
, \5 k: c  K* I$ l0 }read all var xx into xx10;
d10=(xx10-x0)*si*(xx10-x0)`;
print d1 d2 d3 d4 d5 d6 d7 d8 d9 d10;
run;
data md;
$ D8 e! f7 p, P' Minput n d @@;
cards;
' Z* Q  _6 o% ~2 k1  4.0586824
* d7 x3 o5 C1 w8 Q% ?9 C* i; U2  2.1095808
3  2.1074318
2 |  u5 p3 h8 F" `' x( v2 G4 j4  0.6361144
) y) B. y0 a5 ~5  3.2654794
6  0.0079034
7  0.5218616
8  0.6479336
8 D2 [3 S9 h, X. C2 x# ~! z$ f9  2.0590803
! j2 Q" A. R2 A2 I) p10 2.5859323
;
run;
, ?; S! k% V, M8 D8 nproc sort data=md;
7 Z, m0 z8 m6 V. c# Fby d;
run;
proc print data=md;
/ u3 o3 {: y: F8 c7 Q7 C, Arun;
' g$ G) O) ^- l$ E- |proc means  data=md noprint;
var d;
output out=chiqn n=totn;
run;
* R* G) Y" y, k7 N1 b8 qdata chiqq;
; h6 d6 u* R: M5 ^" Sif (_n_=1) then set chiqn;
set md;
novar=2;
" m. m. e% R/ g% h' j5 j! U: Lchisq=cinv(((_n_-0.5)/totn),novar);
prop=0;
$ f1 L) T( L# ]& xd0=cinv(0.5,novar);
6 J$ L+ J; b. T: K- k0 uif d <=d0 then prop=1;
proc univariate data=chiqq;
var prop;
run;
; ]; x* K" f0 r% p9 Z. `proc gplot;
plot d*chisq;
2 Y* [9 r, A; P# Mlabel d='Mahalanobis Distance'
      chisq='Chi-Square Quantile';
; e8 R6 m' M4 G2 \! k8 ?- [          symbol1 v=star;
9 |' J& W2 m/ i: J5 m4 Q0 e1 V          *symbol2 i=join v=+;
run;
+ @3 a3 ?5 M; K& o8 v( m% x4 i) ?: a" p
程序结果：
                                    SAS 系统

1 _( b- _6 j/ {" O( i" y7 d5 H                 x0                   s                  si

( ^* h( B; X8 Z1 Q1 j                7.2      0.97 5.9555556   -1.298889   0.4722924  1.3956265
7 x$ a% o+ u; A                                   -1.298889    0.4395556  1.3956265  6.3991086
" U/ p3 S2 Y1 M8 ^, H2 C
' G1 i0 z* b5 }
% b# \5 M# w  L+ q- q: a6 l        d1               d2            d3                 d4               d5             d6                d7        
4.0586824  2.1095808  2.1074318  0.6361144  3.2654794  0.0079034  0.5218616
/ w' h% @% s& v       d8                d9             d10
3 r# x1 b8 j/ `- @7 m& N5 L; W* Q. g0.6479336  2.0590803  2.5859323
                                    SAS 系统
$ q" z7 Y* z, X7 F7 [
; d( K: P! a' I7 ]; X! C                              Obs     n       d
! O& D$ ]/ G/ N; A- v  h
                                1     6    0.00790
. e' P0 r# a9 O/ [4 W! ?. W                                2     7    0.52186
' J" c( l3 f" f; d7 c1 U; c0 a9 @                                3     4    0.63611
                                4     8    0.64793
. |8 X! R( ~; j$ b- s                                5     9    2.05908
! ~+ s( b9 h. k( C                                6     3    2.10743
                                7     2    2.10958
5 z9 @+ T* u, A1 _5 N0 [                                8    10    2.58593
                                9     5    3.26548
                               10     1    4.05868
7 x% Z" a2 r& W' g6 c5 z2 n                                    SAS 系统

                              UNIVARIATE PROCEDURE
, D: y3 p, B6 n7 P7 ~7 r6 a                                  变量:   prop
- S! H7 l, c# V" O& h
                                       矩

             N                     10    权重总和                10
' j$ L: x& q3 G# [& U0 g% q             均值                 0.4    观测总和                 4
             标准差        0.51639778    方差            0.26666667
             偏度          0.48412292    峰度            -2.2767857
             未校平方和             4    校正平方和             2.4
             变异系数      129.099445    标准误差均值    0.16329932


) }& g* I% _+ K' D                                 基本统计测度

                       位置                    变异性
* t% S# D. s0 J+ e
                 均值     0.400000     标准差          0.51640
                 中位数   0.000000     方差            0.26667
' ~8 E8 s3 C5 y. @) J' N5 ?4 {7 ~7 {                 众数     0.000000     极差            1.00000
& l) _/ T  y- u                                       四分位极差      1.00000


                                位置检验: Mu0=0
/ n2 y5 Y" }6 \3 S; S0 M" W
; r4 I, x, ?+ m$ A: b6 y2 T                  检验      --统计量---    -------P 值-------

                  学生 t    t   2.44949    Pr > |t|    0.0368
5 S# d9 a# O* @, V; _$ }                  符号      M         2    Pr >= |M|   0.1250
                  符号秩    S         5    Pr >= |S|   0.1250


$ r) c( M* Z" k+ S$ k                                分位数（定义 5）

                             分位数          估计值
& k8 h" [8 Z4 ?, Y
; n' k( B9 v; Y3 m5 I: n                             100%  最大值    1
4 \! ~9 ^/ m* I& {, b                             99%                   1
                             95%                   1
9 ]* J, Z; `. A! ~, t$ D# Q& r                             90%                   1
2 f  @% W; M0 I' x7 ?7 r- K, }                             75% Q3             1
3 n- V* Y/ H6 d7 i7 u1 q                             50%  中位数       0
) n( [  O. B1 F3 w2 ?' f/ n' \                             25% Q1              0
                             10%                   0
9 m& p/ J1 f, ~, n7 W, ?3 R                             5%                     0
5 ~+ D! V; n% J                             1%                     0
                             0%  最小值         0


                                    极值观测

                         ---最小值--        ---最大值--

7 H/ P+ H- A' T9 ]% |" x, K                         值     观测        值     观测

                          0       10         0       10
5 t$ c& V  w& j8 J* N) [! L4 n0 s                          0        9         1        1
7 J* g2 N" R. ?" d) O                          0        8         1        2
                          0        7         1        3
: I) j6 }1 S* `9 u; u                          0        6         1        4

# ~  U6 _. }# T

, {& l2 [. @/ u( h! U

厚积薄发

陈华良 发表于 2012-4-17 13:34
; l6 C! Y, s, p- e! _* u' F上面答案不是很正确，下面给出正确答案，以提供交流探讨

所有程序见下文：

6 J" b- A% c" h0 S- ~错误的话，指出错误所在

陈华良

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

没说错误，只是有点答非所问的感觉，正确答案详见5楼

（一路向北）

看不懂，请解析一下

sdccumcm

胖哥

感觉好复杂啊