From a6a1b4f08390c5684e9ad0115eda34ef9a23e22e Mon Sep 17 00:00:00 2001 From: =?UTF-8?q?Olalla=20D=C3=ADaz-Y=C3=A1=C3=B1ez?= Date: Mon, 6 Nov 2017 12:25:41 +0200 Subject: [PATCH] Fix Diam. cut points Fix Diam. cut points according to Brantseg 1967 (Table 2 & 3, pages 703-704), from where I believe the models were obtained. --- R/FuruVol.R | 8 ++++---- 1 file changed, 4 insertions(+), 4 deletions(-) diff --git a/R/FuruVol.R b/R/FuruVol.R index b6adcdc..9b29fab 100644 --- a/R/FuruVol.R +++ b/R/FuruVol.R @@ -16,17 +16,17 @@ function (dbh, trh, bark, enhet) f.dub <- round(dub, 1) } if (bark == "mb") { - i <- f.dmb < 11.1 + i <- f.dmb < 12.1 vol[i] <- 0.6716 + (0.075708 * (dmb[i]^2)) + (0.029679 * (dmb[i]^2) * h[i]) + (0.004341 * dmb[i] * (h[i]^2)) - i <- f.dmb > 11 + i <- f.dmb > 12 vol[i] <- -6.3954 + (0.178053 * (dmb[i]^2)) + (0.03317 * (dmb[i]^2) * h[i]) - (0.003008 * (dmb[i]^2) * b[i]) } else if (bark == "ub") { - i <- f.dub < 11.1 + i <- f.dub < 12.1 vol[i] <- 2.3393 + (0.010045 * (dub[i]^2)) + (0.038834 * (dub[i]^2) * h[i]) + (0.002732 * dub[i] * (h[i]^2)) - i <- f.dub > 11 + i <- f.dub > 12 vol[i] <- -3.5425 + (0.128182 * (dub[i]^2)) + (0.028268 * (dub[i]^2) * h[i]) + (0.008216 * dub[i] * (h[i]^2)) }