library(psych)
library(readxl)
library(ggplot2)
## 
## Attaching package: 'ggplot2'
## The following objects are masked from 'package:psych':
## 
##     %+%, alpha
library(lattice)
library(tidyverse)
## ── Attaching core tidyverse packages ──────────────────────── tidyverse 2.0.0 ──
## ✔ dplyr     1.1.3     ✔ readr     2.1.4
## ✔ forcats   1.0.0     ✔ stringr   1.5.0
## ✔ lubridate 1.9.3     ✔ tibble    3.2.1
## ✔ purrr     1.0.2     ✔ tidyr     1.3.0
## ── Conflicts ────────────────────────────────────────── tidyverse_conflicts() ──
## ✖ ggplot2::%+%()   masks psych::%+%()
## ✖ ggplot2::alpha() masks psych::alpha()
## ✖ dplyr::filter()  masks stats::filter()
## ✖ dplyr::lag()     masks stats::lag()
## ℹ Use the conflicted package (<http://conflicted.r-lib.org/>) to force all conflicts to become errors
library(dplyr)
library(abdiv)
## 
## Attaching package: 'abdiv'
## 
## The following objects are masked from 'package:psych':
## 
##     manhattan, minkowski, shannon
library(vegan)
## Loading required package: permute
## This is vegan 2.6-4
library(codyn)
library(divo)
## Loading required package: cluster
library(fossil)
## Loading required package: sp
## The legacy packages maptools, rgdal, and rgeos, underpinning the sp package,
## which was just loaded, were retired in October 2023.
## Please refer to R-spatial evolution reports for details, especially
## https://r-spatial.org/r/2023/05/15/evolution4.html.
## It may be desirable to make the sf package available;
## package maintainers should consider adding sf to Suggests:.
## Loading required package: maps
## 
## Attaching package: 'maps'
## 
## The following object is masked from 'package:cluster':
## 
##     votes.repub
## 
## The following object is masked from 'package:purrr':
## 
##     map
## 
## Loading required package: shapefiles
## Loading required package: foreign
## 
## Attaching package: 'shapefiles'
## 
## The following objects are masked from 'package:foreign':
## 
##     read.dbf, write.dbf
## 
## The following object is masked from 'package:divo':
## 
##     dp
## 
## 
## Attaching package: 'fossil'
## 
## The following objects are masked from 'package:abdiv':
## 
##     euclidean, jaccard, manhattan, simpson, sorenson
## 
## The following object is masked from 'package:psych':
## 
##     manhattan

CEBADA

Datos y gráficas del seguimiento de parametros de crecimiento del cultivo (muestreo 3)

DMTOTAL3 <- read_excel("cultivo cebada medio m2.xlsx", 
    col_types = c("numeric", "numeric", "numeric", 
        "numeric", "numeric", "numeric", 
        "numeric"))

attach(DMTOTAL3)

# Altura de cebada:
psych::describe(ALTURAPROM0.5M2)
##    vars n  mean   sd median trimmed mad min  max range skew kurtosis   se
## X1    1 8 61.54 6.35   60.7   61.54 6.3  53 71.5  18.5 0.24    -1.55 2.24
describe.by(ALTURAPROM0.5M2, group = PUNTOMUES)
## Warning: describe.by is deprecated.  Please use the describeBy function
## 
##  Descriptive statistics by group 
## group: 1
##    vars n  mean sd median trimmed mad   min   max range skew kurtosis se
## X1    1 1 61.41 NA  61.41   61.41   0 61.41 61.41     0   NA       NA NA
## ------------------------------------------------------------ 
## group: 2
##    vars n mean sd median trimmed mad  min  max range skew kurtosis se
## X1    1 1 64.7 NA   64.7    64.7   0 64.7 64.7     0   NA       NA NA
## ------------------------------------------------------------ 
## group: 3
##    vars n mean sd median trimmed mad min max range skew kurtosis se
## X1    1 1   53 NA     53      53   0  53  53     0   NA       NA NA
## ------------------------------------------------------------ 
## group: 4
##    vars n mean sd median trimmed mad min max range skew kurtosis se
## X1    1 1   60 NA     60      60   0  60  60     0   NA       NA NA
## ------------------------------------------------------------ 
## group: 5
##    vars n mean sd median trimmed mad  min  max range skew kurtosis se
## X1    1 1 71.5 NA   71.5    71.5   0 71.5 71.5     0   NA       NA NA
## ------------------------------------------------------------ 
## group: 6
##    vars n mean sd median trimmed mad  min  max range skew kurtosis se
## X1    1 1 56.2 NA   56.2    56.2   0 56.2 56.2     0   NA       NA NA
## ------------------------------------------------------------ 
## group: 7
##    vars n mean sd median trimmed mad min max range skew kurtosis se
## X1    1 1   57 NA     57      57   0  57  57     0   NA       NA NA
## ------------------------------------------------------------ 
## group: 8
##    vars n mean sd median trimmed mad  min  max range skew kurtosis se
## X1    1 1 68.5 NA   68.5    68.5   0 68.5 68.5     0   NA       NA NA
# Macollas por planta:
psych::describe(MACOLLAS0.5M2)
##    vars n  mean   sd median trimmed  mad min max range skew kurtosis   se
## X1    1 8 14.75 3.65   13.5   14.75 1.48  11  23    12 1.27     0.41 1.29
describe.by(MACOLLAS0.5M2, group = PUNTOMUES)
## Warning: describe.by is deprecated.  Please use the describeBy function
## 
##  Descriptive statistics by group 
## group: 1
##    vars n mean sd median trimmed mad min max range skew kurtosis se
## X1    1 1   13 NA     13      13   0  13  13     0   NA       NA NA
## ------------------------------------------------------------ 
## group: 2
##    vars n mean sd median trimmed mad min max range skew kurtosis se
## X1    1 1   15 NA     15      15   0  15  15     0   NA       NA NA
## ------------------------------------------------------------ 
## group: 3
##    vars n mean sd median trimmed mad min max range skew kurtosis se
## X1    1 1   14 NA     14      14   0  14  14     0   NA       NA NA
## ------------------------------------------------------------ 
## group: 4
##    vars n mean sd median trimmed mad min max range skew kurtosis se
## X1    1 1   23 NA     23      23   0  23  23     0   NA       NA NA
## ------------------------------------------------------------ 
## group: 5
##    vars n mean sd median trimmed mad min max range skew kurtosis se
## X1    1 1   13 NA     13      13   0  13  13     0   NA       NA NA
## ------------------------------------------------------------ 
## group: 6
##    vars n mean sd median trimmed mad min max range skew kurtosis se
## X1    1 1   16 NA     16      16   0  16  16     0   NA       NA NA
## ------------------------------------------------------------ 
## group: 7
##    vars n mean sd median trimmed mad min max range skew kurtosis se
## X1    1 1   13 NA     13      13   0  13  13     0   NA       NA NA
## ------------------------------------------------------------ 
## group: 8
##    vars n mean sd median trimmed mad min max range skew kurtosis se
## X1    1 1   11 NA     11      11   0  11  11     0   NA       NA NA
# densidad de cebada:
psych::describe(DENSIDAD0.5M2)
##    vars n  mean    sd median trimmed   mad min max range skew kurtosis   se
## X1    1 8 70.12 17.13     70   70.12 20.02  46  94    48 0.03    -1.64 6.05
describe.by(DENSIDAD0.5M2, group = PUNTOMUES)
## Warning: describe.by is deprecated.  Please use the describeBy function
## 
##  Descriptive statistics by group 
## group: 1
##    vars n mean sd median trimmed mad min max range skew kurtosis se
## X1    1 1   66 NA     66      66   0  66  66     0   NA       NA NA
## ------------------------------------------------------------ 
## group: 2
##    vars n mean sd median trimmed mad min max range skew kurtosis se
## X1    1 1   52 NA     52      52   0  52  52     0   NA       NA NA
## ------------------------------------------------------------ 
## group: 3
##    vars n mean sd median trimmed mad min max range skew kurtosis se
## X1    1 1   61 NA     61      61   0  61  61     0   NA       NA NA
## ------------------------------------------------------------ 
## group: 4
##    vars n mean sd median trimmed mad min max range skew kurtosis se
## X1    1 1   78 NA     78      78   0  78  78     0   NA       NA NA
## ------------------------------------------------------------ 
## group: 5
##    vars n mean sd median trimmed mad min max range skew kurtosis se
## X1    1 1   94 NA     94      94   0  94  94     0   NA       NA NA
## ------------------------------------------------------------ 
## group: 6
##    vars n mean sd median trimmed mad min max range skew kurtosis se
## X1    1 1   46 NA     46      46   0  46  46     0   NA       NA NA
## ------------------------------------------------------------ 
## group: 7
##    vars n mean sd median trimmed mad min max range skew kurtosis se
## X1    1 1   74 NA     74      74   0  74  74     0   NA       NA NA
## ------------------------------------------------------------ 
## group: 8
##    vars n mean sd median trimmed mad min max range skew kurtosis se
## X1    1 1   90 NA     90      90   0  90  90     0   NA       NA NA
## Gráficos de altura de cebada por metro cuadrado 

par(mfrow=c(1,1))
ggplot(DMTOTAL3, aes(x = PUNTOMUES, y = ALTURAPROM0.5M2, fill = MUESTREO)) + 
  geom_bar(stat = "identity", position = "dodge") +
   guides(fill = guide_legend(title = "Muestreo")) + 
    ggtitle("Altura de cebada por medio metro cuadrado") + 
     theme(panel.border = element_rect(fill = "transparent", color = 4,size = 2))  + 
       theme(plot.background = element_rect(fill = "gray86")) +
     theme(plot.background = element_rect(color = "black", size = 2)) +
  xlab("Grupos") + ylab("Altura de plantas (cm) por medio metro cuadrado") + 
  theme(axis.title.x = element_text(size = 11, color = "black")) + theme(axis.title.y = element_text(size = 13, color = "black")) + theme(axis.text.x = element_text(angle = 0, vjust = 1, hjust = 1))+
  scale_x_continuous(n.breaks = 8)
## Warning: The `size` argument of `element_rect()` is deprecated as of ggplot2 3.4.0.
## ℹ Please use the `linewidth` argument instead.
## This warning is displayed once every 8 hours.
## Call `lifecycle::last_lifecycle_warnings()` to see where this warning was
## generated.

## Gráfico de densidad del cultivo por metro cuadrado  

ggplot(DMTOTAL3, aes(x = PUNTOMUES, y = DENSIDAD0.5M2, fill = MUESTREO)) + 
  geom_bar(stat = "identity", position = "dodge") +
   guides(fill = guide_legend(title = "Muestreo")) + 
    ggtitle("Densisdad de cebada por medio metro cuadrado") + 
     theme(panel.border = element_rect(fill = "transparent", color = 4,size = 2))  + 
       theme(plot.background = element_rect(fill = "gray86")) +
     theme(plot.background = element_rect(color = "black", size = 2)) +
  xlab("Grupos") + ylab("Número de plantas por medio metro cuadrado") + 
  theme(axis.title.x = element_text(size = 11, color = "black")) + theme(axis.title.y = element_text(size = 13, color = "black")) + theme(axis.text.x = element_text(angle = 0, vjust = 1, hjust = 1))  + scale_x_continuous(n.breaks = 8)

## Gráfico de número de macollas por metro cuadrado 

ggplot(DMTOTAL3, aes(x = PUNTOMUES, y = MACOLLAS0.5M2, fill = MUESTREO)) + 
  geom_bar(stat = "identity", position = "dodge") +
   guides(fill = guide_legend(title = "Muestreo")) + 
    ggtitle("Número de macollas de cebada por medio metro cuadrado") + 
     theme(panel.border = element_rect(fill = "transparent", color = 4,size = 2))  + 
       theme(plot.background = element_rect(fill = "gray86")) +
     theme(plot.background = element_rect(color = "black", size = 2)) +
  xlab("Grupos") + ylab("Número de macollas por medio metro cuadrado") + 
  theme(axis.title.x = element_text(size = 11, color = "black")) + theme(axis.title.y = element_text(size = 13, color = "black")) + theme(axis.text.x = element_text(angle = 0, vjust = 1, hjust = 1)) +   scale_x_continuous(n.breaks = 8)

Datos y gráficas del seguimiento de parametros de crecimiento del cultivo (muestreo 1,2,3)

cultivocebada <- read_excel("cebada 1.xlsx", col_types = c("numeric", 
    "numeric", "numeric", "numeric", "numeric", 
    "numeric", "numeric"))
attach(cultivocebada)
## The following objects are masked from DMTOTAL3:
## 
##     DENMALEZAS, MUESTREO, PUNTOMUES
# Altura de cebada:
psych::describe(ALTURAPROMM2)
##    vars  n  mean    sd median trimmed   mad   min  max range skew kurtosis   se
## X1    1 24 39.56 17.42  33.55   38.72 16.23 15.37 71.5 56.13 0.44    -1.35 3.56
describe.by(ALTURAPROMM2, group = PUNTOMUES)
## Warning: describe.by is deprecated.  Please use the describeBy function
## 
##  Descriptive statistics by group 
## group: 1
##    vars n  mean    sd median trimmed   mad  min   max range skew kurtosis  se
## X1    1 3 43.73 16.63  41.39   43.73 19.26 28.4 61.41 33.01 0.14    -2.33 9.6
## ------------------------------------------------------------ 
## group: 2
##    vars n  mean    sd median trimmed  mad   min  max range skew kurtosis    se
## X1    1 3 44.24 20.39  44.09   44.24 29.9 23.92 64.7 40.78 0.01    -2.33 11.77
## ------------------------------------------------------------ 
## group: 3
##    vars n  mean   sd median trimmed mad   min max range skew kurtosis   se
## X1    1 3 38.15 12.9  31.67   38.15 2.8 29.78  53 23.22 0.38    -2.33 7.45
## ------------------------------------------------------------ 
## group: 4
##    vars n  mean    sd median trimmed   mad min max range skew kurtosis    se
## X1    1 3 38.58 19.46  33.75   38.58 17.42  22  60    38 0.23    -2.33 11.23
## ------------------------------------------------------------ 
## group: 5
##    vars n  mean    sd median trimmed   mad   min  max range skew kurtosis   se
## X1    1 3 37.37 29.96  25.25   37.37 14.65 15.37 71.5 56.13 0.34    -2.33 17.3
## ------------------------------------------------------------ 
## group: 6
##    vars n  mean    sd median trimmed   mad   min  max range skew kurtosis    se
## X1    1 3 34.87 19.02  28.77   34.87 13.52 19.65 56.2 36.55 0.29    -2.33 10.98
## ------------------------------------------------------------ 
## group: 7
##    vars n  mean    sd median trimmed   mad  min max range skew kurtosis    se
## X1    1 3 37.58 17.68  33.35   37.58 16.23 22.4  57  34.6 0.23    -2.33 10.21
## ------------------------------------------------------------ 
## group: 8
##    vars n  mean    sd median trimmed   mad  min  max range skew kurtosis   se
## X1    1 3 41.97 23.72   34.6   41.97 17.49 22.8 68.5  45.7 0.28    -2.33 13.7
# Macollas por planta:
psych::describe(MACOLLASM2)
##    vars  n mean   sd median trimmed  mad min max range skew kurtosis   se
## X1    1 24 8.85 4.97   7.25    8.38 4.23 2.5  23  20.5 0.98     0.44 1.01
describe.by(MACOLLASM2, group = PUNTOMUES)
## Warning: describe.by is deprecated.  Please use the describeBy function
## 
##  Descriptive statistics by group 
## group: 1
##    vars n mean   sd median trimmed  mad min max range skew kurtosis   se
## X1    1 3 7.83 4.48    5.5    7.83 0.74   5  13     8 0.38    -2.33 2.59
## ------------------------------------------------------------ 
## group: 2
##    vars n mean   sd median trimmed  mad min max range skew kurtosis   se
## X1    1 3 9.25 4.99   6.75    9.25 1.11   6  15     9 0.38    -2.33 2.88
## ------------------------------------------------------------ 
## group: 3
##    vars n mean   sd median trimmed  mad min max range skew kurtosis   se
## X1    1 3 9.67 3.79      8    9.67 1.48   7  14     7 0.35    -2.33 2.19
## ------------------------------------------------------------ 
## group: 4
##    vars n  mean    sd median trimmed  mad min max range skew kurtosis   se
## X1    1 3 10.42 11.02   5.75   10.42 4.82 2.5  23  20.5 0.35    -2.33 6.36
## ------------------------------------------------------------ 
## group: 5
##    vars n mean   sd median trimmed  mad min max range skew kurtosis   se
## X1    1 3 6.83 5.35      4    6.83 0.74 3.5  13   9.5 0.38    -2.33 3.09
## ------------------------------------------------------------ 
## group: 6
##    vars n mean   sd median trimmed  mad min max range skew kurtosis  se
## X1    1 3  9.7 5.54    7.5     9.7 2.82 5.6  16  10.4 0.33    -2.33 3.2
## ------------------------------------------------------------ 
## group: 7
##    vars n mean   sd median trimmed  mad min max range skew kurtosis   se
## X1    1 3 8.23 4.75    8.2    8.23 6.97 3.5  13   9.5 0.01    -2.33 2.74
## ------------------------------------------------------------ 
## group: 8
##    vars n mean   sd median trimmed  mad min max range  skew kurtosis  se
## X1    1 3  8.9 2.59    9.7     8.9 1.93   6  11     5 -0.28    -2.33 1.5
# densidad de cebada:
psych::describe(DENSIDADM2)
##    vars  n  mean sd median trimmed   mad min max range skew kurtosis   se
## X1    1 24 94.75 31     95    93.7 28.17  43 159   116 0.25    -0.68 6.33
describe.by(DENSIDADM2, group = PUNTOMUES)
## Warning: describe.by is deprecated.  Please use the describeBy function
## 
##  Descriptive statistics by group 
## group: 1
##    vars n mean    sd median trimmed   mad min max range skew kurtosis    se
## X1    1 3  111 46.57    108     111 62.27  66 159    93 0.06    -2.33 26.89
## ------------------------------------------------------------ 
## group: 2
##    vars n mean sd median trimmed mad min max range skew kurtosis    se
## X1    1 3   81 29     81      81  43  52 110    58    0    -2.33 16.74
## ------------------------------------------------------------ 
## group: 3
##    vars n mean    sd median trimmed   mad min max range skew kurtosis    se
## X1    1 3   93 29.46     99      93 29.65  61 119    58 -0.2    -2.33 17.01
## ------------------------------------------------------------ 
## group: 4
##    vars n   mean    sd median trimmed   mad min max range  skew kurtosis   se
## X1    1 3 122.67 39.31    138  122.67 20.76  78 152    74 -0.33    -2.33 22.7
## ------------------------------------------------------------ 
## group: 5
##    vars n   mean    sd median trimmed   mad min max range skew kurtosis   se
## X1    1 3 108.33 16.26    105  108.33 16.31  94 126    32  0.2    -2.33 9.39
## ------------------------------------------------------------ 
## group: 6
##    vars n mean   sd median trimmed  mad min max range skew kurtosis    se
## X1    1 3   57 21.7     46      57 4.45  43  82    39 0.38    -2.33 12.53
## ------------------------------------------------------------ 
## group: 7
##    vars n  mean    sd median trimmed  mad min max range  skew kurtosis   se
## X1    1 3 89.33 13.32     96   89.33 2.97  74  98    24 -0.38    -2.33 7.69
## ------------------------------------------------------------ 
## group: 8
##    vars n  mean    sd median trimmed   mad min max range skew kurtosis    se
## X1    1 3 95.67 20.11     90   95.67 16.31  79 118    39 0.26    -2.33 11.61
## Gráficos de altura de cebada por metro cuadrado 
par(mfrow=c(1,1))
ggplot(cultivocebada, aes(x = PUNTOMUES, y = ALTURAPROMM2, fill = MUESTREO)) + 
  geom_bar(stat = "identity", position = "dodge") +
  scale_color_manual(labels = c("1", "2","3"))+
  theme(legend.key.size =unit(0.5, 'cm'))+
   guides(fill = guide_legend(title = "Muestreo")) + 
    ggtitle("Altura de cebada por metro cuadrado") + 
     theme(panel.border = element_rect(fill = "transparent", color = 4,size = 2))  + 
       theme(plot.background = element_rect(fill = "gray86")) +
     theme(plot.background = element_rect(color = "black", size = 2)) +
  xlab("Grupos") + ylab("Altura de plantas (cm) por metro cuadrado") + 
  theme(axis.title.x = element_text(size = 11, color = "black")) + theme(axis.title.y = element_text(size = 13, color = "black")) + theme(axis.text.x = element_text(angle = 0, vjust = 1, hjust = 1))+
  scale_x_continuous(n.breaks = 8)

## Gráfico de densidad del cultivo por metro cuadrado  

ggplot(cultivocebada, aes(x = PUNTOMUES, y = DENSIDADM2, fill = MUESTREO)) + 
  geom_bar(stat = "identity", position = "dodge") +
   guides(fill = guide_legend(title = "Muestreo")) + 
    ggtitle("Densisdad de cebada por metro cuadrado") + 
     theme(panel.border = element_rect(fill = "transparent", color = 4,size = 2))  + 
       theme(plot.background = element_rect(fill = "gray86")) +
     theme(plot.background = element_rect(color = "black", size = 2)) +
  xlab("Grupos") + ylab("Número de plantas por metro cuadrado") + 
  theme(axis.title.x = element_text(size = 11, color = "black")) + theme(axis.title.y = element_text(size = 13, color = "black")) + theme(axis.text.x = element_text(angle = 0, vjust = 1, hjust = 1))  + scale_x_continuous(n.breaks = 8)

## Gráfico de número de macollas por metro cuadrado 

ggplot(cultivocebada, aes(x = PUNTOMUES, y = MACOLLASM2, fill = MUESTREO)) + 
  geom_bar(stat = "identity", position = "dodge") +
   guides(fill = guide_legend(title = "Muestreo")) + 
    ggtitle("Número de macollas de cebada por metro cuadrado") + 
     theme(panel.border = element_rect(fill = "transparent", color = 4,size = 2))  + 
       theme(plot.background = element_rect(fill = "gray86")) +
     theme(plot.background = element_rect(color = "black", size = 2)) +
  xlab("Grupos") + ylab("Número de macollas por metro cuadrado") + 
  theme(axis.title.x = element_text(size = 11, color = "black")) + theme(axis.title.y = element_text(size = 13, color = "black")) + theme(axis.text.x = element_text(angle = 0, vjust = 1, hjust = 1)) +   scale_x_continuous(n.breaks = 8)

Datos, gráficas e indices del tercer muestreo (arvenses)

MALEZAS3 <- MALEZAS_CEBADA_3 <- read_excel("MALEZAS CEBADA 3.xlsx", 
    col_types = c("numeric", "text", "numeric", 
        "numeric", "numeric", "numeric", 
        "numeric", "numeric", "numeric", 
        "numeric", "numeric", "numeric"))

attach(MALEZAS3)
## The following object is masked from cultivocebada:
## 
##     MUESTREO
## The following object is masked from DMTOTAL3:
## 
##     MUESTREO
# Densidades:
psych::describe(MALEZAS3)
##                     vars  n  mean     sd median trimmed  mad min    max  range
## MUESTREO               1 18  3.00   0.00   3.00    3.00 0.00   3   3.00   0.00
## ESPECIESGL3*           2 18  9.50   5.34   9.50    9.50 6.67   1  18.00  17.00
## DENSIDAD0.5MT3         3 18 28.39 104.16   2.00    4.12 2.97   0 445.00 445.00
## DENSIDADRELATIVA3      4 18  5.56  20.38   0.39    0.81 0.58   0  87.08  87.08
## FRECUENCIATOTAL3       5 18  0.11   0.18   0.06    0.07 0.09   0   0.75   0.75
## FRECUENCIARELATIVA3    6 18  5.55   9.38   3.21    3.83 4.77   0  38.71  38.71
## COBERTURATOTAL3        7 18 11.78  25.11   1.50    7.81 2.22   0  87.00  87.00
## DOMINANCIARELATIVA3    8 18  5.56  11.84   0.70    3.68 1.05   0  41.04  41.04
## DENSIDADPROMT3         9 18  1.78   6.51   0.13    0.26 0.19   0  27.81  27.81
## COBERTURAPROMT3       10 18  0.79   1.71   0.10    0.50 0.15   0   6.21   6.21
## DENSIDADMCUADRADO3    11 18  7.10  26.04   0.50    1.03 0.74   0 111.25 111.25
## IVI3                  12 18 16.66  39.26   4.32    8.32 6.41   0 166.83 166.83
##                     skew kurtosis    se
## MUESTREO             NaN      NaN  0.00
## ESPECIESGL3*        0.00    -1.40  1.26
## DENSIDAD0.5MT3      3.54    11.22 24.55
## DENSIDADRELATIVA3   3.54    11.22  4.80
## FRECUENCIATOTAL3    2.41     5.73  0.04
## FRECUENCIARELATIVA3 2.42     5.79  2.21
## COBERTURATOTAL3     2.16     3.22  5.92
## DOMINANCIARELATIVA3 2.16     3.22  2.79
## DENSIDADPROMT3      3.54    11.22  1.53
## COBERTURAPROMT3     2.26     3.79  0.40
## DENSIDADMCUADRADO3  3.54    11.22  6.14
## IVI3                3.09     8.92  9.25
describe.by(DENSIDAD0.5MT3, group = ESPECIESGL3)
## Warning: describe.by is deprecated.  Please use the describeBy function
## 
##  Descriptive statistics by group 
## group: Amaranthus hybridus
##    vars n mean sd median trimmed mad min max range skew kurtosis se
## X1    1 1    0 NA      0       0   0   0   0     0   NA       NA NA
## ------------------------------------------------------------ 
## group: Brassica sp
##    vars n mean sd median trimmed mad min max range skew kurtosis se
## X1    1 1    2 NA      2       2   0   2   2     0   NA       NA NA
## ------------------------------------------------------------ 
## group: Chenopodium petiolare Kunth
##    vars n mean sd median trimmed mad min max range skew kurtosis se
## X1    1 1    7 NA      7       7   0   7   7     0   NA       NA NA
## ------------------------------------------------------------ 
## group: Cotula australis
##    vars n mean sd median trimmed mad min max range skew kurtosis se
## X1    1 1    0 NA      0       0   0   0   0     0   NA       NA NA
## ------------------------------------------------------------ 
## group: Fuertesimalva sp.
##    vars n mean sd median trimmed mad min max range skew kurtosis se
## X1    1 1    0 NA      0       0   0   0   0     0   NA       NA NA
## ------------------------------------------------------------ 
## group: Fumaria capreolata
##    vars n mean sd median trimmed mad min max range skew kurtosis se
## X1    1 1   23 NA     23      23   0  23  23     0   NA       NA NA
## ------------------------------------------------------------ 
## group: Galinsoga 
## Ciliata
##    vars n mean sd median trimmed mad min max range skew kurtosis se
## X1    1 1    0 NA      0       0   0   0   0     0   NA       NA NA
## ------------------------------------------------------------ 
## group: Lepidium bipinnatifidum
##    vars n mean sd median trimmed mad min max range skew kurtosis se
## X1    1 1    0 NA      0       0   0   0   0     0   NA       NA NA
## ------------------------------------------------------------ 
## group: Lolium multiflorum
##    vars n mean sd median trimmed mad min max range skew kurtosis se
## X1    1 1    0 NA      0       0   0   0   0     0   NA       NA NA
## ------------------------------------------------------------ 
## group: Malva sylvestris
##    vars n mean sd median trimmed mad min max range skew kurtosis se
## X1    1 1    2 NA      2       2   0   2   2     0   NA       NA NA
## ------------------------------------------------------------ 
## group: Oxalis
##    vars n mean sd median trimmed mad min max range skew kurtosis se
## X1    1 1    2 NA      2       2   0   2   2     0   NA       NA NA
## ------------------------------------------------------------ 
## group: Poa annua
##    vars n mean sd median trimmed mad min max range skew kurtosis se
## X1    1 1    0 NA      0       0   0   0   0     0   NA       NA NA
## ------------------------------------------------------------ 
## group: Polygonum aviculare
##    vars n mean sd median trimmed mad min max range skew kurtosis se
## X1    1 1  445 NA    445     445   0 445 445     0   NA       NA NA
## ------------------------------------------------------------ 
## group: Polygonum segetum
##    vars n mean sd median trimmed mad min max range skew kurtosis se
## X1    1 1   14 NA     14      14   0  14  14     0   NA       NA NA
## ------------------------------------------------------------ 
## group: Raphanus
##    vars n mean sd median trimmed mad min max range skew kurtosis se
## X1    1 1   11 NA     11      11   0  11  11     0   NA       NA NA
## ------------------------------------------------------------ 
## group: Rumex
##    vars n mean sd median trimmed mad min max range skew kurtosis se
## X1    1 1    2 NA      2       2   0   2   2     0   NA       NA NA
## ------------------------------------------------------------ 
## group: Sonchus sp.
##    vars n mean sd median trimmed mad min max range skew kurtosis se
## X1    1 1    0 NA      0       0   0   0   0     0   NA       NA NA
## ------------------------------------------------------------ 
## group: Veronica persica
##    vars n mean sd median trimmed mad min max range skew kurtosis se
## X1    1 1    3 NA      3       3   0   3   3     0   NA       NA NA
# densidad relativa: 
psych::describe(DENSIDADRELATIVA3)
##    vars  n mean    sd median trimmed  mad min   max range skew kurtosis  se
## X1    1 18 5.56 20.38   0.39    0.81 0.58   0 87.08 87.08 3.54    11.22 4.8
describe.by(DENSIDADRELATIVA3, group = ESPECIESGL3)
## Warning: describe.by is deprecated.  Please use the describeBy function
## 
##  Descriptive statistics by group 
## group: Amaranthus hybridus
##    vars n mean sd median trimmed mad min max range skew kurtosis se
## X1    1 1    0 NA      0       0   0   0   0     0   NA       NA NA
## ------------------------------------------------------------ 
## group: Brassica sp
##    vars n mean sd median trimmed mad  min  max range skew kurtosis se
## X1    1 1 0.39 NA   0.39    0.39   0 0.39 0.39     0   NA       NA NA
## ------------------------------------------------------------ 
## group: Chenopodium petiolare Kunth
##    vars n mean sd median trimmed mad  min  max range skew kurtosis se
## X1    1 1 1.37 NA   1.37    1.37   0 1.37 1.37     0   NA       NA NA
## ------------------------------------------------------------ 
## group: Cotula australis
##    vars n mean sd median trimmed mad min max range skew kurtosis se
## X1    1 1    0 NA      0       0   0   0   0     0   NA       NA NA
## ------------------------------------------------------------ 
## group: Fuertesimalva sp.
##    vars n mean sd median trimmed mad min max range skew kurtosis se
## X1    1 1    0 NA      0       0   0   0   0     0   NA       NA NA
## ------------------------------------------------------------ 
## group: Fumaria capreolata
##    vars n mean sd median trimmed mad min max range skew kurtosis se
## X1    1 1  4.5 NA    4.5     4.5   0 4.5 4.5     0   NA       NA NA
## ------------------------------------------------------------ 
## group: Galinsoga 
## Ciliata
##    vars n mean sd median trimmed mad min max range skew kurtosis se
## X1    1 1    0 NA      0       0   0   0   0     0   NA       NA NA
## ------------------------------------------------------------ 
## group: Lepidium bipinnatifidum
##    vars n mean sd median trimmed mad min max range skew kurtosis se
## X1    1 1    0 NA      0       0   0   0   0     0   NA       NA NA
## ------------------------------------------------------------ 
## group: Lolium multiflorum
##    vars n mean sd median trimmed mad min max range skew kurtosis se
## X1    1 1    0 NA      0       0   0   0   0     0   NA       NA NA
## ------------------------------------------------------------ 
## group: Malva sylvestris
##    vars n mean sd median trimmed mad  min  max range skew kurtosis se
## X1    1 1 0.39 NA   0.39    0.39   0 0.39 0.39     0   NA       NA NA
## ------------------------------------------------------------ 
## group: Oxalis
##    vars n mean sd median trimmed mad  min  max range skew kurtosis se
## X1    1 1 0.39 NA   0.39    0.39   0 0.39 0.39     0   NA       NA NA
## ------------------------------------------------------------ 
## group: Poa annua
##    vars n mean sd median trimmed mad min max range skew kurtosis se
## X1    1 1    0 NA      0       0   0   0   0     0   NA       NA NA
## ------------------------------------------------------------ 
## group: Polygonum aviculare
##    vars n  mean sd median trimmed mad   min   max range skew kurtosis se
## X1    1 1 87.08 NA  87.08   87.08   0 87.08 87.08     0   NA       NA NA
## ------------------------------------------------------------ 
## group: Polygonum segetum
##    vars n mean sd median trimmed mad  min  max range skew kurtosis se
## X1    1 1 2.74 NA   2.74    2.74   0 2.74 2.74     0   NA       NA NA
## ------------------------------------------------------------ 
## group: Raphanus
##    vars n mean sd median trimmed mad  min  max range skew kurtosis se
## X1    1 1 2.15 NA   2.15    2.15   0 2.15 2.15     0   NA       NA NA
## ------------------------------------------------------------ 
## group: Rumex
##    vars n mean sd median trimmed mad  min  max range skew kurtosis se
## X1    1 1 0.39 NA   0.39    0.39   0 0.39 0.39     0   NA       NA NA
## ------------------------------------------------------------ 
## group: Sonchus sp.
##    vars n mean sd median trimmed mad min max range skew kurtosis se
## X1    1 1    0 NA      0       0   0   0   0     0   NA       NA NA
## ------------------------------------------------------------ 
## group: Veronica persica
##    vars n mean sd median trimmed mad  min  max range skew kurtosis se
## X1    1 1 0.59 NA   0.59    0.59   0 0.59 0.59     0   NA       NA NA
# frecuencia: 
psych::describe(FRECUENCIATOTAL3)
##    vars  n mean   sd median trimmed  mad min  max range skew kurtosis   se
## X1    1 18 0.11 0.18   0.06    0.07 0.09   0 0.75  0.75 2.41     5.73 0.04
describe.by(FRECUENCIATOTAL3, group = ESPECIESGL3)
## Warning: describe.by is deprecated.  Please use the describeBy function
## 
##  Descriptive statistics by group 
## group: Amaranthus hybridus
##    vars n mean sd median trimmed mad min max range skew kurtosis se
## X1    1 1    0 NA      0       0   0   0   0     0   NA       NA NA
## ------------------------------------------------------------ 
## group: Brassica sp
##    vars n mean sd median trimmed mad  min  max range skew kurtosis se
## X1    1 1 0.06 NA   0.06    0.06   0 0.06 0.06     0   NA       NA NA
## ------------------------------------------------------------ 
## group: Chenopodium petiolare Kunth
##    vars n mean sd median trimmed mad  min  max range skew kurtosis se
## X1    1 1 0.19 NA   0.19    0.19   0 0.19 0.19     0   NA       NA NA
## ------------------------------------------------------------ 
## group: Cotula australis
##    vars n mean sd median trimmed mad min max range skew kurtosis se
## X1    1 1    0 NA      0       0   0   0   0     0   NA       NA NA
## ------------------------------------------------------------ 
## group: Fuertesimalva sp.
##    vars n mean sd median trimmed mad min max range skew kurtosis se
## X1    1 1    0 NA      0       0   0   0   0     0   NA       NA NA
## ------------------------------------------------------------ 
## group: Fumaria capreolata
##    vars n mean sd median trimmed mad  min  max range skew kurtosis se
## X1    1 1 0.13 NA   0.13    0.13   0 0.13 0.13     0   NA       NA NA
## ------------------------------------------------------------ 
## group: Galinsoga 
## Ciliata
##    vars n mean sd median trimmed mad min max range skew kurtosis se
## X1    1 1    0 NA      0       0   0   0   0     0   NA       NA NA
## ------------------------------------------------------------ 
## group: Lepidium bipinnatifidum
##    vars n mean sd median trimmed mad min max range skew kurtosis se
## X1    1 1    0 NA      0       0   0   0   0     0   NA       NA NA
## ------------------------------------------------------------ 
## group: Lolium multiflorum
##    vars n mean sd median trimmed mad min max range skew kurtosis se
## X1    1 1    0 NA      0       0   0   0   0     0   NA       NA NA
## ------------------------------------------------------------ 
## group: Malva sylvestris
##    vars n mean sd median trimmed mad  min  max range skew kurtosis se
## X1    1 1 0.06 NA   0.06    0.06   0 0.06 0.06     0   NA       NA NA
## ------------------------------------------------------------ 
## group: Oxalis
##    vars n mean sd median trimmed mad  min  max range skew kurtosis se
## X1    1 1 0.06 NA   0.06    0.06   0 0.06 0.06     0   NA       NA NA
## ------------------------------------------------------------ 
## group: Poa annua
##    vars n mean sd median trimmed mad min max range skew kurtosis se
## X1    1 1    0 NA      0       0   0   0   0     0   NA       NA NA
## ------------------------------------------------------------ 
## group: Polygonum aviculare
##    vars n mean sd median trimmed mad  min  max range skew kurtosis se
## X1    1 1 0.75 NA   0.75    0.75   0 0.75 0.75     0   NA       NA NA
## ------------------------------------------------------------ 
## group: Polygonum segetum
##    vars n mean sd median trimmed mad  min  max range skew kurtosis se
## X1    1 1 0.25 NA   0.25    0.25   0 0.25 0.25     0   NA       NA NA
## ------------------------------------------------------------ 
## group: Raphanus
##    vars n mean sd median trimmed mad  min  max range skew kurtosis se
## X1    1 1 0.25 NA   0.25    0.25   0 0.25 0.25     0   NA       NA NA
## ------------------------------------------------------------ 
## group: Rumex
##    vars n mean sd median trimmed mad  min  max range skew kurtosis se
## X1    1 1 0.06 NA   0.06    0.06   0 0.06 0.06     0   NA       NA NA
## ------------------------------------------------------------ 
## group: Sonchus sp.
##    vars n mean sd median trimmed mad min max range skew kurtosis se
## X1    1 1    0 NA      0       0   0   0   0     0   NA       NA NA
## ------------------------------------------------------------ 
## group: Veronica persica
##    vars n mean sd median trimmed mad  min  max range skew kurtosis se
## X1    1 1 0.13 NA   0.13    0.13   0 0.13 0.13     0   NA       NA NA
# frecuencia relativa:
psych::describe(FRECUENCIARELATIVA3)
##    vars  n mean   sd median trimmed  mad min   max range skew kurtosis   se
## X1    1 18 5.55 9.38   3.21    3.83 4.77   0 38.71 38.71 2.42     5.79 2.21
describe.by(FRECUENCIARELATIVA3, group = ESPECIESGL3)
## Warning: describe.by is deprecated.  Please use the describeBy function
## 
##  Descriptive statistics by group 
## group: Amaranthus hybridus
##    vars n mean sd median trimmed mad min max range skew kurtosis se
## X1    1 1    0 NA      0       0   0   0   0     0   NA       NA NA
## ------------------------------------------------------------ 
## group: Brassica sp
##    vars n mean sd median trimmed mad  min  max range skew kurtosis se
## X1    1 1 3.23 NA   3.23    3.23   0 3.23 3.23     0   NA       NA NA
## ------------------------------------------------------------ 
## group: Chenopodium petiolare Kunth
##    vars n mean sd median trimmed mad  min  max range skew kurtosis se
## X1    1 1 9.68 NA   9.68    9.68   0 9.68 9.68     0   NA       NA NA
## ------------------------------------------------------------ 
## group: Cotula australis
##    vars n mean sd median trimmed mad min max range skew kurtosis se
## X1    1 1    0 NA      0       0   0   0   0     0   NA       NA NA
## ------------------------------------------------------------ 
## group: Fuertesimalva sp.
##    vars n mean sd median trimmed mad min max range skew kurtosis se
## X1    1 1    0 NA      0       0   0   0   0     0   NA       NA NA
## ------------------------------------------------------------ 
## group: Fumaria capreolata
##    vars n mean sd median trimmed mad  min  max range skew kurtosis se
## X1    1 1 6.45 NA   6.45    6.45   0 6.45 6.45     0   NA       NA NA
## ------------------------------------------------------------ 
## group: Galinsoga 
## Ciliata
##    vars n mean sd median trimmed mad min max range skew kurtosis se
## X1    1 1    0 NA      0       0   0   0   0     0   NA       NA NA
## ------------------------------------------------------------ 
## group: Lepidium bipinnatifidum
##    vars n mean sd median trimmed mad min max range skew kurtosis se
## X1    1 1    0 NA      0       0   0   0   0     0   NA       NA NA
## ------------------------------------------------------------ 
## group: Lolium multiflorum
##    vars n mean sd median trimmed mad min max range skew kurtosis se
## X1    1 1    0 NA      0       0   0   0   0     0   NA       NA NA
## ------------------------------------------------------------ 
## group: Malva sylvestris
##    vars n mean sd median trimmed mad  min  max range skew kurtosis se
## X1    1 1 3.23 NA   3.23    3.23   0 3.23 3.23     0   NA       NA NA
## ------------------------------------------------------------ 
## group: Oxalis
##    vars n mean sd median trimmed mad  min  max range skew kurtosis se
## X1    1 1 3.23 NA   3.23    3.23   0 3.23 3.23     0   NA       NA NA
## ------------------------------------------------------------ 
## group: Poa annua
##    vars n mean sd median trimmed mad min max range skew kurtosis se
## X1    1 1    0 NA      0       0   0   0   0     0   NA       NA NA
## ------------------------------------------------------------ 
## group: Polygonum aviculare
##    vars n  mean sd median trimmed mad   min   max range skew kurtosis se
## X1    1 1 38.71 NA  38.71   38.71   0 38.71 38.71     0   NA       NA NA
## ------------------------------------------------------------ 
## group: Polygonum segetum
##    vars n mean sd median trimmed mad  min  max range skew kurtosis se
## X1    1 1 12.9 NA   12.9    12.9   0 12.9 12.9     0   NA       NA NA
## ------------------------------------------------------------ 
## group: Raphanus
##    vars n mean sd median trimmed mad  min  max range skew kurtosis se
## X1    1 1 12.9 NA   12.9    12.9   0 12.9 12.9     0   NA       NA NA
## ------------------------------------------------------------ 
## group: Rumex
##    vars n mean sd median trimmed mad min max range skew kurtosis se
## X1    1 1  3.2 NA    3.2     3.2   0 3.2 3.2     0   NA       NA NA
## ------------------------------------------------------------ 
## group: Sonchus sp.
##    vars n mean sd median trimmed mad min max range skew kurtosis se
## X1    1 1    0 NA      0       0   0   0   0     0   NA       NA NA
## ------------------------------------------------------------ 
## group: Veronica persica
##    vars n mean sd median trimmed mad  min  max range skew kurtosis se
## X1    1 1 6.45 NA   6.45    6.45   0 6.45 6.45     0   NA       NA NA
# Dominancia: 
psych::describe(COBERTURATOTAL3)
##    vars  n  mean    sd median trimmed  mad min max range skew kurtosis   se
## X1    1 18 11.78 25.11    1.5    7.81 2.22   0  87    87 2.16     3.22 5.92
describe.by(COBERTURATOTAL3, group = ESPECIESGL3)
## Warning: describe.by is deprecated.  Please use the describeBy function
## 
##  Descriptive statistics by group 
## group: Amaranthus hybridus
##    vars n mean sd median trimmed mad min max range skew kurtosis se
## X1    1 1    0 NA      0       0   0   0   0     0   NA       NA NA
## ------------------------------------------------------------ 
## group: Brassica sp
##    vars n mean sd median trimmed mad min max range skew kurtosis se
## X1    1 1    2 NA      2       2   0   2   2     0   NA       NA NA
## ------------------------------------------------------------ 
## group: Chenopodium petiolare Kunth
##    vars n mean sd median trimmed mad min max range skew kurtosis se
## X1    1 1    5 NA      5       5   0   5   5     0   NA       NA NA
## ------------------------------------------------------------ 
## group: Cotula australis
##    vars n mean sd median trimmed mad min max range skew kurtosis se
## X1    1 1    0 NA      0       0   0   0   0     0   NA       NA NA
## ------------------------------------------------------------ 
## group: Fuertesimalva sp.
##    vars n mean sd median trimmed mad min max range skew kurtosis se
## X1    1 1    0 NA      0       0   0   0   0     0   NA       NA NA
## ------------------------------------------------------------ 
## group: Fumaria capreolata
##    vars n mean sd median trimmed mad min max range skew kurtosis se
## X1    1 1   71 NA     71      71   0  71  71     0   NA       NA NA
## ------------------------------------------------------------ 
## group: Galinsoga 
## Ciliata
##    vars n mean sd median trimmed mad min max range skew kurtosis se
## X1    1 1    0 NA      0       0   0   0   0     0   NA       NA NA
## ------------------------------------------------------------ 
## group: Lepidium bipinnatifidum
##    vars n mean sd median trimmed mad min max range skew kurtosis se
## X1    1 1    0 NA      0       0   0   0   0     0   NA       NA NA
## ------------------------------------------------------------ 
## group: Lolium multiflorum
##    vars n mean sd median trimmed mad min max range skew kurtosis se
## X1    1 1    0 NA      0       0   0   0   0     0   NA       NA NA
## ------------------------------------------------------------ 
## group: Malva sylvestris
##    vars n mean sd median trimmed mad min max range skew kurtosis se
## X1    1 1    4 NA      4       4   0   4   4     0   NA       NA NA
## ------------------------------------------------------------ 
## group: Oxalis
##    vars n mean sd median trimmed mad min max range skew kurtosis se
## X1    1 1    1 NA      1       1   0   1   1     0   NA       NA NA
## ------------------------------------------------------------ 
## group: Poa annua
##    vars n mean sd median trimmed mad min max range skew kurtosis se
## X1    1 1    0 NA      0       0   0   0   0     0   NA       NA NA
## ------------------------------------------------------------ 
## group: Polygonum aviculare
##    vars n mean sd median trimmed mad min max range skew kurtosis se
## X1    1 1   87 NA     87      87   0  87  87     0   NA       NA NA
## ------------------------------------------------------------ 
## group: Polygonum segetum
##    vars n mean sd median trimmed mad min max range skew kurtosis se
## X1    1 1   14 NA     14      14   0  14  14     0   NA       NA NA
## ------------------------------------------------------------ 
## group: Raphanus
##    vars n mean sd median trimmed mad min max range skew kurtosis se
## X1    1 1   17 NA     17      17   0  17  17     0   NA       NA NA
## ------------------------------------------------------------ 
## group: Rumex
##    vars n mean sd median trimmed mad min max range skew kurtosis se
## X1    1 1    8 NA      8       8   0   8   8     0   NA       NA NA
## ------------------------------------------------------------ 
## group: Sonchus sp.
##    vars n mean sd median trimmed mad min max range skew kurtosis se
## X1    1 1    0 NA      0       0   0   0   0     0   NA       NA NA
## ------------------------------------------------------------ 
## group: Veronica persica
##    vars n mean sd median trimmed mad min max range skew kurtosis se
## X1    1 1    3 NA      3       3   0   3   3     0   NA       NA NA
#Dominancia relativa:
psych::describe(DOMINANCIARELATIVA3)
##    vars  n mean    sd median trimmed  mad min   max range skew kurtosis   se
## X1    1 18 5.56 11.84    0.7    3.68 1.05   0 41.04 41.04 2.16     3.22 2.79
describe.by(DOMINANCIARELATIVA3, group = ESPECIESGL3)
## Warning: describe.by is deprecated.  Please use the describeBy function
## 
##  Descriptive statistics by group 
## group: Amaranthus hybridus
##    vars n mean sd median trimmed mad min max range skew kurtosis se
## X1    1 1    0 NA      0       0   0   0   0     0   NA       NA NA
## ------------------------------------------------------------ 
## group: Brassica sp
##    vars n mean sd median trimmed mad  min  max range skew kurtosis se
## X1    1 1 0.94 NA   0.94    0.94   0 0.94 0.94     0   NA       NA NA
## ------------------------------------------------------------ 
## group: Chenopodium petiolare Kunth
##    vars n mean sd median trimmed mad  min  max range skew kurtosis se
## X1    1 1 2.36 NA   2.36    2.36   0 2.36 2.36     0   NA       NA NA
## ------------------------------------------------------------ 
## group: Cotula australis
##    vars n mean sd median trimmed mad min max range skew kurtosis se
## X1    1 1    0 NA      0       0   0   0   0     0   NA       NA NA
## ------------------------------------------------------------ 
## group: Fuertesimalva sp.
##    vars n mean sd median trimmed mad min max range skew kurtosis se
## X1    1 1    0 NA      0       0   0   0   0     0   NA       NA NA
## ------------------------------------------------------------ 
## group: Fumaria capreolata
##    vars n  mean sd median trimmed mad   min   max range skew kurtosis se
## X1    1 1 33.49 NA  33.49   33.49   0 33.49 33.49     0   NA       NA NA
## ------------------------------------------------------------ 
## group: Galinsoga 
## Ciliata
##    vars n mean sd median trimmed mad min max range skew kurtosis se
## X1    1 1    0 NA      0       0   0   0   0     0   NA       NA NA
## ------------------------------------------------------------ 
## group: Lepidium bipinnatifidum
##    vars n mean sd median trimmed mad min max range skew kurtosis se
## X1    1 1    0 NA      0       0   0   0   0     0   NA       NA NA
## ------------------------------------------------------------ 
## group: Lolium multiflorum
##    vars n mean sd median trimmed mad min max range skew kurtosis se
## X1    1 1    0 NA      0       0   0   0   0     0   NA       NA NA
## ------------------------------------------------------------ 
## group: Malva sylvestris
##    vars n mean sd median trimmed mad  min  max range skew kurtosis se
## X1    1 1 1.89 NA   1.89    1.89   0 1.89 1.89     0   NA       NA NA
## ------------------------------------------------------------ 
## group: Oxalis
##    vars n mean sd median trimmed mad  min  max range skew kurtosis se
## X1    1 1 0.47 NA   0.47    0.47   0 0.47 0.47     0   NA       NA NA
## ------------------------------------------------------------ 
## group: Poa annua
##    vars n mean sd median trimmed mad min max range skew kurtosis se
## X1    1 1    0 NA      0       0   0   0   0     0   NA       NA NA
## ------------------------------------------------------------ 
## group: Polygonum aviculare
##    vars n  mean sd median trimmed mad   min   max range skew kurtosis se
## X1    1 1 41.04 NA  41.04   41.04   0 41.04 41.04     0   NA       NA NA
## ------------------------------------------------------------ 
## group: Polygonum segetum
##    vars n mean sd median trimmed mad min max range skew kurtosis se
## X1    1 1  6.6 NA    6.6     6.6   0 6.6 6.6     0   NA       NA NA
## ------------------------------------------------------------ 
## group: Raphanus
##    vars n mean sd median trimmed mad  min  max range skew kurtosis se
## X1    1 1 8.02 NA   8.02    8.02   0 8.02 8.02     0   NA       NA NA
## ------------------------------------------------------------ 
## group: Rumex
##    vars n mean sd median trimmed mad  min  max range skew kurtosis se
## X1    1 1 3.77 NA   3.77    3.77   0 3.77 3.77     0   NA       NA NA
## ------------------------------------------------------------ 
## group: Sonchus sp.
##    vars n mean sd median trimmed mad min max range skew kurtosis se
## X1    1 1    0 NA      0       0   0   0   0     0   NA       NA NA
## ------------------------------------------------------------ 
## group: Veronica persica
##    vars n mean sd median trimmed mad  min  max range skew kurtosis se
## X1    1 1 1.42 NA   1.42    1.42   0 1.42 1.42     0   NA       NA NA
## Resultados gráficos por especie 
#Densidad por medio metro cuadrado de cada especie en los 4 metros cuadrados en total

ggplot(MALEZAS3, aes(x = ESPECIESGL3, y = DENSIDADMCUADRADO3, fill = ESPECIESGL3)) + 
  geom_bar(stat = "identity") +
  theme(legend.key.size =unit(0.3, 'cm'))+
   guides(fill = guide_legend(title = "ESPECIES")) + 
    ggtitle("Densidad de especies por medio metro cuadrado") + 
     theme(panel.border = element_rect(fill = "transparent", color = 4,size = 2))  + 
       theme(plot.background = element_rect(fill = "gray86")) +
     theme(plot.background = element_rect(color = "black", size = 2)) +
  xlab("Especies") + ylab("Número de individuos") + 
  theme(axis.title.x = element_text(size = 11, color = "black")) + theme(axis.title.y = element_text(size = 13, color = "black")) + theme(axis.text.x = element_text(angle = 60, vjust = 1, hjust = 1))

#Densidad relativa de cada especie en los 4 metros cuadrados en total

ggplot(MALEZAS3, aes(x = ESPECIESGL3, y = DENSIDADRELATIVA3, fill = ESPECIESGL3)) + 
  geom_bar(stat = "identity") +
    theme(legend.key.size =unit(0.3, 'cm'))+
   guides(fill = guide_legend(title = "ESPECIES")) + 
    ggtitle("Densidad relativa de especies en las 16 muestras (%)") + 
     theme(panel.border = element_rect(fill = "transparent", color = 4,size = 2))  + 
       theme(plot.background = element_rect(fill = "gray86")) +
     theme(plot.background = element_rect(color = "black", size = 2)) +
  xlab("Especies") + ylab("Densidad relativa %") + 
  theme(axis.title.x = element_text(size = 13, color = "black")) + theme(axis.title.y = element_text(size = 13, color = "black")) + theme(axis.text.x = element_text(angle = 60, vjust = 1, hjust = 1))

#Frecuencia relativa de cada especie en los 4 metros cuadrados en total

ggplot(MALEZAS3, aes(x = ESPECIESGL3, y = FRECUENCIARELATIVA3, fill = ESPECIESGL3)) + 
  geom_bar(stat = "identity") +
    theme(legend.key.size =unit(0.3, 'cm'))+
   guides(fill = guide_legend(title = "ESPECIES")) + 
    ggtitle("Frecuencia relativa de especies en las 16 muestras (%)") + 
     theme(panel.border = element_rect(fill = "transparent", color = 4,size = 2))  + 
       theme(plot.background = element_rect(fill = "gray86")) +
     theme(plot.background = element_rect(color = "black", size = 2)) +
  xlab("Especies") + ylab("frecuencia relativa %") + 
  theme(axis.title.x = element_text(size = 13, color = "black")) + theme(axis.title.y = element_text(size = 13, color = "black")) + theme(axis.text.x = element_text(angle = 60, vjust = 1, hjust = 1))

#Dominancia relativa de cada especie en los 4 metros cuadrados en total

ggplot(MALEZAS3, aes(x = ESPECIESGL3, y = DOMINANCIARELATIVA3, fill = ESPECIESGL3)) + 
  geom_bar(stat = "identity") +
    theme(legend.key.size =unit(0.3, 'cm'))+
   guides(fill = guide_legend(title = "ESPECIES")) + 
    ggtitle("Dominancia relativa de especies en las 16 muestras (%)") + 
     theme(panel.border = element_rect(fill = "transparent", color = 4,size = 2))  + 
       theme(plot.background = element_rect(fill = "gray86")) +
     theme(plot.background = element_rect(color = "black", size = 2)) +
  xlab("Especies") + ylab("Dominancia relativa %") + 
  theme(axis.title.x = element_text(size = 13, color = "black")) + theme(axis.title.y = element_text(size = 13, color = "black")) + theme(axis.text.x = element_text(angle = 60, vjust = 1, hjust = 1))

#Indice de valor de importancia  cada especie en los 4 metros cuadrados en total

ggplot(MALEZAS3, aes(x = ESPECIESGL3, y = IVI3, fill = ESPECIESGL3)) + 
  geom_bar(stat = "identity") +
    theme(legend.key.size =unit(0.3, 'cm'))+
   guides(fill = guide_legend(title = "ESPECIES")) + 
    ggtitle("Indice de valor de importancia") + 
     theme(panel.border = element_rect(fill = "transparent", color = 4,size = 2))  + 
       theme(plot.background = element_rect(fill = "gray86")) +
     theme(plot.background = element_rect(color = "black", size = 2)) +
  xlab("Especies") + ylab("Valor de importancia") + 
  theme(axis.title.x = element_text(size = 13, color = "black")) + theme(axis.title.y = element_text(size = 13, color = "black")) + theme(axis.text.x = element_text(angle = 60, vjust = 1, hjust = 1))

Gráficos de cobertura y densidad (Cultivo Vs arvenses)

coberturas3 <- cobertura_arv_y_cult <- read_excel("cobertura arv y cult.xlsx", 
    col_types = c("numeric", "numeric", "text", 
        "numeric", "numeric"))

attach(coberturas3)

# densidad y cobertura del cultivo y el total de arvenses por metro cuadrado 

  psych::describe(DENSIDADESCA3)
##    vars  n mean    sd median trimmed   mad min max range skew kurtosis    se
## X1    1 16   67 51.93     63   58.64 23.72  14 237   223    2     4.38 12.98
describe.by(DENSIDADESCA3, group = MUESTRASMC3)
## Warning: describe.by is deprecated.  Please use the describeBy function
## 
##  Descriptive statistics by group 
## group: arvenses
##    vars n  mean    sd median trimmed   mad min max range skew kurtosis    se
## X1    1 8 63.88 73.91     44   63.88 40.77  14 237   223 1.49     0.84 26.13
## ------------------------------------------------------------ 
## group: cultivo
##    vars n  mean    sd median trimmed   mad min max range skew kurtosis   se
## X1    1 8 70.12 17.13     70   70.12 20.02  46  94    48 0.03    -1.64 6.05
psych::describe(COBERTURASCA3)
##    vars  n  mean    sd median trimmed    mad min max range skew kurtosis    se
## X1    1 16 96.94 79.82    117   93.64 105.26   0 240   240 0.09    -1.58 19.95
describe.by(COBERTURASCA3, group = MUESTRASMC3)
## Warning: describe.by is deprecated.  Please use the describeBy function
## 
##  Descriptive statistics by group 
## group: arvenses
##    vars n mean    sd median trimmed   mad min max range skew kurtosis    se
## X1    1 8 26.5 34.07   14.5    26.5 14.08   0 102   102 1.25     0.13 12.05
## ------------------------------------------------------------ 
## group: cultivo
##    vars n   mean    sd median trimmed   mad min max range skew kurtosis    se
## X1    1 8 167.38 33.92    160  167.38 20.02 132 240   108 1.05    -0.16 11.99
## Gráfico de la cobertura del cultivo y arvenses por metro cuadrado 

ggplot(coberturas3, aes(x = PMUEST3, y = COBERTURASCA3, fill = MUESTRASMC3)) + 
  geom_bar(stat = "identity", position = "dodge") +
   guides(fill = guide_legend(title = "Especie")) + 
    ggtitle("Cobertura de cebada y arvenses por medio metro cuadrado") + 
     theme(panel.border = element_rect(fill = "transparent", color = 4,size = 2))  + 
       theme(plot.background = element_rect(fill = "gray86")) +
     theme(plot.background = element_rect(color = "black", size = 2)) +
  xlab("Grupos") + ylab("cobertura por medio metro cuadrado") + 
  theme(axis.title.x = element_text(size = 11, color = "black")) + theme(axis.title.y = element_text(size = 13, color = "black")) + theme(axis.text.x = element_text(angle = 0, vjust = 1, hjust = 1)) +   scale_x_continuous(n.breaks = 8)

## Gráfico de la densidad del cultivo y arvences por metro cuadrado 

ggplot(coberturas3, aes(x = PMUEST3, y = DENSIDADESCA3, fill = MUESTRASMC3)) + 
  geom_bar(stat = "identity", position = "dodge") +
   guides(fill = guide_legend(title = "Especie")) + 
    ggtitle("Densidad de cebada y arvenses por medio metro cuadrado") + 
     theme(panel.border = element_rect(fill = "transparent", color = 4,size = 2))  + 
       theme(plot.background = element_rect(fill = "gray86")) +
     theme(plot.background = element_rect(color = "black", size = 2)) +
  xlab("Grupos") + ylab("densidad por medio metro cuadrado") + 
  theme(axis.title.x = element_text(size = 11, color = "black")) + theme(axis.title.y = element_text(size = 13, color = "black")) + theme(axis.text.x = element_text(angle = 0, vjust = 1, hjust = 1)) +   scale_x_continuous(n.breaks = 8)

INDICE, DENSIDAD, COBERTURA Y FRECUENCIA DE MALEZAS EN LOS 3 MUESTREOS

Último muestreo medio metro cuadrado por grupo

DMTOTAL <- read_excel("MALEZASGRUPALES.xlsx", 
    col_types = c("numeric", "text", "numeric", 
        "numeric", "numeric", "numeric", 
        "numeric", "numeric", "numeric", 
        "numeric", "numeric", "numeric"))
attach(DMTOTAL)
## The following object is masked from MALEZAS3:
## 
##     MUESTREO
## The following object is masked from cultivocebada:
## 
##     MUESTREO
## The following object is masked from DMTOTAL3:
## 
##     MUESTREO
## Resultados gráficos por Muestreo

ggplot(DMTOTAL, aes(x = ESPECIESGL, y = DENSIDADMCUADRADO, fill = MUESTREO)) + 
  geom_bar(stat = "identity") +
   guides(fill = guide_legend(title = "Muestreo")) + 
    ggtitle("Densidad de especies por metro cuadrado") + 
     theme(panel.border = element_rect(fill = "transparent", color = 4,size = 2))  + 
       theme(plot.background = element_rect(fill = "gray86")) +
     theme(plot.background = element_rect(color = "black", size = 2)) +
  xlab("Especies") + ylab("Número de individuos") + 
  theme(axis.title.x = element_text(size = 13, color = "black")) + theme(axis.title.y = element_text(size = 13, color = "black")) + theme(axis.text.x = element_text(angle = 45, vjust = 1, hjust = 1))

#Densidad relativa de cada especie en los 8 metros cuadrados en total

ggplot(DMTOTAL, aes(x = ESPECIESGL, y = DENSIDADRELATIVA, fill = MUESTREO)) + 
  geom_bar(stat = "identity") +
   guides(fill = guide_legend(title = "Muestreo")) + 
    ggtitle("Densidad relativa de especies en las 32 muestras") + 
     theme(panel.border = element_rect(fill = "transparent", color = 4,size = 2))  + 
       theme(plot.background = element_rect(fill = "gray86")) +
     theme(plot.background = element_rect(color = "black", size = 2)) +
  xlab("Especies") + ylab("Densidad relativa %") + 
  theme(axis.title.x = element_text(size = 13, color = "black")) + theme(axis.title.y = element_text(size = 13, color = "black")) + theme(axis.text.x = element_text(angle = 45, vjust = 1, hjust = 1))

#Frecuencia relativa de cada especie en los 8 metros cuadrados en total

ggplot(DMTOTAL, aes(x = ESPECIESGL, y = FRECUENCIARELATIVA, fill = MUESTREO)) + 
  geom_bar(stat = "identity") +
   guides(fill = guide_legend(title = "Muestrep")) + 
    ggtitle("Frecuencia relativa de especies en las 32 muestras") + 
     theme(panel.border = element_rect(fill = "transparent", color = 4,size = 2))  + 
       theme(plot.background = element_rect(fill = "gray86")) +
     theme(plot.background = element_rect(color = "black", size = 2)) +
  xlab("Especies") + ylab("frecuencia relativa %") + 
  theme(axis.title.x = element_text(size = 13, color = "black")) + theme(axis.title.y = element_text(size = 13, color = "black")) + theme(axis.text.x = element_text(angle = 45, vjust = 1, hjust = 1))

#Dominancia relativa de cada especie en los 8 metros cuadrados en total

ggplot(DMTOTAL, aes(x = ESPECIESGL, y = DOMINANCIARELATIVA, fill = MUESTREO)) + 
  geom_bar(stat = "identity") +
   guides(fill = guide_legend(title = "Muestreo")) + 
    ggtitle("Dominancia relativa de especies en las 32 muestras") + 
     theme(panel.border = element_rect(fill = "transparent", color = 4,size = 2))  + 
       theme(plot.background = element_rect(fill = "gray86")) +
     theme(plot.background = element_rect(color = "black", size = 2)) +
  xlab("Especies") + ylab("Dominancia relativa %") + 
  theme(axis.title.x = element_text(size = 13, color = "black")) + theme(axis.title.y = element_text(size = 13, color = "black")) + theme(axis.text.x = element_text(angle = 45, vjust = 1, hjust = 1))

#Dominancia relativa de cada especie en los 8 metros cuadrados en total

ggplot(DMTOTAL, aes(x = ESPECIESGL, y = IVI, fill = MUESTREO)) + 
  geom_bar(stat = "identity") +
   guides(fill = guide_legend(title = "Muestreo")) + 
    ggtitle("Indice de valor de importancia") + 
     theme(panel.border = element_rect(fill = "transparent", color = 4,size = 2))  + 
       theme(plot.background = element_rect(fill = "gray86")) +
     theme(plot.background = element_rect(color = "black", size = 2)) +
  xlab("Especies") + ylab("Valor de importancia") + 
  theme(axis.title.x = element_text(size = 13, color = "black")) + theme(axis.title.y = element_text(size = 13, color = "black")) + theme(axis.text.x = element_text(angle = 45, vjust = 1, hjust = 1))

AVENA

Datos y gráficas del seguimiento de parametros de crecimiento del cultivo de avena (muestreo 3)

variables_avena <- read_excel("variables avena.xlsx", 
    col_types = c("numeric", "numeric", "numeric", 
        "numeric", "numeric", "numeric", 
        "numeric"))

attach(variables_avena)

# Altura de avena:
psych::describe(ALTURAPROMAVM2)
##    vars n  mean   sd median trimmed  mad   min   max range  skew kurtosis   se
## X1    1 8 40.53 4.14  41.15   40.53 4.37 33.88 45.25 11.37 -0.36    -1.59 1.46
# Macollas por planta:
psych::describe(MACOLLASAVM2)
##    vars n  mean   sd median trimmed  mad min max range skew kurtosis   se
## X1    1 8 18.88 4.61   17.5   18.88 4.45  14  27    13 0.54    -1.29 1.63
# densidad de cebada:
psych::describe(DENSIDADAVM2)
##    vars n  mean    sd median trimmed   mad min max range skew kurtosis   se
## X1    1 8 122.5 19.16  122.5   122.5 27.43 102 153    51 0.22    -1.64 6.77
## Gráficos de altura de cebada por metro cuadrado 

par(mfrow=c(1,1))
ggplot(variables_avena, aes(x = PUNTOMUESAV, y = ALTURAPROMAVM2, fill = MUESTREOAV)) + 
  geom_bar(stat = "identity", position = "dodge") +
   guides(fill = guide_legend(title = "Muestreo")) + 
    ggtitle("Altura de avena por metro cuadrado") + 
     theme(panel.border = element_rect(fill = "transparent", color = 4,size = 2))  + 
       theme(plot.background = element_rect(fill = "gray86")) +
     theme(plot.background = element_rect(color = "black", size = 2)) +
  xlab("Grupos") + ylab("Altura de plantas (cm) por metro cuadrado") + 
  theme(axis.title.x = element_text(size = 11, color = "black")) + theme(axis.title.y = element_text(size = 13, color = "black")) + theme(axis.text.x = element_text(angle = 0, vjust = 1, hjust = 1))+
  scale_x_continuous(n.breaks = 8)

## Gráfico de densidad del cultivo por metro cuadrado  

ggplot(variables_avena, aes(x = PUNTOMUESAV, y = DENSIDADAVM2, fill = MUESTREOAV)) + 
  geom_bar(stat = "identity", position = "dodge") +
   guides(fill = guide_legend(title = "Muestreo")) + 
    ggtitle("Densisdad de avena por metro cuadrado") + 
     theme(panel.border = element_rect(fill = "transparent", color = 4,size = 2))  + 
       theme(plot.background = element_rect(fill = "gray86")) +
     theme(plot.background = element_rect(color = "black", size = 2)) +
  xlab("Grupos") + ylab("Número de plantas por metro cuadrado") + 
  theme(axis.title.x = element_text(size = 11, color = "black")) + theme(axis.title.y = element_text(size = 13, color = "black")) + theme(axis.text.x = element_text(angle = 0, vjust = 1, hjust = 1))  + scale_x_continuous(n.breaks = 8)

## Gráfico de número de macollas por metro cuadrado 

ggplot(variables_avena, aes(x = PUNTOMUESAV, y = MACOLLASAVM2, fill = MUESTREOAV)) + 
  geom_bar(stat = "identity", position = "dodge") +
   guides(fill = guide_legend(title = "Muestreo")) + 
    ggtitle("Número de macollas de avena por metro cuadrado") + 
     theme(panel.border = element_rect(fill = "transparent", color = 4,size = 2))  + 
       theme(plot.background = element_rect(fill = "gray86")) +
     theme(plot.background = element_rect(color = "black", size = 2)) +
  xlab("Grupos") + ylab("Número de macollas por metro cuadrado") + 
  theme(axis.title.x = element_text(size = 11, color = "black")) + theme(axis.title.y = element_text(size = 13, color = "black")) + theme(axis.text.x = element_text(angle = 0, vjust = 1, hjust = 1)) +   scale_x_continuous(n.breaks = 8)

Datos, gráficas e indices del tercer muestreo (arvenses)

MALEZAS_AVENA <- read_excel("MALEZAS AVENA.xlsx", 
    col_types = c("numeric", "text", "numeric", 
        "numeric", "numeric", "numeric", 
        "numeric", "numeric", "numeric", 
        "numeric", "numeric", "numeric"))

attach(MALEZAS_AVENA)
## The following object is masked from DMTOTAL:
## 
##     MUESTREO
## The following object is masked from MALEZAS3:
## 
##     MUESTREO
## The following object is masked from cultivocebada:
## 
##     MUESTREO
## The following object is masked from DMTOTAL3:
## 
##     MUESTREO
# Densidades:
psych::describe(MALEZAS_AVENA)
##                       vars  n  mean     sd median trimmed   mad min    max
## MUESTREO                 1 18  3.00   0.00   3.00    3.00  0.00   3   3.00
## ESPECIESAV*              2 18  9.50   5.34   9.50    9.50  6.67   1  18.00
## DENSIDADAVT3             3 18 42.78  83.75  14.00   26.94 16.31   0 339.00
## DENSIDADRELATIVA3AV      4 18  5.56  11.62   1.69    3.25  1.96   0  48.07
## FRECUENCIATOTAL3AV       5 18  0.16   0.22   0.11    0.12  0.12   0   0.94
## FRECUENCIARELATIVA3AV    6 18  5.55   7.83   3.89    4.16  4.12   0  33.33
## COBERTURATOTAL3AV        7 18 40.23 126.55   6.00   11.20  8.90   0 545.00
## DOMINANCIARELATIVA3AV    8 18  5.56  17.48   0.83    1.55  1.23   0  75.26
## DENSIDADPROMT3AV         9 18  1.44   3.01   0.44    0.84  0.51   0  12.47
## COBERTURAPROMT3AV       10 18  1.26   3.95   0.19    0.36  0.28   0  17.03
## DENSIDADMCUADRADO3AV    11 18  5.35  10.47   1.75    3.37  2.04   0  42.38
## IVI3AV                  12 18 16.66  36.23   6.27    8.96  7.83   0 156.66
##                        range skew kurtosis    se
## MUESTREO                0.00  NaN      NaN  0.00
## ESPECIESAV*            17.00 0.00    -1.40  1.26
## DENSIDADAVT3          339.00 2.52     5.85 19.74
## DENSIDADRELATIVA3AV    48.07 2.74     7.03  2.74
## FRECUENCIATOTAL3AV      0.94 2.43     5.95  0.05
## FRECUENCIARELATIVA3AV  33.33 2.43     5.92  1.85
## COBERTURATOTAL3AV     545.00 3.51    11.07 29.83
## DOMINANCIARELATIVA3AV  75.26 3.51    11.07  4.12
## DENSIDADPROMT3AV       12.47 2.74     7.04  0.71
## COBERTURAPROMT3AV      17.03 3.51    11.06  0.93
## DENSIDADMCUADRADO3AV   42.38 2.52     5.85  2.47
## IVI3AV                156.66 3.18     9.41  8.54
# densidad relativa: 
psych::describe(DENSIDADRELATIVA3AV)
##    vars  n mean    sd median trimmed  mad min   max range skew kurtosis   se
## X1    1 18 5.56 11.62   1.69    3.25 1.96   0 48.07 48.07 2.74     7.03 2.74
# frecuencia: 
psych::describe(FRECUENCIATOTAL3AV)
##    vars  n mean   sd median trimmed  mad min  max range skew kurtosis   se
## X1    1 18 0.16 0.22   0.11    0.12 0.12   0 0.94  0.94 2.43     5.95 0.05
# frecuencia relativa:
psych::describe(FRECUENCIARELATIVA3AV)
##    vars  n mean   sd median trimmed  mad min   max range skew kurtosis   se
## X1    1 18 5.55 7.83   3.89    4.16 4.12   0 33.33 33.33 2.43     5.92 1.85
# Dominancia: 
psych::describe(COBERTURATOTAL3AV)
##    vars  n  mean     sd median trimmed mad min max range skew kurtosis    se
## X1    1 18 40.23 126.55      6    11.2 8.9   0 545   545 3.51    11.07 29.83
#Dominancia relativa:
psych::describe(DOMINANCIARELATIVA3AV)
##    vars  n mean    sd median trimmed  mad min   max range skew kurtosis   se
## X1    1 18 5.56 17.48   0.83    1.55 1.23   0 75.26 75.26 3.51    11.07 4.12
## Resultados gráficos por especie 
#Densidad por medio metro cuadrado de cada especie en los 4 metros cuadrados en total

ggplot(MALEZAS_AVENA, aes(x = ESPECIESAV, y = DENSIDADMCUADRADO3AV, fill = ESPECIESAV)) + 
  geom_bar(stat = "identity") +
    theme(legend.key.size =unit(0.3, 'cm'))+
   guides(fill = guide_legend(title = "ESPECIES")) + 
    ggtitle("Densidad de especies por medio metro cuadrado") + 
     theme(panel.border = element_rect(fill = "transparent", color = 4,size = 2))  + 
       theme(plot.background = element_rect(fill = "gray86")) +
     theme(plot.background = element_rect(color = "black", size = 2)) +
  xlab("Especies") + ylab("Número de individuos") + 
  theme(axis.title.x = element_text(size = 11, color = "black")) + theme(axis.title.y = element_text(size = 13, color = "black")) + theme(axis.text.x = element_text(angle = 60, vjust = 1, hjust = 1))

#Densidad relativa de cada especie en los 4 metros cuadrados en total

ggplot(MALEZAS_AVENA, aes(x = ESPECIESAV, y = DENSIDADRELATIVA3AV, fill = ESPECIESAV)) + 
  geom_bar(stat = "identity") +
    theme(legend.key.size =unit(0.3, 'cm'))+
   guides(fill = guide_legend(title = "ESPECIES")) + 
    ggtitle("Densidad relativa de especies en las 32 muestras (%)") + 
     theme(panel.border = element_rect(fill = "transparent", color = 4,size = 2))  + 
       theme(plot.background = element_rect(fill = "gray86")) +
     theme(plot.background = element_rect(color = "black", size = 2)) +
  xlab("Especies") + ylab("Densidad relativa %") + 
  theme(axis.title.x = element_text(size = 13, color = "black")) + theme(axis.title.y = element_text(size = 13, color = "black")) + theme(axis.text.x = element_text(angle = 60, vjust = 1, hjust = 1))

#Frecuencia relativa de cada especie en los 4 metros cuadrados en total

ggplot(MALEZAS_AVENA, aes(x = ESPECIESAV, y = FRECUENCIARELATIVA3AV, fill = ESPECIESAV)) + 
  geom_bar(stat = "identity") +
    theme(legend.key.size =unit(0.3, 'cm'))+
   guides(fill = guide_legend(title = "ESPECIES")) + 
    ggtitle("Frecuencia relativa de especies en las 32 muestras (%)") + 
     theme(panel.border = element_rect(fill = "transparent", color = 4,size = 2))  + 
       theme(plot.background = element_rect(fill = "gray86")) +
     theme(plot.background = element_rect(color = "black", size = 2)) +
  xlab("Especies") + ylab("frecuencia relativa %") + 
  theme(axis.title.x = element_text(size = 13, color = "black")) + theme(axis.title.y = element_text(size = 13, color = "black")) + theme(axis.text.x = element_text(angle = 60, vjust = 1, hjust = 1))

#Dominancia relativa de cada especie en los 4 metros cuadrados en total

ggplot(MALEZAS_AVENA, aes(x = ESPECIESAV, y = DOMINANCIARELATIVA3AV, fill = ESPECIESAV)) + 
  geom_bar(stat = "identity") +
    theme(legend.key.size =unit(0.3, 'cm'))+
   guides(fill = guide_legend(title = "ESPECIES")) + 
    ggtitle("Dominancia relativa de especies en las 32 muestras (%)") + 
     theme(panel.border = element_rect(fill = "transparent", color = 4,size = 2))  + 
       theme(plot.background = element_rect(fill = "gray86")) +
     theme(plot.background = element_rect(color = "black", size = 2)) +
  xlab("Especies") + ylab("Dominancia relativa %") + 
  theme(axis.title.x = element_text(size = 13, color = "black")) + theme(axis.title.y = element_text(size = 13, color = "black")) + theme(axis.text.x = element_text(angle = 60, vjust = 1, hjust = 1))

#Indice de valor de importancia  cada especie en los 8 metros cuadrados en total

ggplot(MALEZAS_AVENA, aes(x = ESPECIESAV, y = IVI3AV, fill = ESPECIESAV)) + 
  geom_bar(stat = "identity") +
    theme(legend.key.size =unit(0.3, 'cm'))+
   guides(fill = guide_legend(title = "ESPECIES")) + 
    ggtitle("Indice de valor de importancia") + 
     theme(panel.border = element_rect(fill = "transparent", color = 4,size = 2))  + 
       theme(plot.background = element_rect(fill = "gray86")) +
     theme(plot.background = element_rect(color = "black", size = 2)) +
  xlab("Especies") + ylab("Valor de importancia") + 
  theme(axis.title.x = element_text(size = 13, color = "black")) + theme(axis.title.y = element_text(size = 13, color = "black")) + theme(axis.text.x = element_text(angle = 60, vjust = 1, hjust = 1))

Gráficos de cobertura y densidad (Cultivo Vs arvenses)

COBERTURA_AVENA_Y_MALEZAS <- read_excel("COBERTURA AVENA Y MALEZAS.xlsx", 
    col_types = c("numeric", "numeric", "text", 
        "numeric", "numeric"))
attach(COBERTURA_AVENA_Y_MALEZAS)

# densidad y cobertura del cultivo y el total de arvenses por metro cuadrado 

psych::describe(DENSIDADESCA3AV)
##    vars  n   mean    sd median trimmed   mad min max range  skew kurtosis    se
## X1    1 16 113.12 41.41  118.5  115.21 28.91  27 170   143 -0.61    -0.54 10.35
psych::describe(COBERTURASCA3AV)
##    vars  n  mean    sd median trimmed   mad  min max range skew kurtosis    se
## X1    1 16 159.7 87.98    138  160.36 99.33 30.2 280 249.8 0.19    -1.58 21.99
describe.by(COBERTURASCA3, group = MUESTRASMC3)
## Warning: describe.by is deprecated.  Please use the describeBy function
## 
##  Descriptive statistics by group 
## group: arvenses
##    vars n mean    sd median trimmed   mad min max range skew kurtosis    se
## X1    1 8 26.5 34.07   14.5    26.5 14.08   0 102   102 1.25     0.13 12.05
## ------------------------------------------------------------ 
## group: cultivo
##    vars n   mean    sd median trimmed   mad min max range skew kurtosis    se
## X1    1 8 167.38 33.92    160  167.38 20.02 132 240   108 1.05    -0.16 11.99
## Gráfico de la cobertura del cultivo y arvenses por metro cuadrado 

ggplot(COBERTURA_AVENA_Y_MALEZAS, aes(x = PMUEST3AV, y = COBERTURASCA3AV, fill = MUESTRASMC3AV)) + 
  geom_bar(stat = "identity", position = "dodge") +
   guides(fill = guide_legend(title = "Especie")) + 
    ggtitle("Cobertura de avena y arvenses por metro cuadrado") + 
     theme(panel.border = element_rect(fill = "transparent", color = 4,size = 2))  + 
       theme(plot.background = element_rect(fill = "gray86")) +
     theme(plot.background = element_rect(color = "black", size = 2)) +
  xlab("Grupos") + ylab("cobertura por metro cuadrado") + 
  theme(axis.title.x = element_text(size = 11, color = "black")) + theme(axis.title.y = element_text(size = 13, color = "black")) + theme(axis.text.x = element_text(angle = 0, vjust = 1, hjust = 1)) +   scale_x_continuous(n.breaks = 8)

## Gráfico de la densidad del cultivo y arvences por metro cuadrado 

ggplot(COBERTURA_AVENA_Y_MALEZAS, aes(x = PMUEST3AV, y = DENSIDADESCA3AV, fill = MUESTRASMC3AV)) + 
  geom_bar(stat = "identity", position = "dodge") +
   guides(fill = guide_legend(title = "Especie")) + 
    ggtitle("Densidad de avena y arvenses por metro cuadrado") + 
     theme(panel.border = element_rect(fill = "transparent", color = 4,size = 2))  + 
       theme(plot.background = element_rect(fill = "gray86")) +
     theme(plot.background = element_rect(color = "black", size = 2)) +
  xlab("Grupos") + ylab("densidad por metro cuadrado") + 
  theme(axis.title.x = element_text(size = 11, color = "black")) + theme(axis.title.y = element_text(size = 13, color = "black")) + theme(axis.text.x = element_text(angle = 0, vjust = 1, hjust = 1)) +   scale_x_continuous(n.breaks = 8)