yield_growth_ptbr.R
In forestmangr: Forest Mensuration and Management

## ----echo = FALSE, message = FALSE--------------------------------------------
knitr::opts_chunk$set(collapse = T, comment = "#>")
knitr::opts_chunk$set(fig.width=7, fig.height=5)
options(tibble.print_min = 6L, tibble.print_max = 6L)
library(forestmangr)
library(dplyr)

## -----------------------------------------------------------------------------
library(forestmangr)
library(dplyr)
data(exfm16)

dados <- exfm16
dados

## -----------------------------------------------------------------------------
index_age <- 64
dados <-  dados %>% 
  nls_table(DH ~ b0 * (1 - exp( -b1 * age )  )^b2, 
            mod_start = c( b0=23, b1=0.03, b2 = 1.3), 
            output = "merge" ) %>% 
  mutate(S = DH *( (  (1- exp( -b1/age ))^b2   ) / 
                     (( 1 - exp(-b1/index_age))^b2 ))  ) %>% 
  select(-b0,-b1,-b2)
head(dados)

## -----------------------------------------------------------------------------
coefs_clutter <- fit_clutter(dados, "age", "DH", "B", "V", "S", "plot")
coefs_clutter

## -----------------------------------------------------------------------------
dados_class <- classify_site(dados, "S", 3, "plot")
head(dados_class)

## -----------------------------------------------------------------------------
dados_est <- est_clutter(dados_class, 20:125,"B", "S", "category_", coefs_clutter) 
dados_est

## -----------------------------------------------------------------------------
est_clutter(dados_class, 20:125,"B", "S", "category_", coefs_clutter,output="plot")