目前,我正试图找到一个拟合模型和r中的实际观测的RMSE系统(RMSEs)和RMSE非系统(RMSEu)。
为了得到RMSE,我使用了以下语句:
sqrt( mean( (df$fitted - df$actual)^2 , na.rm = TRUE ) )
然而,我不确定如何使用r来计算RMSEu和rmse。是否有任何函数/包我可以用来计算这些?
非常感谢!
编辑:
这是RMSEu和rmse的定义
Adrian,在caret
包中有一个RMSE
函数
# Load library
library(caret)
# Check Root Mean Squared Error
RMSE(pred = df$fitted, obs = df$actual)
R包中有几个RMSE函数,但我不知道有一个具有系统和非系统RMSE。可以写我们自己的函数之类的sqrt(1/N * sum((modeled - observed)^2))
同样地。在这里,向量化会派上用场,即(differences)^2
项将产生一个差异平方的向量,可以求和。
找到了相关的库!:-(感谢那些帮助过我的人。对于那些仍在寻找rmse和RMSEu图书馆的人,请查看:https://rdrr.io/github/SticsRPacks/CroPlotR/man/predictor_assessment.html