base.R

setwd("~/Library/CloudStorage/GoogleDrive-icarounam@gmail.com/Mi unidad/Agrosavia/Colaboraciones/Laura/Forestal/Bases")
datos<-read.table("sanfor.csv", header=T, sep=',')
datos$municipio<-as.factor(datos$municipio)
datos$lote<-as.factor(datos$lote)
datos$anof<-as.factor(datos$anof)
datos$bloque<-as.factor(datos$bloque)
datos$forestal<-as.factor(datos$forestal)
datos$planta<-as.factor(datos$planta)
datos$semana<-as.factor(datos$semana)
##Package
library(ggplot2)
library(Rmisc)

## Loading required package: lattice

## Loading required package: plyr

library(dplyr)

## 
## Attaching package: 'dplyr'

## The following objects are masked from 'package:plyr':
## 
##     arrange, count, desc, failwith, id, mutate, rename, summarise,
##     summarize

## The following objects are masked from 'package:stats':
## 
##     filter, lag

## The following objects are masked from 'package:base':
## 
##     intersect, setdiff, setequal, union

library(tidyverse)

## ── Attaching packages ─────────────────────────────────────── tidyverse 1.3.1 ──

## ✓ tibble  3.1.6     ✓ purrr   0.3.4
## ✓ tidyr   1.1.4     ✓ stringr 1.4.0
## ✓ readr   2.1.1     ✓ forcats 0.5.1

## ── Conflicts ────────────────────────────────────────── tidyverse_conflicts() ──
## x dplyr::arrange()   masks plyr::arrange()
## x purrr::compact()   masks plyr::compact()
## x dplyr::count()     masks plyr::count()
## x dplyr::failwith()  masks plyr::failwith()
## x dplyr::filter()    masks stats::filter()
## x dplyr::id()        masks plyr::id()
## x dplyr::lag()       masks stats::lag()
## x dplyr::mutate()    masks plyr::mutate()
## x dplyr::rename()    masks plyr::rename()
## x dplyr::summarise() masks plyr::summarise()
## x dplyr::summarize() masks plyr::summarize()

library(ggpubr)

## 
## Attaching package: 'ggpubr'

## The following object is masked from 'package:plyr':
## 
##     mutate

library(rstatix)

## 
## Attaching package: 'rstatix'

## The following objects are masked from 'package:plyr':
## 
##     desc, mutate

## The following object is masked from 'package:stats':
## 
##     filter

library(emmeans)
library(nlme)

## 
## Attaching package: 'nlme'

## The following object is masked from 'package:dplyr':
## 
##     collapse

## Summarizing data base by year
lista<- datos %>%
  group_by(anof, lote, forestal, planta) %>%
  summarise_all(mean)

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(departamento): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(municipio): argument is not numeric or logical:
## returning NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(semana): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

## Warning in mean.default(bloque): argument is not numeric or logical: returning
## NA

lista %>% print(n=Inf)

## # A tibble: 255 × 11
## # Groups:   anof, lote, forestal [45]
##     anof  lote  forestal planta departamento municipio   ano semana bloque  diam
##     <fct> <fct> <fct>    <fct>         <dbl>     <dbl> <dbl>  <dbl>  <dbl> <dbl>
##   1 2019  Aban… Abarco   1                NA        NA     1     NA     NA  3.32
##   2 2019  Aban… Abarco   2                NA        NA     1     NA     NA  3.23
##   3 2019  Aban… Abarco   3                NA        NA     1     NA     NA  2.53
##   4 2019  Aban… Abarco   4                NA        NA     1     NA     NA  2.93
##   5 2019  Aban… Abarco   5                NA        NA     1     NA     NA  2.94
##   6 2019  Aban… Abarco   6                NA        NA     1     NA     NA  3.94
##   7 2019  Aban… Roble    1                NA        NA     1     NA     NA  6.89
##   8 2019  Aban… Roble    2                NA        NA     1     NA     NA  4.63
##   9 2019  Aban… Roble    3                NA        NA     1     NA     NA  5.07
##  10 2019  Aban… Roble    4                NA        NA     1     NA     NA  7.41
##  11 2019  Aban… Roble    5                NA        NA     1     NA     NA  5.48
##  12 2019  Aban… Termina… 1                NA        NA     1     NA     NA  5.28
##  13 2019  Aban… Termina… 2                NA        NA     1     NA     NA  3.45
##  14 2019  Aban… Termina… 3                NA        NA     1     NA     NA  3.68
##  15 2019  Aban… Termina… 4                NA        NA     1     NA     NA  1.32
##  16 2019  Aban… Termina… 5                NA        NA     1     NA     NA  1.22
##  17 2019  Aban… Termina… 6                NA        NA     1     NA     NA  1.5 
##  18 2019  Caba… Abarco   1                NA        NA     1     NA     NA  2.08
##  19 2019  Caba… Abarco   2                NA        NA     1     NA     NA  2.26
##  20 2019  Caba… Abarco   3                NA        NA     1     NA     NA  2.34
##  21 2019  Caba… Abarco   4                NA        NA     1     NA     NA  4.18
##  22 2019  Caba… Abarco   5                NA        NA     1     NA     NA  3.94
##  23 2019  Caba… Abarco   6                NA        NA     1     NA     NA  2.91
##  24 2019  Caba… Roble    1                NA        NA     1     NA     NA  6.92
##  25 2019  Caba… Roble    2                NA        NA     1     NA     NA  9.06
##  26 2019  Caba… Roble    3                NA        NA     1     NA     NA  6.05
##  27 2019  Caba… Roble    4                NA        NA     1     NA     NA  5.24
##  28 2019  Caba… Roble    5                NA        NA     1     NA     NA  4.51
##  29 2019  Caba… Roble    6                NA        NA     1     NA     NA  6.99
##  30 2019  Caba… Termina… 1                NA        NA     1     NA     NA  3.4 
##  31 2019  Caba… Termina… 2                NA        NA     1     NA     NA  4.2 
##  32 2019  Caba… Termina… 4                NA        NA     1     NA     NA  0.9 
##  33 2019  Liba… Abarco   1                NA        NA     1     NA     NA  4.2 
##  34 2019  Liba… Abarco   2                NA        NA     1     NA     NA  4.01
##  35 2019  Liba… Abarco   3                NA        NA     1     NA     NA  3.26
##  36 2019  Liba… Abarco   4                NA        NA     1     NA     NA  3.34
##  37 2019  Liba… Abarco   5                NA        NA     1     NA     NA  3.54
##  38 2019  Liba… Abarco   6                NA        NA     1     NA     NA  3.62
##  39 2019  Liba… Roble    1                NA        NA     1     NA     NA  3.77
##  40 2019  Liba… Roble    2                NA        NA     1     NA     NA  3.18
##  41 2019  Liba… Roble    4                NA        NA     1     NA     NA  3.69
##  42 2019  Liba… Roble    5                NA        NA     1     NA     NA  4.81
##  43 2019  Liba… Roble    6                NA        NA     1     NA     NA  4.64
##  44 2019  Liba… Termina… 1                NA        NA     1     NA     NA  7.17
##  45 2019  Liba… Termina… 2                NA        NA     1     NA     NA  5.97
##  46 2019  Liba… Termina… 3                NA        NA     1     NA     NA  5.65
##  47 2019  Liba… Termina… 4                NA        NA     1     NA     NA  5.16
##  48 2019  Liba… Termina… 5                NA        NA     1     NA     NA  6.25
##  49 2019  Liba… Termina… 6                NA        NA     1     NA     NA  8.28
##  50 2019  Parc… Abarco   1                NA        NA     1     NA     NA  6.22
##  51 2019  Parc… Abarco   2                NA        NA     1     NA     NA  5.4 
##  52 2019  Parc… Abarco   3                NA        NA     1     NA     NA  5.15
##  53 2019  Parc… Abarco   4                NA        NA     1     NA     NA  8.04
##  54 2019  Parc… Abarco   5                NA        NA     1     NA     NA  6.36
##  55 2019  Parc… Abarco   6                NA        NA     1     NA     NA  5.97
##  56 2019  Parc… Roble    1                NA        NA     1     NA     NA  8.65
##  57 2019  Parc… Roble    2                NA        NA     1     NA     NA  7.72
##  58 2019  Parc… Roble    3                NA        NA     1     NA     NA  7.13
##  59 2019  Parc… Roble    4                NA        NA     1     NA     NA  9.51
##  60 2019  Parc… Roble    5                NA        NA     1     NA     NA 10.7 
##  61 2019  Parc… Roble    6                NA        NA     1     NA     NA  9.76
##  62 2019  Parc… Termina… 1                NA        NA     1     NA     NA 11.7 
##  63 2019  Parc… Termina… 2                NA        NA     1     NA     NA 11.3 
##  64 2019  Parc… Termina… 3                NA        NA     1     NA     NA 12.6 
##  65 2019  Parc… Termina… 4                NA        NA     1     NA     NA 11.1 
##  66 2019  Parc… Termina… 5                NA        NA     1     NA     NA 10.8 
##  67 2019  Parc… Termina… 6                NA        NA     1     NA     NA 11.1 
##  68 2019  Sjose Abarco   1                NA        NA     1     NA     NA  2.51
##  69 2019  Sjose Abarco   2                NA        NA     1     NA     NA  2.51
##  70 2019  Sjose Abarco   3                NA        NA     1     NA     NA  2.75
##  71 2019  Sjose Abarco   4                NA        NA     1     NA     NA  3.56
##  72 2019  Sjose Abarco   5                NA        NA     1     NA     NA  2.77
##  73 2019  Sjose Abarco   6                NA        NA     1     NA     NA  2.54
##  74 2019  Sjose Roble    1                NA        NA     1     NA     NA  6.54
##  75 2019  Sjose Roble    2                NA        NA     1     NA     NA  6.39
##  76 2019  Sjose Roble    3                NA        NA     1     NA     NA  5.34
##  77 2019  Sjose Roble    4                NA        NA     1     NA     NA  5.34
##  78 2019  Sjose Roble    5                NA        NA     1     NA     NA  4.92
##  79 2019  Sjose Roble    6                NA        NA     1     NA     NA  6.02
##  80 2019  Sjose Termina… 1                NA        NA     1     NA     NA  5.79
##  81 2019  Sjose Termina… 2                NA        NA     1     NA     NA  6.44
##  82 2019  Sjose Termina… 3                NA        NA     1     NA     NA  6.17
##  83 2019  Sjose Termina… 4                NA        NA     1     NA     NA  6.58
##  84 2019  Sjose Termina… 5                NA        NA     1     NA     NA  6.44
##  85 2019  Sjose Termina… 6                NA        NA     1     NA     NA  6.89
##  86 2020  Aban… Abarco   1                NA        NA     2     NA     NA  6.25
##  87 2020  Aban… Abarco   2                NA        NA     2     NA     NA  6.4 
##  88 2020  Aban… Abarco   3                NA        NA     2     NA     NA  5.98
##  89 2020  Aban… Abarco   4                NA        NA     2     NA     NA  5.7 
##  90 2020  Aban… Abarco   5                NA        NA     2     NA     NA  5.61
##  91 2020  Aban… Abarco   6                NA        NA     2     NA     NA  7   
##  92 2020  Aban… Roble    1                NA        NA     2     NA     NA  8.76
##  93 2020  Aban… Roble    2                NA        NA     2     NA     NA  7.17
##  94 2020  Aban… Roble    3                NA        NA     2     NA     NA  7.22
##  95 2020  Aban… Roble    4                NA        NA     2     NA     NA 11.8 
##  96 2020  Aban… Roble    5                NA        NA     2     NA     NA 10.9 
##  97 2020  Aban… Termina… 1                NA        NA     2     NA     NA  9.54
##  98 2020  Aban… Termina… 2                NA        NA     2     NA     NA  7.04
##  99 2020  Aban… Termina… 3                NA        NA     2     NA     NA  7.45
## 100 2020  Aban… Termina… 4                NA        NA     2     NA     NA  4.54
## 101 2020  Aban… Termina… 5                NA        NA     2     NA     NA  4.27
## 102 2020  Aban… Termina… 6                NA        NA     2     NA     NA  6.27
## 103 2020  Caba… Abarco   1                NA        NA     2     NA     NA  5.07
## 104 2020  Caba… Abarco   2                NA        NA     2     NA     NA  5.58
## 105 2020  Caba… Abarco   3                NA        NA     2     NA     NA  5.58
## 106 2020  Caba… Abarco   4                NA        NA     2     NA     NA  8.53
## 107 2020  Caba… Abarco   5                NA        NA     2     NA     NA  8.5 
## 108 2020  Caba… Abarco   6                NA        NA     2     NA     NA  7   
## 109 2020  Caba… Roble    1                NA        NA     2     NA     NA  8.48
## 110 2020  Caba… Roble    2                NA        NA     2     NA     NA 11.4 
## 111 2020  Caba… Roble    3                NA        NA     2     NA     NA  6.67
## 112 2020  Caba… Roble    4                NA        NA     2     NA     NA  6.53
## 113 2020  Caba… Roble    5                NA        NA     2     NA     NA 10.3 
## 114 2020  Caba… Roble    6                NA        NA     2     NA     NA  9.02
## 115 2020  Caba… Termina… 1                NA        NA     2     NA     NA  5.53
## 116 2020  Caba… Termina… 2                NA        NA     2     NA     NA  6   
## 117 2020  Caba… Termina… 4                NA        NA     2     NA     NA  5.8 
## 118 2020  Liba… Abarco   1                NA        NA     2     NA     NA  6.18
## 119 2020  Liba… Abarco   2                NA        NA     2     NA     NA  6.47
## 120 2020  Liba… Abarco   3                NA        NA     2     NA     NA  5.68
## 121 2020  Liba… Abarco   4                NA        NA     2     NA     NA  5.62
## 122 2020  Liba… Abarco   5                NA        NA     2     NA     NA  6.28
## 123 2020  Liba… Abarco   6                NA        NA     2     NA     NA  5.62
## 124 2020  Liba… Roble    1                NA        NA     2     NA     NA  5.42
## 125 2020  Liba… Roble    2                NA        NA     2     NA     NA  4.35
## 126 2020  Liba… Roble    4                NA        NA     2     NA     NA  5.32
## 127 2020  Liba… Roble    5                NA        NA     2     NA     NA  5.82
## 128 2020  Liba… Roble    6                NA        NA     2     NA     NA  4.65
## 129 2020  Liba… Termina… 1                NA        NA     2     NA     NA 10.4 
## 130 2020  Liba… Termina… 2                NA        NA     2     NA     NA  9.41
## 131 2020  Liba… Termina… 3                NA        NA     2     NA     NA  9.09
## 132 2020  Liba… Termina… 4                NA        NA     2     NA     NA  8.98
## 133 2020  Liba… Termina… 5                NA        NA     2     NA     NA 11.5 
## 134 2020  Liba… Termina… 6                NA        NA     2     NA     NA 12.4 
## 135 2020  Parc… Abarco   1                NA        NA     2     NA     NA 10.8 
## 136 2020  Parc… Abarco   2                NA        NA     2     NA     NA 10.8 
## 137 2020  Parc… Abarco   3                NA        NA     2     NA     NA 10.4 
## 138 2020  Parc… Abarco   4                NA        NA     2     NA     NA 13.1 
## 139 2020  Parc… Abarco   5                NA        NA     2     NA     NA 11.1 
## 140 2020  Parc… Abarco   6                NA        NA     2     NA     NA 11.5 
## 141 2020  Parc… Roble    1                NA        NA     2     NA     NA 11.4 
## 142 2020  Parc… Roble    2                NA        NA     2     NA     NA 11.0 
## 143 2020  Parc… Roble    3                NA        NA     2     NA     NA  9.23
## 144 2020  Parc… Roble    4                NA        NA     2     NA     NA 11.8 
## 145 2020  Parc… Roble    5                NA        NA     2     NA     NA 14.8 
## 146 2020  Parc… Roble    6                NA        NA     2     NA     NA 12.0 
## 147 2020  Parc… Termina… 1                NA        NA     2     NA     NA 18.2 
## 148 2020  Parc… Termina… 2                NA        NA     2     NA     NA 17.1 
## 149 2020  Parc… Termina… 3                NA        NA     2     NA     NA 17.6 
## 150 2020  Parc… Termina… 4                NA        NA     2     NA     NA 17.9 
## 151 2020  Parc… Termina… 5                NA        NA     2     NA     NA 18.0 
## 152 2020  Parc… Termina… 6                NA        NA     2     NA     NA 16.8 
## 153 2020  Sjose Abarco   1                NA        NA     2     NA     NA  4.43
## 154 2020  Sjose Abarco   2                NA        NA     2     NA     NA  4.75
## 155 2020  Sjose Abarco   3                NA        NA     2     NA     NA  4.94
## 156 2020  Sjose Abarco   4                NA        NA     2     NA     NA  6.04
## 157 2020  Sjose Abarco   5                NA        NA     2     NA     NA  5.67
## 158 2020  Sjose Abarco   6                NA        NA     2     NA     NA  6.4 
## 159 2020  Sjose Roble    1                NA        NA     2     NA     NA  9.75
## 160 2020  Sjose Roble    2                NA        NA     2     NA     NA 10.4 
## 161 2020  Sjose Roble    3                NA        NA     2     NA     NA  9.28
## 162 2020  Sjose Roble    4                NA        NA     2     NA     NA  9.24
## 163 2020  Sjose Roble    5                NA        NA     2     NA     NA  8.89
## 164 2020  Sjose Roble    6                NA        NA     2     NA     NA 10.2 
## 165 2020  Sjose Termina… 1                NA        NA     2     NA     NA 10.7 
## 166 2020  Sjose Termina… 2                NA        NA     2     NA     NA 11.1 
## 167 2020  Sjose Termina… 3                NA        NA     2     NA     NA  9.43
## 168 2020  Sjose Termina… 4                NA        NA     2     NA     NA  9.21
## 169 2020  Sjose Termina… 5                NA        NA     2     NA     NA 10.2 
## 170 2020  Sjose Termina… 6                NA        NA     2     NA     NA  9.72
## 171 2021  Aban… Abarco   1                NA        NA     3     NA     NA  9.12
## 172 2021  Aban… Abarco   2                NA        NA     3     NA     NA  9.13
## 173 2021  Aban… Abarco   3                NA        NA     3     NA     NA  7.7 
## 174 2021  Aban… Abarco   4                NA        NA     3     NA     NA  8.1 
## 175 2021  Aban… Abarco   5                NA        NA     3     NA     NA  7.96
## 176 2021  Aban… Abarco   6                NA        NA     3     NA     NA 10.1 
## 177 2021  Aban… Roble    1                NA        NA     3     NA     NA 12.9 
## 178 2021  Aban… Roble    2                NA        NA     3     NA     NA 10.1 
## 179 2021  Aban… Roble    3                NA        NA     3     NA     NA 11.0 
## 180 2021  Aban… Roble    4                NA        NA     3     NA     NA 17.8 
## 181 2021  Aban… Roble    5                NA        NA     3     NA     NA 14.2 
## 182 2021  Aban… Termina… 1                NA        NA     3     NA     NA 12.9 
## 183 2021  Aban… Termina… 2                NA        NA     3     NA     NA  9   
## 184 2021  Aban… Termina… 3                NA        NA     3     NA     NA 10.0 
## 185 2021  Aban… Termina… 4                NA        NA     3     NA     NA  8.65
## 186 2021  Aban… Termina… 5                NA        NA     3     NA     NA  7.58
## 187 2021  Aban… Termina… 6                NA        NA     3     NA     NA  9.3 
## 188 2021  Caba… Abarco   1                NA        NA     3     NA     NA  5.98
## 189 2021  Caba… Abarco   2                NA        NA     3     NA     NA  6.07
## 190 2021  Caba… Abarco   3                NA        NA     3     NA     NA  6.57
## 191 2021  Caba… Abarco   4                NA        NA     3     NA     NA  9.82
## 192 2021  Caba… Abarco   5                NA        NA     3     NA     NA 10.1 
## 193 2021  Caba… Abarco   6                NA        NA     3     NA     NA  9.55
## 194 2021  Caba… Roble    1                NA        NA     3     NA     NA 15.2 
## 195 2021  Caba… Roble    2                NA        NA     3     NA     NA 19.6 
## 196 2021  Caba… Roble    3                NA        NA     3     NA     NA 14.9 
## 197 2021  Caba… Roble    4                NA        NA     3     NA     NA 17.4 
## 198 2021  Caba… Roble    5                NA        NA     3     NA     NA 18.3 
## 199 2021  Caba… Roble    6                NA        NA     3     NA     NA 16.1 
## 200 2021  Caba… Termina… 1                NA        NA     3     NA     NA  8.53
## 201 2021  Caba… Termina… 2                NA        NA     3     NA     NA  7.3 
## 202 2021  Caba… Termina… 4                NA        NA     3     NA     NA  9.17
## 203 2021  Liba… Abarco   1                NA        NA     3     NA     NA  6.93
## 204 2021  Liba… Abarco   2                NA        NA     3     NA     NA  7.57
## 205 2021  Liba… Abarco   3                NA        NA     3     NA     NA  8   
## 206 2021  Liba… Abarco   4                NA        NA     3     NA     NA  6.87
## 207 2021  Liba… Abarco   5                NA        NA     3     NA     NA  7.8 
## 208 2021  Liba… Abarco   6                NA        NA     3     NA     NA  7   
## 209 2021  Liba… Roble    1                NA        NA     3     NA     NA  7.45
## 210 2021  Liba… Roble    2                NA        NA     3     NA     NA  5.05
## 211 2021  Liba… Roble    4                NA        NA     3     NA     NA  7.15
## 212 2021  Liba… Roble    5                NA        NA     3     NA     NA  6.4 
## 213 2021  Liba… Roble    6                NA        NA     3     NA     NA  5.45
## 214 2021  Liba… Termina… 1                NA        NA     3     NA     NA 14.4 
## 215 2021  Liba… Termina… 2                NA        NA     3     NA     NA 12.9 
## 216 2021  Liba… Termina… 3                NA        NA     3     NA     NA 11.9 
## 217 2021  Liba… Termina… 4                NA        NA     3     NA     NA 12.4 
## 218 2021  Liba… Termina… 5                NA        NA     3     NA     NA 15.7 
## 219 2021  Liba… Termina… 6                NA        NA     3     NA     NA 20.0 
## 220 2021  Parc… Abarco   1                NA        NA     3     NA     NA 15.4 
## 221 2021  Parc… Abarco   2                NA        NA     3     NA     NA 16.1 
## 222 2021  Parc… Abarco   3                NA        NA     3     NA     NA 16.0 
## 223 2021  Parc… Abarco   4                NA        NA     3     NA     NA 18.0 
## 224 2021  Parc… Abarco   5                NA        NA     3     NA     NA 15.5 
## 225 2021  Parc… Abarco   6                NA        NA     3     NA     NA 16.3 
## 226 2021  Parc… Roble    1                NA        NA     3     NA     NA 14.8 
## 227 2021  Parc… Roble    2                NA        NA     3     NA     NA 14.7 
## 228 2021  Parc… Roble    3                NA        NA     3     NA     NA 12.9 
## 229 2021  Parc… Roble    4                NA        NA     3     NA     NA 16.1 
## 230 2021  Parc… Roble    5                NA        NA     3     NA     NA 20.2 
## 231 2021  Parc… Roble    6                NA        NA     3     NA     NA 16.1 
## 232 2021  Parc… Termina… 1                NA        NA     3     NA     NA 22.6 
## 233 2021  Parc… Termina… 2                NA        NA     3     NA     NA 20.9 
## 234 2021  Parc… Termina… 3                NA        NA     3     NA     NA 23.5 
## 235 2021  Parc… Termina… 4                NA        NA     3     NA     NA 22.0 
## 236 2021  Parc… Termina… 5                NA        NA     3     NA     NA 22.9 
## 237 2021  Parc… Termina… 6                NA        NA     3     NA     NA 20.5 
## 238 2021  Sjose Abarco   1                NA        NA     3     NA     NA  6.94
## 239 2021  Sjose Abarco   2                NA        NA     3     NA     NA  7.09
## 240 2021  Sjose Abarco   3                NA        NA     3     NA     NA  7.42
## 241 2021  Sjose Abarco   4                NA        NA     3     NA     NA  9.21
## 242 2021  Sjose Abarco   5                NA        NA     3     NA     NA  9.12
## 243 2021  Sjose Abarco   6                NA        NA     3     NA     NA  7.83
## 244 2021  Sjose Roble    1                NA        NA     3     NA     NA 14.0 
## 245 2021  Sjose Roble    2                NA        NA     3     NA     NA 15.3 
## 246 2021  Sjose Roble    3                NA        NA     3     NA     NA 13.9 
## 247 2021  Sjose Roble    4                NA        NA     3     NA     NA 14.2 
## 248 2021  Sjose Roble    5                NA        NA     3     NA     NA 13.5 
## 249 2021  Sjose Roble    6                NA        NA     3     NA     NA 15.8 
## 250 2021  Sjose Termina… 1                NA        NA     3     NA     NA 14.4 
## 251 2021  Sjose Termina… 2                NA        NA     3     NA     NA 15.4 
## 252 2021  Sjose Termina… 3                NA        NA     3     NA     NA 11.5 
## 253 2021  Sjose Termina… 4                NA        NA     3     NA     NA 11.9 
## 254 2021  Sjose Termina… 5                NA        NA     3     NA     NA 14.1 
## 255 2021  Sjose Termina… 6                NA        NA     3     NA     NA 13.4 
## # … with 1 more variable: alt <dbl>

lista2 <- lista                                   # Replicate example data
lista2 <- lista2 %>%                              # Create numbering variable
  group_by(anof) %>%
  mutate(id = row_number())
lista2 %>% print(n=Inf)      

## # A tibble: 255 × 12
## # Groups:   anof [3]
##     anof  lote  forestal planta departamento municipio   ano semana bloque  diam
##     <fct> <fct> <fct>    <fct>         <dbl>     <dbl> <dbl>  <dbl>  <dbl> <dbl>
##   1 2019  Aban… Abarco   1                NA        NA     1     NA     NA  3.32
##   2 2019  Aban… Abarco   2                NA        NA     1     NA     NA  3.23
##   3 2019  Aban… Abarco   3                NA        NA     1     NA     NA  2.53
##   4 2019  Aban… Abarco   4                NA        NA     1     NA     NA  2.93
##   5 2019  Aban… Abarco   5                NA        NA     1     NA     NA  2.94
##   6 2019  Aban… Abarco   6                NA        NA     1     NA     NA  3.94
##   7 2019  Aban… Roble    1                NA        NA     1     NA     NA  6.89
##   8 2019  Aban… Roble    2                NA        NA     1     NA     NA  4.63
##   9 2019  Aban… Roble    3                NA        NA     1     NA     NA  5.07
##  10 2019  Aban… Roble    4                NA        NA     1     NA     NA  7.41
##  11 2019  Aban… Roble    5                NA        NA     1     NA     NA  5.48
##  12 2019  Aban… Termina… 1                NA        NA     1     NA     NA  5.28
##  13 2019  Aban… Termina… 2                NA        NA     1     NA     NA  3.45
##  14 2019  Aban… Termina… 3                NA        NA     1     NA     NA  3.68
##  15 2019  Aban… Termina… 4                NA        NA     1     NA     NA  1.32
##  16 2019  Aban… Termina… 5                NA        NA     1     NA     NA  1.22
##  17 2019  Aban… Termina… 6                NA        NA     1     NA     NA  1.5 
##  18 2019  Caba… Abarco   1                NA        NA     1     NA     NA  2.08
##  19 2019  Caba… Abarco   2                NA        NA     1     NA     NA  2.26
##  20 2019  Caba… Abarco   3                NA        NA     1     NA     NA  2.34
##  21 2019  Caba… Abarco   4                NA        NA     1     NA     NA  4.18
##  22 2019  Caba… Abarco   5                NA        NA     1     NA     NA  3.94
##  23 2019  Caba… Abarco   6                NA        NA     1     NA     NA  2.91
##  24 2019  Caba… Roble    1                NA        NA     1     NA     NA  6.92
##  25 2019  Caba… Roble    2                NA        NA     1     NA     NA  9.06
##  26 2019  Caba… Roble    3                NA        NA     1     NA     NA  6.05
##  27 2019  Caba… Roble    4                NA        NA     1     NA     NA  5.24
##  28 2019  Caba… Roble    5                NA        NA     1     NA     NA  4.51
##  29 2019  Caba… Roble    6                NA        NA     1     NA     NA  6.99
##  30 2019  Caba… Termina… 1                NA        NA     1     NA     NA  3.4 
##  31 2019  Caba… Termina… 2                NA        NA     1     NA     NA  4.2 
##  32 2019  Caba… Termina… 4                NA        NA     1     NA     NA  0.9 
##  33 2019  Liba… Abarco   1                NA        NA     1     NA     NA  4.2 
##  34 2019  Liba… Abarco   2                NA        NA     1     NA     NA  4.01
##  35 2019  Liba… Abarco   3                NA        NA     1     NA     NA  3.26
##  36 2019  Liba… Abarco   4                NA        NA     1     NA     NA  3.34
##  37 2019  Liba… Abarco   5                NA        NA     1     NA     NA  3.54
##  38 2019  Liba… Abarco   6                NA        NA     1     NA     NA  3.62
##  39 2019  Liba… Roble    1                NA        NA     1     NA     NA  3.77
##  40 2019  Liba… Roble    2                NA        NA     1     NA     NA  3.18
##  41 2019  Liba… Roble    4                NA        NA     1     NA     NA  3.69
##  42 2019  Liba… Roble    5                NA        NA     1     NA     NA  4.81
##  43 2019  Liba… Roble    6                NA        NA     1     NA     NA  4.64
##  44 2019  Liba… Termina… 1                NA        NA     1     NA     NA  7.17
##  45 2019  Liba… Termina… 2                NA        NA     1     NA     NA  5.97
##  46 2019  Liba… Termina… 3                NA        NA     1     NA     NA  5.65
##  47 2019  Liba… Termina… 4                NA        NA     1     NA     NA  5.16
##  48 2019  Liba… Termina… 5                NA        NA     1     NA     NA  6.25
##  49 2019  Liba… Termina… 6                NA        NA     1     NA     NA  8.28
##  50 2019  Parc… Abarco   1                NA        NA     1     NA     NA  6.22
##  51 2019  Parc… Abarco   2                NA        NA     1     NA     NA  5.4 
##  52 2019  Parc… Abarco   3                NA        NA     1     NA     NA  5.15
##  53 2019  Parc… Abarco   4                NA        NA     1     NA     NA  8.04
##  54 2019  Parc… Abarco   5                NA        NA     1     NA     NA  6.36
##  55 2019  Parc… Abarco   6                NA        NA     1     NA     NA  5.97
##  56 2019  Parc… Roble    1                NA        NA     1     NA     NA  8.65
##  57 2019  Parc… Roble    2                NA        NA     1     NA     NA  7.72
##  58 2019  Parc… Roble    3                NA        NA     1     NA     NA  7.13
##  59 2019  Parc… Roble    4                NA        NA     1     NA     NA  9.51
##  60 2019  Parc… Roble    5                NA        NA     1     NA     NA 10.7 
##  61 2019  Parc… Roble    6                NA        NA     1     NA     NA  9.76
##  62 2019  Parc… Termina… 1                NA        NA     1     NA     NA 11.7 
##  63 2019  Parc… Termina… 2                NA        NA     1     NA     NA 11.3 
##  64 2019  Parc… Termina… 3                NA        NA     1     NA     NA 12.6 
##  65 2019  Parc… Termina… 4                NA        NA     1     NA     NA 11.1 
##  66 2019  Parc… Termina… 5                NA        NA     1     NA     NA 10.8 
##  67 2019  Parc… Termina… 6                NA        NA     1     NA     NA 11.1 
##  68 2019  Sjose Abarco   1                NA        NA     1     NA     NA  2.51
##  69 2019  Sjose Abarco   2                NA        NA     1     NA     NA  2.51
##  70 2019  Sjose Abarco   3                NA        NA     1     NA     NA  2.75
##  71 2019  Sjose Abarco   4                NA        NA     1     NA     NA  3.56
##  72 2019  Sjose Abarco   5                NA        NA     1     NA     NA  2.77
##  73 2019  Sjose Abarco   6                NA        NA     1     NA     NA  2.54
##  74 2019  Sjose Roble    1                NA        NA     1     NA     NA  6.54
##  75 2019  Sjose Roble    2                NA        NA     1     NA     NA  6.39
##  76 2019  Sjose Roble    3                NA        NA     1     NA     NA  5.34
##  77 2019  Sjose Roble    4                NA        NA     1     NA     NA  5.34
##  78 2019  Sjose Roble    5                NA        NA     1     NA     NA  4.92
##  79 2019  Sjose Roble    6                NA        NA     1     NA     NA  6.02
##  80 2019  Sjose Termina… 1                NA        NA     1     NA     NA  5.79
##  81 2019  Sjose Termina… 2                NA        NA     1     NA     NA  6.44
##  82 2019  Sjose Termina… 3                NA        NA     1     NA     NA  6.17
##  83 2019  Sjose Termina… 4                NA        NA     1     NA     NA  6.58
##  84 2019  Sjose Termina… 5                NA        NA     1     NA     NA  6.44
##  85 2019  Sjose Termina… 6                NA        NA     1     NA     NA  6.89
##  86 2020  Aban… Abarco   1                NA        NA     2     NA     NA  6.25
##  87 2020  Aban… Abarco   2                NA        NA     2     NA     NA  6.4 
##  88 2020  Aban… Abarco   3                NA        NA     2     NA     NA  5.98
##  89 2020  Aban… Abarco   4                NA        NA     2     NA     NA  5.7 
##  90 2020  Aban… Abarco   5                NA        NA     2     NA     NA  5.61
##  91 2020  Aban… Abarco   6                NA        NA     2     NA     NA  7   
##  92 2020  Aban… Roble    1                NA        NA     2     NA     NA  8.76
##  93 2020  Aban… Roble    2                NA        NA     2     NA     NA  7.17
##  94 2020  Aban… Roble    3                NA        NA     2     NA     NA  7.22
##  95 2020  Aban… Roble    4                NA        NA     2     NA     NA 11.8 
##  96 2020  Aban… Roble    5                NA        NA     2     NA     NA 10.9 
##  97 2020  Aban… Termina… 1                NA        NA     2     NA     NA  9.54
##  98 2020  Aban… Termina… 2                NA        NA     2     NA     NA  7.04
##  99 2020  Aban… Termina… 3                NA        NA     2     NA     NA  7.45
## 100 2020  Aban… Termina… 4                NA        NA     2     NA     NA  4.54
## 101 2020  Aban… Termina… 5                NA        NA     2     NA     NA  4.27
## 102 2020  Aban… Termina… 6                NA        NA     2     NA     NA  6.27
## 103 2020  Caba… Abarco   1                NA        NA     2     NA     NA  5.07
## 104 2020  Caba… Abarco   2                NA        NA     2     NA     NA  5.58
## 105 2020  Caba… Abarco   3                NA        NA     2     NA     NA  5.58
## 106 2020  Caba… Abarco   4                NA        NA     2     NA     NA  8.53
## 107 2020  Caba… Abarco   5                NA        NA     2     NA     NA  8.5 
## 108 2020  Caba… Abarco   6                NA        NA     2     NA     NA  7   
## 109 2020  Caba… Roble    1                NA        NA     2     NA     NA  8.48
## 110 2020  Caba… Roble    2                NA        NA     2     NA     NA 11.4 
## 111 2020  Caba… Roble    3                NA        NA     2     NA     NA  6.67
## 112 2020  Caba… Roble    4                NA        NA     2     NA     NA  6.53
## 113 2020  Caba… Roble    5                NA        NA     2     NA     NA 10.3 
## 114 2020  Caba… Roble    6                NA        NA     2     NA     NA  9.02
## 115 2020  Caba… Termina… 1                NA        NA     2     NA     NA  5.53
## 116 2020  Caba… Termina… 2                NA        NA     2     NA     NA  6   
## 117 2020  Caba… Termina… 4                NA        NA     2     NA     NA  5.8 
## 118 2020  Liba… Abarco   1                NA        NA     2     NA     NA  6.18
## 119 2020  Liba… Abarco   2                NA        NA     2     NA     NA  6.47
## 120 2020  Liba… Abarco   3                NA        NA     2     NA     NA  5.68
## 121 2020  Liba… Abarco   4                NA        NA     2     NA     NA  5.62
## 122 2020  Liba… Abarco   5                NA        NA     2     NA     NA  6.28
## 123 2020  Liba… Abarco   6                NA        NA     2     NA     NA  5.62
## 124 2020  Liba… Roble    1                NA        NA     2     NA     NA  5.42
## 125 2020  Liba… Roble    2                NA        NA     2     NA     NA  4.35
## 126 2020  Liba… Roble    4                NA        NA     2     NA     NA  5.32
## 127 2020  Liba… Roble    5                NA        NA     2     NA     NA  5.82
## 128 2020  Liba… Roble    6                NA        NA     2     NA     NA  4.65
## 129 2020  Liba… Termina… 1                NA        NA     2     NA     NA 10.4 
## 130 2020  Liba… Termina… 2                NA        NA     2     NA     NA  9.41
## 131 2020  Liba… Termina… 3                NA        NA     2     NA     NA  9.09
## 132 2020  Liba… Termina… 4                NA        NA     2     NA     NA  8.98
## 133 2020  Liba… Termina… 5                NA        NA     2     NA     NA 11.5 
## 134 2020  Liba… Termina… 6                NA        NA     2     NA     NA 12.4 
## 135 2020  Parc… Abarco   1                NA        NA     2     NA     NA 10.8 
## 136 2020  Parc… Abarco   2                NA        NA     2     NA     NA 10.8 
## 137 2020  Parc… Abarco   3                NA        NA     2     NA     NA 10.4 
## 138 2020  Parc… Abarco   4                NA        NA     2     NA     NA 13.1 
## 139 2020  Parc… Abarco   5                NA        NA     2     NA     NA 11.1 
## 140 2020  Parc… Abarco   6                NA        NA     2     NA     NA 11.5 
## 141 2020  Parc… Roble    1                NA        NA     2     NA     NA 11.4 
## 142 2020  Parc… Roble    2                NA        NA     2     NA     NA 11.0 
## 143 2020  Parc… Roble    3                NA        NA     2     NA     NA  9.23
## 144 2020  Parc… Roble    4                NA        NA     2     NA     NA 11.8 
## 145 2020  Parc… Roble    5                NA        NA     2     NA     NA 14.8 
## 146 2020  Parc… Roble    6                NA        NA     2     NA     NA 12.0 
## 147 2020  Parc… Termina… 1                NA        NA     2     NA     NA 18.2 
## 148 2020  Parc… Termina… 2                NA        NA     2     NA     NA 17.1 
## 149 2020  Parc… Termina… 3                NA        NA     2     NA     NA 17.6 
## 150 2020  Parc… Termina… 4                NA        NA     2     NA     NA 17.9 
## 151 2020  Parc… Termina… 5                NA        NA     2     NA     NA 18.0 
## 152 2020  Parc… Termina… 6                NA        NA     2     NA     NA 16.8 
## 153 2020  Sjose Abarco   1                NA        NA     2     NA     NA  4.43
## 154 2020  Sjose Abarco   2                NA        NA     2     NA     NA  4.75
## 155 2020  Sjose Abarco   3                NA        NA     2     NA     NA  4.94
## 156 2020  Sjose Abarco   4                NA        NA     2     NA     NA  6.04
## 157 2020  Sjose Abarco   5                NA        NA     2     NA     NA  5.67
## 158 2020  Sjose Abarco   6                NA        NA     2     NA     NA  6.4 
## 159 2020  Sjose Roble    1                NA        NA     2     NA     NA  9.75
## 160 2020  Sjose Roble    2                NA        NA     2     NA     NA 10.4 
## 161 2020  Sjose Roble    3                NA        NA     2     NA     NA  9.28
## 162 2020  Sjose Roble    4                NA        NA     2     NA     NA  9.24
## 163 2020  Sjose Roble    5                NA        NA     2     NA     NA  8.89
## 164 2020  Sjose Roble    6                NA        NA     2     NA     NA 10.2 
## 165 2020  Sjose Termina… 1                NA        NA     2     NA     NA 10.7 
## 166 2020  Sjose Termina… 2                NA        NA     2     NA     NA 11.1 
## 167 2020  Sjose Termina… 3                NA        NA     2     NA     NA  9.43
## 168 2020  Sjose Termina… 4                NA        NA     2     NA     NA  9.21
## 169 2020  Sjose Termina… 5                NA        NA     2     NA     NA 10.2 
## 170 2020  Sjose Termina… 6                NA        NA     2     NA     NA  9.72
## 171 2021  Aban… Abarco   1                NA        NA     3     NA     NA  9.12
## 172 2021  Aban… Abarco   2                NA        NA     3     NA     NA  9.13
## 173 2021  Aban… Abarco   3                NA        NA     3     NA     NA  7.7 
## 174 2021  Aban… Abarco   4                NA        NA     3     NA     NA  8.1 
## 175 2021  Aban… Abarco   5                NA        NA     3     NA     NA  7.96
## 176 2021  Aban… Abarco   6                NA        NA     3     NA     NA 10.1 
## 177 2021  Aban… Roble    1                NA        NA     3     NA     NA 12.9 
## 178 2021  Aban… Roble    2                NA        NA     3     NA     NA 10.1 
## 179 2021  Aban… Roble    3                NA        NA     3     NA     NA 11.0 
## 180 2021  Aban… Roble    4                NA        NA     3     NA     NA 17.8 
## 181 2021  Aban… Roble    5                NA        NA     3     NA     NA 14.2 
## 182 2021  Aban… Termina… 1                NA        NA     3     NA     NA 12.9 
## 183 2021  Aban… Termina… 2                NA        NA     3     NA     NA  9   
## 184 2021  Aban… Termina… 3                NA        NA     3     NA     NA 10.0 
## 185 2021  Aban… Termina… 4                NA        NA     3     NA     NA  8.65
## 186 2021  Aban… Termina… 5                NA        NA     3     NA     NA  7.58
## 187 2021  Aban… Termina… 6                NA        NA     3     NA     NA  9.3 
## 188 2021  Caba… Abarco   1                NA        NA     3     NA     NA  5.98
## 189 2021  Caba… Abarco   2                NA        NA     3     NA     NA  6.07
## 190 2021  Caba… Abarco   3                NA        NA     3     NA     NA  6.57
## 191 2021  Caba… Abarco   4                NA        NA     3     NA     NA  9.82
## 192 2021  Caba… Abarco   5                NA        NA     3     NA     NA 10.1 
## 193 2021  Caba… Abarco   6                NA        NA     3     NA     NA  9.55
## 194 2021  Caba… Roble    1                NA        NA     3     NA     NA 15.2 
## 195 2021  Caba… Roble    2                NA        NA     3     NA     NA 19.6 
## 196 2021  Caba… Roble    3                NA        NA     3     NA     NA 14.9 
## 197 2021  Caba… Roble    4                NA        NA     3     NA     NA 17.4 
## 198 2021  Caba… Roble    5                NA        NA     3     NA     NA 18.3 
## 199 2021  Caba… Roble    6                NA        NA     3     NA     NA 16.1 
## 200 2021  Caba… Termina… 1                NA        NA     3     NA     NA  8.53
## 201 2021  Caba… Termina… 2                NA        NA     3     NA     NA  7.3 
## 202 2021  Caba… Termina… 4                NA        NA     3     NA     NA  9.17
## 203 2021  Liba… Abarco   1                NA        NA     3     NA     NA  6.93
## 204 2021  Liba… Abarco   2                NA        NA     3     NA     NA  7.57
## 205 2021  Liba… Abarco   3                NA        NA     3     NA     NA  8   
## 206 2021  Liba… Abarco   4                NA        NA     3     NA     NA  6.87
## 207 2021  Liba… Abarco   5                NA        NA     3     NA     NA  7.8 
## 208 2021  Liba… Abarco   6                NA        NA     3     NA     NA  7   
## 209 2021  Liba… Roble    1                NA        NA     3     NA     NA  7.45
## 210 2021  Liba… Roble    2                NA        NA     3     NA     NA  5.05
## 211 2021  Liba… Roble    4                NA        NA     3     NA     NA  7.15
## 212 2021  Liba… Roble    5                NA        NA     3     NA     NA  6.4 
## 213 2021  Liba… Roble    6                NA        NA     3     NA     NA  5.45
## 214 2021  Liba… Termina… 1                NA        NA     3     NA     NA 14.4 
## 215 2021  Liba… Termina… 2                NA        NA     3     NA     NA 12.9 
## 216 2021  Liba… Termina… 3                NA        NA     3     NA     NA 11.9 
## 217 2021  Liba… Termina… 4                NA        NA     3     NA     NA 12.4 
## 218 2021  Liba… Termina… 5                NA        NA     3     NA     NA 15.7 
## 219 2021  Liba… Termina… 6                NA        NA     3     NA     NA 20.0 
## 220 2021  Parc… Abarco   1                NA        NA     3     NA     NA 15.4 
## 221 2021  Parc… Abarco   2                NA        NA     3     NA     NA 16.1 
## 222 2021  Parc… Abarco   3                NA        NA     3     NA     NA 16.0 
## 223 2021  Parc… Abarco   4                NA        NA     3     NA     NA 18.0 
## 224 2021  Parc… Abarco   5                NA        NA     3     NA     NA 15.5 
## 225 2021  Parc… Abarco   6                NA        NA     3     NA     NA 16.3 
## 226 2021  Parc… Roble    1                NA        NA     3     NA     NA 14.8 
## 227 2021  Parc… Roble    2                NA        NA     3     NA     NA 14.7 
## 228 2021  Parc… Roble    3                NA        NA     3     NA     NA 12.9 
## 229 2021  Parc… Roble    4                NA        NA     3     NA     NA 16.1 
## 230 2021  Parc… Roble    5                NA        NA     3     NA     NA 20.2 
## 231 2021  Parc… Roble    6                NA        NA     3     NA     NA 16.1 
## 232 2021  Parc… Termina… 1                NA        NA     3     NA     NA 22.6 
## 233 2021  Parc… Termina… 2                NA        NA     3     NA     NA 20.9 
## 234 2021  Parc… Termina… 3                NA        NA     3     NA     NA 23.5 
## 235 2021  Parc… Termina… 4                NA        NA     3     NA     NA 22.0 
## 236 2021  Parc… Termina… 5                NA        NA     3     NA     NA 22.9 
## 237 2021  Parc… Termina… 6                NA        NA     3     NA     NA 20.5 
## 238 2021  Sjose Abarco   1                NA        NA     3     NA     NA  6.94
## 239 2021  Sjose Abarco   2                NA        NA     3     NA     NA  7.09
## 240 2021  Sjose Abarco   3                NA        NA     3     NA     NA  7.42
## 241 2021  Sjose Abarco   4                NA        NA     3     NA     NA  9.21
## 242 2021  Sjose Abarco   5                NA        NA     3     NA     NA  9.12
## 243 2021  Sjose Abarco   6                NA        NA     3     NA     NA  7.83
## 244 2021  Sjose Roble    1                NA        NA     3     NA     NA 14.0 
## 245 2021  Sjose Roble    2                NA        NA     3     NA     NA 15.3 
## 246 2021  Sjose Roble    3                NA        NA     3     NA     NA 13.9 
## 247 2021  Sjose Roble    4                NA        NA     3     NA     NA 14.2 
## 248 2021  Sjose Roble    5                NA        NA     3     NA     NA 13.5 
## 249 2021  Sjose Roble    6                NA        NA     3     NA     NA 15.8 
## 250 2021  Sjose Termina… 1                NA        NA     3     NA     NA 14.4 
## 251 2021  Sjose Termina… 2                NA        NA     3     NA     NA 15.4 
## 252 2021  Sjose Termina… 3                NA        NA     3     NA     NA 11.5 
## 253 2021  Sjose Termina… 4                NA        NA     3     NA     NA 11.9 
## 254 2021  Sjose Termina… 5                NA        NA     3     NA     NA 14.1 
## 255 2021  Sjose Termina… 6                NA        NA     3     NA     NA 13.4 
## # … with 2 more variables: alt <dbl>, id <int>

write.csv(lista2, "~/Library/CloudStorage/GoogleDrive-icarounam@gmail.com/Mi unidad/Agrosavia/Colaboraciones/Laura/Forestal/basan.csv")
##Summary statistics
# Diameter
summ<-lista2 %>%
  group_by(lote, forestal, anof) %>%
  get_summary_stats(diam, type = "mean_sd")
summ %>% as_tibble() %>% print(n=Inf)

## # A tibble: 45 × 7
##    anof  lote    forestal   variable     n  mean    sd
##    <fct> <fct>   <fct>      <chr>    <dbl> <dbl> <dbl>
##  1 2019  Abandon Abarco     diam         6  3.15 0.478
##  2 2020  Abandon Abarco     diam         6  6.16 0.514
##  3 2021  Abandon Abarco     diam         6  8.68 0.907
##  4 2019  Abandon Roble      diam         5  5.90 1.2  
##  5 2020  Abandon Roble      diam         5  9.17 2.12 
##  6 2021  Abandon Roble      diam         5 13.2  3.05 
##  7 2019  Abandon Terminalia diam         6  2.74 1.66 
##  8 2020  Abandon Terminalia diam         6  6.52 1.96 
##  9 2021  Abandon Terminalia diam         6  9.58 1.83 
## 10 2019  Cabana  Abarco     diam         6  2.95 0.906
## 11 2020  Cabana  Abarco     diam         6  6.71 1.54 
## 12 2021  Cabana  Abarco     diam         6  8.01 2.00 
## 13 2019  Cabana  Roble      diam         6  6.46 1.60 
## 14 2020  Cabana  Roble      diam         6  8.74 1.94 
## 15 2021  Cabana  Roble      diam         6 16.9  1.87 
## 16 2019  Cabana  Terminalia diam         3  2.83 1.72 
## 17 2020  Cabana  Terminalia diam         3  5.78 0.234
## 18 2021  Cabana  Terminalia diam         3  8.33 0.949
## 19 2019  Libano  Abarco     diam         6  3.66 0.372
## 20 2020  Libano  Abarco     diam         6  5.97 0.38 
## 21 2021  Libano  Abarco     diam         6  7.36 0.49 
## 22 2019  Libano  Roble      diam         5  4.02 0.688
## 23 2020  Libano  Roble      diam         5  5.12 0.601
## 24 2021  Libano  Roble      diam         5  6.3  1.04 
## 25 2019  Libano  Terminalia diam         6  6.41 1.13 
## 26 2020  Libano  Terminalia diam         6 10.3  1.39 
## 27 2021  Libano  Terminalia diam         6 14.5  2.99 
## 28 2019  Parcela Abarco     diam         6  6.19 1.02 
## 29 2020  Parcela Abarco     diam         6 11.3  0.937
## 30 2021  Parcela Abarco     diam         6 16.2  0.95 
## 31 2019  Parcela Roble      diam         6  8.92 1.35 
## 32 2020  Parcela Roble      diam         6 11.7  1.83 
## 33 2021  Parcela Roble      diam         6 15.8  2.47 
## 34 2019  Parcela Terminalia diam         6 11.4  0.657
## 35 2020  Parcela Terminalia diam         6 17.6  0.544
## 36 2021  Parcela Terminalia diam         6 22.1  1.17 
## 37 2019  Sjose   Abarco     diam         6  2.77 0.401
## 38 2020  Sjose   Abarco     diam         6  5.37 0.78 
## 39 2021  Sjose   Abarco     diam         6  7.94 1.00 
## 40 2019  Sjose   Roble      diam         6  5.76 0.654
## 41 2020  Sjose   Roble      diam         6  9.62 0.594
## 42 2021  Sjose   Roble      diam         6 14.5  0.891
## 43 2019  Sjose   Terminalia diam         6  6.39 0.374
## 44 2020  Sjose   Terminalia diam         6 10.1  0.747
## 45 2021  Sjose   Terminalia diam         6 13.4  1.49

# Altura
summ<-lista2 %>%
  group_by(lote, forestal, anof) %>%
  get_summary_stats(alt, type = "mean_sd")
summ %>% as_tibble() %>% print(n=Inf)

## # A tibble: 45 × 7
##    anof  lote    forestal   variable     n  mean    sd
##    <fct> <fct>   <fct>      <chr>    <dbl> <dbl> <dbl>
##  1 2019  Abandon Abarco     alt          6  2.13 0.386
##  2 2020  Abandon Abarco     alt          6  3.77 0.537
##  3 2021  Abandon Abarco     alt          6  5.27 0.496
##  4 2019  Abandon Roble      alt          5  2.82 0.615
##  5 2020  Abandon Roble      alt          5  4.13 0.866
##  6 2021  Abandon Roble      alt          5  5.22 0.843
##  7 2019  Abandon Terminalia alt          6  1.32 0.388
##  8 2020  Abandon Terminalia alt          6  2.76 1.29 
##  9 2021  Abandon Terminalia alt          6  3.37 0.743
## 10 2019  Cabana  Abarco     alt          6  2.13 0.646
## 11 2020  Cabana  Abarco     alt          6  3.54 0.951
## 12 2021  Cabana  Abarco     alt          6  4.72 1.27 
## 13 2019  Cabana  Roble      alt          6  3.43 0.832
## 14 2020  Cabana  Roble      alt          6  5.79 0.865
## 15 2021  Cabana  Roble      alt          6  7.24 0.982
## 16 2019  Cabana  Terminalia alt          3  1.62 0.652
## 17 2020  Cabana  Terminalia alt          3  2.76 0.424
## 18 2021  Cabana  Terminalia alt          3  4.19 1.11 
## 19 2019  Libano  Abarco     alt          6  2.89 0.221
## 20 2020  Libano  Abarco     alt          6  4.29 0.372
## 21 2021  Libano  Abarco     alt          6  5.75 0.443
## 22 2019  Libano  Roble      alt          5  2.08 0.36 
## 23 2020  Libano  Roble      alt          5  2.82 0.233
## 24 2021  Libano  Roble      alt          5  3.60 1.09 
## 25 2019  Libano  Terminalia alt          6  2.73 0.711
## 26 2020  Libano  Terminalia alt          6  4.66 1.25 
## 27 2021  Libano  Terminalia alt          6  6.50 1.99 
## 28 2019  Parcela Abarco     alt          6  3.56 0.219
## 29 2020  Parcela Abarco     alt          6  6.44 0.386
## 30 2021  Parcela Abarco     alt          6  8.25 0.352
## 31 2019  Parcela Roble      alt          6  3.75 0.557
## 32 2020  Parcela Roble      alt          6  4.76 0.879
## 33 2021  Parcela Roble      alt          6  5.84 0.93 
## 34 2019  Parcela Terminalia alt          6  4.97 0.329
## 35 2020  Parcela Terminalia alt          6  8.86 0.733
## 36 2021  Parcela Terminalia alt          6 11.8  0.642
## 37 2019  Sjose   Abarco     alt          6  2.14 0.265
## 38 2020  Sjose   Abarco     alt          6  3.77 0.459
## 39 2021  Sjose   Abarco     alt          6  5.34 0.551
## 40 2019  Sjose   Roble      alt          6  3.05 0.505
## 41 2020  Sjose   Roble      alt          6  4.93 0.445
## 42 2021  Sjose   Roble      alt          6  6.31 0.568
## 43 2019  Sjose   Terminalia alt          6  3.02 0.059
## 44 2020  Sjose   Terminalia alt          6  4.32 0.292
## 45 2021  Sjose   Terminalia alt          6  6.28 0.611

##Visualization
# Diámetro
bxp <- ggboxplot(
  lista2, x = "lote", y = "diam",
  color = "anof", palette = "jco",
  facet.by =  "forestal"
)
bxp

#Altura
bxp <- ggboxplot(
  lista2, x = "lote", y = "alt",
  color = "anof", palette = "jco",
  facet.by =  "forestal"
)
bxp

##Check assumptions
##Outliers
# Diámetro
lista2 %>%
  group_by(lote, forestal, anof) %>%
  identify_outliers(diam) %>% print(n=Inf)

## # A tibble: 11 × 14
##    anof  lote   forestal planta departamento municipio   ano semana bloque  diam
##    <fct> <fct>  <fct>    <fct>         <dbl>     <dbl> <dbl>  <dbl>  <dbl> <dbl>
##  1 2019  Aband… Abarco   6                NA        NA     1     NA     NA  3.94
##  2 2021  Aband… Termina… 1                NA        NA     3     NA     NA 12.9 
##  3 2021  Libano Termina… 6                NA        NA     3     NA     NA 20.0 
##  4 2019  Parce… Abarco   4                NA        NA     1     NA     NA  8.04
##  5 2020  Parce… Abarco   4                NA        NA     2     NA     NA 13.1 
##  6 2021  Parce… Abarco   4                NA        NA     3     NA     NA 18.0 
##  7 2020  Parce… Roble    3                NA        NA     2     NA     NA  9.23
##  8 2020  Parce… Roble    5                NA        NA     2     NA     NA 14.8 
##  9 2021  Parce… Roble    5                NA        NA     3     NA     NA 20.2 
## 10 2019  Parce… Termina… 3                NA        NA     1     NA     NA 12.6 
## 11 2019  Sjose  Abarco   4                NA        NA     1     NA     NA  3.56
## # … with 4 more variables: alt <dbl>, id <int>, is.outlier <lgl>,
## #   is.extreme <lgl>

# Altura
lista2 %>%
  group_by(lote, forestal, anof) %>%
  identify_outliers(alt) %>% print(n=Inf)

## # A tibble: 12 × 14
##    anof  lote   forestal planta departamento municipio   ano semana bloque  diam
##    <fct> <fct>  <fct>    <fct>         <dbl>     <dbl> <dbl>  <dbl>  <dbl> <dbl>
##  1 2021  Aband… Abarco   6                NA        NA     3     NA     NA 10.1 
##  2 2019  Aband… Termina… 1                NA        NA     1     NA     NA  5.28
##  3 2020  Libano Roble    6                NA        NA     2     NA     NA  4.65
##  4 2021  Libano Roble    1                NA        NA     3     NA     NA  7.45
##  5 2021  Libano Roble    6                NA        NA     3     NA     NA  5.45
##  6 2019  Libano Termina… 6                NA        NA     1     NA     NA  8.28
##  7 2021  Libano Termina… 6                NA        NA     3     NA     NA 20.0 
##  8 2019  Parce… Abarco   4                NA        NA     1     NA     NA  8.04
##  9 2021  Parce… Abarco   4                NA        NA     3     NA     NA 18.0 
## 10 2019  Parce… Termina… 3                NA        NA     1     NA     NA 12.6 
## 11 2019  Sjose  Abarco   4                NA        NA     1     NA     NA  3.56
## 12 2019  Sjose  Roble    1                NA        NA     1     NA     NA  6.54
## # … with 4 more variables: alt <dbl>, id <int>, is.outlier <lgl>,
## #   is.extreme <lgl>

##Normality assumption
##Compute Shapiro-Wilk test for each combinations of factor levels:
#Diámetro
norm<-lista2 %>%
  group_by(lote, forestal, anof) %>%
  shapiro_test(diam)
norm %>% as_tibble() %>% print(n=Inf)

## # A tibble: 45 × 6
##    anof  lote    forestal   variable statistic      p
##    <fct> <fct>   <fct>      <chr>        <dbl>  <dbl>
##  1 2019  Abandon Abarco     diam         0.954 0.775 
##  2 2020  Abandon Abarco     diam         0.941 0.670 
##  3 2021  Abandon Abarco     diam         0.910 0.438 
##  4 2019  Abandon Roble      diam         0.912 0.477 
##  5 2020  Abandon Roble      diam         0.882 0.319 
##  6 2021  Abandon Roble      diam         0.943 0.688 
##  7 2019  Abandon Terminalia diam         0.866 0.212 
##  8 2020  Abandon Terminalia diam         0.946 0.709 
##  9 2021  Abandon Terminalia diam         0.888 0.310 
## 10 2019  Cabana  Abarco     diam         0.854 0.170 
## 11 2020  Cabana  Abarco     diam         0.849 0.155 
## 12 2021  Cabana  Abarco     diam         0.785 0.0427
## 13 2019  Cabana  Roble      diam         0.960 0.821 
## 14 2020  Cabana  Roble      diam         0.936 0.627 
## 15 2021  Cabana  Roble      diam         0.946 0.704 
## 16 2019  Cabana  Terminalia diam         0.919 0.448 
## 17 2020  Cabana  Terminalia diam         0.993 0.843 
## 18 2021  Cabana  Terminalia diam         0.967 0.650 
## 19 2019  Libano  Abarco     diam         0.924 0.531 
## 20 2020  Libano  Abarco     diam         0.835 0.117 
## 21 2021  Libano  Abarco     diam         0.870 0.225 
## 22 2019  Libano  Roble      diam         0.913 0.483 
## 23 2020  Libano  Roble      diam         0.943 0.685 
## 24 2021  Libano  Roble      diam         0.927 0.574 
## 25 2019  Libano  Terminalia diam         0.938 0.644 
## 26 2020  Libano  Terminalia diam         0.891 0.324 
## 27 2021  Libano  Terminalia diam         0.863 0.200 
## 28 2019  Parcela Abarco     diam         0.888 0.308 
## 29 2020  Parcela Abarco     diam         0.838 0.127 
## 30 2021  Parcela Abarco     diam         0.811 0.0735
## 31 2019  Parcela Roble      diam         0.971 0.900 
## 32 2020  Parcela Roble      diam         0.931 0.591 
## 33 2021  Parcela Roble      diam         0.900 0.376 
## 34 2019  Parcela Terminalia diam         0.866 0.212 
## 35 2020  Parcela Terminalia diam         0.908 0.426 
## 36 2021  Parcela Terminalia diam         0.939 0.651 
## 37 2019  Sjose   Abarco     diam         0.730 0.0125
## 38 2020  Sjose   Abarco     diam         0.939 0.652 
## 39 2021  Sjose   Abarco     diam         0.851 0.160 
## 40 2019  Sjose   Roble      diam         0.916 0.477 
## 41 2020  Sjose   Roble      diam         0.948 0.722 
## 42 2021  Sjose   Roble      diam         0.876 0.250 
## 43 2019  Sjose   Terminalia diam         0.965 0.855 
## 44 2020  Sjose   Terminalia diam         0.949 0.730 
## 45 2021  Sjose   Terminalia diam         0.947 0.713

#Altura
norm<-lista2 %>%
  group_by(lote, forestal, anof) %>%
  shapiro_test(alt)
norm %>% as_tibble() %>% print(n=Inf)

## # A tibble: 45 × 6
##    anof  lote    forestal   variable statistic      p
##    <fct> <fct>   <fct>      <chr>        <dbl>  <dbl>
##  1 2019  Abandon Abarco     alt          0.947 0.715 
##  2 2020  Abandon Abarco     alt          0.914 0.466 
##  3 2021  Abandon Abarco     alt          0.826 0.0989
##  4 2019  Abandon Roble      alt          0.976 0.912 
##  5 2020  Abandon Roble      alt          0.912 0.481 
##  6 2021  Abandon Roble      alt          0.847 0.184 
##  7 2019  Abandon Terminalia alt          0.911 0.442 
##  8 2020  Abandon Terminalia alt          0.858 0.183 
##  9 2021  Abandon Terminalia alt          0.898 0.363 
## 10 2019  Cabana  Abarco     alt          0.946 0.707 
## 11 2020  Cabana  Abarco     alt          0.946 0.706 
## 12 2021  Cabana  Abarco     alt          0.900 0.374 
## 13 2019  Cabana  Roble      alt          0.890 0.319 
## 14 2020  Cabana  Roble      alt          0.917 0.485 
## 15 2021  Cabana  Roble      alt          0.983 0.965 
## 16 2019  Cabana  Terminalia alt          0.75  0     
## 17 2020  Cabana  Terminalia alt          0.952 0.580 
## 18 2021  Cabana  Terminalia alt          0.999 0.931 
## 19 2019  Libano  Abarco     alt          0.914 0.462 
## 20 2020  Libano  Abarco     alt          0.918 0.489 
## 21 2021  Libano  Abarco     alt          0.977 0.937 
## 22 2019  Libano  Roble      alt          0.967 0.853 
## 23 2020  Libano  Roble      alt          0.948 0.725 
## 24 2021  Libano  Roble      alt          0.856 0.213 
## 25 2019  Libano  Terminalia alt          0.815 0.0798
## 26 2020  Libano  Terminalia alt          0.819 0.0860
## 27 2021  Libano  Terminalia alt          0.857 0.178 
## 28 2019  Parcela Abarco     alt          0.827 0.102 
## 29 2020  Parcela Abarco     alt          0.906 0.410 
## 30 2021  Parcela Abarco     alt          0.918 0.488 
## 31 2019  Parcela Roble      alt          0.916 0.474 
## 32 2020  Parcela Roble      alt          0.976 0.930 
## 33 2021  Parcela Roble      alt          0.936 0.628 
## 34 2019  Parcela Terminalia alt          0.810 0.0723
## 35 2020  Parcela Terminalia alt          0.906 0.412 
## 36 2021  Parcela Terminalia alt          0.979 0.948 
## 37 2019  Sjose   Abarco     alt          0.799 0.0578
## 38 2020  Sjose   Abarco     alt          0.928 0.563 
## 39 2021  Sjose   Abarco     alt          0.861 0.193 
## 40 2019  Sjose   Roble      alt          0.941 0.670 
## 41 2020  Sjose   Roble      alt          0.846 0.147 
## 42 2021  Sjose   Roble      alt          0.877 0.255 
## 43 2019  Sjose   Terminalia alt          0.932 0.594 
## 44 2020  Sjose   Terminalia alt          0.958 0.804 
## 45 2021  Sjose   Terminalia alt          0.971 0.902

##Create QQ plot for each cell of design:
#Diámetro
ggqqplot(lista2, "diam", ggtheme = theme_bw()) +
  facet_grid(anof~ lote*forestal, labeller = "label_both")

#Altura
ggqqplot(lista2, "alt", ggtheme = theme_bw()) +
  facet_grid(anof~ lote*forestal, labeller = "label_both")

##Homogneity of variance assumption
##Compute the Levene’s test at each level of the within-subjects factor, here time variable:
#Diámetro
lev<-lista2 %>%
  group_by(anof) %>%
  levene_test(diam ~ lote*forestal)
lev %>% as_tibble() %>% print(n=Inf)

## # A tibble: 3 × 5
##   anof    df1   df2 statistic       p
##   <fct> <int> <int>     <dbl>   <dbl>
## 1 2019     14    70      2.16 0.0184 
## 2 2020     14    70      2.57 0.00494
## 3 2021     14    70      1.67 0.0830

#Altura
lev<-lista2 %>%
  group_by(anof) %>%
  levene_test(alt ~ lote*forestal)
lev %>% as_tibble() %>% print(n=Inf)

## # A tibble: 3 × 5
##   anof    df1   df2 statistic      p
##   <fct> <int> <int>     <dbl>  <dbl>
## 1 2019     14    70      1.41 0.170 
## 2 2020     14    70      1.94 0.0368
## 3 2021     14    70      1.30 0.231

##Computation
lista3<-as.data.frame(lista2)
#Diámetro
res.aov.diam <- anova_test(
  data = lista3, dv = diam, wid = id,
  within = anof, between = c(lote, forestal)
)
get_anova_table(res.aov.diam)

## ANOVA Table (type III tests)
## 
##               Effect   DFn    DFd        F        p p<.05   ges
## 1               lote  4.00  70.00   86.397 2.53e-26     * 0.796
## 2           forestal  2.00  70.00   52.295 1.28e-14     * 0.542
## 3               anof  1.44 101.01 1675.943 7.54e-72     * 0.832
## 4      lote:forestal  8.00  70.00   19.943 2.61e-15     * 0.644
## 5          lote:anof  5.77 101.01   19.571 1.01e-14     * 0.188
## 6      forestal:anof  2.89 101.01   17.772 3.99e-09     * 0.095
## 7 lote:forestal:anof 11.54 101.01   15.520 1.12e-17     * 0.269

#Altura
res.aov.alt <- anova_test(
  data = lista3, dv = alt, wid = id,
  within = anof, between = c(lote, forestal)
)
get_anova_table(res.aov.alt)

## ANOVA Table (type III tests)
## 
##               Effect   DFn    DFd        F        p p<.05   ges
## 1               lote  4.00  70.00   55.635 5.05e-21     * 0.721
## 2           forestal  2.00  70.00    1.836 1.67e-01       0.041
## 3               anof  1.76 123.33 1347.158 1.59e-81     * 0.784
## 4      lote:forestal  8.00  70.00   22.145 2.14e-16     * 0.672
## 5          lote:anof  7.05 123.33   18.483 9.95e-17     * 0.166
## 6      forestal:anof  3.52 123.33   12.883 3.89e-08     * 0.065
## 7 lote:forestal:anof 14.09 123.33   15.145 9.95e-21     * 0.246

## Contrastes para lote * forestal
attach(lista3)
lista3$id<-as.factor(lista3$id)
#Diámetro
het.diam <- gls(diam ~ anof*forestal*lote, 
            data=lista3, corr=corAR1(, form= ~ 1 | id), 
            weights=varIdent(form = ~ 1 | anof), 
            na.action=na.exclude)
anova(het.diam)

## Denom. DF: 210 
##                    numDF   F-value p-value
## (Intercept)            1 2416.6372  <.0001
## anof                   2 1316.9912  <.0001
## forestal               2   54.3852  <.0001
## lote                   4   65.1838  <.0001
## anof:forestal          4   28.2908  <.0001
## anof:lote              8   15.1951  <.0001
## forestal:lote          8   12.8156  <.0001
## anof:forestal:lote    16   13.3001  <.0001

# Comparando diámetros de forestales para cada lote
contrast.diam <- emmeans(het.diam, ~forestal|lote)

## NOTE: Results may be misleading due to involvement in interactions

plot(contrast.diam, comparisons = TRUE, xlab ="Diámetro")

cld_lote_forestal.diam <-multcomp::cld(contrast.diam, alpha = 0.05, Letters = LETTERS, reversed=T)
cld_lote_forestal.diam

## lote = Abandon:
##  forestal   emmean    SE   df lower.CL upper.CL .group
##  Roble        9.42 0.540 94.5     8.35    10.49  A    
##  Terminalia   6.28 0.493 83.6     5.30     7.26   B   
##  Abarco       5.99 0.493 83.6     5.01     6.97   B   
## 
## lote = Cabana:
##  forestal   emmean    SE   df lower.CL upper.CL .group
##  Roble       10.71 0.493 84.2     9.73    11.69  A    
##  Abarco       5.89 0.493 84.2     4.91     6.87   B   
##  Terminalia   5.65 0.697 82.6     4.26     7.04   B   
## 
## lote = Libano:
##  forestal   emmean    SE   df lower.CL upper.CL .group
##  Terminalia  10.42 0.493 83.6     9.44    11.40  A    
##  Abarco       5.66 0.493 83.2     4.68     6.65   B   
##  Roble        5.14 0.540 94.1     4.07     6.22   B   
## 
## lote = Parcela:
##  forestal   emmean    SE   df lower.CL upper.CL .group
##  Terminalia  17.03 0.493 84.6    16.05    18.01  A    
##  Roble       12.15 0.493 84.6    11.17    13.13   B   
##  Abarco      11.24 0.493 84.4    10.26    12.22   B   
## 
## lote = Sjose:
##  forestal   emmean    SE   df lower.CL upper.CL .group
##  Terminalia   9.96 0.493 84.2     8.98    10.95  A    
##  Roble        9.95 0.493 84.5     8.97    10.93  A    
##  Abarco       5.36 0.493 84.6     4.38     6.34   B   
## 
## Results are averaged over the levels of: anof 
## Degrees-of-freedom method: satterthwaite 
## Confidence level used: 0.95 
## P value adjustment: tukey method for comparing a family of 3 estimates 
## significance level used: alpha = 0.05 
## NOTE: Compact letter displays can be misleading
##       because they show NON-findings rather than findings.
##       Consider using 'pairs()', 'pwpp()', or 'pwpm()' instead.

# Comparando diámetros de cada forestal a lo largo de los lotes
contrast.diam <- emmeans(het.diam, ~lote|forestal)

## NOTE: Results may be misleading due to involvement in interactions

plot(contrast.diam, comparisons = TRUE, xlab ="Diámetro")

cld_lote_forestal2.diam <-multcomp::cld(contrast.diam, alpha = 0.05, Letters = LETTERS, reversed=T)
cld_lote_forestal2.diam

## forestal = Abarco:
##  lote    emmean    SE   df lower.CL upper.CL .group
##  Parcela  11.24 0.493 84.4    10.26    12.22  A    
##  Abandon   5.99 0.493 83.6     5.01     6.97   B   
##  Cabana    5.89 0.493 84.2     4.91     6.87   B   
##  Libano    5.66 0.493 83.2     4.68     6.65   B   
##  Sjose     5.36 0.493 84.6     4.38     6.34   B   
## 
## forestal = Roble:
##  lote    emmean    SE   df lower.CL upper.CL .group
##  Parcela  12.15 0.493 84.6    11.17    13.13  A    
##  Cabana   10.71 0.493 84.2     9.73    11.69  AB   
##  Sjose     9.95 0.493 84.5     8.97    10.93   B   
##  Abandon   9.42 0.540 94.5     8.35    10.49   B   
##  Libano    5.14 0.540 94.1     4.07     6.22    C  
## 
## forestal = Terminalia:
##  lote    emmean    SE   df lower.CL upper.CL .group
##  Parcela  17.03 0.493 84.6    16.05    18.01  A    
##  Libano   10.42 0.493 83.6     9.44    11.40   B   
##  Sjose     9.96 0.493 84.2     8.98    10.95   B   
##  Abandon   6.28 0.493 83.6     5.30     7.26    C  
##  Cabana    5.65 0.697 82.6     4.26     7.04    C  
## 
## Results are averaged over the levels of: anof 
## Degrees-of-freedom method: satterthwaite 
## Confidence level used: 0.95 
## P value adjustment: tukey method for comparing a family of 5 estimates 
## significance level used: alpha = 0.05 
## NOTE: Compact letter displays can be misleading
##       because they show NON-findings rather than findings.
##       Consider using 'pairs()', 'pwpp()', or 'pwpm()' instead.

#Altura
het.alt <- gls(alt ~ anof*forestal*lote, 
                data=lista3, corr=corAR1(, form= ~ 1 | id), 
                weights=varIdent(form = ~ 1 | anof), 
                na.action=na.exclude)
anova(het.alt)

## Denom. DF: 210 
##                    numDF   F-value p-value
## (Intercept)            1 2194.6587  <.0001
## anof                   2  959.5158  <.0001
## forestal               2    9.3945   1e-04
## lote                   4   33.5380  <.0001
## anof:forestal          4    9.5238  <.0001
## anof:lote              8   13.3822  <.0001
## forestal:lote          8    9.2985  <.0001
## anof:forestal:lote    16   11.5432  <.0001

# Comparando alturas de forestales para cada lote
contrast.alt <- emmeans(het.alt, ~forestal|lote)

## NOTE: Results may be misleading due to involvement in interactions

plot(contrast.alt, comparisons = TRUE, xlab ="Diámetro")

cld_lote_forestal.alt <-multcomp::cld(contrast.alt, alpha = 0.05, Letters = LETTERS, reversed=T)
cld_lote_forestal.alt

## lote = Abandon:
##  forestal   emmean    SE   df lower.CL upper.CL .group
##  Roble        4.06 0.294 87.3     3.47     4.64  A    
##  Abarco       3.72 0.269 80.5     3.19     4.26  A    
##  Terminalia   2.48 0.269 81.0     1.95     3.02   B   
## 
## lote = Cabana:
##  forestal   emmean    SE   df lower.CL upper.CL .group
##  Roble        5.48 0.269 81.0     4.95     6.02  A    
##  Abarco       3.47 0.269 81.3     2.93     4.00   B   
##  Terminalia   2.86 0.380 81.4     2.10     3.61   B   
## 
## lote = Libano:
##  forestal   emmean    SE   df lower.CL upper.CL .group
##  Terminalia   4.63 0.269 81.1     4.09     5.16  A    
##  Abarco       4.31 0.269 80.0     3.77     4.84  A    
##  Roble        2.83 0.294 87.0     2.25     3.42   B   
## 
## lote = Parcela:
##  forestal   emmean    SE   df lower.CL upper.CL .group
##  Terminalia   8.56 0.269 81.5     8.02     9.09  A    
##  Abarco       6.09 0.269 81.4     5.55     6.62   B   
##  Roble        4.78 0.269 81.5     4.25     5.32    C  
## 
## lote = Sjose:
##  forestal   emmean    SE   df lower.CL upper.CL .group
##  Roble        4.76 0.269 81.4     4.23     5.30  A    
##  Terminalia   4.54 0.269 81.2     4.00     5.07  AB   
##  Abarco       3.75 0.269 81.5     3.21     4.28   B   
## 
## Results are averaged over the levels of: anof 
## Degrees-of-freedom method: satterthwaite 
## Confidence level used: 0.95 
## P value adjustment: tukey method for comparing a family of 3 estimates 
## significance level used: alpha = 0.05 
## NOTE: Compact letter displays can be misleading
##       because they show NON-findings rather than findings.
##       Consider using 'pairs()', 'pwpp()', or 'pwpm()' instead.

# Comparando alturas de cada forestal a lo largo de los lotes
contrast.alt <- emmeans(het.alt, ~lote|forestal)

## NOTE: Results may be misleading due to involvement in interactions

plot(contrast.alt, comparisons = TRUE, xlab ="Altura")

cld_lote_forestal.alt2 <-multcomp::cld(contrast.alt, alpha = 0.05, Letters = LETTERS, reversed=T)
cld_lote_forestal.alt2

## forestal = Abarco:
##  lote    emmean    SE   df lower.CL upper.CL .group
##  Parcela   6.09 0.269 81.4     5.55     6.62  A    
##  Libano    4.31 0.269 80.0     3.77     4.84   B   
##  Sjose     3.75 0.269 81.5     3.21     4.28   B   
##  Abandon   3.72 0.269 80.5     3.19     4.26   B   
##  Cabana    3.47 0.269 81.3     2.93     4.00   B   
## 
## forestal = Roble:
##  lote    emmean    SE   df lower.CL upper.CL .group
##  Cabana    5.48 0.269 81.0     4.95     6.02  A    
##  Parcela   4.78 0.269 81.5     4.25     5.32  AB   
##  Sjose     4.76 0.269 81.4     4.23     5.30  AB   
##  Abandon   4.06 0.294 87.3     3.47     4.64   B   
##  Libano    2.83 0.294 87.0     2.25     3.42    C  
## 
## forestal = Terminalia:
##  lote    emmean    SE   df lower.CL upper.CL .group
##  Parcela   8.56 0.269 81.5     8.02     9.09  A    
##  Libano    4.63 0.269 81.1     4.09     5.16   B   
##  Sjose     4.54 0.269 81.2     4.00     5.07   B   
##  Cabana    2.86 0.380 81.4     2.10     3.61    C  
##  Abandon   2.48 0.269 81.0     1.95     3.02    C  
## 
## Results are averaged over the levels of: anof 
## Degrees-of-freedom method: satterthwaite 
## Confidence level used: 0.95 
## P value adjustment: tukey method for comparing a family of 5 estimates 
## significance level used: alpha = 0.05 
## NOTE: Compact letter displays can be misleading
##       because they show NON-findings rather than findings.
##       Consider using 'pairs()', 'pwpp()', or 'pwpm()' instead.