сообщение об ошибке в команде "sampleMcm c" в пакете HMS C в R - PullRequest
0 голосов
/ 06 февраля 2020

Я получаю сообщение об ошибке, которое не могу понять при запуске приведенного ниже кода в R

У меня есть следующий пример набора данных графиков x численности видов

dput(fern[1:5,c(1,3:5)])

structure(c("plot1", "plot2", "plot3", "plot4", "plot5", " 1", 
" 7", " 8", " 10", " 0", "5", "3", "2", "1", "4", " 18", " 13", 
" 0", " 9", " 10"), .Dim = 5:4, .Dimnames = list(NULL, c("plot", 
"Actin_penn", "Actin_subt", "Adian_adia")))

и окружающей среды data

dput(env2[1:5,-c(2)])

structure(c("plot1", "plot2", "plot3", "plot4", "plot5", "2902", 
"3047", "2797", "2752", "2832", "1303", "1355", "1252", "1230", 
"1272", "468", "488", "445", "444", "449", "2958", "2962", "3874", 
"3874", "3782", "127", "129", "145", "147", "148"), .Dim = 5:6, .Dimnames = list(
    NULL, c("plot", "Band4median15", "Band5median15", "Band7median15", 
    "Chelsa450bio12", "SRTM90")))

Сначала я запускаю модель

XData=as.data.frame(env2)
m = Hmsc(Y = fern,
         XData = XData,
         distr = "lognormal poisson")

и, кажется, в порядке, по крайней мере, никаких сообщений об ошибках.

Затем я пробую образец

thin = 1
samples = 500
nChains = 1
set.seed(1)
m = sampleMcmc(m, samples = samples, thin = thin,
               adaptNf = rep(ceiling(0.4*samples*thin),1),
               transient = ceiling(0.5*samples*thin),
               nChains = nChains, nParallel = nChains)

И получите

Error in checkForRemoteErrors(val) : 
  4 nodes produced errors; first error: non-numeric argument to binary operator

Есть идеи, что не так?

Добро пожаловать на сайт PullRequest, где вы можете задавать вопросы и получать ответы от других членов сообщества.
...