source("wk8funcs.r") x <- c(0.13, 0.21, 0.29, 0.43, 1.32, 1.41, 1.57, 1.77, 4.13, 4.39) fit <- vector("list",4) for (i in 1:length(fit)) { set.seed(i) fit[[i]] <- mix.em (x, 3, 100) cat("\n\nRESPONSIBILITIES FOR FIT",i,"\n\n") print(round(fit[[i]]$r,3)) }