#initialize beta to a random number of non-zero random value
β = array(0, dim=c(p,m,k))
for (j in 1:p)
{
nonZeroCount = sample(1:m, 1)
#initialize beta to a random number of non-zero random value
β = array(0, dim=c(p,m,k))
for (j in 1:p)
{
nonZeroCount = sample(1:m, 1)