Skip to content

Instantly share code, notes, and snippets.

Francisco Rodriguez-Sanchez Pakillo

Block or report user

Report or block Pakillo

Hide content and notifications from this user.

Learn more about blocking users

Contact Support about this user’s behavior.

Learn more about reporting abuse

Report abuse
View GitHub Profile
@Pakillo
Pakillo / cc_paleotemp.R
Created Sep 3, 2018
Figure of Cenozoic (last 66 million years) and future global temperatures expected with climate change
View cc_paleotemp.R
library(ggplot2)
library(cowplot)
library(latex2exp)
#### PALEOTEMPERATURES ####
## Data from Hansen: http://www.columbia.edu/~jeh1/mailings/2012/20120508_ClimateSensitivity.pdf
## Available here: http://www.columbia.edu/~mhs119/Sensitivity+SL+CO2/Table.txt
@Pakillo
Pakillo / repromap.R
Last active Mar 21, 2018
another reproducible map with Rstats
View repromap.R
library(sf)
library(ggspatial)
library(ggthemes)
coast <- rnaturalearth::ne_coastline(scale = "small", returnclass = "sp")
forests <- st_read("data/forestmap/cfc_gen.shp")
forest.redux <- rmapshaper::ms_simplify(forests)
forest.redux.sp <- as(forest.redux, "Spatial")
@Pakillo
Pakillo / Worldclim_Precipitation_Uncertainty.md
Last active May 16, 2017
Uncertainty in global climate databases: comparing Worldclim with Iberian climate atlas
View Worldclim_Precipitation_Uncertainty.md

Precipitation: Worldclim vs Iberian Climatic Atlas

Francisco Rodriguez-Sanchez 2017-05-16

library(knitr)
opts_chunk$set(cache = TRUE, message = FALSE)
library(raster)
@Pakillo
Pakillo / fig_label.R
Created Apr 10, 2017
R function to draw figure labels in base plots
View fig_label.R
## from: https://logfc.wordpress.com/2017/03/15/adding-figure-labels-a-b-c-in-the-top-left-corner-of-the-plotting-region/
fig_label <- function(text, region="figure", pos="topleft", cex=NULL, ...) {
region <- match.arg(region, c("figure", "plot", "device"))
pos <- match.arg(pos, c("topleft", "top", "topright",
"left", "center", "right",
"bottomleft", "bottom", "bottomright"))
if(region %in% c("figure", "device")) {
@Pakillo
Pakillo / simulate.gamm.R
Created Mar 3, 2017 — forked from gavinsimpson/simulate.gamm.R
S3 method for simulate() for "gamm" objects from package mgcv
View simulate.gamm.R
`simulate.gamm` <- function(object, nsim = 1, seed = NULL, newdata,
freq = FALSE, unconditional = FALSE, ...) {
if (!exists(".Random.seed", envir = .GlobalEnv, inherits = FALSE))
runif(1)
if (is.null(seed))
RNGstate <- get(".Random.seed", envir = .GlobalEnv)
else {
R.seed <- get(".Random.seed", envir = .GlobalEnv)
set.seed(seed)
RNGstate <- structure(seed, kind = as.list(RNGkind()))
@Pakillo
Pakillo / ShinyPerspectivePlot.R
Created Dec 19, 2016 — forked from leeper/ShinyPerspectivePlot.R
Shiny App for 3D Interaction plots
View ShinyPerspectivePlot.R
library("shiny")
library("shinyjs")
ui <- shinyUI(pageWithSidebar(
titlePanel("", "3D Perspective Plot for Interaction Effects"),
sidebarPanel(
tabsetPanel(
tabPanel("Data",
@Pakillo
Pakillo / beamer_two_col.Rmd
Created Nov 7, 2016 — forked from jhollist/beamer_two_col.Rmd
This .Rmd shows how to use two (or more) columns in a Beamer Presentation from RStudio
View beamer_two_col.Rmd
---
title: "Untitled"
author: "Jeff W. Hollister"
date: "2/2/2016"
output: beamer_presentation
---
```{r setup, include=FALSE}
knitr::opts_chunk$set(echo = FALSE)
@Pakillo
Pakillo / gist:d2ee89f78a5a85dc4b1593e7bf82afe9
Last active Nov 4, 2016
ploting prior and posterior Bayes
View gist:d2ee89f78a5a85dc4b1593e7bf82afe9
- https://mvuorre.github.io/r/ggplot2-visualizing-prior-posterior/
- https://users.soe.ucsc.edu/~draper/beta-distributions-r.txt
- http://shinyapps.org/apps/RGraphCompendium/index.php#prior-and-posterior
- https://alexanderetz.com/2015/07/25/understanding-bayes-updating-priors-via-the-likelihood/
- http://doingbayesiandataanalysis.blogspot.com.es/2015/10/posterior-predicted-distribution-for.html
- http://www.stat.cmu.edu/~hseltman/rube/rubeManual.pdf
- http://stackoverflow.com/questions/19476744/how-to-plot-bayesian-prior-and-posterior-distributions-in-one-panel-using-r
@Pakillo
Pakillo / themes.R
Last active Sep 6, 2018
ggplot_themes
View themes.R
### My favourite ggplot2 themes, always accessible ###
## just source this gist from any project ##
#####################################################################################
# Noam's default ggplot2 styling #####
# Source: https://github.com/noamross/noamtools/blob/master/R/theme_nr.R
theme_nr <-theme(text=element_text(family="Helvetica", size=14),
panel.grid.major.x=element_blank(),
You can’t perform that action at this time.