Plot op error

[BoburaII]

Hi,
I’m trying to plot op from a PM_result object (res), using: res$op$plot(outeq = 2)
I get the error: Error in summary(mod)$coefficients[2, 2] : subscript out of bounds

This didn’t happen previously with the same model and data. The model has 3 output equations, plotting op with outeq = 1 or outeq = 3 is fine. I'm using Pmetrics 2.1.1, R4.3.1, and recently updated macOS sequoia.
I appreciate any help with this.

Bob

[mhovd]

Thanks for raising this issue - this definitely sounds like a bug in how the coefficients for the regression line are calculated.

While waiting for us to patch the bug, you can also access the underlying data directly with res$op$data, which can be then used in your favorite plotting library, e.g. ggplot2 or plotly. Just remember to filter for your desired outeq, icen, and pred.type.

Let me know if you need any help with the code to make the figures!

[BoburaII]

Thank you for the response. On inspection of the PM_op data, I've noticed that pop and post predictions aren't being generated for outeq = 2. I've double checked the model file and all seems ok. Previously, the same model and datafile did run successfully. I'm unsure what's causing it to break now..

[mhovd]

Thanks for that information. Can you verify that you 1) have an output equation for outeq = 2 in your model file and 2) have at least one individual with outeq = 2 in your data?