我想编写一个带有 NumericMatrix 参数的 Rcpp 函数。如果任何矩阵元素为 NA,则返回 true,否则返回 false。我尝试在所有列上循环 is_na 但我正在寻找一种更简洁的方法。我也关心速度。
bool check(NumericMatrix M){
n=M.ncol();
for(int i=0; i < n; i ++){
if(is_na( M(_,i) ){ return T;}
}
return F;
}
Run Code Online (Sandbox Code Playgroud) 我使用 tidymodels 框架拟合多项式逻辑回归模型来预测鸢尾花数据集中的物种。
library(tidymodels)
iris.lr = multinom_reg(
mode="classification",
penalty=NULL,
mixture=NULL
) %>%
set_engine("glmnet")
iris.fit = iris.lr %>%
fit(Species ~. , data = iris)
Run Code Online (Sandbox Code Playgroud)
然后我想查看模型的系数并写出公式。我的理解是我应该从 iris.fit 获取这个。
iris.fit 的输出有一个 100 行表,其中包含 Df、%Dev 、Lambda。iris 数据集只有 4 个预测变量。如何将此输出转换为系数?