I have a large behavioural data set, and I would like to measure pairwise correlations between several of my outcome variables, and binned categories of outcome variables. However, I have multiple 'subjects' (animal nests), with highly variable numbers of fixed observation periods between each (from 2-30+). What is the best method to test correlation, that properly models the distribution (typically Poisson or negative binomial) and the mixed design?
I have found the Rmcorr R package, but its solution to normal data is to bootstraping. Is there a rank-based method? Should I reshape my data and use a glmm with my outcome variables as predictor and predicted?