Assistance with Predictions Using PLNmodels 1.2.1

[ivangalvan]

Hello,

I am encountering differences in predictions when using the predict() function compared to the myPLN$fitted object. Here is a small reproducible example:

library(PLNmodels)
data(trichoptera)
trichoptera <- prepare_data(trichoptera$Abundance, trichoptera$Covariate)

myPLN <- PLN(Abundance ~ 1, data = trichoptera)

new_predictions = predict(myPLN,newdata = trichoptera)

plot(myPLN$fitted,new_predictions)

I would like to make predictions on both training and test sets and ensure that I am correctly using the predict() function for the test set. Could you please assist me with this issue? I am currently using PLNmodels_1.2.1.

Thank you for your help.

Best regards,
Ivan