Datos

Los datos se encuentran en el archivo ‘ProdPrim.csv’, y se cargan a R usando la función read_csv.

library(readr)
ProdPrim <- read_csv("ProdPrim.csv", na = "NA")
head(ProdPrim)
## # A tibble: 6 x 28
##   site_name site_id country landmark state    lat    long elevation map_mm
##   <chr>     <chr>   <chr>   <chr>    <chr>  <dbl>   <dbl>     <dbl>  <dbl>
## 1 Badkhyz   BDK     Turkme… Serheta… Mary   35.7    62          700   292.
## 2 Beacon H… BCN     UK      Portsmo… West…  50.9    -0.85       205   858.
## 3 Bridger   BRD     USA     Bozeman  Mont…  45.8  -111.        2340   382 
## 4 Calabozo  CLB     Venezu… Calabozo Guár…   8.93  -67.4         98  1257 
## 5 Canas     CNS     Costa … Canas    Guan…  10.4   -85.1         45  1539.
## 6 Charlevi… CHR1    Austra… Charlev… Quee… -26.4   146.         304   483.
## # … with 19 more variables: meantmin_C <dbl>, meantmax_C <dbl>,
## #   biome <chr>, soil_type <chr>, ecoregion <chr>, species <chr>,
## #   vegetation <chr>, plot_mgmt <chr>, mgmt_hist <chr>, agb_gm2 <dbl>,
## #   bgb_gm2 <dbl>, ANPP_gm2yr <dbl>, BNPP_gm2yr <dbl>, start_date <dbl>,
## #   end_date <dbl>, s_interval <chr>, s_period <chr>, m_periods <chr>,
## #   reference <chr>

Gráficas de relación de productividad primaria con atributos geográficos

Gráficas de ANPP versus latitud y elevación.

Latitud

library(ggplot2)
#gráfica de puntos
NPP_lat <- ggplot(data=ProdPrim, aes(x=lat, y=ANPP_gm2yr)) +
  geom_point(pch=19, color="green", size=2) +
  labs(x="Latitud", y="ANPP, g/m^2/yr")
NPP_lat

#gráfica con tendencia
NPP_lat <- ggplot(data=ProdPrim, aes(x=lat, y=ANPP_gm2yr)) +
  geom_point(pch=19, color="green", size=2) +
  geom_smooth(stat = "smooth", method = "auto") +
  labs(x="Latitud", y="ANPP, g/m^2/yr")
NPP_lat

Prueba de hipótesis - Regresión polinomial

library(sjPlot)
## Learn more about sjPlot with 'browseVignettes("sjPlot")'.
regpoly <- lm(ANPP_gm2yr ~ lat + I(lat^2), data=ProdPrim)
tab_model(regpoly)
  ANPP gm 2 yr
Predictors Estimates CI p
(Intercept) 989.40 867.34 – 1111.45 <0.001
lat 6.84 2.00 – 11.67 0.006
lat^2 -0.31 -0.41 – -0.22 <0.001
Observations 122
R2 / R2 adjusted 0.303 / 0.291

Elevación

library(ggplot2)
#gráfica de puntos
NPP_elev <- ggplot(data=ProdPrim, aes(x=elevation, y=ANPP_gm2yr)) +
  geom_point(pch=19, color="green", size=2) +
  labs(x="Elevación, m.s.n.m.", y="ANPP, g/m^2/yr")
NPP_elev

#gráfica con tendencia
NPP_elev <- ggplot(data=ProdPrim, aes(x=elevation, y=ANPP_gm2yr)) +
  geom_point(pch=19, color="green", size=2) +
  geom_smooth(stat = "smooth", method = "auto") +
  labs(x="Elevación, m.s.n.m.", y="ANPP, g/m^2/yr")
NPP_elev