Predicting with missing values in non-Gaussian case about kfas

[helske]

Can you give a reproducible example of this? The line 329 is related to variance computation based on the importance sampling, there is likely some issue in computing the weights in the previous stages, meaning that importance sampling likely does not work in your case for some reason (poor model specification, approximation does not work, numerical issues, ...), leading for example infinite weights.

For these cases there probably should be some checks in place after the importance sampling and more informative error message.

[jhal324]

Yes, sorry for the delay - I have been away. I've attached a snippet of the problem (I don't have permission to share fully).

predictProblem.zip

[helske]

Thanks. There are couple of issues in your example:

First, the argument for the new data is called newdata, whereas in your script you use argument testData which is ignored as there is not such argument for the predict method.

Second, as you have NA values in the Z array, it is impossible to predict the corresponding values of y (as y_t = link_function(Z_t*alpha_t)). So that causes the issue where there are missing values in the importance samples, which then cause issues in the later computations.

While you cannot get predictions for those y values, you can still get predictions for others, if you just set NA values in Z for example to zero. That shouldn't affect the predictions of other observations as corresponding values of y are still NA.

I now added the NAOK argument anyway so things should work now as expected.