绘制套索L1约束

时间:2018-07-06 17:50:44

标签: r linear-regression absolute-value

在套索回归中,使用L1约束:

我正在尝试使用R绘制约束。一个示例如下所示: enter image description here

这是我编写的简单R代码:

beta= seq(-1, 1, length=100)
lambda=2
penalty=lambda*abs(beta)
plot(penalty, type="l") 

它仅绘制情节的下部。我怎么也可以在一个图中绘制上部? 这是情节:

enter image description here

1 个答案:

答案 0 :(得分:1)

要绘制abs(b1) + abs(b2) = 1,可以绘制b2 = 1 - abs(b1),然后绘制负数(因为abs表示其中任一个都是解决方案)

b1 <- seq(-1, 1, .01)
b2 <- 1 - abs(b1)
plot(b1, b2, ylim = c(-1, 1), type = 'l')
b2 <- -b2
lines(b1, b2)