Я пишу GAM, используя пакет mgcv, который предсказывает численность норы и распределение двух разных видов на острове, используя данные, полученные во время экскурсии, и изображения, полученные со спутника Sentinel.Было обследовано 101 участок.922 нор, относящихся к виду 1, были зарегистрированы на 66 участках, а 29 нор, относящихся к виду 2, были зарегистрированы на 8 участках.
Я использовал отрицательное биномиальное распределение для вида 1, так как при использовании распределения Пуассона получилась модель болееразошлись.Максимальная модель была:
gam(Species_1 ~ s(x, y, bs="ts") +
Sentinel2_band_1 + Sentinel2_band_2 + Sentinel2_band_3 + Sentinel2_band_4 + Sentinel2_band_5 +
Sentinel2_band_6 + Sentinel2_band_7 + Sentinel2_band_8 + Sentinel2_band_9 + Sentinel2_band_10 +
I(Sentinel2_band_1^2) + I(Sentinel2_band_2^2) + I(Sentinel2_band_3^2) + I(Sentinel2_band_4^2) + I(Sentinel2_band_5^2) +
I(Sentinel2_band_6^2) + I(Sentinel2_band_7^2) + I(Sentinel2_band_8^2) + I(Sentinel2_band_9^2) + I(Sentinel2_band_10^2) +
aspect + elevation + slope +
I(aspect^2) + I(elevation^2) + I(slope^2) +
aspect:elevation + aspect:slope + elevation:slope,
data = dat,
family = nb(1))
Процесс выбора модели привел к модели, которая дает приемлемые результаты.
Когда я запускаю ту же модель, используя вид 2 в качестве переменной отклика, я получаюследующее сообщение об ошибке:
Warning message:
In newton(lsp = lsp, X = G$X, y = G$y, Eb = G$Eb, UrS = G$UrS, L = G$L, :
Fitting terminated with step failure - check results carefully
Диагностические графики также выглядят довольно хитроумно:
Мое предположение о проблеме, с которой я столкнулсяпроисходит из-за гораздо меньшего размера выборки для вида 2.
Есть идеи, что я могу сделать, чтобы решить эту проблему?