Giải thích đầu ra .L & .Q từ GLM nhị thức âm với dữ liệu phân loại


14

Tôi vừa chạy GLM nhị thức âm và đây là đầu ra:

Call:
glm.nb(formula = small ~ method + site + depth, data = size.dat, 
    init.theta = 1.080668549, link = log)

Deviance Residuals: 
    Min       1Q   Median       3Q      Max  
-2.2452  -0.9973  -0.3028   0.3864   1.8727  

Coefficients:
            Estimate Std. Error z value Pr(>|z|)    
(Intercept)   1.6954     0.1152  14.720  < 2e-16 ***
method.L     -0.6828     0.1637  -4.171 3.04e-05 ***
site.L        0.9952     0.2050   4.854 1.21e-06 ***
site.Q       -0.4634     0.1941  -2.387    0.017 *  
depth.L       0.8951     0.1988   4.502 6.74e-06 ***
depth.Q       0.2060     0.1984   1.038    0.299    
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1   1

(Dispersion parameter for Negative Binomial(1.0807) family taken to be 1)

    Null deviance: 185.1  on 89  degrees of freedom
Residual deviance: 100.8  on 84  degrees of freedom
AIC: 518.24

Dự đoán của tôi là tất cả các loại. Đây có phải là lý do tại sao tôi nhận được .L.Q. Tôi đoán họ đại diện cho các danh mục khác nhau nhưng có ai biết mã tôi có thể sử dụng để gắn nhãn cho chúng trước khi chạy GLM để chúng hiển thị dưới dạng các danh mục khác nhau không?



Tôi thấy mã đó khá khó hiểu là có mã đơn giản hơn không? Các danh mục rất đơn giản: phương pháp 1 và phương pháp 2, trang web 1, trang web 2 và trang web 3 và độ sâu là 5, 10 và 15
Vivienne

Câu trả lời:


12

Các biến của bạn không chỉ được mã hóa thành các yếu tố (để biến chúng thành phân loại), chúng được mã hóa thành các yếu tố theo thứ tự . Sau đó, theo mặc định, R phù hợp với một loạt các hàm đa thức cho các mức của biến. Đầu tiên là tuyến tính ( .L), thứ hai là bậc hai ( .Q), thứ ba (nếu bạn có đủ cấp độ) sẽ là khối, v.v. R sẽ phù hợp với một hàm đa thức ít hơn số cấp trong biến của bạn. Ví dụ: nếu bạn chỉ có hai cấp độ, chỉ xu hướng tuyến tính sẽ phù hợp. Hơn nữa, các cơ sở đa thức được sử dụng là trực giao. (Đối với những gì đáng giá, không có gì trong số này là dành riêng cho R phiên bản hoặc các mô hình nhị thức âm tính, tất cả các phần mềm và loại mô hình hồi quy sẽ làm như vậy.)


Tập trung đặc biệt vào R, nếu bạn muốn biến của bạn được mã hóa như ra lệnh hoặc có thứ tự, bạn sẽ sử dụng yếu tố? :

my.variable <- factor(my.variable, ordered=TRUE)   # an ordered factor
my.variable <- factor(my.variable, ordered=FALSE)  # an unordered factor

Oh tôi đã đặt hàng chúng bây giờ nó có nhãn, cảm ơn bạn rất nhiều gung!
Vivienne
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.