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ê] Kiểm định ý nghĩa hệ số hồi quy bội

    Sau khi ước lượng mô hình hồi quy bội y^=b0+b1x1+b2x2\hat y = b_0 + b_1x_1 + b_2x_2y^​=b0​+b1​x1​+b2​x2​ bằng phương pháp bình phương nhỏ nhất (xem hệ phương trình chuẩn tắc ở bài "Hồi quy tuyến tính bội với hai biến độc lập"), ta cần kiểm định ý nghĩa thống kê của từng hệ số bằng thống kê t:

    tj=bjse(bj),se(bj)=σ^2⋅[(XTX)−1]jj,σ^2=SSEn−pt_j = \frac{b_j}{se(b_j)}, \qquad se(b_j) = \sqrt{\hat\sigma^2 \cdot \big[(X^TX)^{-1}\big]_{jj}}, \qquad \hat\sigma^2 = \frac{SSE}{n-p}tj​=se(bj​)bj​​,se(bj​)=σ^2⋅[(XTX)−1]jj​​,σ^2=n−pSSE​

    với j=0,1,2j=0,1,2j=0,1,2; p=3p=3p=3 là số tham số của mô hình (kể cả hệ số chặn); SSE=∑i=1n(yi−y^i)2SSE = \sum_{i=1}^n (y_i-\hat y_i)^2SSE=∑i=1n​(yi​−y^​i​)2 là tổng bình phương phần dư; [(XTX)−1]jj\big[(X^TX)^{-1}\big]_{jj}[(XTX)−1]jj​ là phần tử trên đường chéo thứ jjj của ma trận nghịch đảo (XTX)−1(X^TX)^{-1}(XTX)−1 (đánh số từ 0, ứng với b0,b1,b2b_0,b_1,b_2b0​,b1​,b2​).

    Cho nnn quan sát (x1i,x2i,yi)(x_{1i},x_{2i},y_i)(x1i​,x2i​,yi​), hãy tính bjb_jbj​, se(bj)se(b_j)se(bj​) và tjt_jtj​ cho j=0,1,2j=0,1,2j=0,1,2.

    Ví dụ (minh họa cách trình bày, không phải kết quả chính xác của test): dòng đầu ứng với b0b_0b0​: 1.7500 0.0527 33.2039, nghĩa là b0=1.7500b_0=1.7500b0​=1.7500, se(b0)=0.0527se(b_0)=0.0527se(b0​)=0.0527, t0=33.2039t_0=33.2039t0​=33.2039.

    • Định dạng đầu vào:
      • Dòng 1: số nguyên nnn (5≤n≤10005 \le n \le 10005≤n≤1000) — đảm bảo n≥4n \ge 4n≥4 để n−p≥2n-p \ge 2n−p≥2.
      • nnn dòng tiếp theo, mỗi dòng gồm ba số thực x1i x2i yix_{1i}\ x_{2i}\ y_ix1i​ x2i​ yi​.
      • Đảm bảo XTXX^TXXTX khả nghịch và SSE>0SSE > 0SSE>0 (mô hình không khớp hoàn hảo 100%, để se(bj)≠0se(b_j) \ne 0se(bj​)=0).
    • Định dạng đầu ra:

      In ra đúng 3 dòng, dòng thứ j+1j+1j+1 (j=0,1,2j=0,1,2j=0,1,2) gồm ba số bj se(bj) tjb_j\ se(b_j)\ t_jbj​ se(bj​) tj​ (cách nhau một khoảng trắng), mỗi số làm tròn đến 4 chữ số thập phân. Thứ tự dòng: b0b_0b0​ rồi b1b_1b1​ rồi b2b_2b2​.

    Ví dụ:

    Đầu vào:

    5
    1 1 6.2
    2 1 7.8
    1 2 9.1
    3 2 12.3
    2 3 13.5
    

    Đầu ra:

    1.7500 0.0527 33.2039
    1.6056 0.0232 69.0842
    2.8556 0.0232 122.8695
    

    Đầu vào:

    6
    0 0 1.2
    1 0 3.1
    0 1 2.2
    1 1 4.3
    2 1 6.1
    1 2 5.4
    

    Đầu ra:

    1.1636 0.0617 18.8726
    1.9318 0.0519 37.2273
    1.1318 0.0519 21.8108
    

    Đang tải editor...