I’m working with a dataset that has missing values, and I am using multiple imputation with the mice package in R. I then fit a Generalized Additive Model (GAM) using the mgcv package on the imputed datasets. I would like to plot the overall fit of the model across all imputations, but I’m not sure how to combine the plots from each imputed dataset.

How do you extract the pooled estimated degrees of freedom (edf) from a pooled generalized additive model ?