取得各成分
上一題配適完成的模型已存於物件 fit_with_covariance。本題的目標是擷取兩個叢集的平均數與共變異數矩陣。
本練習屬於課程
R 中的混合模型
練習說明
- 針對每個叢集，使用
parameters()函式並搭配對應的成分編號，擷取平均數與共變異數矩陣。將結果儲存到物件comp_number。 comp_number的前兩個值對應到平均數。將它們儲存到物件mean_comp_number。comp_number的後四個值對應到共變異數矩陣。把這些值組成 2x2 矩陣。
動手互動練習
試著完成這個範例程式碼，體驗一下這個練習。
# Get the parameters
comp_1 <- ___(___, component = 1)
comp_2 <- parameters(fit_with_covariance, component = ___)
# The means
mean_comp_1 <- ___[1:2]
mean_comp_1
mean_comp_2 <- comp_2[___]
mean_comp_2
# The covariance matrices
covariance_comp_1 <- matrix(comp_1[___], nrow = 2)
covariance_comp_1
covariance_comp_2 <- matrix(comp_2[___], nrow = ___)
covariance_comp_2