# NOT RUN {
library(repmis)
source_data("https://github.com/jimclarkatduke/gjam/blob/master/forestTraits.RData?raw=True")
xdata <- forestTraits$xdata
plotByTree <- gjamReZero(forestTraits$treesDeZero) # re-zero
traitTypes <- forestTraits$traitTypes
specByTrait <- forestTraits$specByTrait
tmp <- gjamSpec2Trait(pbys = plotByTree, sbyt = specByTrait,
tTypes = traitTypes)
tTypes <- tmp$traitTypes
traity <- tmp$plotByCWM
censor <- tmp$censor
ml <- list(ng=2000, burnin=500, typeNames = tTypes, censor = censor)
out <- gjam(~ temp + stdage + deficit, xdata, ydata = traity, modelList = ml)
gjamPlot( output = out )
# }
Run the code above in your browser using DataLab