Thêm hệ số để có được hiệu ứng tương tác - phải làm gì với SE?


13

Tôi có một hồi quy đa biến, bao gồm các tương tác. Ví dụ, để có được ước tính về hiệu quả điều trị cho nhóm ngũ phân vị nghèo nhất, tôi cần thêm các hệ số từ bộ hồi quy điều trị vào hệ số từ biến tương tác (tương tác giữa điều trị và nhóm ngũ phân vị 1). Khi thêm hai hệ số từ hồi quy, làm thế nào để có được các lỗi tiêu chuẩn? Có thể thêm các lỗi tiêu chuẩn từ hai hệ số không? Còn chỉ số t thì sao? Có thể thêm những thứ này là tốt? Tôi đoán là không nhưng tôi không thể tìm thấy bất kỳ hướng dẫn nào về điều này.

Cảm ơn rất nhiều trước sự giúp đỡ của bạn!


Điều này thực sự hữu ích! Tôi đang tìm cách làm một cái gì đó tương tự trong R, nhưng tôi có kích thước mẫu hơi khác nhau giữa các nhóm. Tôi vẫn có thể sử dụng cùng một phương trình để kết hợp hai lỗi để đưa cho tôi Std mới. lỗi? Cảm ơn bạn trước bất kỳ sự giúp đỡ nào Crystal
Crystal

1
Hey @Crystal - chào mừng đến với trang web! Đây là một câu hỏi hay, nhưng bạn nên đặt nó như một câu hỏi mới (nút "Đặt câu hỏi" ở trên cùng bên phải). Ngay bây giờ, bạn đã gửi nó dưới dạng "câu trả lời" cho câu hỏi cũ này. Nếu bạn chỉ sao chép và dán URL của câu hỏi này vào câu hỏi mới của mình, tất cả chúng ta sẽ hiểu bạn đang nói về điều gì.
Matt Parker

Câu trả lời:


10

Tôi nghĩ đây là biểu thức cho :SEbnew

SE12+SE22+2Cov(b1,b2)

Bạn có thể làm việc với sai số chuẩn mới này để tìm số liệu thống kê thử nghiệm mới của bạn để thử nghiệm Ho:β=0


Xin chào Sarah, bạn nên đóng câu hỏi này nếu bạn nghĩ rằng nó đã được trả lời.
suncoolsu

Xin chào - Cảm ơn một lần nữa cho câu trả lời của bạn. Tôi quên đề cập rằng tôi đang sử dụng Stata. Khi tôi thêm hai hệ số với nhau (sử dụng đầu ra từ Stata), tôi cũng có thể thêm các lỗi tiêu chuẩn không? Nếu vậy, thì tôi sẽ có thể có được các lỗi tiêu chuẩn bằng cách chia tổng các hệ số cho tổng các lỗi tiêu chuẩn. Bạn có đồng ý không? Cảm ơn một lần nữa.
Sarah

Sarah, ở Stata, sử dụng chức năng 'lincom'. Giả sử bạn có các biến var1 và var2 và muốn thêm 3 lần hệ số trên var1 và 2 lần hệ số trên var2. Nhập 'lincom 3 * var1 + 2 * var2'. Điều này đưa ra lỗi tiêu chuẩn và khoảng tin cậy cho ước tính này.
Charlie

5

Tôi giả sử bạn có nghĩa là 'hồi quy đa biến', không phải 'đa biến'. "Đa biến" đề cập đến việc có nhiều biến phụ thuộc.

Nó không được coi là thực hành thống kê có thể chấp nhận được để có một công cụ dự đoán liên tục và cắt nó thành các khoảng. Điều này sẽ dẫn đến nhiễu loạn dư và sẽ làm cho các tương tác có ý nghĩa sai lệch vì một số tương tác chỉ có thể phản ánh sự thiếu phù hợp (ở đây, thiếu) một số tác động chính. Có rất nhiều biến thể không giải thích được trong các nhóm tinh túy bên ngoài. Thêm vào đó, thực sự không thể giải thích chính xác các "hiệu ứng tinh túy".

Để so sánh lợi ích, dễ nhất là hình dung chúng là sự khác biệt trong các giá trị dự đoán. Dưới đây là một ví dụ sử dụng rmsgói R.

require(rms)
f <- ols(y ~ x1 + rcs(x2,3)*treat)  # or lrm, cph, psm, Rq, Gls, Glm, ...
# This model allows nonlinearity in x2 and interaction between x2 and treat.
# x2 is modeled as two separate restricted cubic spline functions with 3
# knots or join points in common (one function for the reference treatment
# and one function for the difference in curves between the 2 treatments)
contrast(f, list(treat='B', x2=c(.2, .4)),
            list(treat='A', x2=c(.2, .4)))
# Provides a comparison of treatments at 2 values of x2
anova(f) # provides 2 d.f. interaction test and test of whether treatment
# is effective at ANY value of x2 (combined treat main effect + treat x x2
# interaction - this has 3 d.f. here)

2

RRβRV^R'V^Rχr2r là số lượng hàng trong ma trận của bạn (giả sử rằng các hàng độc lập tuyến tính).


Cảm ơn. Tôi sẽ đặt ra một câu hỏi khác vì tôi không phải là chuyên gia thống kê và không chắc câu hỏi của tôi đã rõ ràng.
Sarah
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.