- for (dk_index in 1:length(d_K))
- {
- d = d_K[[dk_index]][1]
- K = d_K[[dk_index]][2]
- #NOTE: sometimes large errors if pr is not balanced enough (e.g. random);
- # 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")
- .computeMuCheckDiag(io$X, io$Y, K, jd_method="uwedge", β_ref)
- #TODO: some issues with jedi method (singular system)
- #.computeMuCheckDiag(io$X, io$Y, K, jd_method="jedi", β_ref)
- }
+ for (dk_index in 1:length(d_K))
+ {
+ d <- d_K[[dk_index]][1]
+ K <- d_K[[dk_index]][2]
+ #NOTE: sometimes large errors if pr is not balanced enough (e.g. random);
+ # 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")
+ .computeMuCheckDiag(io$X, io$Y, K, jd_method="uwedge", β_ref)
+ #TODO: some issues with jedi method (singular system)
+ #.computeMuCheckDiag(io$X, io$Y, K, jd_method="jedi", β_ref)
+ }