Hello, I am looking for some guidance on how to calculate p-values for mixed-effects multinomial logistic regression using the npmlt() function from the mixcat package in R. I have fitted a model using this function and obtained the estimates and standard errors for each parameter. However, I am not sure how to derive the p-values from these values. I have tried to compute the Z values by dividing the estimates by the standard errors, and then compare them with the critical values from a standard normal distribution. Is this a valid method? I have yet found any documentation or examples on how to calculate p-values using npmlt() in R. I would appreciate any discussion or suggestion on this topic. Thank you very much.