Mô hình sô-cô-la với tất cả hệ số là ngẫu nhiên

Giờ đây, khi đã lưu mã hóa hiệu ứng trong dữ liệu chocolate, bạn sẽ fit một mô hình với giả định tất cả các hệ số tuân theo phân phối chuẩn. Để làm được điều đó, chúng ta cần tạo vector rpar để truyền vào mlogit(). Phần này hơi khó nên mình đã viết sẵn mã; bạn hãy chạy để xem nó hoạt động thế nào. Sau đó, bạn sẽ tự viết lời gọi mlogit().