Datos, gráficas e indices del segundo muestreo (arvenses)
library(psych)
library(readxl)
library(ggplot2)
##
## Attaching package: 'ggplot2'
## The following objects are masked from 'package:psych':
##
## %+%, alpha
library(lattice)
library(tidyverse)
## -- Attaching packages --------------------------------------- tidyverse 1.3.1 --
## v tibble 3.1.6 v dplyr 1.0.8
## v tidyr 1.2.0 v stringr 1.4.0
## v readr 2.1.2 v forcats 0.5.1
## v purrr 0.3.4
## -- Conflicts ------------------------------------------ tidyverse_conflicts() --
## x ggplot2::%+%() masks psych::%+%()
## x ggplot2::alpha() masks psych::alpha()
## x dplyr::filter() masks stats::filter()
## x dplyr::lag() masks stats::lag()
library(dplyr)
library(abdiv)
##
## Attaching package: 'abdiv'
## The following objects are masked from 'package:psych':
##
## manhattan, minkowski, shannon
## Loading required package: permute
## This is vegan 2.6-4
library(codyn)
library(divo)
## 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':
##
## manhattan
DMTOTAL2 <- read_excel("C:\\Users\\andre\\OneDrive\\Escritorio\\SEGUNDO CÓDIGO MIP.xlsx",
col_types = c("numeric", "text", "numeric",
"numeric", "numeric", "numeric",
"numeric", "numeric", "numeric",
"numeric", "numeric", "numeric"))
attach(DMTOTAL2)
# Densidades:
psych::describe(DENSIDADMT2)
## vars n mean sd median trimmed mad min max range skew kurtosis se
## X1 1 19 66.11 109.15 12 54.35 11.86 0 332 332 1.5 0.59 25.04
describe.by(DENSIDADMT2, group = ESPECIESGL2)
## 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 12 NA 12 12 0 12 12 0 NA NA NA
## ------------------------------------------------------------
## group: Brassica sp
## 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: Chenopodium petiolare Kunth
## 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: Cotula australis
## 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: Fuertesimalva sp.
## vars n mean sd median trimmed mad min max range skew kurtosis se
## X1 1 1 6 NA 6 6 0 6 6 0 NA NA NA
## ------------------------------------------------------------
## group: Fumaria capreolata
## vars n mean sd median trimmed mad min max range skew kurtosis se
## X1 1 1 20 NA 20 20 0 20 20 0 NA NA NA
## ------------------------------------------------------------
## group: Galinsoga
## Ciliata
## 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: Lepidium bipinnatifidum
## vars n mean sd median trimmed mad min max range skew kurtosis se
## X1 1 1 10 NA 10 10 0 10 10 0 NA NA NA
## ------------------------------------------------------------
## group: Lolium multiflorum
## vars n mean sd median trimmed mad min max range skew kurtosis se
## X1 1 1 39 NA 39 39 0 39 39 0 NA NA NA
## ------------------------------------------------------------
## group: Malva sylvestris
## vars n mean sd median trimmed mad min max range skew kurtosis se
## X1 1 1 6 NA 6 6 0 6 6 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: Poace sp
## 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: Polygonum aviculare
## vars n mean sd median trimmed mad min max range skew kurtosis se
## X1 1 1 278 NA 278 278 0 278 278 0 NA NA NA
## ------------------------------------------------------------
## group: Polygonum segetum
## vars n mean sd median trimmed mad min max range skew kurtosis se
## X1 1 1 332 NA 332 332 0 332 332 0 NA NA NA
## ------------------------------------------------------------
## group: Polygonum sp
## vars n mean sd median trimmed mad min max range skew kurtosis se
## X1 1 1 290 NA 290 290 0 290 290 0 NA NA NA
## ------------------------------------------------------------
## group: Raphanus
## vars n mean sd median trimmed mad min max range skew kurtosis se
## X1 1 1 116 NA 116 116 0 116 116 0 NA NA NA
## ------------------------------------------------------------
## group: Rumex
## 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: Sonchus sp.
## 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: Veronica persica
## vars n mean sd median trimmed mad min max range skew kurtosis se
## X1 1 1 98 NA 98 98 0 98 98 0 NA NA NA
# densidad relativa:
psych::describe(DENSIDADRELATIVA2)
## vars n mean sd median trimmed mad min max range skew kurtosis se
## X1 1 19 5.26 8.69 0.96 4.33 0.95 0 26.43 26.43 1.5 0.59 1.99
describe.by(DENSIDADRELATIVA2, group = ESPECIESGL2)
## 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.96 NA 0.96 0.96 0 0.96 0.96 0 NA NA NA
## ------------------------------------------------------------
## group: Brassica sp
## vars n mean sd median trimmed mad min max range skew kurtosis se
## X1 1 1 0.56 NA 0.56 0.56 0 0.56 0.56 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.32 NA 0.32 0.32 0 0.32 0.32 0 NA NA NA
## ------------------------------------------------------------
## group: Cotula australis
## vars n mean sd median trimmed mad min max range skew kurtosis se
## X1 1 1 1.35 NA 1.35 1.35 0 1.35 1.35 0 NA NA NA
## ------------------------------------------------------------
## group: Fuertesimalva sp.
## vars n mean sd median trimmed mad min max range skew kurtosis se
## X1 1 1 0.48 NA 0.48 0.48 0 0.48 0.48 0 NA NA NA
## ------------------------------------------------------------
## group: Fumaria capreolata
## vars n mean sd median trimmed mad min max range skew kurtosis se
## X1 1 1 1.59 NA 1.59 1.59 0 1.59 1.59 0 NA NA NA
## ------------------------------------------------------------
## group: Galinsoga
## Ciliata
## vars n mean sd median trimmed mad min max range skew kurtosis se
## X1 1 1 0.08 NA 0.08 0.08 0 0.08 0.08 0 NA NA NA
## ------------------------------------------------------------
## group: Lepidium bipinnatifidum
## vars n mean sd median trimmed mad min max range skew kurtosis se
## X1 1 1 0.8 NA 0.8 0.8 0 0.8 0.8 0 NA NA NA
## ------------------------------------------------------------
## group: Lolium multiflorum
## vars n mean sd median trimmed mad min max range skew kurtosis se
## X1 1 1 3.11 NA 3.11 3.11 0 3.11 3.11 0 NA NA NA
## ------------------------------------------------------------
## group: Malva sylvestris
## vars n mean sd median trimmed mad min max range skew kurtosis se
## X1 1 1 0.48 NA 0.48 0.48 0 0.48 0.48 0 NA NA NA
## ------------------------------------------------------------
## group: Oxalis
## vars n mean sd median trimmed mad min max range skew kurtosis se
## X1 1 1 0.16 NA 0.16 0.16 0 0.16 0.16 0 NA NA NA
## ------------------------------------------------------------
## group: Poace sp
## vars n mean sd median trimmed mad min max range skew kurtosis se
## X1 1 1 0.32 NA 0.32 0.32 0 0.32 0.32 0 NA NA NA
## ------------------------------------------------------------
## group: Polygonum aviculare
## vars n mean sd median trimmed mad min max range skew kurtosis se
## X1 1 1 22.13 NA 22.13 22.13 0 22.13 22.13 0 NA NA NA
## ------------------------------------------------------------
## group: Polygonum segetum
## vars n mean sd median trimmed mad min max range skew kurtosis se
## X1 1 1 26.43 NA 26.43 26.43 0 26.43 26.43 0 NA NA NA
## ------------------------------------------------------------
## group: Polygonum sp
## vars n mean sd median trimmed mad min max range skew kurtosis se
## X1 1 1 23.09 NA 23.09 23.09 0 23.09 23.09 0 NA NA NA
## ------------------------------------------------------------
## group: Raphanus
## vars n mean sd median trimmed mad min max range skew kurtosis se
## X1 1 1 9.24 NA 9.24 9.24 0 9.24 9.24 0 NA NA NA
## ------------------------------------------------------------
## group: Rumex
## 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: Sonchus sp.
## vars n mean sd median trimmed mad min max range skew kurtosis se
## X1 1 1 1.11 NA 1.11 1.11 0 1.11 1.11 0 NA NA NA
## ------------------------------------------------------------
## group: Veronica persica
## vars n mean sd median trimmed mad min max range skew kurtosis se
## X1 1 1 7.8 NA 7.8 7.8 0 7.8 7.8 0 NA NA NA
# frecuencia:
psych::describe(FRECUENCIATOTAL2)
## vars n mean sd median trimmed mad min max range skew kurtosis se
## X1 1 19 0.14 0.15 0.09 0.13 0.09 0 0.56 0.56 1.29 0.85 0.03
describe.by(FRECUENCIATOTAL2, group = ESPECIESGL2)
## 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.06 NA 0.06 0.06 0 0.06 0.06 0 NA NA NA
## ------------------------------------------------------------
## group: Brassica sp
## vars n mean sd median trimmed mad min max range skew kurtosis se
## X1 1 1 0.09 NA 0.09 0.09 0 0.09 0.09 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.06 NA 0.06 0.06 0 0.06 0.06 0 NA NA NA
## ------------------------------------------------------------
## group: Cotula australis
## 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: Fuertesimalva 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: Fumaria capreolata
## 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: Galinsoga
## Ciliata
## vars n mean sd median trimmed mad min max range skew kurtosis se
## X1 1 1 0.03 NA 0.03 0.03 0 0.03 0.03 0 NA NA NA
## ------------------------------------------------------------
## group: Lepidium bipinnatifidum
## vars n mean sd median trimmed mad min max range skew kurtosis se
## X1 1 1 0.03 NA 0.03 0.03 0 0.03 0.03 0 NA NA NA
## ------------------------------------------------------------
## group: Lolium multiflorum
## vars n mean sd median trimmed mad min max range skew kurtosis se
## X1 1 1 0.09 NA 0.09 0.09 0 0.09 0.09 0 NA NA NA
## ------------------------------------------------------------
## group: Malva sylvestris
## vars n mean sd median trimmed mad min max range skew kurtosis se
## X1 1 1 0.03 NA 0.03 0.03 0 0.03 0.03 0 NA NA NA
## ------------------------------------------------------------
## group: Oxalis
## vars n mean sd median trimmed mad min max range skew kurtosis se
## X1 1 1 0.03 NA 0.03 0.03 0 0.03 0.03 0 NA NA NA
## ------------------------------------------------------------
## group: Poace sp
## vars n mean sd median trimmed mad min max range skew kurtosis se
## X1 1 1 0.03 NA 0.03 0.03 0 0.03 0.03 0 NA NA NA
## ------------------------------------------------------------
## group: Polygonum aviculare
## vars n mean sd median trimmed mad min max range skew kurtosis se
## X1 1 1 0.34 NA 0.34 0.34 0 0.34 0.34 0 NA NA NA
## ------------------------------------------------------------
## group: Polygonum segetum
## vars n mean sd median trimmed mad min max range skew kurtosis se
## X1 1 1 0.31 NA 0.31 0.31 0 0.31 0.31 0 NA NA NA
## ------------------------------------------------------------
## group: Polygonum sp
## vars n mean sd median trimmed mad min max range skew kurtosis se
## X1 1 1 0.16 NA 0.16 0.16 0 0.16 0.16 0 NA NA NA
## ------------------------------------------------------------
## group: Raphanus
## vars n mean sd median trimmed mad min max range skew kurtosis se
## X1 1 1 0.56 NA 0.56 0.56 0 0.56 0.56 0 NA NA NA
## ------------------------------------------------------------
## group: Rumex
## 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: Sonchus sp.
## 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: Veronica persica
## vars n mean sd median trimmed mad min max range skew kurtosis se
## X1 1 1 0.34 NA 0.34 0.34 0 0.34 0.34 0 NA NA NA
# frecuencia relativa:
psych::describe(FRECUENCIARELATIVA2)
## vars n mean sd median trimmed mad min max range skew kurtosis se
## X1 1 19 5.26 5.4 3.41 4.68 3.37 0 20.45 20.45 1.31 0.88 1.24
describe.by(FRECUENCIARELATIVA2, group = ESPECIESGL2)
## 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 2.27 NA 2.27 2.27 0 2.27 2.27 0 NA NA NA
## ------------------------------------------------------------
## group: Brassica sp
## vars n mean sd median trimmed mad min max range skew kurtosis se
## X1 1 1 3.41 NA 3.41 3.41 0 3.41 3.41 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.27 NA 2.27 2.27 0 2.27 2.27 0 NA NA NA
## ------------------------------------------------------------
## group: Cotula australis
## vars n mean sd median trimmed mad min max range skew kurtosis se
## X1 1 1 4.55 NA 4.55 4.55 0 4.55 4.55 0 NA NA NA
## ------------------------------------------------------------
## group: Fuertesimalva sp.
## vars n mean sd median trimmed mad min max range skew kurtosis se
## X1 1 1 2.27 NA 2.27 2.27 0 2.27 2.27 0 NA NA NA
## ------------------------------------------------------------
## group: Fumaria capreolata
## vars n mean sd median trimmed mad min max range skew kurtosis se
## X1 1 1 6.82 NA 6.82 6.82 0 6.82 6.82 0 NA NA NA
## ------------------------------------------------------------
## group: Galinsoga
## Ciliata
## vars n mean sd median trimmed mad min max range skew kurtosis se
## X1 1 1 1.14 NA 1.14 1.14 0 1.14 1.14 0 NA NA NA
## ------------------------------------------------------------
## group: Lepidium bipinnatifidum
## vars n mean sd median trimmed mad min max range skew kurtosis se
## X1 1 1 1.14 NA 1.14 1.14 0 1.14 1.14 0 NA NA NA
## ------------------------------------------------------------
## group: Lolium multiflorum
## vars n mean sd median trimmed mad min max range skew kurtosis se
## X1 1 1 3.41 NA 3.41 3.41 0 3.41 3.41 0 NA NA NA
## ------------------------------------------------------------
## group: Malva sylvestris
## vars n mean sd median trimmed mad min max range skew kurtosis se
## X1 1 1 1.14 NA 1.14 1.14 0 1.14 1.14 0 NA NA NA
## ------------------------------------------------------------
## group: Oxalis
## vars n mean sd median trimmed mad min max range skew kurtosis se
## X1 1 1 1.14 NA 1.14 1.14 0 1.14 1.14 0 NA NA NA
## ------------------------------------------------------------
## group: Poace sp
## vars n mean sd median trimmed mad min max range skew kurtosis se
## X1 1 1 1.14 NA 1.14 1.14 0 1.14 1.14 0 NA NA NA
## ------------------------------------------------------------
## group: Polygonum aviculare
## vars n mean sd median trimmed mad min max range skew kurtosis se
## X1 1 1 12.5 NA 12.5 12.5 0 12.5 12.5 0 NA NA NA
## ------------------------------------------------------------
## group: Polygonum segetum
## vars n mean sd median trimmed mad min max range skew kurtosis se
## X1 1 1 11.36 NA 11.36 11.36 0 11.36 11.36 0 NA NA NA
## ------------------------------------------------------------
## group: Polygonum sp
## vars n mean sd median trimmed mad min max range skew kurtosis se
## X1 1 1 5.68 NA 5.68 5.68 0 5.68 5.68 0 NA NA NA
## ------------------------------------------------------------
## group: Raphanus
## vars n mean sd median trimmed mad min max range skew kurtosis se
## X1 1 1 20.45 NA 20.45 20.45 0 20.45 20.45 0 NA NA NA
## ------------------------------------------------------------
## group: Rumex
## 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: Sonchus sp.
## vars n mean sd median trimmed mad min max range skew kurtosis se
## X1 1 1 6.82 NA 6.82 6.82 0 6.82 6.82 0 NA NA NA
## ------------------------------------------------------------
## group: Veronica persica
## vars n mean sd median trimmed mad min max range skew kurtosis se
## X1 1 1 12.5 NA 12.5 12.5 0 12.5 12.5 0 NA NA NA
# Dominancia:
psych::describe(COBERTURATOTAL2)
## vars n mean sd median trimmed mad min max range skew kurtosis se
## X1 1 19 19.56 28.16 7 16.09 7.41 0 98.2 98.2 1.39 0.77 6.46
describe.by(COBERTURATOTAL2, group = ESPECIESGL2)
## 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 4 NA 4 4 0 4 4 0 NA NA NA
## ------------------------------------------------------------
## group: Brassica sp
## vars n mean sd median trimmed mad min max range skew kurtosis se
## X1 1 1 2.7 NA 2.7 2.7 0 2.7 2.7 0 NA NA NA
## ------------------------------------------------------------
## group: Chenopodium petiolare Kunth
## 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: Cotula australis
## 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: Fuertesimalva sp.
## 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: Fumaria capreolata
## vars n mean sd median trimmed mad min max range skew kurtosis se
## X1 1 1 10.3 NA 10.3 10.3 0 10.3 10.3 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 8 NA 8 8 0 8 8 0 NA NA NA
## ------------------------------------------------------------
## group: Lolium multiflorum
## 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: 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 0 NA 0 0 0 0 0 0 NA NA NA
## ------------------------------------------------------------
## group: Poace 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: Polygonum aviculare
## vars n mean sd median trimmed mad min max range skew kurtosis se
## X1 1 1 58 NA 58 58 0 58 58 0 NA NA NA
## ------------------------------------------------------------
## group: Polygonum segetum
## vars n mean sd median trimmed mad min max range skew kurtosis se
## X1 1 1 58 NA 58 58 0 58 58 0 NA NA NA
## ------------------------------------------------------------
## group: Polygonum sp
## vars n mean sd median trimmed mad min max range skew kurtosis se
## X1 1 1 48 NA 48 48 0 48 48 0 NA NA NA
## ------------------------------------------------------------
## group: Raphanus
## vars n mean sd median trimmed mad min max range skew kurtosis se
## X1 1 1 98.2 NA 98.2 98.2 0 98.2 98.2 0 NA NA NA
## ------------------------------------------------------------
## group: Rumex
## 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: Sonchus sp.
## 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: Veronica persica
## vars n mean sd median trimmed mad min max range skew kurtosis se
## X1 1 1 49.5 NA 49.5 49.5 0 49.5 49.5 0 NA NA NA
#Dominancia relativa:
psych::describe(DOMINANCIARELATIVA2)
## vars n mean sd median trimmed mad min max range skew kurtosis se
## X1 1 19 5.26 7.58 1.88 4.33 1.99 0 26.42 26.42 1.4 0.77 1.74
describe.by(DOMINANCIARELATIVA2, group = ESPECIESGL2)
## 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 1.08 NA 1.08 1.08 0 1.08 1.08 0 NA NA NA
## ------------------------------------------------------------
## group: Brassica sp
## vars n mean sd median trimmed mad min max range skew kurtosis se
## X1 1 1 0.73 NA 0.73 0.73 0 0.73 0.73 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.15 NA 2.15 2.15 0 2.15 2.15 0 NA NA NA
## ------------------------------------------------------------
## group: Cotula australis
## vars n mean sd median trimmed mad min max range skew kurtosis se
## X1 1 1 1.88 NA 1.88 1.88 0 1.88 1.88 0 NA NA NA
## ------------------------------------------------------------
## group: Fuertesimalva sp.
## vars n mean sd median trimmed mad min max range skew kurtosis se
## X1 1 1 1.35 NA 1.35 1.35 0 1.35 1.35 0 NA NA NA
## ------------------------------------------------------------
## group: Fumaria capreolata
## vars n mean sd median trimmed mad min max range skew kurtosis se
## X1 1 1 2.77 NA 2.77 2.77 0 2.77 2.77 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 2.15 NA 2.15 2.15 0 2.15 2.15 0 NA NA NA
## ------------------------------------------------------------
## group: Lolium multiflorum
## 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: Malva sylvestris
## vars n mean sd median trimmed mad min max range skew kurtosis se
## X1 1 1 0.54 NA 0.54 0.54 0 0.54 0.54 0 NA NA NA
## ------------------------------------------------------------
## group: Oxalis
## 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: Poace 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: Polygonum aviculare
## vars n mean sd median trimmed mad min max range skew kurtosis se
## X1 1 1 15.6 NA 15.6 15.6 0 15.6 15.6 0 NA NA NA
## ------------------------------------------------------------
## group: Polygonum segetum
## vars n mean sd median trimmed mad min max range skew kurtosis se
## X1 1 1 15.6 NA 15.6 15.6 0 15.6 15.6 0 NA NA NA
## ------------------------------------------------------------
## group: Polygonum sp
## vars n mean sd median trimmed mad min max range skew kurtosis se
## X1 1 1 12.91 NA 12.91 12.91 0 12.91 12.91 0 NA NA NA
## ------------------------------------------------------------
## group: Raphanus
## vars n mean sd median trimmed mad min max range skew kurtosis se
## X1 1 1 26.42 NA 26.42 26.42 0 26.42 26.42 0 NA NA NA
## ------------------------------------------------------------
## group: Rumex
## 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: Sonchus sp.
## vars n mean sd median trimmed mad min max range skew kurtosis se
## X1 1 1 1.35 NA 1.35 1.35 0 1.35 1.35 0 NA NA NA
## ------------------------------------------------------------
## group: Veronica persica
## vars n mean sd median trimmed mad min max range skew kurtosis se
## X1 1 1 13.32 NA 13.32 13.32 0 13.32 13.32 0 NA NA NA
## Resultados gráficos por especie
#Densidad por metro cuadrado de cada especie en los 8 metros cuadrados en total
ggplot(DMTOTAL2, aes(x = ESPECIESGL2, y = DENSIDADMCUADRADO2, fill = ESPECIESGL2)) +
geom_bar(stat = "identity") +
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(DMTOTAL2, aes(x = ESPECIESGL2, y = DENSIDADRELATIVA2, fill = ESPECIESGL2)) +
geom_bar(stat = "identity") +
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(DMTOTAL2, aes(x = ESPECIESGL2, y = FRECUENCIARELATIVA2, fill = ESPECIESGL2)) +
geom_bar(stat = "identity") +
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 8 metros cuadrados en total
ggplot(DMTOTAL2, aes(x = ESPECIESGL2, y = DOMINANCIARELATIVA2, fill = ESPECIESGL2)) +
geom_bar(stat = "identity") +
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(DMTOTAL2, aes(x = ESPECIESGL2, y = IVI2, fill = ESPECIESGL2)) +
geom_bar(stat = "identity") +
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))
#Índice de Shannon
diversity(DMTOTAL2$DENSIDADMT2, index = "shannon")
## [1] 1.940511
#Índice de Simpson
diversity(DMTOTAL2$DENSIDADRELATIVA2, index = "simpson")
## [1] 0.8114902
#Índice de Pielou
pielou_e(DMTOTAL2$DENSIDADRELATIVA2)
## [1] 0.6715454
#Índice de Margalef
margalef(DMTOTAL2$DENSIDADRELATIVA2)
## [1] 3.691423
soronsen <- li(DMTOTAL2, CI = 0.95, resample = 100, graph = FALSE, csv_output = FALSE,
PlugIn = FALSE, size = 1, saveBootstrap = FALSE)
## Warning in li(DMTOTAL2, CI = 0.95, resample = 100, graph = FALSE, csv_output =
## FALSE, : NAs introducidos por coerción
## $Mean
## MUESTREO ESPECIESGL2 DENSIDADMT2 DENSIDADRELATIVA2
## MUESTREO 1.0000000 0.0000000 0.0000000 0.0000000
## ESPECIESGL2 1.0000000 1.0000000 0.0000000 0.0000000
## DENSIDADMT2 0.9729730 0.9729730 1.0000000 0.0000000
## DENSIDADRELATIVA2 0.6428571 0.6428571 0.6666667 1.0000000
## FRECUENCIATOTAL2 0.0000000 0.0000000 0.0000000 0.0000000
## FRECUENCIARELATIVA2 0.9729730 0.9729730 1.0000000 0.6666667
## COBERTURATOTAL2 0.8823529 0.8823529 0.9090909 0.7500000
## DOMINANCIARELATIVA2 0.8125000 0.8125000 0.8387097 0.8181818
## DENSIDADPROMT2 0.4800000 0.4800000 0.5000000 0.8000000
## COBERTURAPROMT2 0.4166667 0.4166667 0.4347826 0.7142857
## DENSIDADMCUADRADO2 0.7333333 0.7333333 0.7586207 0.9000000
## IVI2 0.9729730 0.9729730 1.0000000 0.6666667
## FRECUENCIATOTAL2 FRECUENCIARELATIVA2 COBERTURATOTAL2
## MUESTREO 0 0.0000000 0.0000000
## ESPECIESGL2 0 0.0000000 0.0000000
## DENSIDADMT2 0 0.0000000 0.0000000
## DENSIDADRELATIVA2 0 0.0000000 0.0000000
## FRECUENCIATOTAL2 1 0.0000000 0.0000000
## FRECUENCIARELATIVA2 0 1.0000000 0.0000000
## COBERTURATOTAL2 0 0.9090909 1.0000000
## DOMINANCIARELATIVA2 0 0.8387097 0.9285714
## DENSIDADPROMT2 0 0.5000000 0.5714286
## COBERTURAPROMT2 0 0.4347826 0.5000000
## DENSIDADMCUADRADO2 0 0.7586207 0.8461538
## IVI2 0 1.0000000 0.9090909
## DOMINANCIARELATIVA2 DENSIDADPROMT2 COBERTURAPROMT2
## MUESTREO 0.0000000 0.0000000 0.0000000
## ESPECIESGL2 0.0000000 0.0000000 0.0000000
## DENSIDADMT2 0.0000000 0.0000000 0.0000000
## DENSIDADRELATIVA2 0.0000000 0.0000000 0.0000000
## FRECUENCIATOTAL2 0.0000000 0.0000000 0.0000000
## FRECUENCIARELATIVA2 0.0000000 0.0000000 0.0000000
## COBERTURATOTAL2 0.0000000 0.0000000 0.0000000
## DOMINANCIARELATIVA2 1.0000000 0.0000000 0.0000000
## DENSIDADPROMT2 0.6315789 1.0000000 0.0000000
## COBERTURAPROMT2 0.5555556 0.9090909 1.0000000
## DENSIDADMCUADRADO2 0.9166667 0.7058824 0.6250000
## IVI2 0.8387097 0.5000000 0.4347826
## DENSIDADMCUADRADO2 IVI2
## MUESTREO 0.0000000 0
## ESPECIESGL2 0.0000000 0
## DENSIDADMT2 0.0000000 0
## DENSIDADRELATIVA2 0.0000000 0
## FRECUENCIATOTAL2 0.0000000 0
## FRECUENCIARELATIVA2 0.0000000 0
## COBERTURATOTAL2 0.0000000 0
## DOMINANCIARELATIVA2 0.0000000 0
## DENSIDADPROMT2 0.0000000 0
## COBERTURAPROMT2 0.0000000 0
## DENSIDADMCUADRADO2 1.0000000 0
## IVI2 0.7586207 1
##
## $Lower.Quantile
## MUESTREO ESPECIESGL2 DENSIDADMT2 DENSIDADRELATIVA2
## MUESTREO 1.0000000 0.0000000 0.0000000 0.0000000
## ESPECIESGL2 1.0000000 1.0000000 0.0000000 0.0000000
## DENSIDADMT2 0.9729730 0.9729730 1.0000000 0.0000000
## DENSIDADRELATIVA2 0.6428571 0.6428571 0.6666667 1.0000000
## FRECUENCIATOTAL2 0.0000000 0.0000000 0.0000000 0.0000000
## FRECUENCIARELATIVA2 0.9729730 0.9729730 1.0000000 0.6666667
## COBERTURATOTAL2 0.8823529 0.8823529 0.9090909 0.7500000
## DOMINANCIARELATIVA2 0.8125000 0.8125000 0.8387097 0.8181818
## DENSIDADPROMT2 0.4800000 0.4800000 0.5000000 0.8000000
## COBERTURAPROMT2 0.4166667 0.4166667 0.4347826 0.7142857
## DENSIDADMCUADRADO2 0.7333333 0.7333333 0.7586207 0.9000000
## IVI2 0.9729730 0.9729730 1.0000000 0.6666667
## FRECUENCIATOTAL2 FRECUENCIARELATIVA2 COBERTURATOTAL2
## MUESTREO 0 0.0000000 0.0000000
## ESPECIESGL2 0 0.0000000 0.0000000
## DENSIDADMT2 0 0.0000000 0.0000000
## DENSIDADRELATIVA2 0 0.0000000 0.0000000
## FRECUENCIATOTAL2 1 0.0000000 0.0000000
## FRECUENCIARELATIVA2 0 1.0000000 0.0000000
## COBERTURATOTAL2 0 0.9090909 1.0000000
## DOMINANCIARELATIVA2 0 0.8387097 0.9285714
## DENSIDADPROMT2 0 0.5000000 0.5714286
## COBERTURAPROMT2 0 0.4347826 0.5000000
## DENSIDADMCUADRADO2 0 0.7586207 0.8461538
## IVI2 0 1.0000000 0.9090909
## DOMINANCIARELATIVA2 DENSIDADPROMT2 COBERTURAPROMT2
## MUESTREO 0.0000000 0.0000000 0.0000000
## ESPECIESGL2 0.0000000 0.0000000 0.0000000
## DENSIDADMT2 0.0000000 0.0000000 0.0000000
## DENSIDADRELATIVA2 0.0000000 0.0000000 0.0000000
## FRECUENCIATOTAL2 0.0000000 0.0000000 0.0000000
## FRECUENCIARELATIVA2 0.0000000 0.0000000 0.0000000
## COBERTURATOTAL2 0.0000000 0.0000000 0.0000000
## DOMINANCIARELATIVA2 1.0000000 0.0000000 0.0000000
## DENSIDADPROMT2 0.6315789 1.0000000 0.0000000
## COBERTURAPROMT2 0.5555556 0.9090909 1.0000000
## DENSIDADMCUADRADO2 0.9166667 0.7058824 0.6250000
## IVI2 0.8387097 0.5000000 0.4347826
## DENSIDADMCUADRADO2 IVI2
## MUESTREO 0.0000000 0
## ESPECIESGL2 0.0000000 0
## DENSIDADMT2 0.0000000 0
## DENSIDADRELATIVA2 0.0000000 0
## FRECUENCIATOTAL2 0.0000000 0
## FRECUENCIARELATIVA2 0.0000000 0
## COBERTURATOTAL2 0.0000000 0
## DOMINANCIARELATIVA2 0.0000000 0
## DENSIDADPROMT2 0.0000000 0
## COBERTURAPROMT2 0.0000000 0
## DENSIDADMCUADRADO2 1.0000000 0
## IVI2 0.7586207 1
##
## $Upper.Quantile
## MUESTREO ESPECIESGL2 DENSIDADMT2 DENSIDADRELATIVA2
## MUESTREO 1.0000000 0.0000000 0.0000000 0.0000000
## ESPECIESGL2 1.0000000 1.0000000 0.0000000 0.0000000
## DENSIDADMT2 0.9729730 0.9729730 1.0000000 0.0000000
## DENSIDADRELATIVA2 0.6428571 0.6428571 0.6666667 1.0000000
## FRECUENCIATOTAL2 0.0000000 0.0000000 0.0000000 0.0000000
## FRECUENCIARELATIVA2 0.9729730 0.9729730 1.0000000 0.6666667
## COBERTURATOTAL2 0.8823529 0.8823529 0.9090909 0.7500000
## DOMINANCIARELATIVA2 0.8125000 0.8125000 0.8387097 0.8181818
## DENSIDADPROMT2 0.4800000 0.4800000 0.5000000 0.8000000
## COBERTURAPROMT2 0.4166667 0.4166667 0.4347826 0.7142857
## DENSIDADMCUADRADO2 0.7333333 0.7333333 0.7586207 0.9000000
## IVI2 0.9729730 0.9729730 1.0000000 0.6666667
## FRECUENCIATOTAL2 FRECUENCIARELATIVA2 COBERTURATOTAL2
## MUESTREO 0 0.0000000 0.0000000
## ESPECIESGL2 0 0.0000000 0.0000000
## DENSIDADMT2 0 0.0000000 0.0000000
## DENSIDADRELATIVA2 0 0.0000000 0.0000000
## FRECUENCIATOTAL2 1 0.0000000 0.0000000
## FRECUENCIARELATIVA2 0 1.0000000 0.0000000
## COBERTURATOTAL2 0 0.9090909 1.0000000
## DOMINANCIARELATIVA2 0 0.8387097 0.9285714
## DENSIDADPROMT2 0 0.5000000 0.5714286
## COBERTURAPROMT2 0 0.4347826 0.5000000
## DENSIDADMCUADRADO2 0 0.7586207 0.8461538
## IVI2 0 1.0000000 0.9090909
## DOMINANCIARELATIVA2 DENSIDADPROMT2 COBERTURAPROMT2
## MUESTREO 0.0000000 0.0000000 0.0000000
## ESPECIESGL2 0.0000000 0.0000000 0.0000000
## DENSIDADMT2 0.0000000 0.0000000 0.0000000
## DENSIDADRELATIVA2 0.0000000 0.0000000 0.0000000
## FRECUENCIATOTAL2 0.0000000 0.0000000 0.0000000
## FRECUENCIARELATIVA2 0.0000000 0.0000000 0.0000000
## COBERTURATOTAL2 0.0000000 0.0000000 0.0000000
## DOMINANCIARELATIVA2 1.0000000 0.0000000 0.0000000
## DENSIDADPROMT2 0.6315789 1.0000000 0.0000000
## COBERTURAPROMT2 0.5555556 0.9090909 1.0000000
## DENSIDADMCUADRADO2 0.9166667 0.7058824 0.6250000
## IVI2 0.8387097 0.5000000 0.4347826
## DENSIDADMCUADRADO2 IVI2
## MUESTREO 0.0000000 0
## ESPECIESGL2 0.0000000 0
## DENSIDADMT2 0.0000000 0
## DENSIDADRELATIVA2 0.0000000 0
## FRECUENCIATOTAL2 0.0000000 0
## FRECUENCIARELATIVA2 0.0000000 0
## COBERTURATOTAL2 0.0000000 0
## DOMINANCIARELATIVA2 0.0000000 0
## DENSIDADPROMT2 0.0000000 0
## COBERTURAPROMT2 0.0000000 0
## DENSIDADMCUADRADO2 1.0000000 0
## IVI2 0.7586207 1
vegdist(DMTOTAL2$DENSIDADRELATIVA2, method = "jaccard", diag = F)
## Warning in vegdist(DMTOTAL2$DENSIDADRELATIVA2, method = "jaccard", diag = F): you have empty rows: their dissimilarities may be
## meaningless in method "jaccard"
## 1 2 3 4 5 6 7
## 2 0.85714286
## 3 0.49549550 0.92792793
## 4 0.97469498 0.99638500 0.94984184
## 5 0.97881196 0.99697314 0.95800227 0.16269391
## 6 0.97574708 0.99653530 0.95192724 0.04157644 0.12637155
## 7 0.93939394 0.99134199 0.87987013 0.58246724 0.65039728 0.59982676
## 8 0.92820513 0.98974359 0.85769231 0.64753728 0.70488082 0.66219142 0.15584416
## 9 0.42857143 0.75000000 0.71171171 0.98553999 0.98789255 0.98614119 0.96536797
## 10 0.42857143 0.75000000 0.71171171 0.98553999 0.98789255 0.98614119 0.96536797
## 11 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000
## 12 0.58518519 0.94074074 0.17777778 0.93899684 0.94892168 0.94153313 0.85389610
## 13 0.81993569 0.97427653 0.64308682 0.85946679 0.88233068 0.86530966 0.66341991
## 14 0.41666667 0.91666667 0.13513514 0.95661997 0.96367764 0.95842356 0.89610390
## 15 0.64779874 0.94968553 0.30188679 0.92815183 0.93984109 0.93113902 0.82792208
## 16 0.71428571 0.50000000 0.85585586 0.99277000 0.99394627 0.99307059 0.98268398
## 17 0.14285714 0.83333333 0.56756757 0.97830999 0.98183882 0.97921178 0.94805195
## 18 0.30000000 0.90000000 0.27927928 0.96384998 0.96973137 0.96535297 0.91341991
## 19 0.14285714 0.83333333 0.56756757 0.97830999 0.98183882 0.97921178 0.94805195
## 8 9 10 11 12 13 14
## 2
## 3
## 4
## 5
## 6
## 7
## 8
## 9 0.95897436
## 10 0.95897436 0.00000000
## 11 1.00000000 1.00000000 1.00000000
## 12 0.82692308 0.76296296 0.76296296 1.00000000
## 13 0.60128205 0.89710611 0.89710611 1.00000000 0.56591640
## 14 0.87692308 0.66666667 0.66666667 1.00000000 0.28888889 0.69131833
## 15 0.79615385 0.79874214 0.79874214 1.00000000 0.15094340 0.48874598 0.39622642
## 16 0.97948718 0.50000000 0.50000000 1.00000000 0.88148148 0.94855305 0.83333333
## 17 0.93846154 0.33333333 0.33333333 1.00000000 0.64444444 0.84565916 0.50000000
## 18 0.89743590 0.60000000 0.60000000 1.00000000 0.40740741 0.74276527 0.16666667
## 19 0.93846154 0.33333333 0.33333333 1.00000000 0.64444444 0.84565916 0.50000000
## 15 16 17 18
## 2
## 3
## 4
## 5
## 6
## 7
## 8
## 9
## 10
## 11
## 12
## 13
## 14
## 15
## 16 0.89937107
## 17 0.69811321 0.66666667
## 18 0.49685535 0.80000000 0.40000000
## 19 0.69811321 0.66666667 0.00000000 0.40000000
Datos y gráficas del seguimiento de parametros de crecimiento del
cultivo (muestreo 1 y 2)
## Parámetros del cultivo
cultivocebada <- read_excel("C:\\Users\\andre\\OneDrive\\Escritorio\\cebada 1.xlsx",
col_types = c("text", "numeric", "numeric",
"numeric", "numeric", "numeric", "numeric"))
attach(cultivocebada)
## The following object is masked from DMTOTAL2:
##
## MUESTREO
# Altura de cebada:
psych::describe(ALTURAPROMM2)
## vars n mean sd median trimmed mad min max range skew kurtosis se
## X1 1 16 28.57 7.76 28.59 28.41 8.12 15.37 44.09 28.72 0.33 -0.77 1.94
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 2 34.89 9.19 34.89 34.89 9.63 28.4 41.39 12.99 0 -2.75 6.5
## ------------------------------------------------------------
## group: 2
## vars n mean sd median trimmed mad min max range skew kurtosis
## X1 1 2 34.01 14.26 34.01 34.01 14.95 23.92 44.09 20.17 0 -2.75
## se
## X1 10.09
## ------------------------------------------------------------
## group: 3
## vars n mean sd median trimmed mad min max range skew kurtosis se
## X1 1 2 30.73 1.34 30.73 30.73 1.4 29.78 31.67 1.89 0 -2.75 0.95
## ------------------------------------------------------------
## group: 4
## vars n mean sd median trimmed mad min max range skew kurtosis se
## X1 1 2 27.88 8.31 27.88 27.88 8.71 22 33.75 11.75 0 -2.75 5.87
## ------------------------------------------------------------
## group: 5
## vars n mean sd median trimmed mad min max range skew kurtosis se
## X1 1 2 20.31 6.99 20.31 20.31 7.32 15.37 25.25 9.88 0 -2.75 4.94
## ------------------------------------------------------------
## group: 6
## vars n mean sd median trimmed mad min max range skew kurtosis se
## X1 1 2 24.21 6.45 24.21 24.21 6.76 19.65 28.77 9.12 0 -2.75 4.56
## ------------------------------------------------------------
## group: 7
## vars n mean sd median trimmed mad min max range skew kurtosis se
## X1 1 2 27.88 7.74 27.88 27.88 8.12 22.4 33.35 10.95 0 -2.75 5.48
## ------------------------------------------------------------
## group: 8
## vars n mean sd median trimmed mad min max range skew kurtosis se
## X1 1 2 28.7 8.34 28.7 28.7 8.75 22.8 34.6 11.8 0 -2.75 5.9
# Macollas por planta:
psych::describe(MACOLLASM2)
## vars n mean sd median trimmed mad min max range skew kurtosis se
## X1 1 16 5.91 1.94 5.88 5.88 2.04 2.5 9.7 7.2 0.06 -0.89 0.48
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 2 5.25 0.35 5.25 5.25 0.37 5 5.5 0.5 0 -2.75 0.25
## ------------------------------------------------------------
## group: 2
## vars n mean sd median trimmed mad min max range skew kurtosis se
## X1 1 2 6.38 0.53 6.38 6.38 0.56 6 6.75 0.75 0 -2.75 0.37
## ------------------------------------------------------------
## group: 3
## vars n mean sd median trimmed mad min max range skew kurtosis se
## X1 1 2 7.5 0.71 7.5 7.5 0.74 7 8 1 0 -2.75 0.5
## ------------------------------------------------------------
## group: 4
## vars n mean sd median trimmed mad min max range skew kurtosis se
## X1 1 2 4.12 2.3 4.12 4.12 2.41 2.5 5.75 3.25 0 -2.75 1.62
## ------------------------------------------------------------
## group: 5
## vars n mean sd median trimmed mad min max range skew kurtosis se
## X1 1 2 3.75 0.35 3.75 3.75 0.37 3.5 4 0.5 0 -2.75 0.25
## ------------------------------------------------------------
## group: 6
## vars n mean sd median trimmed mad min max range skew kurtosis se
## X1 1 2 6.55 1.34 6.55 6.55 1.41 5.6 7.5 1.9 0 -2.75 0.95
## ------------------------------------------------------------
## group: 7
## vars n mean sd median trimmed mad min max range skew kurtosis se
## X1 1 2 5.85 3.32 5.85 5.85 3.48 3.5 8.2 4.7 0 -2.75 2.35
## ------------------------------------------------------------
## group: 8
## vars n mean sd median trimmed mad min max range skew kurtosis se
## X1 1 2 7.85 2.62 7.85 7.85 2.74 6 9.7 3.7 0 -2.75 1.85
# densidad de cebada:
psych::describe(DENSIDADM2)
## vars n mean sd median trimmed mad min max range skew kurtosis se
## X1 1 16 107.06 29.19 106.5 107.93 23.72 43 159 116 -0.13 -0.35 7.3
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 2 133.5 36.06 133.5 133.5 37.81 108 159 51 0 -2.75 25.5
## ------------------------------------------------------------
## group: 2
## vars n mean sd median trimmed mad min max range skew kurtosis se
## X1 1 2 95.5 20.51 95.5 95.5 21.5 81 110 29 0 -2.75 14.5
## ------------------------------------------------------------
## group: 3
## vars n mean sd median trimmed mad min max range skew kurtosis se
## X1 1 2 109 14.14 109 109 14.83 99 119 20 0 -2.75 10
## ------------------------------------------------------------
## group: 4
## vars n mean sd median trimmed mad min max range skew kurtosis se
## X1 1 2 145 9.9 145 145 10.38 138 152 14 0 -2.75 7
## ------------------------------------------------------------
## group: 5
## vars n mean sd median trimmed mad min max range skew kurtosis se
## X1 1 2 115.5 14.85 115.5 115.5 15.57 105 126 21 0 -2.75 10.5
## ------------------------------------------------------------
## group: 6
## vars n mean sd median trimmed mad min max range skew kurtosis se
## X1 1 2 62.5 27.58 62.5 62.5 28.91 43 82 39 0 -2.75 19.5
## ------------------------------------------------------------
## group: 7
## vars n mean sd median trimmed mad min max range skew kurtosis se
## X1 1 2 97 1.41 97 97 1.48 96 98 2 0 -2.75 1
## ------------------------------------------------------------
## group: 8
## vars n mean sd median trimmed mad min max range skew kurtosis se
## X1 1 2 98.5 27.58 98.5 98.5 28.91 79 118 39 0 -2.75 19.5
## 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") +
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)
Gráficos de cobertura y densidad (Cultivo Vs arvenses)
#Cobertura cultivo Vs cobertura de las arvenses
coberturas <- read_excel("C:\\Users\\andre\\OneDrive\\Escritorio\\Coberturas metro cuadrado.xlsx",
col_types = c("numeric", "numeric", "text",
"numeric", "numeric"))
attach(coberturas)
# densidad y cobertura del cultivo y el total de arvenses por metro cuadrado
psych::describe(DENSIDADESCA)
## vars n mean sd median trimmed mad min max range skew kurtosis se
## X1 1 16 136.06 66.41 112.5 129.07 39.29 79 291 212 1.35 0.54 16.6
describe.by(DENSIDADESCA, group = MUESTRASMC)
## 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 157 87.83 121.5 157 61.53 79 291 212 0.59 -1.58 31.05
## ------------------------------------------------------------
## group: cultivo
## vars n mean sd median trimmed mad min max range skew kurtosis se
## X1 1 8 115.12 27.11 107.5 115.12 16.31 82 159 77 0.54 -1.39 9.59
psych::describe(COBERTURASCA)
## vars n mean sd median trimmed mad min max range skew kurtosis
## X1 1 16 166.48 131.83 147.5 160.84 177.17 10 402 392 0.2 -1.63
## se
## X1 32.96
describe.by(COBERTURASCA, group = MUESTRASMC)
## 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 46.46 30.26 34.5 46.46 27.65 10 89 79 0.27 -1.88 10.7
## ------------------------------------------------------------
## group: cultivo
## vars n mean sd median trimmed mad min max range skew kurtosis se
## X1 1 8 286.5 58.3 277 286.5 31.13 206 402 196 0.66 -0.56 20.61
## Gráfico de la cobertura del cultivo y arvenses por metro cuadrado
ggplot(coberturas, aes(x = PMUEST, y = COBERTURASCA, fill = MUESTRASMC)) +
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(coberturas, aes(x = PMUEST, y = DENSIDADESCA, fill = MUESTRASMC)) +
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 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)
DMTOTAL <- read_excel("C:\\Users\\andre\\OneDrive\\Escritorio\\MALEZASGRUPALES.xlsx",
col_types = c("text", "text", "numeric",
"numeric", "numeric", "numeric",
"numeric", "numeric", "numeric",
"numeric", "numeric", "numeric"))
attach(DMTOTAL)
## The following object is masked from cultivocebada:
##
## MUESTREO
## The following object is masked from DMTOTAL2:
##
## 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))
