Misure di incrementi non disponibili

increment_borings %>% filter(is.na(IncR5anni))
## # A tibble: 3 x 10
##   id_ads id_fusto  d130 IncR5anni IncR_periodoVariabile ColonneExtra   cld
##    <chr>    <int> <int>     <dbl>                 <dbl>        <lgl> <int>
## 1    G-A       20    19        NA                    NA           NA    20
## 2    P-G       31    21        NA                    NA           NA    20
## 3    P-H       15    38        NA                    NA           NA    40
## # ... with 3 more variables: `misure in mm` <chr>,
## #   `da-sottoc-a-5-anni` <int>, `da-sottoc-a-CIRCA-40-anni` <int>
df <- increment_borings %>% 
  filter(!is.na(IncR5anni)) %>%
  select(1:5) %>%
  left_join(select(stands, id_ads, eta2017, etaIniziale_periodoVariabile)) %>%
  mutate(IncR_su5anni =  IncR5anni / 5,
         delta_eta = eta2017 - etaIniziale_periodoVariabile,
         IncR_suPeriodoVariabile = IncR_periodoVariabile / delta_eta,
         delta_eta = as.factor(delta_eta))
## Joining, by = "id_ads"

Confronto tra i due rilievi incrementali

library(ggplot2)
ggplot(df, aes(IncR_su5anni, IncR_suPeriodoVariabile, color = delta_eta)) + 
  geom_point(alpha = .3) +
  geom_smooth(linetype="dashed") +
  geom_smooth(method = "lm")
## `geom_smooth()` using method = 'loess'
confronto grafico

confronto grafico 

lm(data=df, IncR_suPeriodoVariabile~IncR_su5anni * delta_eta) %>%
  summary()
## 
## Call:
## lm(formula = IncR_suPeriodoVariabile ~ IncR_su5anni * delta_eta, 
##     data = df)
## 
## Residuals:
##     Min      1Q  Median      3Q     Max 
## -0.9628 -0.1868 -0.0225  0.1685  1.6901 
## 
## Coefficients:
##                           Estimate Std. Error t value Pr(>|t|)    
## (Intercept)              -0.017603   0.081948  -0.215 0.830010    
## IncR_su5anni              0.903986   0.057204  15.803  < 2e-16 ***
## delta_eta7                0.487063   0.141704   3.437 0.000638 ***
## delta_eta8                0.236847   0.145816   1.624 0.104964    
## delta_eta12               0.393200   0.089949   4.371 1.51e-05 ***
## delta_eta15               0.295956   0.096443   3.069 0.002270 ** 
## IncR_su5anni:delta_eta7  -0.005364   0.070609  -0.076 0.939479    
## IncR_su5anni:delta_eta8   0.031936   0.072878   0.438 0.661433    
## IncR_su5anni:delta_eta12  0.013947   0.059332   0.235 0.814251    
## IncR_su5anni:delta_eta15  0.083174   0.062936   1.322 0.186935    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.3139 on 486 degrees of freedom
## Multiple R-squared:  0.9348, Adjusted R-squared:  0.9335 
## F-statistic: 773.6 on 9 and 486 DF,  p-value: < 2.2e-16

Conclusioni

La lunghezza del periodo considerato influenza la stima dell’incremento annuo in modo limitato:
* non modifica significativamente la pendenza della retta
* modifica significativamente la quota, per alcuni valori di delta_eta, ma in misura molto piccola