一个预测器案例的 GAMens 包答案

【问题标题】：GAMens package for one predictor case一个预测器案例的 GAMens 包
【发布时间】：2012-03-16 16:05:06
【问题描述】：

我正在尝试将 GAMens 包用于非常简单的人工数据集的“一维数据集”（一个预测变量和一个响应）：

library(GAMens)
trainingdata <- read.table("D:/data.txt", header = TRUE, sep="\t")
trainingdata.GAMrsm <- GAMens(class~., trainingdata, 1,  autoform=TRUE, iter=10, rsm=FALSE)

不幸的是我得到了：

[.data.frame(data, , as.character(formula[[2]])) 中的错误：选择了未定义的列

也许我遗漏了一些东西，但是应该有可能让它适用于一维情况（“维度”更大的一个工作正常）。谢谢。

编辑：

这是生成人工数据的方式：

n.site <- 150

    X1<- sort(runif(n = n.site, min = -1, max =1))

    xb <- 0.0 + 3.0*X1 

    occ.prob <- 1/(1+exp(-xb))

    true.presence <- rbinom(n = n.site, size = 1, prob = occ.prob)

    data <- data.frame(X1, true.presence)

    write.matrix(data, file = "data.txt", sep = "\t")

您也可以在此处下载数据：

Data

【问题讨论】：

你能用我们都能得到的数据集（或合成示例）复制这个失败吗？可能来自 mlbench 的 Ionosphere 数据，如 help(GAMens) 中所示的示例中使用的那样？还是其中的一部分？
请参阅编辑。谢谢。

标签： r

【解决方案1】：

您的data 输入没有名为“类”的列，因此您使用的公式class ~ . 无效。当然，错误消息不是很有帮助...

这似乎有效（但速度很慢，10 分钟后仍在运行）：

trainingdata.GAMrsm

【讨论】：

是的，对不起。我在本地数据副本 (data.txt) 中将 true.presence 重命名为 class。我知道 autoform=TRUE 但 true.presence ~ X1 不限制“模型空间”吗？我想用true.presence~。如果一个有多个预测器，这将起作用。谢谢。

【解决方案2】：

该问题已在新版本的软件包中得到解决。请通过 CRAN 存储库下载 GAMens 包的 1.1.2 版本。

最好的问候

【讨论】：