correlation.mata

来自「是一个经济学管理应用软件 很难找的 但是经济学学生又必须用到」· MATA 代码 · 共 29 行

*! version 1.0.0  15oct2004
version 9.0
mata:

real matrix correlation(real matrix X, real colvector w)
{
	real rowvector	CP
	real rowvector	means
	real matrix	res
	real scalar	i, j
	real scalar	n 

	CP = quadcross(w,0, X,1)
	n  = cols(CP)
	means = CP[|1\n-1|] :/ CP[n]
	res = crossdev(X,0,means, w, X,0,means)

	for (i=1; i<=rows(res); i++) {
		res[i,i] = sqrt(res[i,i])
		for (j=1; j<i; j++) {
			res[i,j] = res[j,i] = res[i,j]/(res[i,i]*res[j,j])
		}
	}
	for (i=1; i<=rows(res); i++) res[i,i] = 1
	return(res)
}

end