library(sitree)
result.sitree <- sitree (tree.df = tr,stand.df = fl,
functions = list(fn.growth ='grow.dbhinc.hgtinc',
fn.mort ='mort.B2007',
fn.recr ='recr.BBG2008',
fn.management ='management.prob',
fn.tree.removal ='mng.tree.removal',
fn.modif = NULL,
fn.prep.common.vars ='prep.common.vars.fun'),
n.periods = 5,
period.length = 5,mng.options = NA,
print.comments = FALSE,
fn.dbh.inc = "dbhi.BN2009",
fn.hgt.inc = "height.korf",
species.spruce = c(1, 2, 3),
species.pine = c(10, 11, 20, 21, 29),species.harw = c(30, 31),
fun.final.felling = "harv.prob",
fun.thinning = "thin.prob",per.vol.harv = 0.83)
PlotDataToLong(result.sitree$plot.data)
Run the code above in your browser using DataLab