search

ZELLMECHANIK-DRESDEN / ShapeOut

Shape-Out has been superseded by Shape-Out 2.
https://shapeout.readthedocs.io
GNU General Public License v2.0
5 stars 5 forks source link

IndexError in lin_mix_mod.py line 511 #248

Closed paulmueller closed 5 years ago

paulmueller commented 5 years ago

This error happens when performing an lmm analysis for volume in ShapeOut for some sample combinations (for some others it works).

Error_lmm_volume-1

The data are located at MPL guck\Members\Marta\Vimentin\lmm_analysis\Wt_VimKO_Filtered\Wt_VimKO.zmso

paulmueller commented 5 years ago

@maikherbig Maik, did see this type of error before? It looks like R does not return the correct number of coefficients. The line affected is: https://github.com/ZELLMECHANIK-DRESDEN/ShapeOut/blob/0.9.5/shapeout/lin_mix_mod.py#L511

Could this be a nan-value issue?

paulmueller commented 5 years ago

Subscribe @MartaUrb

maikherbig commented 5 years ago

Yes, this is a nan issue. For example look at \20170825_Marta_DC_VimKO_d11\02_VimKO_cnt\Online\M2_0.320000ul_s.tdms ->Cell track->time: the first entry is 119, the second 166 Then look into the contours file and you will see that it starts like this: Contour in frame 166 ( 157 , 32 ) ( 156 , 33 ) ... Apparently, the first conotur (time 119) is missing. Hence, no volume or principal inertia ratio can be calculated. Five of Martas files are affected. One can circument the issue in ShapeOut by using a filter for volume from -999999999 to 9999999. Then, the filter throws out the nan-event and lmm runs nicely :)

paulmueller commented 5 years ago

OK, so a fix in lin_mix_mod.py would be to remove all nan/inf values in the arrays given to linmixmod(xs=..)?

MartaUrb commented 5 years ago

Thank you guys for stepping in! And especially thank you Maik for the workaround! You are the best!

On Sat, 3 Aug 2019, 23:23 Paul Müller notifications@github.com wrote:

OK, so a fix in lin_mix_mod.py would be to remove all nan/inf values in the arrays given to linmixmod in xs?

— You are receiving this because you were mentioned. Reply to this email directly, view it on GitHub https://github.com/ZELLMECHANIK-DRESDEN/ShapeOut/issues/248?email_source=notifications&email_token=ADXH5S2NJEHP263KMJGGEHTQCXZNDA5CNFSM4II6HMQKYY3PNVWWK3TUL52HS4DFVREXG43VMVBW63LNMVXHJKTDN5WW2ZLOORPWSZGOD3PWASY#issuecomment-517955659, or mute the thread https://github.com/notifications/unsubscribe-auth/ADXH5S32ZJHDUVHF2YDW2X3QCXZNDANCNFSM4II6HMQA .

paulmueller commented 5 years ago

Thank you for your help. This will be fixed in the next release of Shape-Out.