含交互項的模型配適
在影片中,我示範了如何在對角共變異數矩陣(也就是不含交互項)的情況下,配適二元高斯混合模型。這裡的目標是重現相同的估計,但這次要把交互項納入考量。
本練習屬於課程
R 中的混合模型
練習說明
使用 flexmix() 函式,對 gender 資料集中的 Weight 與 BMI 兩個變數配適二元高斯混合模型,設定為 2 個叢集,並包含交互項。
動手互動練習
試著完成這個範例程式碼,體驗一下這個練習。
set.seed(1313)
fit_with_covariance <- ___(cbind(___,___) ~ ___,
data = gender,
k = ___,
model = FLXMCmvnorm(diag = ___),
control = list(tolerance = 1e-15, iter.max = 1000))