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, Bestof = 3,
df1.nl = c(Zoraspin=2.5, 3), Crow1positive = TRUE)
index = 1:ncol(p1@y)
lvplot(p1, lcol=index, pcol=index, y=TRUE, las=1)
Run the code above in your browser using DataLab