能够


11

如果 β=一个[RG一世ñβÿ-Xβ22+λβ1个, 能够 β2 增加时间 λ 增加?

我认为这是可能的。虽然β1个 何时不增加 λ增加(我的证明),β2可以增加。下图显示了一种可能性。什么时候λ 增加,如果 β 从(线性)行进 P, 然后 β2 增加而 β1个减少。但是我不知道如何构造一个具体的例子(即Xÿ),这样 β证明这种行为。有任何想法吗?谢谢。

在此处输入图片说明

Answers:


10

答案是肯定的,并且您在 2 在那里。

查找向量规范的等价定义。你会发现

X2X1个ñX2
哪里 ñ 是向量的维数 X。因此,有一些回旋的余地2 规范,相比 1个 规范。

实际上,您要解决的问题可以表述为:

d 这样

X+d2>X2
同时
X+d1个<X1个

平方第一个不等式,扩大并看到

2一世X一世d一世>-一世d一世2
并假设 X一世0X一世+d一世0,我们从第二个不等式中得出
一世d一世<0。
任何 d 满足这些约束将增加 2 规范,同时减少 1个 规范。

在您的示例中 d[-0.40.3]ŤX:=P[0.50.6]Ť

一世d一世-0.1<0
2一世P一世d一世-0.04>-0.25-一世d一世2

但是,它与 Xÿ
ziyuang,2015年

3

感谢@TommyL的回答,但他的回答并不直接涉及到 Xÿ。我自己以某种方式“解决”了这个问题。首先,什么时候λ 增加, β2 当每个 β一世单调减少。这发生在X 是正交的,我们有

β一世=s一世Gñβ一世大号小号β一世大号小号-λ+

从几何上讲,在这种情况下 β 垂直于 1个 规范,所以 β2 不能增加。

实际上,Hastie等人。在论文中提到正向逐步回归和单调套索,轮廓路径的单调性的充要条件:

在此处输入图片说明

在本文的第6节中,他们基于分段线性基函数构造了一个人工数据集,该数据集违反了上述条件,显示出非单调性。但是,如果运气好的话,我们还可以创建一个随机数据集,以相似的方式展示其行为,但方式更简单。这是我的R代码:

library(glmnet)
set.seed(0)
N <- 10
p <- 15
x1 <- rnorm(N)
X <- mat.or.vec(N, p)
X[, 1] <- x1
for (i in 2:p) {X[, i] <- x1 + rnorm(N, sd=0.2)}
beta <- rnorm(p, sd=10)
y <- X %*% beta + rnorm(N, sd=0.01)
model <- glmnet(X, y, family="gaussian", alpha=1, intercept=FALSE)

我故意让 X的列高度相关(与正常情况相去甚远), β具有大量正面和负面条目。这是β的配置文件(不足为奇,只有5个变量被激活):

在此处输入图片说明

和之间的关系 λβ2

在此处输入图片说明

所以我们可以看到一段时间 λβ2 随着增加 λ 增加。

By using our site, you acknowledge that you have read and understand our Cookie Policy and Privacy Policy.
Licensed under cc by-sa 3.0 with attribution required.