我的代码如下:
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和天数),所以这就是我现在所能做的。
也许会为你敲钟