# NOT RUN {
library(dplyr)
library(minpack.lm)
library(timbeR)
tree_scaling <- tree_scaling %>%
mutate(did = di/dbh,
hih = hi/h)
poli5 <- lm(did~hih+I(hih^2)+I(hih^3)+I(hih^4)+I(hih^5),tree_scaling)
coef_poli <- coef(poli5)
dbh <- 25
h <- 20
hi <- 15
poly5_hi(dbh, h, hi, coef_poli)
# }
Run the code above in your browser using DataLab