Tình huống
Tôi có một tập dữ liệu với một biến phụ thuộc và một biến độc lập . Tôi muốn điều chỉnh hồi quy tuyến tính từng phần liên tục với điểm dừng cố định / đã biết xảy ra tại (a_ {1}, a_ {2}, \ ldots, a_ {k}) . Các breakpoins được biết mà không có sự không chắc chắn, vì vậy tôi không muốn ước tính chúng. Sau đó, tôi phù hợp với hồi quy (OLS) có dạng y_ {i} = \ beta_ {0} + \ beta_ {1} x_ {i} + \ beta_ {2} \ operatorname {max} (x_ {i} -a_ { 1}, 0) + \ beta_ {3} \ operatorname {max} (x_ {i} -a_ {2}, 0) + \ ldots + \ beta_ {k + 1} \ operatorname {max} (x_ {i} - a_ {k}, 0) + \ epsilon_ {i} Đây là một ví dụ trongx k ( một 1 , một 2 , ... , một k ) y i = β 0 + β 1 x i + β 2 max ( x i - một 1 , 0 ) + β 3 max ( x i - một 2 , 0 ) + ... + β k + 1 max ( x
R
set.seed(123)
x <- c(1:10, 13:22)
y <- numeric(20)
y[1:10] <- 20:11 + rnorm(10, 0, 1.5)
y[11:20] <- seq(11, 15, len=10) + rnorm(10, 0, 2)
Giả sử rằng điểm dừng xảy ra ở :
mod <- lm(y~x+I(pmax(x-9.6, 0)))
summary(mod)
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) 21.7057 1.1726 18.511 1.06e-12 ***
x -1.1003 0.1788 -6.155 1.06e-05 ***
I(pmax(x - 9.6, 0)) 1.3760 0.2688 5.120 8.54e-05 ***
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
Phần chặn và độ dốc của hai đoạn lần lượt là: và cho phần thứ nhất và và cho phần thứ hai.
Câu hỏi
- Làm thế nào để dễ dàng tính toán đánh chặn và độ dốc của từng đoạn? Mô hình có thể được sửa lại để làm điều này trong một phép tính không?
- Làm thế nào để tính sai số chuẩn của từng độ dốc của từng đoạn?
- Làm thế nào để kiểm tra xem hai sườn liền kề có cùng độ dốc hay không (nghĩa là điểm dừng có thể được bỏ qua không)?
x
vàI(pmax(x-9.6,0))
, điều đó có đúng không?