Hello,

We've been going through the code to understand the potential variability within the model, and spotted that the model$logResponsibility is assigned with model$responsibility in this line - is that a typo, or is there any reason to do specifically that?

if (!is.matrix(model$responsibility) ) {
  model$responsibility <- t ( as.matrix(model$responsibility)  )
  model$logResponsibility <- t ( as.matrix(model$responsibility)  )
}

Thank you in advance.

Best regards,
Nadezda