Skip to content

Instantly share code, notes, and snippets.

@ito4303
ito4303 / ggcube_map_animation.R
Last active June 20, 2026 02:23
Generate 3D animated map using ggcube
# 3D animated map from DEM
library(terra)
library(ggplot2)
library(geometry)
library(ggcube)
# Load DEM data from GeoTIFF
# 国土地理院基盤地図情報数値標高モデル(DEM10B)を使用
geotiff_file <- "data/543654.tif"
@ito4303
ito4303 / plot_dem_using_ggcube.R
Created June 2, 2026 10:15
Plot DEM data using ggcube
# Plot DEM data using ggcube
library(terra)
library(ggplot2)
library(geometry)
library(ggcube)
# Load DEM data from GeoTIFF
# 国土地理院基盤地図情報数値標高モデル(DEM5A)を使用
geotiff_file <- "data/543654.tif"
# Mexican hat using ggcube
#
# https://tororolab.hatenablog.com/entry/2024/12/16/225348
library(ggplot2)
library(ggcube)
fun <- function(x, y) {
r <- sqrt(x^2 + y^2)
1.5 * cos(r) - 0.75 * cos(3 * r)
@ito4303
ito4303 / ggplot2_aes_string.R
Last active October 16, 2023 10:46
Use string to specify variables in aes() function of ggplot2
library(ggplot2)
library(palmerpenguins)
point_plot <- function(d, x, y) {
ggplot(d) +
geom_point(aes(x = .data[[x]], y = .data[[y]]))
}
(var_names <- colnames(penguins))
@ito4303
ito4303 / D-H.R
Last active July 22, 2023 23:24
Estimation of diameter-height curves
# 『たのしいベイズモデリング2』第5章のRコード
# CmdStanR版
# R 4.3.1 Mac, Stan 2.23.2, CmdStanR 0.5.3 で動作確認
library(readxl)
library(dplyr)
library(cmdstanr)
options(mc.cores = parallel::detectCores())
# データファイルのURL
@ito4303
ito4303 / tmap_draw_sapporo.R
Created July 16, 2023 02:43
draw Sapporo city map using the tmap package
library(sf)
library(tmap)
# 参考にしたところ
# https://qiita.com/ocean_f/items/700aff67f3e35266b0fe
# 国土数値情報からデータをダウンロード
# https://nlftp.mlit.go.jp/ksj/gml/datalist/KsjTmplt-N03-v3_1.html
datafile <- file.path("data", "N03-20230101_01_GML",
"N03-23_01_230101.geojson")
@ito4303
ito4303 / zib_regression.Rmd
Last active February 16, 2023 08:31
Zero-inflated beta regression
---
title: "Zero-inflated beta regression"
output: html_notebook
---
## Setup
```{r setup}
library(ggplot2)
library(zoib)
@ito4303
ito4303 / beta_regression.Rmd
Created February 15, 2023 10:32
Beta regression using R and Stan
---
title: "Beta regression"
output: html_notebook
---
## Setup
```{r setup}
library(ggplot2)
library(betareg)
@ito4303
ito4303 / lrtest.Rmd
Created December 22, 2022 06:50
Likelihood ratio test in R (using base functions or lrtest function)
---
title: "R Notebook"
output: html_notebook
---
## An example of likelihood ratio tests in R
```{r setup}
library(lmtest)
library(palmerpenguins)
@ito4303
ito4303 / glht.Rmd
Created December 22, 2022 05:19
An example to use glht function in the multcomp package
---
title: "R Notebook"
output: html_notebook
---
## An example to use glht function in the multcomp package
```{r setup}
library(multcomp)
library(palmerpenguins)