Cách thích hợp nhất để biến đổi tỷ lệ khi chúng là một biến độc lập là gì?


12

Tôi nghĩ rằng tôi đã hiểu vấn đề này, nhưng bây giờ tôi không chắc chắn và tôi muốn kiểm tra với những người khác trước khi tôi tiến hành.

Tôi có hai biến, XY. Ylà một tỷ lệ, và nó không bị giới hạn bởi 0 và 1 và thường được phân phối bình thường. Xlà một tỷ lệ và được giới hạn bởi 0 và 1 (nó chạy từ 0,0 đến 0,6). Khi tôi chạy hồi quy tuyến tính Y ~ Xvà tôi phát hiện ra điều đó XYcó liên quan tuyến tính đáng kể. Càng xa càng tốt.

Nhưng sau đó tôi điều tra thêm và tôi bắt đầu nghĩ rằng có lẽ XYmối quan hệ của họ có thể có nhiều đường cong hơn là tuyến tính. Đối với tôi, nó trông giống như mối quan hệ của XYcó thể gần gũi hơn với Y ~ log(X), Y ~ sqrt(X)hay Y ~ X + X^2, hoặc một cái gì đó như thế. Tôi có những lý do thực nghiệm để cho rằng mối quan hệ có thể là đường cong, nhưng không phải là lý do để cho rằng bất kỳ mối quan hệ phi tuyến tính nào có thể tốt hơn bất kỳ mối quan hệ nào khác.

Tôi có một vài câu hỏi liên quan từ đây. Đầu tiên, Xbiến của tôi có bốn giá trị: 0, 0,2, 0,4 và 0,6. Khi tôi đăng nhập hoặc chuyển đổi căn bậc hai các dữ liệu này, khoảng cách giữa các giá trị này sẽ biến dạng để các giá trị 0 nằm cách xa tất cả các dữ liệu khác. Vì thiếu một cách tốt hơn để hỏi, đây có phải là điều tôi muốn? Tôi cho rằng không phải như vậy, vì tôi nhận được kết quả rất khác nhau tùy thuộc vào mức độ biến dạng mà tôi chấp nhận. Nếu đây không phải là điều tôi muốn, tôi nên tránh bằng cách nào?

Thứ hai, để chuyển đổi nhật ký các dữ liệu này, tôi phải thêm một số tiền vào mỗi Xgiá trị vì bạn không thể lấy nhật ký bằng 0. Khi tôi thêm một lượng rất nhỏ, giả sử 0,001, tôi bị biến dạng rất đáng kể. Khi tôi thêm một số tiền lớn hơn, nói 1, tôi nhận được rất ít biến dạng. Có một số tiền "chính xác" để thêm vào một Xbiến? Hoặc không phù hợp để thêm bất cứ điều gì vào một Xbiến thay vì chọn một phép biến đổi thay thế (ví dụ: khối lập phương) hoặc mô hình (ví dụ hồi quy logistic)?

Những gì tôi đã có thể tìm ra ở đó về vấn đề này khiến tôi cảm thấy mình nên bước đi cẩn thận. Đối với người dùng R, mã này sẽ tạo một số dữ liệu có cấu trúc tương tự như của tôi.

X = rep(c(0, 0.2,0.4,0.6), each = 20)
Y1 = runif(20, 6, 10)
Y2 = runif(20, 6, 9.5)
Y3 = runif(20, 6, 9)
Y4 = runif(20, 6, 8.5)
Y = c(Y4, Y3, Y2, Y1)
plot(Y~X)

Bạn nói rằng Y là một tỷ lệ, nhưng trong dữ liệu của bạn, nó nằm trong khoảng từ 6 đến 10?

Vâng, tôi đã sửa lỗi này ở trên - đó là tỷ lệ, không phải tỷ lệ.
Bajcz

Câu trả lời:


13

Câu hỏi chính về chuyển đổi tỷ lệ (tôi sẽ sử dụng làm biểu tượng, tương tự nhưng không giống với ký hiệu của bạn) cho phép một số nhận xét chung.x

Trong phần tiếp theo tôi cho rằng động cơ chính của việc chuyển đổi tỷ lệ là hiệp phương sai (dự đoán, biến độc lập) là cải thiện sự gần đúng với tuyến tính của mối quan hệ, hoặc nếu trong chế độ thăm dò để có được ý tưởng rõ ràng hơn về hình dạng hoặc thực sự tồn tại bất kỳ mối quan hệ. Như thường lệ, liệu một hiệp phương sai (ví dụ) được phân phối bình thường không quan trọng như vậy. (Tỷ lệ là một mối quan hệ không quá xa của các biến chỉ báo với các giá trị không bao giờ có thể được phân phối bình thường và tỷ lệ cũng nhất thiết phải bị ràng buộc.)0,1

Nếu các tỷ lệ có thể đạt được các số 0 chính xác hoặc chính xác, thì điều cần thiết là một phép biến đổi được xác định cho các giới hạn đó, quy định rõ ràng về , vì log 0 không xác định. Ngoài ra, một hình dạng cụ thể lý tưởng đòi hỏi một số biện minh thực tế (khoa học, thực tế), nhưng thiếu đi theo một số phân tích đơn giản rằng log ( x + c ) rất nhạy cảm với giá trị của c , như bạn gợi ý. logxlog0log(x+c)c

Điều này dễ thấy hơn một chút với logarit đến cơ sở , vì vậy tạm thời hãy xem xét c = 10 k để log 10 ( x + 10 k ) ánh xạ x = 0 đến k .10c=10klog10(x+10k)x=0k

Do đó ánh xạ x = 0 đến 0x = 1 đến khoảng 0,01 , trong khi k = - 3 , c = 0,001 ánh xạ x = 0 đến - 3x = 1 chỉ một smidgen nhiều hơn 0 .k=0,c=1x=00x=10.301k=3,c=0.001x=03x=10

Tương tự, bất cứ điều gì có nghĩa là 0 được ánh xạ tới các giới hạn tương tự, trong khi với xấp xỉ ngày càng tốt x = 1 được ánh xạ thành 0 .k=6,9,0x=10

Vì vậy, giới hạn dưới được kéo dài ra bên ngoài với các hằng số nhỏ hơn và nhỏ hơn , trong khi giới hạn trên vẫn giữ nguyên. Do đó, các phép biến đổi như vậy có thể kéo dài phần dưới của phạm vi và thậm chí tạo ra các ngoại lệ từ các giá trị rất nhỏ tại hoặc gần 0 .c0

Đơn giản, mọi người gợi ý điều này có lẽ tưởng tượng rằng (bây giờ với bất kỳ cơ sở nào bạn thích) nên hoạt động rất giống với log x cho small c , điều này rõ ràng đúng với x lớn , nhưng hoàn toàn không đúng với x nhỏ . Mặt khác, độ dốc lớn hơn và dốc hơn của log x là hàm của xx 0 có thể cắn ở đây rất khó.log(x+c)logxcxxlogxxx0

Có vẻ tốt hơn là tập trung vào các biến đổi thay đổi dần dần gần và (vì các lý do khác, nhưng có liên quan) cũng gần x = 1 .x=0x=1

Căn bậc hai và căn bậc ba và các lũy thừa khác được xác định hoàn toàn tốt cho x = 0 , 1 và thường giúp khi cần kéo dài các giá trị gần 0 . Nhưng những biến đổi này là nổi tiếng và tôi tập trung ở đây nhiều hơn vào một khả năng khác.xpx=0,10

Họ các quyền lực gấp được phổ biến bởi JW Tukey ( Phân tích dữ liệu khám phá , đọc, MA: Addison-Wesley, 1977) là một khả năng, và là . Mặc dù không có ép buộc phải chọn quyền hạn cho phép tên gợi nhiều liên tưởng đơn giản, những lựa chọn p = 1 / 2 (gấp root) và p = 1 / 3 (căn bậc ba gấp) dường như các thành viên hữu ích nhất của gia đình này.xp(1x)pp=1/2p=1/3

Gia đình giống như quen thuộc chuyển đổi logit và thực sự logit là một trường hợp hạn chế như p có xu hướng 0 . Một sự khác biệt quan trọng là sức mạnh gấp được định nghĩa cho x = 0 , 1p 0 .logit x=logxlog(1x)p0x=0,1p0

Các quyền hạn được gấp lại, bao gồm cả logit, xử lý các trường hợp cực đoan gần 1 đối xứng và vẽ đồ thị như các đường cong sigmoid nghịch đảo (một số biểu đồ bên dưới) trộn lẫn hành vi cộng gộp và nhân, lặp lại định tính thường xuyên (nếu không phải là vật lý, sinh học, kinh tế, bất cứ điều gì) sự thật cho hiện tượng cơ bản đó01

  • sự khác biệt từ đến 0,02 có thể là một "vấn đề lớn" (chắc chắn, x thay đổi chỉ 0,01 , nhưng nó cũng tăng gấp đôi)0.010.02x0.01

  • sự khác biệt từ đến 0,99 cũng có thể là một "vấn đề lớn" (chắc chắn, x chỉ thay đổi 0,01 , nhưng "phân số không có" 1 - x cũng giảm một nửa)0.980.99x0.011x

  • sự khác biệt từ đến 0,51 có thể là "thỏa thuận nhỏ hơn" (chắc chắn, x cũng thay đổi 0,01 , nhưng thay đổi tỷ lệ nhỏ hơn nhiều)0.500.51x0.01

Điều này có lẽ dễ nghĩ nhất khi một số động lực cơ bản được tưởng tượng: phần ngày càng tăng của những người biết chữ cần một cú hích lớn để đi, tăng tốc và sau đó chậm lại khi nó tiếp cận với sự không biết chữ của phổ cập. Vì vậy, đường cong trong thời gian có thể giống như một logistic tăng hoặc giảm. Việc tỷ lệ 1 được tiếp cận chậm hơn hoặc chậm hơn tự nhiên là một trong một số động lực cho logit và các mô hình tương tự cho các phản ứng tỷ lệ; mặc dù chúng tôi ở đây tập trung vào các đồng biến tỷ lệ, sigmoids cũng có thể hữu ích ở đây.01

Các quyền hạn được gấp lại như gốc gấp hoặc gốc khối không phải là sigmoid mạnh như logit, nhưng một giá trị có giá trị ở đây là chúng được xác định trực tiếp và dễ dàng mà không cần fud, kyd hay nud cho .x=0,1

Chuyển sang tập dữ liệu giả nhưng có vẻ thực tế của bạn (mà tôi đã nhập vào phần mềm yêu thích của riêng tôi, nhưng phân tích đơn giản ở bất cứ thứ gì tử tế), hóa ra không có biến đổi nào trong số này thực sự giúp ích cả. Nhưng việc vẽ đồ thị dữ liệu đưa ra một cảnh báo rõ ràng rằng ngay cả là một phép biến đổi mạnh mẽ, cũng có thể được nhìn thấy bằng cách vẽ trực tiếp.log(x+0.001)

Hai điểm chính tôi muốn làm là

  1. thường được đề xuất và thường được coi là vô hại, là một biến đổi nguy hiểm trừ khi được hiểu và thường không phù hợp bất cứ khi nào nó mở rộng phân phối cho x nhỏ(trừ khi đây thực sự là hành vi mong muốn).log(x+c)x

  2. Đối với dữ liệu ví dụ của bạn, không có chuyển đổi nào tôi đã cố gắng để giúp đỡ.

Đồng thời, các khả năng khác là hết. (Đáng chú ý, tôi đã không thử căn bậc hai hoặc căn bậc hai, và nhấn mạnh rằng trong nhiều vấn đề khác, đó có thể là ứng cử viên rõ ràng và nghiêm trọng.)

Nhóm đồ thị đầu tiên chỉ đơn giản hiển thị một số phép biến đổi ứng viên cho các tỷ lệ có thể đạt được cả 1 . (Tôi đã sử dụng logarit tự nhiên, nhưng hình dạng không phụ thuộc vào cơ sở được chọn).01

nhập mô tả hình ảnh ở đây

Nhóm biểu đồ thứ hai cho thấy không có biến đổi nào giúp ích nhiều cho dữ liệu mẫu. (Để so sánh, hồi quy đơn giản trên dữ liệu gốc mang lại %, RMSE = 0,994 .)R2=3.7=0.994

nhập mô tả hình ảnh ở đây

Câu đố nhỏ. của bạn được cho là một tỷ lệ, nhưng giá trị của nó là khoảng 6 đến 10 .y610

EDIT: Dữ liệu gốc có thể được vẽ ở đây vì OP đã đăng dữ liệu ngắn gọn, nhưng sau đó đã xóa chúng.

Các chủ đề khác ở đây sử dụng quyền hạn gấp bao gồm

Chuyển đổi dữ liệu tỷ lệ: khi căn bậc hai arcsin không đủ

Hồi quy: Scatterplot với R bình phương thấp và giá trị p cao

Vẽ một tập dữ liệu rất sai lệch


Câu trả lời tuyệt vời và rất kỹ lưỡng. Tôi nghĩ tôi nên nói rằng Ytỷ lệ của tôi chứ không phải là tỷ lệ, có lẽ là một sự khác biệt khá lớn, vì vậy thật tốt khi bạn chỉ ra.
Bajcz

Tỷ lệ bị giới hạn khi tôi xác định chúng. Cảm ơn đã làm rõ, điều này không tạo ra bất kỳ sự khác biệt nào đối với phân tích của tôi (đó là lý do tại sao tôi gắn nhãn cho nó một chi tiết nhỏ).
Nick Cox

2
Nhận xét thêm: Về nguyên tắc, bạn có thể kiểm tra độ cong v.v ... bằng cách sử dụng spline hoặc máy làm mịn, nhưng chỉ với 4 cấp độ khác nhau của công cụ dự đoán là không dễ dàng. Tôi sẽ xem xét hồi quy lượng tử cho dữ liệu của bạn.
Nick Cox

x2x310x=0,1
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.