In turf-reserves/mpa.tools: Tools to simulate and analyze no-take Marine Reserves

suppressPackageStartupMessages({
  library(tidyverse)
  library(magrittr)
  })

jc <- read.csv("../MAREAmanuscript/Data/natividad_fish.csv") %>%
  MPAtools::richness(location = "Isla Natividad") %>%
  mutate(Post = ifelse(Ano <= 2006, 0, 1),
         Ano = as.factor(Ano)) %>%
  lm(Indicador ~ Ano + Zona + Zona * Post + Temperatura + Profundidad + Visibilidad, data = .)

source('./R/robust_se.R')
source("./R/reg_table.R")

reg_table(model = jc,
          title = "Tabla de regresion para riqueza especifica",
          dep.var.labels = "Riqueza (Spp/transecto)") %>% 
  cat()

turf-reserves/mpa.tools documentation built on June 1, 2019, 2:55 a.m.