first_order_latent_A =~ ag_item1 + ag_item2 + ag_item3 + ag_item4
first_order_latent_B =~ ag_item5 + ag_item6 + ag_item7 + ag_item8
second_order_latent =~ first_order_latent_A + first_order_latent_B "
Now, I am using 'reliability()' from semTools to get the reliability indices for the first-order-latents and I know that I have to use 'reliabilityL2()' for the second-order-latents.
However, I get an error with the following code that says 'Error in ve[[i]][secondFactor, secondFactor, drop = FALSE] : subscript out of bounds'
Is it just not possible to calculate reliability indices with 'reliabilityL2()' for multilevel models? Is there another option how I can arrive at those indices?
Thanks for any advice,