Can


11

Nếu , có thể β * 2 tăng khiλtăng?β=argminβyXβ22+λβ1β2λ

Tôi nghĩ rằng điều này là có thể. Mặc dù không tăng khi λ tăng (tôi bằng chứng ), β * 2 có thể tăng. Hình dưới đây cho thấy một khả năng. Khi λ tăng, nếu β * hành trình (tuyến tính) từ P đến Q , sau đó β * 2 tăng trong khi β * 1 giảm. Nhưng tôi không biết cách xây dựng một ví dụ cụ thể (nghĩa là xây dựng Xyβ1λβ2λβPQβ2β1Xy), để hồ sơ của thể hiện hành vi này. Có ý kiến ​​gì không? Cảm ơn bạn.β

nhập mô tả hình ảnh ở đây

Câu trả lời:


10

Câu trả lời là có, và bạn có một bằng chứng đồ họa trong ngay tại đó.2

x2x1nx2,
nx21

Trong thực tế, vấn đề bạn muốn giải quyết có thể được nêu là:

d

x+d2>x2
x+d1<x1.

2ixidi>idi2
xi0xi+di0
idi<0.
d21

d[0.4,0.3]Tx:=P[0.5,0.6]T

idi0.1<0,
2iPidi0.04>0.25idi2.

Xy

3

Xyλβ2βiX

βi=sign(βiLS)(βiLSλ)+

β1β2

Trên thực tế, Hastie et al. được đề cập trong bài báo Chuyển tiếp hồi quy theo giai đoạn và Lasso đơn điệu , một điều kiện cần và đủ về tính đơn điệu của các đường dẫn hồ sơ:

nhập mô tả hình ảnh ở đây

Trong Phần 6 của bài báo, họ đã xây dựng một bộ dữ liệu nhân tạo dựa trên các hàm cơ sở tuyến tính từng phần vi phạm điều kiện trên, cho thấy tính không đơn điệu. Nhưng nếu chúng ta gặp may mắn, chúng ta cũng có thể tạo một tập dữ liệu ngẫu nhiên thể hiện hành vi tương tự nhưng theo cách đơn giản hơn. Đây là mã R của tôi:

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ββ

nhập mô tả hình ảnh ở đây

λβ2

nhập mô tả hình ảnh ở đây

λβ2λ

Khi sử dụng trang web của chúng tôi, bạn xác nhận rằng bạn đã đọc và hiểu Chính sách cookieChính sách bảo mật của chúng tôi.
Licensed under cc by-sa 3.0 with attribution required.