相关系数计算.cpp

2个关于概率统计的C程序

#include "math.h"
   ////////////////////////////////////////////////////////////////
   //   n整数,数据的组数;                                        //
   //   aa实型数组,存放给定的n组数据中Xi的值;                    //
   //   bb实型数组,存放给定的n组数据中yi的值;                    //
   //   函数返回相关系数,计算结果返回序列X,Y的相关系数           //
   //   调用方式:  lirskn(n,aa,bb)                               //
   ////////////////////////////////////////////////////////////////
float lirskn(int n,float aa[],float bb[])
{  float sx,sy,sxx,syy,sxy,rr;
   int i;
   
   sx= 0.0;  sy= 0.0;
   sxx=0.0;  syy=0.0;
   sxy=0.0;
   for (i=0;i<n;i++)
      {  sx=sx+aa[i];
         sy=sy+bb[i];
         sxy=sxy+aa[i]*bb[i];
         sxx=sxx+aa[i]*aa[i];
         syy=syy+bb[i]*bb[i];
         }
   rr=sqrt((sxx-sx*sx/n)*(syy-sy*sy/n));
   rr=(sxy-sx*sy/n)/rr;
   return rr;
   }