I'm trying to estimate profiles using tidyLPA, which I have used in the past. Code is below and very straightforward. However, for model_1_class_3 I'm getting the warning "$ operator is invalid for atomic vectors" and for model_1_class_5 I'm getting the warning "subscript out of bounds." Not sure what's causing this given I'm passing three columns of ints into estimate_profiles, but seems like maybe something in how my dataframe is set up? Any advice would be great!
lpa <- data %>% select(CNITot,ICUTot,BESTot) %>% single_imputation(missForest) %>% estimate_profiles(1:5,package = "mclust") %>% compare_solutions(statistics = c("AIC","BIC","BLRT_val"))