这是R中的代码：

marsModel <- earth(eval(parse(text=paste(ResponseVariable,"~."))), data = data) #build model
ev <- evimp (marsModel)

响应变量不言自明，数据包含 178 个特征和 1 个响应变量。每个特征有 56 个观测值。

'ev' 最终从 178 个特征中选择了 10 个特征，以及它们的 gcv（广义交叉验证）、rss（残差平方和）和 nsubsets（变量出现的子集模型的数量）得分。