--- /dev/null
+%call generateIO with default parameters (random means, covariances = identity, equirepartition)
+function[X,Y,Z] = generateIOdefault(n, p, m, k)
+
+ rangeX = 100;
+ meanX = rangeX * (1 - 2*rand(k, p));
+ covX = zeros(p,p,k);
+ covY = zeros(m,m,k);
+ for r=1:k
+ covX(:,:,r) = eye(p);
+ covY(:,:,r) = eye(m);
+ end
+ pi = (1/k) * ones(1,k);
+
+ %initialize beta to a random number of non-zero random value
+ beta = zeros(p,m,k);
+ for j=1:p
+ nonZeroCount = ceil(m*rand(1));
+ beta(j,1:nonZeroCount,:) = rand(nonZeroCount, k);
+ end
+
+ [X,Y,Z] = generateIO(meanX, covX, covY, pi, beta, n);
+
+end