In turfeffect/MPAtools: Tools to simulate and analyze no-take Marine Reserves

title: `r params$title`

knitr::opts_chunk$set(echo = F, fig.height = 3, warning = F, message = F, results = 'asis')

suppressPackageStartupMessages({
  library(ggplot2)
  library(dplyr)
  library(tidyr)
  library(MPAtools)
  library(reshape)
  library(stargazer)
  library(ggExtra)
})

Este documento fue generado por el paquete MPAtools, uno de los entregables del grupo TURFeffect. Los resultados presentados son una version preliminar del software que se esta desarrollando, y no recomendamos la toma de decisiones basado en la informacion aqui presentada

Acentos y caracteres especiales omitidos

Resumen

Leyenda

peces <-  params$peces
socioeco <- params$socioeco
invertebrados <-  params$invertebrados
pesca <-  params$pesca
comunidad <-  params$comunidad
reserva <-  params$reserva
control <-  params$control

Dp <- summary(turfeffect(density(peces, comunidad), reserva, control))

Sl <- summary(lm(Landings~Year, socioeco))


##

summary <-list(
  Bio = list(
    P = list(
      Dp = paste(
        score(x = data.frame(est = coefficients(Dp)[7], p = coefficients(Dp)[28])), ".jpg", sep = ""
      )
    )
  ),
  Soc = list(
    L = paste(
      score(x = data.frame(est = coefficients(Sl)[1], p = coefficients(Sl)[7])), ".jpg", sep = ""
    )
  )
)

Cat | Indicador | Valor
-------|-------------------------|---------------------------------------------- Bio | Densidad Peces | r knitr::include_graphics(summary$Bio$P$Dp) Soc | Arribos de Langosta | r knitr::include_graphics(summary$Soc$L)

Ecologicos

Peces

Densidad

Np <- density(data = peces, location = comunidad)

mpa_plot3(Np, reserve = reserva, control = control, y.lab = "Densidad (org / m^2)")

stargazer(turfeffect(Np, reserve = reserva, control = control), dep.var.labels = "Densidad (org / transecto)", type = "html", dep.var.caption = "", report = "vc*", single.row = T, omit.stat = c("adj.rsq", "n"), digits = 2, df = F, notes = "+p < 0.1, ++p<0.05, +++p<0.001", notes.append = FALSE, star.char = "+")

Socioeconomicos

Arribos

Arribos de langosta

ggplot(socioeco, aes(x = Year, y = Landings)) +
  geom_point(color = "red") +
  geom_line() +
  theme_bw()

stargazer(lm(Landings~Year, socioeco), dep.var.labels = "Densidad (org / transecto)", type = "html", dep.var.caption = "", report = "vc*", single.row = T, omit.stat = c("adj.rsq", "n"), digits = 2, df = F, notes = "+p < 0.1, ++p<0.05, +++p<0.001", notes.append = FALSE, star.char = "+")