您的訂閱是我製作影片的動力
訂閱點這裡~
影片程式碼
#延續上一支影片
gmm = GMM(iris, 4, dist_mode = "eucl_dist", em_iter = 10) # EM調參
gmm_out <- as.data.frame(gmm$Log_likelihood) #大好
final <- cbind(iris,gmm_out)
library(dplyr)
final <- final %>%
mutate(最大=pmax(V1,V2,V3,V4)) %>%
mutate(分群=ifelse(最大==V1,"c1",
ifelse(最大==V2,"c2",
ifelse(最大==V3,"c3","c4"))))
若內容有誤,還請留言指正,謝謝您的指教