Skip to content

Instantly share code, notes, and snippets.

View agila5's full-sized avatar

Andrea Gilardi agila5

72 followers · 138 following
View GitHub Profile
@agila5
agila5 / gist:1880071f01228cf4ea3686677e1882e5
Created April 12, 2019 10:03
``` r
# packages
library(dplyr)
#>
#> Attaching package: 'dplyr'
#> The following objects are masked from 'package:stats':
#>
#> filter, lag
#> The following objects are masked from 'package:base':
#>
@agila5
agila5 / inla.md
Created February 24, 2020 17:47 
# packages
library(INLA)
#> Loading required package: Matrix
#> Loading required package: sp
#> Loading required package: parallel
#> This is INLA_19.09.03 built 2020-02-20 16:12:21 UTC.
#> See www.r-inla.org/contact-us for how to get help.
#> To enable PARDISO sparse library; see inla.pardiso()
@agila5
agila5 / inla_zeroinflated0spatial.md
Last active March 2, 2020 11:29 
# packages
library(INLA)
#> Loading required package: Matrix
#> Loading required package: sp
#> Loading required package: parallel
#> Loading required package: foreach
#> This is INLA_20.02.19 built 2020-02-19 06:22:22 UTC.
#> See www.r-inla.org/contact-us for how to get help.
library(igraph)
@agila5
agila5 / reprex_for_som.md
Created March 14, 2020 16:34 
# packages
library(sf)
#> Linking to GEOS 3.6.1, GDAL 2.2.3, PROJ 4.9.3
library(osmdata)
#> Data (c) OpenStreetMap contributors, ODbL 1.0. https://www.openstreetmap.org/copyright
library(tmap)
#> Warning: replacing previous import 'sf::st_make_valid' by
#> 'lwgeom::st_make_valid' when loading 'tmap'
tmap_mode("view")
@agila5
agila5 / countreg.md
Created April 14, 2020 11:54 
# packages
remotes::install_cran("countreg", repos = "http://R-Forge.R-project.org")
#> Skipping install of 'countreg' from a cran remote, the SHA1 (0.2-1) has not changed since last install.
#>   Use `force = TRUE` to force installation
library(countreg)
#> Loading required package: MASS

# Examples
dhpois(x = 0, lambda = 1, pi = 0.5, log = TRUE)
@agila5
agila5 / PITs.md
Created April 14, 2020 20:18 
# packages
library(countreg) # for CDF of Hurdle Poisson distribution
#> Loading required package: MASS

# Take nonrandomized PIT function from Czado et al. (2009)
czado_pit <- function(x, Px, Px1, n.bins=10, y.max=2.75, my.title="PIT Histogram")
{
  a.mat <- matrix(0,n.bins,length(x))
  k.vec <- pmax(ceiling(n.bins*Px1),1)
@agila5
agila5 / inla_nbinomial.md
Created April 20, 2020 10:16 
# packages
library(INLA)
#> Loading required package: Matrix
#> Loading required package: sp
#> Loading required package: parallel
#> Loading required package: foreach
#> This is INLA_20.03.09 built 2020-03-09 09:12:35 UTC.
#> See www.r-inla.org/contact-us for how to get help.
@agila5
agila5 / inla_make_lincomb.md
Last active June 15, 2020 07:22 
# packages
library(INLA)
#> Loading required package: Matrix
#> Loading required package: sp
#> Loading required package: parallel
#> Loading required package: foreach
#> This is INLA_20.03.09 built 2020-03-09 09:12:35 UTC.
#> See www.r-inla.org/contact-us for how to get help.
inla.version()
@agila5
agila5 / example_sfnetworks.md
Created July 12, 2020 16:46 
# packages
library(sf)
#> Linking to GEOS 3.8.0, GDAL 3.0.4, PROJ 6.3.1
library(sfnetworks)
library(tidygraph)
#> 
#> Attaching package: 'tidygraph'
#> The following object is masked from 'package:stats':
#>
@agila5
agila5 / example_W.md
Created July 29, 2020 15:31 
# packages
library(spdep)
#> Loading required package: sp
#> Loading required package: spData
#> Loading required package: sf
#> Linking to GEOS 3.8.0, GDAL 3.0.4, PROJ 6.3.1
library(rgdal)
#> rgdal: version: 1.5-12, (SVN revision 1018)
#> Geospatial Data Abstraction Library extensions to R successfully loaded