for(ind in 1:nind) { ## nind = nrow(d$X)
      for(yr in 1:nyear) {
        logit(phi[ind,yr]) <-
          phi.sp[species[ind]] + ## effect of species
            phi.year[yr] + ## effect of year
               phi.pc1[species[ind]]*pc1[ind] + ## Effect of morphology, but only the PCA scores
               phi.pc2[species[ind]]*pc2[ind]
of individuals captured)
        } ## (yr in 1:nyear)

PC1_2[pc1[ind],pc2[ind]]
PC1_2[,] ~ dnorm(0, 0.01) # I want it to be a fixed factor

for(ind in 1:nind) { ## nind = nrow(d$X)
      for(yr in 1:nyear) {
        logit(phi[ind,yr]) <-
          phi.sp[species[ind]] + ## effect of species
            phi.year[yr] + ## effect of year
               phi.pc1[species[ind]]*pc1[ind] + ## Effect of morphology, but only the PCA scores
               phi.pc2[species[ind]]*pc2[ind] +
PC1_2[pc1[ind],pc2[ind]]
of individuals captured)
        } ## (yr in 1:year)