shfujim
P=matrix(c(0.0529,0.00621,0.00621,0.0225),nrow=2,ncol=2,byrow=T)
P
G=matrix(c(0.02116,0.004632,0.02277,0.004632,0.0099,0.01068),nrow=2,ncol=3,byrow=T)
G
C=matrix(c(0.02116,0.004632,0.02277,0.004632,0.0099,0.01068,0.02277,0.01068,0.05),nrow=3,ncol=3,byrow=T)
C
V=matrix(c(370,-20,-50),nrow=3,ncol=1,byrow=T)
b=solve(P)%*%G%*%V
b
rHI=sqrt((t(b)%*%G%*%V)/(t(V)%*%C%*%V))
rHI
Sh=(t(b)%*%G%*%V)/sqrt(t(b)%*%P%*%b)
Sh
Sg=(t(b)%*%G)/sqrt(t(b)%*%G%*%V)
误于(t(b) %*% G)/sqrt(t(b) %*% G %*% V) : 非整合陈列
麻烦大虾帮忙看看，先谢谢了^_^
