- html - 出于某种原因,IE8 对我的 Sass 文件中继承的 html5 CSS 不友好?
- JMeter 在响应断言中使用 span 标签的问题
- html - 在 :hover and :active? 上具有不同效果的 CSS 动画
- html - 相对于居中的 html 内容固定的 CSS 重复背景?
我想在弹性网回归模型中强制一些自变量。使用时出现错误惩罚因子在 cva.glmnet() 函数中。任何建议和意见表示赞赏。
# --------------- import data
# outcome: am
# ten predictors: mpg+cyl+disp+hp+drat+wt+qsec+vs+gear+carb
mydata<-mtcars
mydata$am<-factor(mydata$am,levels=c(0,1),labels=c("no","yes"))
mydata$gear<-factor(mydata$gear,levels=c(3,4,5),labels=c("level3","level4","level5"))
str(mydata)
'data.frame': 32 obs. of 11 variables:
$ mpg : num 21 21 22.8 21.4 18.7 18.1 14.3 24.4 22.8 19.2 ...
$ cyl : num 6 6 4 6 8 6 8 4 4 6 ...
$ disp: num 160 160 108 258 360 ...
$ hp : num 110 110 93 110 175 105 245 62 95 123 ...
$ drat: num 3.9 3.9 3.85 3.08 3.15 2.76 3.21 3.69 3.92 3.92 ...
$ wt : num 2.62 2.88 2.32 3.21 3.44 ...
$ qsec: num 16.5 17 18.6 19.4 17 ...
$ vs : num 0 0 1 1 0 1 0 1 1 1 ...
$ am : Factor w/ 2 levels "no","yes": 2 2 2 1 1 1 1 1 1 1 ...
$ gear: Factor w/ 3 levels "level3","level4",..: 2 2 2 1 1 1 1 2 2 2 ...
$ carb: num 4 4 1 1 2 1 4 2 2 4 ...
# --------------- do elastic net cross-validation for alpha and lambda simultaneously
# works
library("glmnetUtils")
set.seed(12345)
cvfit<-cva.glmnet(am~mpg+cyl+disp+hp+drat+wt+qsec+vs+gear+carb,
family="binomial",
alpha=seq(from=0,to=1,by=0.05),
nfolds=10,
data=mydata)
# --------------- force three variables of "cyl", "disp", "hp" in the final model
# does not work
set.seed(12345)
cvfit<-cva.glmnet(am~mpg+cyl+disp+hp+drat+wt+qsec+vs+gear+carb,
family="binomial",
penalty.factor=c(0,0,0,1,1,1,1,1,1,1),
alpha=seq(from=0,to=1,by=0.05),
nfolds=10,
data=mydata)
Error in approx(lambda, seq(lambda), sfrac) :
need at least two non-NA values to interpolate
最佳答案
此外,由于齿轮是一个因素,因此您需要在惩罚因子中再增加 2 个。和 cva.glmnet()
保留一个因子的所有虚拟变量,而不是省略一个,这样正则化将缩小对整体均值的影响。如果省略 1 个级别,则将剩余系数缩小为 0 会强制将预测推向基线级别(请参阅下面作者的评论)。
在任何情况下,您都可以检查您拥有的第一个示例:
set.seed(12345)
cvfit<-cva.glmnet(am~mpg+cyl+disp+hp+drat+wt+qsec+vs+gear+carb,
family="binomial",
alpha=seq(from=0,to=1,by=0.05),
nfolds=10,
data=mydata)
dim(cvfit$modlist[[1]]$glmnet.fit$beta)
[1] 12 100
所以你需要一个 12 的向量作为惩罚因子。如果我运行下面的,它的工作原理:
cvfit<-cva.glmnet(am~mpg+cyl+disp+hp+drat+wt+qsec+vs+gear+carb,
family="binomial",
penalty.factor=c(0,0,0,1,1,1,1,1,1,1,1,1),
alpha=seq(from=0,to=1,by=0.05),
nfolds=10,
data=mydata)
您可以看到前 3 个变量永远不会为零:
[,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [,11] [,12] [,13]
mpg 0 0 0 0 0 0 0 0 0 0 0 0 0
cyl 0 0 0 0 0 0 0 0 0 0 0 0 0
disp 0 0 0 0 0 0 0 0 0 0 0 0 0
hp 0 6 7 7 7 8 9 9 10 11 12 14 15
drat 0 7 7 7 8 8 8 9 10 10 11 13 14
wt 0 13 14 15 16 17 18 19 20 21 22 23 25
qsec 0 1 1 1 1 1 1 1 1 1 1 1 1
vs 0 11 12 13 14 15 15 16 17 17 18 19 19
gearlevel3 0 4 4 5 5 5 5 5 5 6 6 6 7
gearlevel4 0 50 49 50 51 52 53 55 58 61 65 70 77
gearlevel5 0 1 1 1 1 1 1 1 1 1 1 1 1
carb 0 12 17 23 29 37 46 62 65 70 75 80 88
[,14] [,15] [,16] [,17] [,18] [,19] [,20] [,21]
mpg 0 0 0 0 0 0 0 0
cyl 0 0 0 0 0 0 0 0
disp 0 0 0 0 0 0 0 0
hp 17 20 22 26 32 40 62 73
drat 16 19 23 27 33 41 58 73
wt 26 28 30 33 38 48 78 73
qsec 1 1 1 1 1 1 1 2
vs 20 22 23 25 29 35 58 73
gearlevel3 7 8 8 9 11 14 37 73
gearlevel4 84 87 85 84 82 80 78 73
gearlevel5 1 1 1 1 1 1 1 1
carb 88 87 85 84 82 80 78 73
关于r - 如何强制 cva.glmnet() 在弹性网络回归模型中包含特定变量?,我们在Stack Overflow上找到一个类似的问题: https://stackoverflow.com/questions/64509382/
这个问题在这里已经有了答案: 关闭 10 年前。 Possible Duplicate: How to nest OR statements in JavaScript? 有没有办法做到这一点:
在 JavaScript 中有没有办法让一个变量总是等于一个变量?喜欢var1 = var2但是当var2更新,也是var1 . 例子 var var1 = document.getElementBy
我正在努力理解这代表什么 var1 = var2 == var3 我的猜测是这等同于: if (var2 == var3): var1 = var2 最佳答案 赋值 var1 = var2
这个问题已经有答案了: What does the PHP error message "Notice: Use of undefined constant" mean? (2 个回答) 已关闭 8
我在临时表中有几条记录,我想从每条记录中获取一个值并将其添加到一个变量中,例如 color | caption -------------------------------- re
如何将字符串转为变量(字符串变量--> $variable)? 或者用逗号分隔的变量列表然后转换为实际变量。 我有 2 个文件: 列名文件 行文件 我需要根据字符串匹配行文件中的整行,并根据列名文件命
我有一个我无法解决的基本 php 问题,我也想了解为什么! $upperValueCB = 10; $passNodeMatrixSource = 'CB'; $topValue= '$uppe
这可能吗? php $variable = $variable1 || $variable2? 如果 $variable1 为空则使用 $variable2 是否存在类似的东西? 最佳答案 PHP 5
在 Perl 5.20 中,for 循环似乎能够修改模块作用域的变量,但不能修改父作用域中的词法变量。 #!/usr/bin/env perl use strict; use warnings; ou
为什么这不起作用: var variable; variable = variable.concat(variable2); $('#lunk').append(variable) 我无法弄清楚这一点
根据我的理解,在32位机器上,指针的sizeof是32位(4字节),而在64位机器上,它是8字节。无论它们指向什么数据类型,它们都有固定的大小。我的计算机在 64 位上运行,但是当我打印包含 * 的大
例如: int a = 10; a += 1.5; 这运行得很完美,但是 a = a+1.5; 此作业表示类型不匹配:无法从 double 转换为 int。所以我的问题是:+= 运算符 和= 运算符
您好,我写了这个 MySQL 存储过程,但我一直收到这个语法错误 #1064 - You have an error in your SQL syntax; check the manual that
我试图在我的场景中显示特定的奖牌,这取决于你的高分是基于关卡的目标。 // Get Medal Colour if levelHighscore goalScore { sc
我必须维护相当古老的 Visual C++ 源代码的大型代码库。我发现代码如下: bIsOk = !!m_ptr->isOpen(some Parameters) bIsOk的数据类型是bool,is
我有一个从 MySQL 数据库中提取的动态产品列表。在 list 上有一个立即联系 按钮,我正在使用一个 jquery Modal 脚本,它会弹出一个表单。 我的问题是尝试将产品信息变量传递给该弹出窗
这个问题在这里已经有了答案: 关闭 10 年前。 Possible Duplicate: What is the difference between (type)value and type(va
jQuery Core Style Guidelines建议两种不同的方法来检查变量是否已定义。 全局变量:typeof variable === "undefined" 局部变量:variable
这个问题已经有答案了: 已关闭11 年前。 Possible Duplicate: “Variable” Variables in Javascript? 我想肯定有一种方法可以在 JavaScrip
在语句中使用多重赋值有什么优点或缺点吗?在简单的例子中 var1 = var2 = true; 赋值是从右到左的(我相信 C# 中的所有赋值都是如此,而且可能是 Java,尽管我没有检查后者)。但是,
我是一名优秀的程序员,十分优秀!