我试图使用R中的“ mda”程序包运行混合物鉴别分析。但是,当我遵循手册并运行以下脚本时。我意识到结果不是唯一的:
library(mda)
data(iris)
irisfit <- mda(Species ~ ., data = iris)
irisfit
结果有时看起来像这样:
致电: mda(公式=物种〜。,数据=虹膜)
维度:4
组间差异百分比说明: v1 v2 v3 v4 95.06 97.78 99.59 100.00
自由度(每维度):5
训练错误分类错误:0.01333(N = 150)
距离:13.302
有时是这样的:
致电: mda(公式=物种〜。,数据=虹膜)
维度:5
组间差异百分比说明: v1 v2 v3 v4 v5 96.14 98.47 99.89 100.00 100.00
自由度(每维度):5
训练分类错误:0.02(N = 150)
距离:15.249
我在不同的工作站和不同的R版本上尝试过,但是情况保持不变。因此,这不应该是系统的问题。
据我了解,MDA的执行方式应产生独特的结果。有人可能会想到这是怎么发生的吗?