C 代码实现：计算相关系数

 

资料链接

简单相关性分析(两个连续型变量) - 知乎 (zhihu.com)

参考公式：

//计算相关系数
//A		自变量A数组
//B		自变量B数组
//Length 	数组长度
double Corrcoef(float32_t *A,float32_t *B, u16 Length) 
{
	double sumA, sumB, aveA, aveB;
	u16 i;
	//求和
	sumA=0;
	for(i=0;i<Length;i++)
	{
		sumA+=A[i];
	}
	sumB=0;
	for(i=0;i<Length;i++)
	{
		sumB+=B[i];
	}	
 
	//求平均值
	aveA = sumA / double(Length);
	aveB = sumB / double(Length);
 
	//计算相关系数
	double R1=0, R2=0, R3=0;
	for (long i = 0; i < Length; i++)
	{
		R1 += (A[i] - aveA) * (B[i] - aveB);
		R2 += pow((A[i] - aveA), 2);
		R3 += pow((B[i] - aveB), 2);
	}
 
	return (R1 / sqrt(R2*R3));
}