Optimizing planarflexors fails with LaiArnold2017 model update

[GEMcConnochie]

Hi there,
Thank you for sharing your code. I am trying to run the tool using this model: https://simtk.org/frs/?group_id=1302

I am not able to optimise the gastroc, soleus and plantar flexors with the following error:

Do you know why this occurs?
Thanks

[modenaxe]

Hello, the error is about indexing, so I would recommend first checking that you are looping starting from 0, if you are going through OpenSim sets.
If that's not the issue, could you please try this version of the MATLAB tool? https://github.com/modenaxe/MuscleParamOptimizer/tree/updOS4 and let me know if it works?
Cheers,
Luca

[GEMcConnochie]

Hi Luca, Thanks for getting back to me, I am not sure what you mean by looping starting from 0, I am just using the example script and replacing the ref and target models. I tried the other tool and had the same error. I think the issue is the matrix returned as Mus_ref from sampleMuscleQuantities(osimModel_ref,curr_mus,'all',N_eval) is size 100x5, whilst sampleMuscleQuantities(osimModel_targ,curr_mus_scaled,'MTL',N_eval) returns MTL_targ as a 100x1 double. Hence it cannot be indexed by the okList as size 1000x1. Cheers, Grace

…

________________________________ From: Luca Modenese ***@***.***> Sent: Friday, 2 September 2022 10:35 AM To: modenaxe/MuscleParamOptimizer ***@***.***> Cc: Grace Eva McConnochie (Student) ***@***.***>; Author ***@***.***> Subject: Re: [modenaxe/MuscleParamOptimizer] Optimizing planarflexors fails with LaiArnold2017 model update (Issue #21) CAUTION: External email. Only click on links or open attachments from trusted senders.

________________________________ Hello, the error is about indexing, so I would recommend first checking that you are looping starting from 0, if you are going through OpenSim sets. If that's not the issue, could you please try this version of the MATLAB tool? https://github.com/modenaxe/MuscleParamOptimizer/tree/updOS4<https://github.com/modenaxe/MuscleParamOptimizer/tree/updOS4> and let me know if it works? Cheers, Luca — Reply to this email directly, view it on GitHub<#21 (comment)>, or unsubscribe<https://github.com/notifications/unsubscribe-auth/AVJ4KJPOL6OBC26ST2SGBTLV4FHG7ANCNFSM6AAAAAAQCAX3XY>. You are receiving this because you authored the thread.Message ID: ***@***.***>

[modenaxe]

It could also be that the tool has troubles navigating through the knee, since some dof are coupled with flexion.
Could you verify if the issue is present also with models with a simpler knee, e.g. gait2392? At least we will know if the issue is due to the tool or to the application to the Lai model specifically.
Luca

[GEMcConnochie]

Yes I think it is an issue with the model - I tried with 2 different gait2392 models and didn't have any errors.

[modenaxe]

Hi @GraceEMc , any update about this issue?

[GEMcConnochie]

Hi @modenaxe, thanks for following up on this.
I have only be able to continue using the Lai 2017 model by avoiding adjusting the muscles: ["tibpost","gaslat","gasmed","perbrev","perlong","soleus","tibant"] which all fail. So it likely is an issue with the knee flexion. If you have any potential workarounds that would allow me to also optimize these muscles please let me know.