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()
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.