X <- as.matrix(read.table("Cortest.txt"))	#多変量データの読み込みXとする
Y <- cor(X)				#相関行列の計算結果をYとする
Y					#Yを表示
invY <- solve(Y)			#Yの逆行列をinvYとする
O <- diag(1/sqrt(diag(invY)))	#O:invYの対角成分の平方根の逆数の対角行列
P <- -O %*% invY %*% O		#偏相関係数行列をPとする
P					#Pを表示