Các lượt nộp
    Danh sách bài
    Trang chủ
    Báo lỗi

    solution

    Đề bài: [Xác suất - Thống kê] Hồi quy tuyến tính có trọng số

    Trong thực tế, các quan sát (xi,yi)(x_i,y_i)(xi​,yi​), i=1,…,ni=1,\dots,ni=1,…,n có thể có độ tin cậy khác nhau, thể hiện qua trọng số wi≥0w_i \ge 0wi​≥0 (quan sát đáng tin hơn có trọng số lớn hơn). Hồi quy tuyến tính có trọng số (weighted least squares) tìm b0,b1b_0,b_1b0​,b1​ cực tiểu hoá ∑iwi(yi−b0−b1xi)2\sum_i w_i (y_i-b_0-b_1x_i)^2∑i​wi​(yi​−b0​−b1​xi​)2, cho bởi:

    xˉw=∑wixi∑wi,yˉw=∑wiyi∑wi\bar x_w=\dfrac{\sum w_ix_i}{\sum w_i},\quad \bar y_w=\dfrac{\sum w_iy_i}{\sum w_i}xˉw​=∑wi​∑wi​xi​​,yˉ​w​=∑wi​∑wi​yi​​ b1=∑wi(xi−xˉw)(yi−yˉw)∑wi(xi−xˉw)2,b0=yˉw−b1xˉwb_1=\dfrac{\sum w_i(x_i-\bar x_w)(y_i-\bar y_w)}{\sum w_i(x_i-\bar x_w)^2},\qquad b_0=\bar y_w-b_1\bar x_wb1​=∑wi​(xi​−xˉw​)2∑wi​(xi​−xˉw​)(yi​−yˉ​w​)​,b0​=yˉ​w​−b1​xˉw​

    Nếu n<2n<2n<2, tổng trọng số ∑wi≤0\sum w_i \le 0∑wi​≤0, hoặc mẫu số ∑wi(xi−xˉw)2=0\sum w_i(x_i-\bar x_w)^2 = 0∑wi​(xi​−xˉw​)2=0 (không xác định được b1b_1b1​), chỉ in KHONG_XAC_DINH.

    Ví dụ: khi tất cả trọng số bằng nhau, kết quả trùng với hồi quy bình phương nhỏ nhất thông thường: với w=(1,1,1,1,1)w=(1,1,1,1,1)w=(1,1,1,1,1), x=(1,2,3,4,5)x=(1,2,3,4,5)x=(1,2,3,4,5), y=(2,4,5,4,5)y=(2,4,5,4,5)y=(2,4,5,4,5) ta được b0=2.2000b_0=2.2000b0​=2.2000, b1=0.6000b_1=0.6000b1​=0.6000.

    • Định dạng đầu vào:

      Dòng 1: số nguyên nnn (n≥1n \ge 1n≥1). Dòng 2: nnn số thực w1 … wnw_1\ \dots\ w_nw1​ … wn​ (trọng số, không âm). Dòng 3: nnn số thực x1 … xnx_1\ \dots\ x_nx1​ … xn​. Dòng 4: nnn số thực y1 … yny_1\ \dots\ y_ny1​ … yn​.

    • Định dạng đầu ra:

      Nếu không xác định được: in đúng 1 dòng KHONG_XAC_DINH. Ngược lại in 1 dòng: b0b_0b0​ và b1b_1b1​, mỗi giá trị làm tròn 4 chữ số thập phân, cách nhau 1 dấu cách.

    Ví dụ:

    Đầu vào:

    5
    1 1 1 1 1
    1 2 3 4 5
    2 4 5 4 5

    Đầu ra:

    2.2000 0.6000
    

    Đầu vào:

    5
    1 1 1 1 10
    1 2 3 4 5
    2 4 5 4 5

    Đầu ra:

    2.3125 0.5437
    

    Đang tải editor...