X-Git-Url: https://git.auder.net/?a=blobdiff_plain;f=R%2FgenerateIOdefault.R;h=b0d748ae4e5955c54ccb4a03e6819ae92b299463;hb=c3bc47052f3ccb659659c59a82e9a99ea842398d;hp=3613f2b9b58998386bdb175f3d2d44b0e31729b2;hpb=0b216f854a21821f9be375d07c2932b31e227e78;p=valse.git diff --git a/R/generateIOdefault.R b/R/generateIOdefault.R index 3613f2b..b0d748a 100644 --- a/R/generateIOdefault.R +++ b/R/generateIOdefault.R @@ -8,22 +8,22 @@ #----------------------------------------------------------------------- generateIOdefault = function(n, p, m, k) { - covX = diag(p) - covY = array(0, dim=c(m,m,k)) - for(r in 1:k) - { - covY[,,r] = diag(m) - } - - pi = rep(1./k,k) - - beta = array(0, dim=c(p,m,k)) - for(j in 1:p) - { - nonZeroCount = ceiling(m * runif(1)) - beta[j,1:nonZeroCount,] = matrix(runif(nonZeroCount*k), ncol=k) - } - - sample_IO = generateIO(covX, covY, pi, beta, n) - return (list(X=sample_IO$X,Y=sample_IO$Y)) + covX = diag(p) + covY = array(0, dim=c(m,m,k)) + for(r in 1:k) + { + covY[,,r] = diag(m) + } + + pi = rep(1./k,k) + + beta = array(0, dim=c(p,m,k)) + for(j in 1:p) + { + nonZeroCount = ceiling(m * runif(1)) + beta[j,1:nonZeroCount,] = matrix(runif(nonZeroCount*k), ncol=k) + } + + sample_IO = generateIO(covX, covY, pi, beta, n) + return (list(X=sample_IO$X,Y=sample_IO$Y)) }