- io = generateSampleIO(n, p, βs_ref[,,i], rep(0,K), model)
- μ = computeMu(io$X, io$Y, list(K=K))
- μ_aligned = alignMatrices(list(μ), ref=μ_ref, ls_mode="exact")[[1]]
+ io <- generateSampleIO(n, p, βs_ref[,,i], rep(0,K), link)
+ μ <- computeMu(io$X, io$Y, list(K=K))
+ μ_aligned <- alignMatrices(list(μ), ref=μ_ref, ls_mode="exact")[[1]]