R/vp_rv_Wi.R

In ClustBlock: Clustering of Datasets

.vp_rv_Wi=function(Wi){
  # Data size (n, all the var)
  # Blocks vector given the nb of var of each block

  #initialisation
  n=dim(Wi)[1]
  nblo=dim(Wi)[3]

  # RV matrix:
  RV=matrix(0,nblo,nblo)
  diag(RV)=rep(1,nblo)
  if(nblo>1)
  {
    for (i in 1:(nblo-1)) {
      for (j in (i+1):nblo) {
        RV[i,j]=sum(diag(crossprod(Wi[,,i], Wi[,,j])))
        RV[j,i]=RV[i,j]
      } }
  }


  #eigenvalues of RV matrix
  ressvd=svd(RV)
  lambda=ressvd$d
  return(lambda[-nblo])
}