class.plot {forestat}

Calculate the site classes based on stand height growth

Description

class.plot adds new variables: the original height classes and the adjusted height classes. And the existing variables are retained.

Usage

class.plot(
  data,
  model = "Richards",
  interval = 5,
  number = 5,
  maxiter = 1000,
  H_start = c(a = 20, b = 0.05, c = 1),
  BA_start = c(a = 80, b = 1e-04, c = 8, d = 0.1),
  Bio_start = c(a = 450, b = 1e-04, c = 12, d = 0.1)
)

Arguments

data

A data.frame data in which at least four columns are required as input: ID, code, AGE, H.

model

Type of model used for building the H-model (stand height model), options are 'Logistic', 'Richards', 'Korf', 'Gompertz', 'Weibull', or 'Schumacher'.

interval

The initial stand age interval for height classes.

number

The maximum number of initial height classes.

maxiter

The maximum number of iterations to fit the H-model.

H_start

The initial parameters for fitting the H-model, the default value is c(a=20,b=0.05,c=1.0).

BA_start

The initial parameters for fitting the BA-model, the default value is c(a = 80, b = 0.0001, c = 8, d = 0.1).

Bio_start

The initial parameters for fitting the Bio-model, the default value is c(a=450, b=0.0001, c=12, d=0.1).

Details

Input takes a data.frame with three variables ID, AGE, H and returns height classes of every sample (rows in the data.frame).

Value

A data of forestData class with output values, models and model parameters.

Examples

# Load sample data
data("forestData")

# Build a model based on the forestData and return a forestData class object
forestData <- class.plot(forestData,model="Richards",
                         interval=5,number=5,maxiter=1000,
                         H_start=c(a=20,b=0.05,c=1.0))