potential.productivity {forestat}

Calculate the potential productivity.

Description

potential.productivity calculate the potential productivity of stand based on model parameters(obtained from the parameterOutput function).

Usage

potential.productivity(
  forestData,
  code = 1,
  age.min = 5,
  age.max = 150,
  left = 0.05,
  right = 100,
  e = 1e-05,
  maxiter = 50
)

Arguments

forestData

A forestData class data

code

Codes for forest types.

age.min

The minimum age of the stand.

age.max

The maximum age of the stand.

left

Solving for the left boundary of the potential productivity.

right

Solving for the right boundary of the potential productivity.

e

Accuracy parameters for solving the stand density index according to Newton's iterative method.

maxiter

Maximum number of iterations parameter for solving the stand density index according to Newton's iteration method.

Details

potential.productivity takes data_BA,data_V parameters as required inputs.

Value

A forestData class in which a data.frame with potential productivity parameters is added.

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))

# Calculate the potential productivity of the forestData object
forestData <- potential.productivity(forestData,code=1,
                                     age.min=5,age.max=150,
                                     left=0.05,right=100,
                                     e=1e-05,maxiter=50)