####################SimSIP statistics
SimSIP<-function(r1,r2)
{
  hF<-1/((r1)^0.5)
  hG<-1/((r2)^0.5)
  Texs<-hF%*%hG
  return(Texs)
}