NHMM模型日数据运行出现NaN值报错的解决方法咨询

运行代码 #

# devtools::install_github("cran/NHMM")
library(NHMM)
gfs.nhmm <- NHMM(y = obs_daily_train, X = atmvar_mme, K = 6, 
  iters = 1000,  burnin = 100, emdist = "gamma", nmix = 1, 
  delta = TRUE, outdir = paste(outdir_MME, "/", sep = ""), 
  ypred = ypred.gfs, Xp = xp.mme)

报错信息 #

运行至8%时触发如下错误：

in if (sum(pro) == 0) { : missing value where TRUE/FALSE needed
Called from: Cgetz(z, QQ, denzity, subseqy)

问题排查细节 #

• 调试发现Pro向量（长度为6）的取值为：NaN 2.554379e-127 NaN NaN NaN 0.000000e+00
• 此前使用3小时分辨率数据运行NHMM模型可成功执行，切换为日分辨率数据后出现该问题
• 已确认输入数据中无NA、NaN或无穷值，也无字符类型数据

需求 #

寻求解决该问题的方法，确保代码无报错完成运行。

内容的提问来源于stack exchange，提问作者Sanghita Basu