我的代码如下:

lat <- 24.419
lon <- 54.502
days <- 1
Tmax <- 44.3
Tmin <- 26.0
tal <- 0.72
BCb_guess <- 0.13
epsilon <- 0.5
rad_mea <- 254
sirad::bccal(lat, days, rad_mea, extraT = NULL, Tmax, Tmin, tal)


当我从bccal()包执行函数sirad时,出现以下错误:

Error in prettyNum(.Internal(format(x, trim, digits, nsmall, width, 3L,  :
  invalid 'trim' argument


我对R完全陌生,将非常感谢您对解决问题的任何帮助,谢谢

最佳答案

我已经调试了。
bccal函数在这里定义。

https://github.com/cran/sirad/blob/master/R/bccal.R

问题出在第16行。

dtempM <- mean(as.numeric(aggregate(Zdtemp, by = format(time(Zdtemp),
    "%m"), FUN = mean, na.rm = TRUE)), na.rm = T)


特别是在这里

format(time(Zdtemp), "%m")


我真的不知道你想要实现什么(不要知道关于bccal的线索,我只能说它取决于Tmax,Tmin和天数),所以这就是我现在所能做的。
也许会为你敲钟

07-28 13:18