久久91这里精品国产2020,人妻无码精品久久亚瑟影视 ,亚洲国产精品无码久久久久久曰

如何一遍計(jì)算相關(guān)系數(shù)

一般在計(jì)算相關(guān)系數(shù)時(shí)都要遍歷4次數(shù)據(jù)，這在數(shù)據(jù)比較多，需要讀寫文件時(shí)，顯然是費(fèi)時(shí)的。這里給出Single Pass版本的計(jì)算方法,但是其中的變化依據(jù)還沒有搞懂，希望各位指教:
偽代碼:

sum_sq_x = 0
2

sum_sq_y = 0
3

sum_coproduct = 0
4

mean_x = x[1]
5

mean_y = y[1]
6

for i in 2 to N:
7

sweep = (i - 1.0) / i
8

delta_x = x[i] - mean_x
9

delta_y = y[i] - mean_y
10

sum_sq_x += delta_x * delta_x * sweep
11

sum_sq_y += delta_y * delta_y * sweep
12

sum_coproduct += delta_x * delta_y * sweep
13

mean_x += delta_x / i
14

mean_y += delta_y / i
15

pop_sd_x = sqrt( sum_sq_x / N )
16

pop_sd_y = sqrt( sum_sq_y / N )
17

cov_x_y = sum_coproduct / N
18

correlation = cov_x_y / (pop_sd_x * pop_sd_y)
19

下面是我用C改寫的，這份改寫同時(shí)放在了Wiki上面:

#include "stdlib.h"
2

#include "math.h"
3

void SinglePassCorrelation(FILE *file,int N)
4

{
5

double *x=new double[N+1];
6

double *y=new double[N+1];
7

int i;
8

for(i=1;i<=N;i++)
9

{
10

fscanf(file,"%lf %lf",&x[i],&y[i]);
11

printf("%10.6f %10.6f\n",x[i],y[i]);
12

}
13

double sum_sq_x=0;
14

double sum_sq_y=0;
15

double sum_coproduct=0;
16

double mean_x=x[1];
17

double mean_y=y[1];
18

double sweep=0;
19

double delta_x;
20

double delta_y;
21

double pop_sd_x;
22

double pop_sd_y;
23

double cov_x_y;
24

double correlation;
25

for(i=2;i<=N;i++)
26

{
27

sweep=(i-1.0)/i;
28

delta_x=x[i]-mean_x;
29

delta_y=y[i]-mean_y;
30

sum_sq_x+=delta_x*delta_y*sweep;
31

sum_sq_y+=delta_y*delta_y*sweep;
32

sum_coproduct+=delta_x*delta_y*sweep;
33

mean_x+=delta_x/i;
34

mean_y+=delta_y/i;
35

}//end for i
36

pop_sd_x=sqrt(sum_sq_x/N);
37

pop_sd_y=sqrt(sum_sq_y/N);
38

cov_x_y=sum_coproduct/N;
39

correlation=cov_x_y/(pop_sd_x*pop_sd_y);
40

}
41

Posted on 2008-11-28 20:13 鄒敏閱讀(2114) 評(píng)論(0) 編輯收藏引用

只有注冊(cè)用戶登錄后才能發(fā)表評(píng)論。
【推薦】100%開源！大型工業(yè)跨平臺(tái)軟件C++源碼提供，建模，組態(tài)！



網(wǎng)站導(dǎo)航: 博客園 IT新聞 BlogJava 博問 Chat2DB 管理

2025年7月

日

一

二

三

四

五

六

常用鏈接

留言簿(1)

隨筆檔案

文章檔案

記

如何一遍計(jì)算相關(guān)系數(shù)

常用鏈接

留言簿(1)

隨筆檔案

文章檔案

面試題

搜索

最新評(píng)論

閱讀排行榜

評(píng)論排行榜