vech.mata

是一个经济学管理应用软件 很难找的 但是经济学学生又必须用到

*! version 1.0.0  15oct2004
version 9.0
mata:

transmorphic colvector vech(transmorphic matrix x)
{
	transmorphic matrix	res
	real scalar		n, k, j

	if ((n=rows(x)) != cols(x)) _error(3205) // not square

	res = J((n*(n+1))/2, 1, missingof(x))
	if (n>0) {
		for (k=j=1; j<=n; j++) {
			res[|k\(k+n-j)|] = x[|(j,j)\(n,j)|]
			k = k + n-j + 1
		}
	}
	return(res)
}

end