3

R 3 天大，无法弄清楚我做错了什么。我正在尝试将一些具有两种交互方式的列发送到 glmnet cox 模型中。我有一些名为 dtable 的 data.frame()

编辑以使代码可重现

xs<-c("Col1", "Col2", "Col3")
v<-c(1, 2, 3, 4, 5, 6, 7, 8, 9, 10, NA, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, NA, 25, 26, 27, 28, 29, 30)
df<-data.frame(matrix(v,ncol=3))
dm<-as.matrix(df)
dm<-matrix(dm[complete.cases(dm)], ncol=3)
colnames(dm)<-xs
dfdata<-data.frame(dm)
f<-as.formula(time~.*.)
xmatrix<-model.matrix(f, dfdata)[,-1]

当我运行它时，我得到了错误

Error in model.frame.default(object, data, xlev = xlev) :    
object is not a matrix

提前致谢

4

3 回答 3

4

这是因为公式：时间~.*。在 data.frame() 中，没有时间列供公式跨越所有内容。

于 2016-11-18T18:19:19.583 回答
0

我认为不是这样：

dm<-matrix(dm[complete.cases(dm)], ncol=3)

你要这个：

dm <- as.matrix(dm[complete.cases(dm), ])

不要自吹自擂（太多），而是考虑使用glmnetUtils来简化管理数据帧、公式和模型矩阵的任务。

于 2016-11-18T17:15:49.527 回答
-2

我怀疑该错误准确地说明了问题所在。您传递给的对象model.matrix()不是矩阵。结果是class(data)什么？可能是一个data.frame。

尝试添加as.matrix()到model.matrix()对data.

另外两个注意事项 - 不要调用您的 data.frames data。还使用可重现的示例在此处提出问题。你会得到更好的回应。

于 2016-11-18T17:04:30.610 回答