Yes!! Got likelihood function working (without data), at least for a small SpeciesXSite system. Need to add in data part and phylogenetic transformation of normally distributed Latent Niche Traits (in the prior) by applying the Cholesky decomposition of the phylogenetic covariance matrix.