Closed martin-gmx closed 4 years ago
Just thought of a situation where your solution won't work for shrinkage. If we have an omega that is nonzero but fixed, we will certainly still have shrinkage. Your solution to the fixing problem ignores that parameter. I think we can only ignore if the OMEGA is zero. Otherwise we need to account for it.
You're right. Good that you spotted this. :-)
Without notfixed in shrinkage.R problems inverting matrix. Without setting temp_1 in mf3.R crash with notfixed_bpop == 0.