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 bội với hai biến độc lập

    Cho nnn quan sát (x1,i,x2,i,yi)(x_{1,i}, x_{2,i}, y_i)(x1,i​,x2,i​,yi​). Xét mô hình hồi quy tuyến tính bội

    y=β0+β1x1+β2x2.y = \beta_0 + \beta_1 x_1 + \beta_2 x_2.y=β0​+β1​x1​+β2​x2​.

    Các hệ số β0,β1,β2\beta_0, \beta_1, \beta_2β0​,β1​,β2​ được ước lượng bằng bình phương tối thiểu, tức là nghiệm của hệ phương trình chuẩn (normal equations) XTXβ=XTyX^T X \boldsymbol{\beta} = X^T \mathbf{y}XTXβ=XTy, cụ thể:

    {nβ0+β1∑x1,i+β2∑x2,i=∑yiβ0∑x1,i+β1∑x1,i2+β2∑x1,ix2,i=∑x1,iyiβ0∑x2,i+β1∑x1,ix2,i+β2∑x2,i2=∑x2,iyi\begin{cases} n\beta_0 + \beta_1\sum x_{1,i} + \beta_2 \sum x_{2,i} = \sum y_i \\ \beta_0 \sum x_{1,i} + \beta_1 \sum x_{1,i}^2 + \beta_2 \sum x_{1,i}x_{2,i} = \sum x_{1,i} y_i \\ \beta_0 \sum x_{2,i} + \beta_1 \sum x_{1,i}x_{2,i} + \beta_2 \sum x_{2,i}^2 = \sum x_{2,i} y_i \end{cases}⎩⎨⎧​nβ0​+β1​∑x1,i​+β2​∑x2,i​=∑yi​β0​∑x1,i​+β1​∑x1,i2​+β2​∑x1,i​x2,i​=∑x1,i​yi​β0​∑x2,i​+β1​∑x1,i​x2,i​+β2​∑x2,i2​=∑x2,i​yi​​

    Đề bài đảm bảo hệ có nghiệm duy nhất (ma trận XTXX^TXXTX khả nghịch).

    Ví dụ: Với dữ liệu khớp đúng mô hình y=1+2x1+3x2y = 1 + 2x_1 + 3x_2y=1+2x1​+3x2​ tại 4 điểm (1,1,6),(2,1,8),(1,2,9),(2,2,11)(1,1,6),(2,1,8),(1,2,9),(2,2,11)(1,1,6),(2,1,8),(1,2,9),(2,2,11), kết quả là β0=1.0000,β1=2.0000,β2=3.0000\beta_0=1.0000, \beta_1=2.0000, \beta_2=3.0000β0​=1.0000,β1​=2.0000,β2​=3.0000.

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

      Dòng 1: số nguyên nnn (3≤n≤10003 \le n \le 10003≤n≤1000). nnn dòng tiếp theo, mỗi dòng gồm 3 số thực x1,i x2,i yix_{1,i}\ x_{2,i}\ y_ix1,i​ x2,i​ yi​ cách nhau bởi khoảng trắng.

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

      In ra 3 số β0,β1,β2\beta_0, \beta_1, \beta_2β0​,β1​,β2​ trên cùng một dòng, cách nhau bởi khoảng trắng, mỗi số làm tròn đến 4 chữ số thập phân, theo đúng thứ tự trên.

    Ví dụ:

    Đầu vào:

    3
    1 0 1
    0 1 2
    1 1 4
    

    Đầu ra:

    -1.0000 2.0000 3.0000
    

    Đầu vào:

    4
    1 1 6
    2 1 8
    1 2 9
    2 2 11
    

    Đầu ra:

    1.0000 2.0000 3.0000
    

    Đang tải editor...