# same note for β. However we could be more random than that...
β_ref <- rbind(diag(K),matrix(0,nrow=d-K,ncol=K))
io <- generateSampleIO(n, p=rep(1/K,K-1), β=β_ref, rep(0,K), link="logit")
# same note for β. However we could be more random than that...
β_ref <- rbind(diag(K),matrix(0,nrow=d-K,ncol=K))
io <- generateSampleIO(n, p=rep(1/K,K-1), β=β_ref, rep(0,K), link="logit")