R/initial_cor3.R

In CorReg: Linear Regression Based on Linear Structure Between Variables

# ' calcul d'une structure initiale bas?e sur les corr?lations
# ' correl=cor(X_appr)
initial_cor3 <- function(correl = correl, p2max = Inf, rmax = Inf) {
  p = ncol(correl)
  Z = Matrix(0, nrow = p, ncol = p)
  qui1 = which(matrix(runif(p^2), ncol = p) < correl)
  qui1 = arrayInd(sample(qui1), c(p, p)) # ordre al?atoire
  for (i in 1:nrow(qui1)) {
    if (qui1[i, 1] != qui1[i, 2] & (sum(Z[, qui1[i, 1]]) == 0 & sum(Z[qui1[i, 2], ]) == 0) & sum(Z[, qui1[i, 2]]) < rmax & (length(which(colSums(Z) > 0)) < p2max | (length(which(colSums(Z) > 0)) == p2max & sum(Z[, qui1[i, 2]]) > 0))) {
      Z[qui1[i, 1], qui1[i, 2]] = 1
    }
  }
  return(Z)
}