Pacoche_Analisis_2020
Diego J. Lizcano, Enrique de la Montaña, et al.
24/4/2020
R Markdown
This is an R Markdown document. Markdown is a simple formatting syntax for authoring HTML, PDF, and MS Word documents. For more details on using R Markdown see http://rmarkdown.rstudio.com.
When you click the Knit button a document will be generated that includes both content as well as the output of any embedded R code chunks within the document. You can embed an R code chunk like this:
library(xtable)
library(lubridate)##
## Attaching package: 'lubridate'## The following object is masked from 'package:base':
##
## datelibrary(dplyr)## Warning: package 'dplyr' was built under R version 3.6.3##
## Attaching package: 'dplyr'## The following objects are masked from 'package:lubridate':
##
## intersect, setdiff, union## The following objects are masked from 'package:stats':
##
## filter, lag## The following objects are masked from 'package:base':
##
## intersect, setdiff, setequal, unionlibrary(maptools)## Warning: package 'maptools' was built under R version 3.6.2## Loading required package: sp## Warning: package 'sp' was built under R version 3.6.2## Checking rgeos availability: TRUE##
## Attaching package: 'maptools'## The following object is masked from 'package:xtable':
##
## label# Load 'rgdal' package, which is used to read/write shapefiles and rasters
library(rgdal)## rgdal: version: 1.4-8, (SVN revision 845)
## Geospatial Data Abstraction Library extensions to R successfully loaded
## Loaded GDAL runtime: GDAL 2.2.3, released 2017/11/20
## Path to GDAL shared files: C:/Users/usuario/Documents/R/win-library/3.6/rgdal/gdal
## GDAL binary built with GEOS: TRUE
## Loaded PROJ.4 runtime: Rel. 4.9.3, 15 August 2016, [PJ_VERSION: 493]
## Path to PROJ.4 shared files: C:/Users/usuario/Documents/R/win-library/3.6/rgdal/proj
## Linking to sp version: 1.3-2# source("C:/Users/diego.lizcano/Box Sync/CodigoR/Toshiba/Pacoche_unmakerd/code/TEAM_code.R")
# source("C:/Users/diego.lizcano/Box Sync/CodigoR/Toshiba/Pacoche_unmakerd/code/calendar.R")
source("D:/BoxFiles/Box Sync/CodigoR/Toshiba/Pacoche_unmakerd/code/TEAM_code.R")## Loading required package: reshape2## Loading required package: plyr## Warning: package 'plyr' was built under R version 3.6.2## ------------------------------------------------------------------------------## You have loaded plyr after dplyr - this is likely to cause problems.
## If you need functions from both plyr and dplyr, please load plyr first, then dplyr:
## library(plyr); library(dplyr)## ------------------------------------------------------------------------------##
## Attaching package: 'plyr'## The following objects are masked from 'package:dplyr':
##
## arrange, count, desc, failwith, id, mutate, rename, summarise,
## summarize## The following object is masked from 'package:lubridate':
##
## here## Loading required package: ggplot2########## Call saved and fixed data set arrays 1 a 4
### pacoche.raw object
load("D:/BoxFiles/Box Sync/CodigoR/Toshiba/Pacoche_unmakerd/data/Pacoche_arrays_1_a_4.RData")Especies registradas
Las especies registradas en Pacoche fueron
Numero_de_registros naive_occu especie
-------------------- ----------- --------------------------
1040 1.0000000
100 0.1612903 Bos primigenius
93 0.3870968 Cuniculus paca
52 0.0806452 Capra aegagrus
41 0.2096774 Dasypus novemcinctus
39 0.3709677 Eira barbara
38 0.1451613 Equus ferus
36 0.1451613 Tinamus major
32 0.1129032 Rattus rattus
30 0.1935484 Homo sapiens
27 0.2258065 Leptotila verreauxi
24 0.1612903 Odocoileus virginianus
23 0.2096774 Ortalis erythroptera
17 0.1774194 Tamandua mexicana
17 0.1774194 Procyon cancrivorus
16 0.0967742 Sylvilagus brasiliensis
11 0.0806452 Canis lupus
11 0.0967742 Dasyprocta punctata
10 0.0967742 Leopardus wiedii
10 0.0322581 Equus africanus
9 0.1129032 Leopardus pardalis
7 0.1129032 Puma yauguaroundi
7 0.0645161 Sciurus stramineus
4 0.0161290 Zenaida auriculata
4 0.0322581 Gallus gallus
3 0.0161290 Momotus momota
3 0.0161290 Alouatta palliata
3 0.0483871 Ortalis vetula
3 0.0483871 Buteogallus urubitinga
3 0.0483871 Didelphis marsupialis
3 0.0161290 Cebus capucinus
3 0.0322581 Felis silvestris
2 0.0322581 Galictis vittata
2 0.0161290 Coragyps atratus
2 0.0161290 Pipistrellus pipistrellus
1 0.0161290 Damophila julie
1 0.0161290 Carduelis siemiradzkii
1 0.0161290 Sus scrofa
1 0.0161290 Columbina cruziana
1 0.0161290 Artibeus jamaicensis Covariables
Inicialmente se probó con cinco covariables para ajustar los modelos de ocupación. Estas cinco covariables fueron: Altitud (elev), Pendiente (slope) y Distancia a la carretera pavimentada (dis_rd) como covariables geográficas y altura del dosel, cobertura del dosel y área basal como covariables que se midieron momento de retirar las cámaras. Estas medidas se tomaron usando la metodología del cuadrante centrado en un punto.
La Altitud se obtuvo de una imagen SRTM del repositorio de [CGIAR](http://srtm.csi.cgiar.org/). La pendiente se infirió a partir de la altitud y la distancia a las carreteras se obtuvo de un mapa del Ministerio del Medio Ambiente del Ecuador.
<img src="analisis_2020_files/figure-html/covs-1.png" width="768" />
```
## [1] 45
```
```
## [1] 17
```
<img src="analisis_2020_files/figure-html/covs-2.png" width="768" /><img src="analisis_2020_files/figure-html/covs-3.png" width="768" />Función de modelos
Modelos de ocupacion por especie
La Guanta (Cuniculus paca)
Matriz de datos colapsada a 15 dias
### Selección de Modelos
[1] "Cuniculus paca"
Evidence ratio between models 'p(basal_a)psi(in_out)' and 'p(elev)psi(in_out)':
1.96 analisis Cuniculus paca p(basal_a) psi(in_out) sp=3
## Loading required package: unmarked## Loading required package: lattice## Loading required package: parallel## Loading required package: Rcpp##
## Attaching package: 'unmarked'## The following object is masked from 'package:sp':
##
## coordinatesEvidence ratio between models ‘p(basal_a)psi(in_out)’ and ‘p(elev)psi(in_out)’: 1.92
mod.avg.pred uncond.se lower.CL upper.CL[1,] 0.5759323 0.09076592 0.39596662 0.7377875 [2,] 0.2090277 0.10954273 0.06727373 0.4919019 Predicted SE lower upper 1 0.3960222 0.05278911 0.2984619 0.5026246 2 0.3960222 0.05278911 0.2984619 0.5026246 
El Armadillo (Dasypus novemcinctus)
Matriz de datos colapsada a 15 dias
Selección de Modelos
[1] “Dasypus novemcinctus”
Evidence ratio between models ‘p(basal_a)psi(canopy_h)’ and ‘p(elev)psi(in_out)’: 1.01
analisis Dasypus novemcinctus sp=14
El tigrillo (Eira barbara)
Matriz de datos colapsada a 15 dias
### Selección de Modelos
[1] "Eira barbara"
Evidence ratio between models 'p(elev)psi(canopy_c)' and 'p(.)psi(.)':
1.07 analisis Eira barbara p(elev)psi(slope) sp=6
El Venado de cola blanca (Odocoileus virginianus)
Matriz de datos colapsada a 15 dias
Selección de Modelos
[1] "Odocoileus virginianus"
Evidence ratio between models 'p(elev)psi(in_out)' and 'p(basal_a)psi(in_out)':
5.42 analisis Odocoileus virginianus p(elev)psi(in_out) sp=18
El mapache (Procyon cancrivorus)
Matriz de datos colapsada a 15 dias
Selección de Modelos Procyon cancrivorus
[1] "Procyon cancrivorus"
Evidence ratio between models 'p(basal_a)psi(canopy_c)' and 'p(elev)psi(canopy_c)':
1.14 analisis Procyon cancrivorus p(basal_a)psi(canopy_c) sp=12
El oso hormigero (Tamandua mexicana)
Matriz de datos colapsada a 15 dias
Selección de Modelos Tamandua mexicana
[1] "Tamandua mexicana"
Evidence ratio between models 'p(elev)psi(basal_a)' and 'p(elev)psi(canopy_h)':
2.72 analisis Tamandua mexicana p(.)psi(.) sp=11
El conejo (Sylvilagus brasiliensis)
Matriz de datos colapsada a 15 dias
Selección de Modelos Sylvilagus brasiliensis
[1] "Sylvilagus brasiliensis"
Evidence ratio between models 'p(.)psi(.)' and 'p(basal_a)psi(canopy_h)':
1.76 analisis Sylvilagus brasiliensis p(.)psi(.) sp=22
El margay (Leopardus wiedii)
Matriz de datos colapsada a 15 dias
Selección de Modelos Leopardus wiedii
[1] "Leopardus wiedii"
Evidence ratio between models 'p(elev)psi(canopy_c)' and 'p(elev)psi(dist_def)':
15.15 analisis Leopardus wiedii p(elev)psi(canopy_c) sp=10
El nieque (Dasyprocta punctata)
Matriz de datos colapsada a 15 dias
Selección de Modelos Dasyprocta punctata
[1] "Dasyprocta punctata"
Evidence ratio between models 'p(basal_a)psi(canopy_h)' and 'p(elev)psi(dist_rd)':
2.13 analisis Dasyprocta punctata p(basal_a)psi(dist_rd) sp=26
Evidence ratio between models ‘p(basal_a)psi(canopy_h)’ and ‘p(elev)psi(dist_rd)’: 2.15
##
## Attaching package: 'MuMIn'## The following objects are masked from 'package:AICcmodavg':
##
## AICc, DIC, importanceEl tigrillo (Leopardus pardalis)
Matriz de datos colapsada a 15 dias
Selección de Modelos Leopardus pardalis
[1] "Leopardus pardalis"
Evidence ratio between models 'p(elev)psi(canopy_c)' and 'p(.)psi(.)':
1.1 analisis Leopardus pardalis p(.)psi(.) sp=13
Ardilla (Sciurus stramineus)
Matriz de datos colapsada a 15 dias
Selección de Modelos Sciurus stramineus
[1] "Sciurus stramineus"
Evidence ratio between models 'p(.)psi(.)' and 'p(basal_a)psi(in_out)':
2.23