R/generateurZ.R

In CorReg: Linear Regression Based on Linear Structure Between Variables

# ' generateur de Z aleatoire
# MH seulement mais all?g? pour ?tre plus rapide
generateurZ <- function(p = p, maxit = NULL) {
  if (is.null(maxit)) {
    # maxit=round(p^2/4)
    maxit = floor(runif(1, min = 0, max = p^2))
  }
  Z = Matrix(0, ncol = p, nrow = p)
  candidats = (1:p^2)[-seq(1, p^2, p + 1)] # tout sauf la diagonale
  # on initialise
  for (i in 1:maxit) {
    qui = arrayInd(sample(candidats, size = 1), c(p, p)) # ordre al?atoire des candidats
    i = qui[1]
    j = qui[2]
    Z[, i] = Z[i, j] * Z[, i] # annule si ? 0 (donc si on le passe ? 1) neutre sinon
    Z[j, ] = Z[i, j] * Z[j, ] # annule si ? 0 (donc si on le passe ? 1) neutre sinon
    Z[i, j] = 1 - Z[i, j] # on modifie le point candidat
  }
  return(Z)
}