如何计算Iris数据集的真值和估算值之间的RMSE。
我在虹膜数据集中生成10%的NA值。 后来我使用Amelia包来估算和预测这些值。
现在我需要在原始虹膜数据集和模拟虹膜数据集之间计算RMSE
这是我的代码:
library(Amelia)
#load data
data("iris")
#produce 10% missing data
iris.mis <- prodNA(iris, noNA = 0.1)
#summary(iris.mis)
#run amelia
#m – Refers to 2 imputed data sets
#noms – keep nominal variables here
amelia_imp <- amelia(iris.mis, m=2, noms = "Species")
#check imputed outputs
amelia_imp$imputations[[1]]
amelia_imp$imputations[[2]]
谢谢。