Các mô hình dọc trong R và WINBUGS hoặc JAGS


8

Tôi đã thử sử dụng R để phù hợp với một số mô hình dọc, chủ yếu là thông qua lmernlmecác gói. Tuy nhiên, dường như nhiều mô hình tiêu chuẩn đang thiếu, chẳng hạn như mô hình chống phụ thuộc hoặc mô hình phân tích nhân tố cho ma trận hiệp phương sai. Những mô hình này có sẵn trong SAS.

Bất cứ ai sẽ giới thiệu các gói khác cho công việc trong R? Tôi không thực sự quan tâm nếu tôi sử dụng để làm việc trong một thế giới thường xuyên hoặc bayes miễn là tôi có sự linh hoạt mô hình hơn. Tôi cũng sẽ thích làm điều đó trong WINBUGS / JAGS.


1
Khó để khuyên - có lẽ bạn nên đăng một liên kết hoặc chỉ định chi tiết những mô hình bạn muốn sử dụng. Bạn có thể làm bất cứ điều gì trong WINBUGS / JAGS vì bạn tự viết mô hình, vì vậy nếu bạn biết chính xác những gì bạn muốn, bạn có thể làm điều đó.
Tò mò

1
Ngoài ra, nếu bạn đang tìm kiếm sự linh hoạt của mô hình, thì chắc chắn hãy tìm kiếm WINBUGS / JAGS. Trong R, bạn chỉ có thể làm những gì họ đã chuẩn bị cho bạn, nhưng trong bayesian bạn hoàn toàn tự do làm bất cứ điều gì (với giá của thời gian học, và bạn cần phải hiểu rất rõ các mô hình).
Tò mò

Tôi thấy bạn cũng đã hỏi câu hỏi này tại r-sig-mixed-models@r-project.org... tốt nhất là không đăng bài chéo, hoặc ít nhất là nói rõ rằng bạn đang đăng chéo. (Tôi nghĩ r-sig-mixed-modelscó thể sẽ có kết quả hơn.)
Ben Bolker

Tôi có một vài bộ dữ liệu và tôi muốn thử nghiệm với các cấu trúc hiệp phương sai khác nhau. Tôi không có một mô hình cụ thể, duy nhất nào tôi muốn phù hợp nhưng thay vào đó tôi muốn thử hầu hết các mô hình cơ bản, chẳng hạn như các mô hình có sẵn trong SAS support.sas.com/rnd/app/ con / glimmix.pdf

1
Nhân tiện, điều này có thể tốt hơn cho trang web Stack Exchange thống kê, vì đây là về việc tìm kiếm các gói hữu ích cho mục đích thống kê cụ thể này hơn là tạo mã mới trong R. Như vậy, có thể có nhiều đối tượng phù hợp hơn trên trang web đó .
Lặp lại

Câu trả lời:



0

Tôi không chắc ý của bạn là gì khi R không có "mô hình phân tích nhân tố cho ma trận hiệp phương sai" - bạn có thể làm rõ những gì bạn muốn sao chép từ SAS không? Theo hiểu biết của tôi thì điều này là khả thi với rất nhiều gói khác nhau trong R.

Về các mô hình chống phụ thuộc, có một cuốn sách về chính chủ đề này có liên quan đến mã R và các ví dụ, tại trang web của tác giả đầu tiên .

Tôi không chắc liệu WinBUGS có mang lại cho bạn bất kỳ may mắn nào không, nhưng tôi sẽ bắt đầu với sách giáo khoa đã nói ở trên - nó dường như có thẩm quyền trên các mô hình chống phụ thuộc. :)


1
Có lẽ OP muốn các cấu trúc hiệp phương sai ở đây: support.sas.com/documentation/cdl/en/statug/63033/HTML/default/iêu ... những nguyên tắc này có thể được xây dựng thông qua corClassescơ sở trong nlme, nhưng nó không tầm thường và như theo như tôi biết chưa được thực hiện ...
Ben Bolker

Ý tôi là, SAS có nhiều lựa chọn hơn cho các cấu trúc hiệp phương sai so với nlme trong R, và tất cả chúng đều trưởng thành; so sánh: support.sas.com/rnd/app/ con / glimmix.pdf

Tôi cũng quen thuộc với cuốn sách được đề cập nhưng với tôi nó giống như một cuốn sách lý thuyết với một số mã R nhưng không có gói R trưởng thành.

Không có ý xúc phạm, nhưng không chính xác khi nói "các tùy chọn cho cấu trúc hiệp phương sai" - Tôi cho rằng bạn muốn nói rằng SAS có nhiều tùy chọn hơn để chỉ định cấu trúc hiệp phương sai trong việc ước tính các tham số của một số loại mô hình tham số nhất định, đúng không? Nhưng chúng tôi không đạt được nhiều tiến bộ ở đây: bạn muốn đạt được điều gì? Tôi đã đọc lướt tài liệu GLMMIX (khá tồi tệ từ góc độ thống kê - rất nhiều bóng bẩy, không nhiều toán học). Bạn thực sự muốn làm gì? (còn tiếp)
Iterator

0

Tôi tin rằng, với một đường cong học tập nhẹ, bạn có thể sử dụng một trong các gói SEM trong R: lavaan, OpenMX, hoặc sem. Tôi chỉ tìm hiểu về SEM và các gói này, nhưng theo tôi lavaancó cú pháp công thức giống như mô hình hóa khác ( lm, lmer) trong R và SEM cho phép bạn thực hiện nhiều điều với cấu trúc hiệp phương sai của mình.

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.