判别分析程序

cuijxin

大家好 工作需要 我在做一个气象统计分析与预报方法的 算法 参考的教材是《气象统计分析与预报方法》这本书 是黄嘉佑编著的 参考里面的第三章 判别分析 费歇尔判别准则 给出的判别方法 但是不是很明白 能不能请大家给我讲解一下啊 我看到论坛里有人发 “气象常用程序包-相当给力”的帖子 里面的第三章 多元统计分析里有Fortran代码 但是读来很迷糊 因为刚开始接触 所以不熟 希望大家给予帮助

mofangbao

这怎么讲解。。。上课么。。。

cuijxin

那段Fortran程序 调试就不能通过 有很多错误 我改了点代码 调试通过了 但是输出的结果还是有问题  

下面是我调试时输出的结果：

总的均值=   -8.10  -17.36   16.47    7.53

第 1组的均值=  -14.43  -17.34   12.71   31.14
第 2组的均值=     .80  -17.43   17.50     .00
第 3组的均值=   -6.65  -17.33   20.17  -15.00

总的离差阵 T
29652.28   898.14   654.50  9566.50
   898.14   215.46   193.57  1976.03
   654.50   193.57   558.24   283.76
  9566.50  1976.03   283.76 28070.24

组内的离差阵 S
29042.47   900.98   419.29 11076.63
   900.98   215.44   193.70  1974.84
   419.29   193.70   373.26  1435.29
11076.63  1974.84  1435.29 20894.86

组间离差阵 B
   609.81    -2.84   235.21 -1510.13
    -2.84      .02     -.13     1.19
   235.21     -.13   184.97 -1151.52
-1510.13     1.19 -1151.52  7175.38

     特征值
  -.1651E-15   .2540E-15

     特征向量
   .5868E-03  -.2102E-03
   .5215E-01   .1919E+00
   .5066E-02  -.7823E-01
   .9279E-03  -.1263E-01
NP1= 1
     X2检验
X2= 17.6860
X22= 15.1617
    样品与判别函数中心的距离
kind= 1  .000
kind=

=============================================
这里是文档里给出的结果：

总的均值=   -8.10  -17.36   16.47    7.53

第 1组的均值=  -14.43  -17.34   12.71   31.14
第 2组的均值=     .80  -17.43   17.50     .00
第 3组的均值=   -6.65  -17.33   20.17  -15.00

总的离差阵T
29652.28   898.14   654.50  9566.50
   898.14   215.46   193.57  1976.03
   654.50   193.57   558.24   283.76
  9566.50  1976.03   283.76 28070.24

组内的离差阵S
29042.47   900.98   419.29 11076.63
   900.98   215.44   193.70  1974.84
   419.29   193.70   373.26  1435.29
11076.63  1974.84  1435.29 20894.86

组间离差阵B
   609.81    -2.84   235.21 -1510.13
    -2.84      .02     -.13     1.19
   235.21     -.13   184.97 -1151.52
-1510.13     1.19 -1151.52  7175.38

     特征值
   .3116E+01   .1239E-01

     特征(列)向量
   .2638E-02  -.6104E-02
   .1451E+00   .2850E-01
   .1264E-01   .6449E-02
  -.1828E-01  -.2678E-03
NP1= 1
      χ2检验
X2= 17.6861
X22= 15.1617
    样品与判别函数重心的距离
序号       1     2     3     4     5    6     7     8     9    10    11   12
组别=1   .009  1.497  .800   .002  .297  .739  1.471  .108  .084   .002  .250  .480
组别=2   .566   .320  .056   .497  .013  .041   .308  .109  .135   .492  .025  .001
组别=3  1.111   .070  .004  1.012  .172  .010   .064  .398  .448  1.006  .211  .071  
最后分组  1     3    3(×)   1     2     3     3     1     1   1(×)   2     2
序号      13     14   15    16     17  
组别=1   .016   .129  .008  1.047  1.317
组别=2   .283  1.035  .323   .133   .240
组别=3   .695  1.739  .757   .004   .035
最后分组  1     1     1     3     3

cuijxin

我开始的时候 也只是想按照Fortran代码 写出C++代码 但是这个Fortran程序好像本身就有问题 至少文档中给出的程序是有问题的 我就只能去看书 看书里给出的计算方法 但是书里有好多数学知识 我都忘了 有人给我再上一课 我是求之不得的啊

mofangbao

cuijxin 发表于 2012-4-9 11:29

                               
登录/注册后可看大图

我开始的时候 也只是想按照Fortran代码 写出C++代码 但是这个Fortran程序好像本身就有问题 至少文档中给出的 ...

看看统计学的课件 可能比书本容易懂

topmad

如果专门为了做出统计结果 还是用sas spss 或者R这些专门的统计软件吧  毕竟要的是结果，把主要精力放在对结果的分析上

cuijxin

topmad 发表于 2012-4-9 13:27

                               
登录/注册后可看大图

如果专门为了做出统计结果 还是用sas spss 或者R这些专门的统计软件吧  毕竟要的是结果，把主要精力放在对结 ...

很遗憾 不单是要做出结果 而且还要掌握计算过程

topmad

cuijxin 发表于 2012-4-9 14:37

                               
登录/注册后可看大图

很遗憾 不单是要做出结果 而且还要掌握计算过程

那就先把所有的计算过程都搞清楚，然后根算法看代码 并找出其中的可能的错误

cuijxin

topmad 发表于 2012-4-9 16:49

                               
登录/注册后可看大图

那就先把所有的计算过程都搞清楚，然后根算法看代码 并找出其中的可能的错误

大家一起帮我研究啊 我个人能力有限呐
这是书上写的关于多级判别计算步骤
第一步：选择适当因子，并根据预报量类别确定不同类别的样本，计算各组因子的平均值和总平均值；
第二步：计算总离差交叉积阵T，组内离差交叉积阵W及组间离差交叉积阵B；
第三步：求W（-1）B的特征向量，得V阵；
第四步：对判别函数进行显著性检验以便确定选取多少个判别函数构成判别空间；
第五步：计算个样品点与各组重心距离并进行分类判别
第六步：计算每一样品和各组判别函数重心的距离。例如计算第一个样品与第一组判别函数重心距离

现在前三步都得出了结果

beesfly

帮你顶顶  偶不会这个东西