## Not run:
# hspider[, 1:6] <- scale(hspider[, 1:6]) # Standardized environmental vars
# set.seed(123)
# p1 <- cao(cbind(Pardlugu, Pardmont, Pardnigr, Pardpull, Zoraspin) ~
# WaterCon + BareSand + FallTwig + CoveMoss + CoveHerb + ReflLux,
# family = poissonff, data = hspider, Rank = 1, df1.nl =
# c(Zoraspin = 2.5, 3), Bestof = 3, Crow1positive = TRUE)
#
# var(latvar(p1)) # Scaled to unit variance # Scaled to unit variance
# c(latvar(p1)) # Estimated site scores
# ## End(Not run)
Run the code above in your browser using DataLab