error in PLN network model postTreatment

PLN-team / PLNmodels

A collection of Poisson lognormal models for multivariate count data analysis

https://pln-team.github.io/PLNmodels

GNU General Public License v3.0

54 stars 18 forks source link

error in PLN network model postTreatment #138

Open maddyduran opened 1 month ago

maddyduran commented 1 month ago

Hello,

I get this error when trying to fit a PLN network model with the latest version. This model with this data worked on an older version.

Error in `colnames<-`(`*tmp*`, value = colnames(responses)) : 
  attempt to set 'colnames' on an object with less than two dimensions

It appears that the issue is because there is an issue with the last of the models in the self$models list while going through postTreatment:

mahendra-mariadassou commented 1 month ago

Hello,

Thank you for reporting this bug. This seems to arise from automatic conversion of one-column matrices to vectors but I'm not quite sure where it happens nor why it happened only for a given value of $\lambda$.

Is there any chance you could share the dataset (or a subset that still generates the bug and/or by email if it's private) and exact call used so we can reproduce and fix it ?

Thanks