some IMM examples fit the wrong model

[venpopov]

bmm_examples/Example5_Oberauer2017_IMMabc.R - it fits the '3p' model, not the abc bmm_examples/Example5_Oberauer2017_IMMbsc.R - it fits the '3p' model, not the bsc bmm_examples/Example5_Oberauer2017_IMMfull.R - it fits the '3p' model, not the full

brms_examples/Example5_Oberauer2017_IMMabc.R - it uses bmm code, not brms, and also fits the '3p' model

[venpopov]

Actually, I see that you did it in separate sections, and all files first fit the 3p model. But there are also errors in the filenames for the output files.

bmm_examples/Example5_Oberauer2017_IMMabc.R

line 1 - the description only mentions the Bays model line 72 - filename is for the abc model, but the fit code is for the 3p model

bmm_examples/Example5_Oberauer2017_IMMbsc.R

line 232 - it checks if (!file.exists(here("output","fit_E5_OL2017_IMMabc.RData"))), but the fit, the saved output and the loaded output is the bsc model

bmm_examples/Example5_Oberauer2017_IMMfull.R

line 234-263 - three different model files are specified for the check if exists, for the output, or for the load.

My suggestions:

• fix the filepaths in all scripts. To be safe, specify the filename as a separate variable and then just use that variable in all spots
• delete or change brms_examples/Example5_Oberauer2017_IMMabc.R to reflect the brms implementation
• I think it's confusing to have both the '3p' model and the other models in each file. Especially since they are not directly compared anywhere in the code. Extract the '3p' fit to a separate file, and remove that part from the imm scripts. Also update the description at the top of each file to say that the imm is being fit and which version it is.

After that, everything is ready I think.

[GidonFrischkorn]

I have updated the Scripts as you proposed and set up a new script for fitting the 3-parameter Mixture Model for Example 5.