ANÁLISIS MALEZAS Y PLAGAS EN CULTIVO DE CEBADA, AVENA Y FEIJOA
Grupo 6, manejos integrados
2023-2S
##
## Attaching package: 'ggplot2'## The following objects are masked from 'package:psych':
##
## %+%, alpha## ── 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##
## Attaching package: 'abdiv'
##
## The following objects are masked from 'package:psych':
##
## manhattan, minkowski, shannon## Loading required package: permute
## This is vegan 2.6-4## Loading required package: cluster## Loading required package: sp
## 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':
##
## manhattanCEBADA
Datos y gráficas del seguimiento de parametros de crecimiento del cultivo (muestreo 4)
DMTOTAL4 <- read_excel("PARAMETROSCEBADA02NOV.xlsx",
col_types = c("numeric", "numeric", "numeric",
"numeric", "numeric", "numeric",
"numeric"))
attach(DMTOTAL4)
# Altura de cebada:
psych::describe(ALTURAPROM2NOV)## vars n mean sd median trimmed mad min max range skew kurtosis se
## X1 1 8 88.62 15.03 84 88.62 8.15 71 122 51 1.14 0.29 5.31## vars n mean sd median trimmed mad min max range skew kurtosis se
## X1 1 8 41.25 12.85 37.5 41.25 6.67 30 71 41 1.41 0.7 4.54## vars n mean sd median trimmed mad min max range skew kurtosis se
## X1 1 8 85.62 39.64 88.5 85.62 24.46 0 134 134 -0.98 0.03 14.02## Gráficos de altura de cebada por metro cuadrado
par(mfrow=c(1,1))
ggplot(DMTOTAL4, aes(x = PUNTOMUES, y = ALTURAPROM2NOV, fill = MUESTREO)) +
geom_bar(stat = "identity", position = "dodge") +
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)## 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(DMTOTAL4, aes(x = PUNTOMUES, y = DENSIDAD.2NOV, 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(DMTOTAL4, aes(x = PUNTOMUES, y = TALLOSR2NOV, fill = MUESTREO)) +
geom_bar(stat = "identity", position = "dodge") +
guides(fill = guide_legend(title = "Muestreo")) +
ggtitle("Número de tallos reproductivos 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 tallos reprodcutivos 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 y 4)
cultivocebada4 <- read_excel("cebada 1.xlsx", col_types = c("numeric",
"numeric", "numeric", "numeric", "numeric",
"numeric", "numeric"))
attach(cultivocebada4)## The following objects are masked from DMTOTAL4:
##
## MUESTREO, PUNTOMUES## vars n mean sd median trimmed mad min max range skew kurtosis se
## X1 1 32 51.83 27.24 48.55 49.83 31.58 15.37 122 106.63 0.55 -0.66 4.82## vars n mean sd median trimmed mad min max range skew kurtosis se
## X1 1 32 16.95 16.08 10.35 14.46 7.12 2.5 71 68.5 1.5 1.81 2.84## vars n mean sd median trimmed mad min max range skew kurtosis se
## X1 1 32 92.47 32.93 93 93.12 25.95 0 159 159 -0.34 0.46 5.82## Gráficos de altura de cebada por metro cuadrado
par(mfrow=c(1,1))
ggplot(cultivocebada4, aes(x = PUNTOMUES, y = ALTURAPROMM2, fill = as.factor(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(cultivocebada4, aes(x = PUNTOMUES, y = DENSIDADM2, fill = as.factor(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(cultivocebada4, aes(x = PUNTOMUES, y = MACOLLASM2, fill = as.factor(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 cuarto muestreo (arvenses)
MALEZAS4 <- read_excel("MALEZAS CEBADA 4 NOV2.xlsx",
col_types = c("numeric", "text", "numeric",
"numeric", "numeric", "numeric",
"numeric", "numeric", "numeric",
"numeric", "numeric", "numeric"))
attach(MALEZAS4)## The following object is masked from cultivocebada4:
##
## MUESTREO## The following object is masked from DMTOTAL4:
##
## MUESTREO## vars n mean sd median trimmed mad min max range
## MUESTREO 1 18 4.00 0.00 4.00 4.00 0.00 4 4.00 0.00
## ESPECIESGL4* 2 18 9.50 5.34 9.50 9.50 6.67 1 18.00 17.00
## DENSIDAD4 3 18 35.39 107.44 4.00 11.06 5.93 0 460.00 460.00
## DENSIDADRELATIVA4 4 18 5.56 16.87 0.63 1.74 0.93 0 72.21 72.21
## FRECUENCIATOTAL4 5 18 0.10 0.12 0.06 0.08 0.09 0 0.44 0.44
## FRECUENCIARELATIVA4 6 18 5.56 6.71 3.51 4.71 5.20 0 24.56 24.56
## COBERTURATOTAL4 7 18 31.48 83.15 1.56 13.75 2.32 0 346.60 346.60
## DOMINANCIARELATIVA4 8 18 5.56 14.68 0.28 2.43 0.41 0 61.17 61.17
## DENSIDADPROMT4 9 18 1.11 3.36 0.13 0.35 0.19 0 14.38 14.38
## COBERTURAPROMT4 10 18 0.98 2.60 0.05 0.43 0.07 0 10.83 10.83
## DENSIDADMCUADRADO3 11 18 8.85 26.86 1.00 2.77 1.48 0 115.00 115.00
## IVI4 12 18 16.67 37.14 5.14 8.88 7.61 0 157.94 157.94
## skew kurtosis se
## MUESTREO NaN NaN 0.00
## ESPECIESGL4* 0.00 -1.40 1.26
## DENSIDAD4 3.41 10.56 25.32
## DENSIDADRELATIVA4 3.41 10.56 3.98
## FRECUENCIATOTAL4 1.43 1.31 0.03
## FRECUENCIARELATIVA4 1.44 1.33 1.58
## COBERTURATOTAL4 3.05 8.52 19.60
## DOMINANCIARELATIVA4 3.05 8.52 3.46
## DENSIDADPROMT4 3.41 10.56 0.79
## COBERTURAPROMT4 3.05 8.53 0.61
## DENSIDADMCUADRADO3 3.41 10.56 6.33
## IVI4 3.06 8.68 8.75## Resultados gráficos por especie
#Densidad por metro cuadrado de cada especie en los 8 metros cuadrados en total
ggplot(MALEZAS4, aes(x = ESPECIESGL4, y = DENSIDADMCUADRADO3, fill = ESPECIESGL4)) +
geom_bar(stat = "identity") +
theme(legend.key.size =unit(0.3, 'cm'))+
guides(fill = guide_legend(title = "ESPECIES")) +
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 = 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 8 metros cuadrados en total
ggplot(MALEZAS4, aes(x = ESPECIESGL4, y = DENSIDADRELATIVA4, fill = ESPECIESGL4)) +
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 8 metros cuadrados en total
ggplot(MALEZAS4, aes(x = ESPECIESGL4, y = FRECUENCIARELATIVA4, fill = ESPECIESGL4)) +
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(MALEZAS4, aes(x = ESPECIESGL4, y = DOMINANCIARELATIVA4, fill = ESPECIESGL4)) +
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 4 metros cuadrados en total
ggplot(MALEZAS4, aes(x = ESPECIESGL4, y = IVI4, fill = ESPECIESGL4)) +
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)
coberturas4 <- read_excel("COBERTYDENSI_NOV2CEBADA.xlsx",
col_types = c("numeric", "numeric", "text",
"numeric", "numeric"))
attach(coberturas4)
# densidad y cobertura del cultivo y el total de arvenses por metro cuadrado
psych::describe(coberturas4)## vars n mean sd median trimmed mad min max range skew
## MUESTREOO 1 16 3.00 0.00 3.0 3.00 0.00 3 3 0 NaN
## PMUEST3 2 16 4.50 2.37 4.5 4.50 2.97 1 8 7 0.00
## MUESTRASMC3* 3 16 1.50 0.52 1.5 1.50 0.74 1 2 1 0.00
## DENSIDADESCA3 4 16 82.62 52.85 82.5 82.21 59.30 0 171 171 -0.08
## COBERTURASCA3 5 16 182.23 142.97 132.5 173.69 161.60 4 480 476 0.55
## kurtosis se
## MUESTREOO NaN 0.00
## PMUEST3 -1.45 0.59
## MUESTRASMC3* -2.12 0.13
## DENSIDADESCA3 -1.23 13.21
## COBERTURASCA3 -0.96 35.74## Gráfico de la cobertura del cultivo y arvenses por metro cuadrado
ggplot(coberturas4, 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 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(coberturas4, 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 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 4
MUESTREOS
DMTOTAL4 <- read_excel("MALEZASGRUPALES.xlsx",
col_types = c("numeric", "text", "numeric",
"numeric", "numeric", "numeric",
"numeric", "numeric", "numeric",
"numeric", "numeric", "numeric"))
attach(DMTOTAL4)## The following object is masked from MALEZAS4:
##
## MUESTREO## The following object is masked from cultivocebada4:
##
## MUESTREO## The following object is masked from DMTOTAL4 (pos = 6):
##
## MUESTREO## Resultados gráficos por Muestreo
ggplot(DMTOTAL4, aes(x = ESPECIESGL, y = DENSIDADMCUADRADO, fill = as.factor(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(DMTOTAL4, aes(x = ESPECIESGL, y = DENSIDADRELATIVA, fill = as.factor(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(DMTOTAL4, aes(x = ESPECIESGL, y = FRECUENCIARELATIVA, fill = as.factor(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(DMTOTAL4, aes(x = ESPECIESGL, y = DOMINANCIARELATIVA, fill = as.factor(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(DMTOTAL4, aes(x = ESPECIESGL, y = IVI, fill = as.factor(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 4)
variables_avena2 <- read_excel("AVENA 2NOV.xlsx",
col_types = c("numeric", "numeric", "numeric",
"numeric", "numeric", "numeric"))
attach(variables_avena2)
# Estadisticos variables avena
psych::describe(variables_avena2)## vars n mean sd median trimmed mad min max range skew
## MUESTREOAV 1 8 4.00 0.00 4.00 4.00 0.00 4.0 4.00 0.00 NaN
## PUNTOMUESAV 2 8 4.50 2.45 4.50 4.50 2.97 1.0 8.00 7.00 0.00
## NTALLOSR 3 8 21.25 10.96 21.50 21.25 5.93 0.0 39.00 39.00 -0.36
## HMACOLLASAV 4 8 15.38 1.77 15.00 15.38 2.22 13.0 18.00 5.00 0.18
## ALTPROMAV 5 8 108.44 26.19 100.84 108.44 8.50 89.2 170.74 81.54 1.59
## DENSIAVE 6 8 56.25 40.43 59.00 56.25 43.74 0.0 108.00 108.00 -0.25
## kurtosis se
## MUESTREOAV NaN 0.00
## PUNTOMUESAV -1.65 0.87
## NTALLOSR -0.36 3.88
## HMACOLLASAV -1.71 0.62
## ALTPROMAV 1.07 9.26
## DENSIAVE -1.57 14.29## Gráficos de altura de cebada 0.75 metros cuadrados
par(mfrow=c(1,1))
ggplot(variables_avena2, aes(x = PUNTOMUESAV, y = ALTPROMAV, fill = MUESTREOAV)) +
geom_bar(stat = "identity", position = "dodge") +
guides(fill = guide_legend(title = "Muestreo")) +
ggtitle("Altura de avena por 0.75 metros cuadrados") +
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 0.75 metros cuadrados") +
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 0.75 metros cuadrados
ggplot(variables_avena2, aes(x = PUNTOMUESAV, y = DENSIAVE, fill = MUESTREOAV)) +
geom_bar(stat = "identity", position = "dodge") +
guides(fill = guide_legend(title = "Muestreo")) +
ggtitle("Densisdad de avena por 0.75 metros cuadrados") +
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 0.75 metros cuadrados") +
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 0.75 metros cuadrados
ggplot(variables_avena2, aes(x = PUNTOMUESAV, y = NTALLOSR, fill = MUESTREOAV)) +
geom_bar(stat = "identity", position = "dodge") +
guides(fill = guide_legend(title = "Muestreo")) +
ggtitle("Número de macollas de avena por 0.75 metros cuadrados") +
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 cuarto muestreo (arvenses)
MALEZAS_AVENA3 <- read_excel("MALEZAS_2NOV_AVENA.xlsx",
col_types = c("text", "numeric", "numeric",
"numeric", "numeric", "numeric",
"numeric", "numeric", "numeric",
"numeric"))
attach(MALEZAS_AVENA3)
# Estadisticos malezas en cultivo de avena:
psych::describe(MALEZAS_AVENA3)## vars n mean sd median trimmed mad min max range
## ESPECIESAV2NOV* 1 21 11.00 6.20 11.00 11.00 7.41 1 21.00 20.00
## DENSIDADAVT4 2 21 20.71 28.74 4.00 15.76 5.93 0 91.00 91.00
## DENSIDADRELATIVA4AV 3 21 4.77 6.61 0.90 3.63 1.33 0 20.90 20.90
## FRECUENCIATOTAL4AV 4 21 0.11 0.10 0.08 0.10 0.12 0 0.33 0.33
## FRECUENCIARELATIVA4AV 5 21 4.76 4.63 3.80 4.32 5.63 0 15.10 15.10
## COBERTURATOTAL4AV 6 21 14.51 24.95 2.00 9.06 2.97 0 96.00 96.00
## DOMINANCIARELATIVA4AV 7 21 4.76 8.18 0.70 2.97 1.04 0 31.50 31.50
## DENSIDADPROMT4AV 8 21 0.88 1.20 0.17 0.68 0.25 0 3.79 3.79
## COBERTURAPROMT4AV 9 21 0.62 1.04 0.09 0.40 0.13 0 4.00 4.00
## IVI4AV 10 21 14.29 17.47 5.70 11.62 8.45 0 54.20 54.20
## skew kurtosis se
## ESPECIESAV2NOV* 0.00 -1.37 1.35
## DENSIDADAVT4 1.20 -0.05 6.27
## DENSIDADRELATIVA4AV 1.20 -0.07 1.44
## FRECUENCIATOTAL4AV 0.54 -1.08 0.02
## FRECUENCIARELATIVA4AV 0.57 -0.98 1.01
## COBERTURATOTAL4AV 1.91 2.98 5.44
## DOMINANCIARELATIVA4AV 1.91 3.00 1.79
## DENSIDADPROMT4AV 1.15 -0.16 0.26
## COBERTURAPROMT4AV 1.87 2.86 0.23
## IVI4AV 1.14 -0.19 3.81## Resultados gráficos por especie
#Densidad 0.75 metros cuadrados de cada especie en los 6 metros cuadrados en total
ggplot(MALEZAS_AVENA3, aes(x = ESPECIESAV2NOV, y = DENSIDADAVT4, fill = ESPECIESAV2NOV)) +
geom_bar(stat = "identity") +
theme(legend.key.size =unit(0.3, 'cm'))+
guides(fill = guide_legend(title = "ESPECIES")) +
ggtitle("Densidad de especies 0.75 metros cuadrados") +
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 6 metros cuadrados en total
ggplot(MALEZAS_AVENA3, aes(x = ESPECIESAV2NOV, y = DENSIDADRELATIVA4AV, fill = ESPECIESAV2NOV)) +
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 24 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 6 metros cuadrados en total
ggplot(MALEZAS_AVENA3, aes(x = ESPECIESAV2NOV, y = FRECUENCIARELATIVA4AV, fill = ESPECIESAV2NOV)) +
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 24 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_AVENA3, aes(x = ESPECIESAV2NOV, y = DOMINANCIARELATIVA4AV, fill = ESPECIESAV2NOV)) +
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 24 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_AVENA3, aes(x = ESPECIESAV2NOV, y = IVI4AV, fill = ESPECIESAV2NOV)) +
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_MALEZAS2 <- read_excel("COBERT Y DENSI_2NOV_AVENA.xlsx",
col_types = c("numeric", "numeric", "text",
"numeric", "numeric"))
attach(COBERTURA_AVENA_Y_MALEZAS2)
# densidad y cobertura del cultivo y el total de arvenses por 0.75 metros cuadrados
psych::describe(COBERTURA_AVENA_Y_MALEZAS2)## vars n mean sd median trimmed mad min max range skew
## MUESTREOOAV 1 16 2.00 0.00 2.0 2.00 0.00 2 2 0 NaN
## PMUEST3AV 2 16 4.50 2.37 4.5 4.50 2.97 1 8 7 0.00
## MUESTRASMC3AV* 3 16 1.50 0.52 1.5 1.50 0.74 1 2 1 0.00
## COBERTURASCA3AV 4 16 145.49 117.62 174.5 144.84 163.83 0 300 300 -0.05
## DENSIAVENA 5 16 55.31 46.21 50.0 50.57 40.03 0 177 177 0.94
## kurtosis se
## MUESTREOOAV NaN 0.00
## PMUEST3AV -1.45 0.59
## MUESTRASMC3AV* -2.12 0.13
## COBERTURASCA3AV -1.86 29.41
## DENSIAVENA 0.56 11.55## Gráfico de la cobertura del cultivo y arvenses por metro cuadrado
ggplot(COBERTURA_AVENA_Y_MALEZAS2, 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 0.75 metros cuadrados") +
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 0.75 metros cuadrados") +
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_MALEZAS2, aes(x = PMUEST3AV, y = DENSIAVENA, fill = MUESTRASMC3AV)) +
geom_bar(stat = "identity", position = "dodge") +
guides(fill = guide_legend(title = "Especie")) +
ggtitle("Densidad de avena y arvenses 0.75 metros cuadrados") +
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 0.75 metros cuadrados") +
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)