Rで解析：月の満ち欠けを表現するのに便利！！「gggibbous」パッケージ

月の満ち欠けを表現するのに便利なパッケージの紹介です。「ggplot2」パッケージと組み合わせて利用できるので大変便利です。

パッケージバージョンは0.1.1。実行コマンドはwindows 11のR version 4.1.2で確認しています。

パッケージのインストール

下記、コマンドを実行してください。

#&#12497;&#12483;&#12465;&#12540;&#12472;&#12398;&#12452;&#12531;&#12473;&#12488;&#12540;&#12523;
install.packages("gggibbous")

実行コマンド

詳細はコメント、パッケージのヘルプを確認してください。

#&#12497;&#12483;&#12465;&#12540;&#12472;&#12398;&#35501;&#12415;&#36796;&#12415;
library("gggibbous")

###&#12487;&#12540;&#12479;&#20363;&#12398;&#20316;&#25104;#####
#tidyverse&#12497;&#12483;&#12465;&#12540;&#12472;&#12364;&#12394;&#12369;&#12428;&#12400;&#12452;&#12531;&#12473;&#12488;&#12540;&#12523;
if(!require("tidyverse", quietly = TRUE)){
  install.packages("tidyverse");require("tidyverse")
}
#&#12487;&#12540;&#12479;&#20363;&#12398;&#20869;&#23481;
Moon <- seq(-3, 3, by = 0.3) #&#12467;&#12469;&#12452;&#12531;&#12459;&#12540;&#12502;&#12398;&#31684;&#22258;
MoonData <- tibble(x = seq(Moon),
                   y = cos(Moon/2),
                   #&#26376;&#12398;&#28288;&#12385;&#27424;&#12369;&#29575;:1&#12391;&#28288;&#26376;,0&#12391;&#26032;&#26376;
                   ratio = cos(Moon/2)) %>%
  #&#26376;&#12398;&#28288;&#12385;&#12398;&#26041;&#21521;:FALSE;&#24038;&#21491;,TRUE;&#21491;&#24038;
  mutate(right = rep(c(FALSE, TRUE), each = 11,
                     length = length(Moon)))
########

#&#26376;&#12398;&#28288;&#12385;&#27424;&#12369;&#12434;&#34920;&#29694;:geom_moon&#12467;&#12510;&#12531;&#12489;
#geom_moon&#12467;&#12510;&#12531;&#12489;&#12434;&#35079;&#25968;&#20351;&#12358;&#12392;&#34920;&#29694;&#12364;&#24195;&#12364;&#12426;&#12414;&#12377;
ggplot(MoonData, aes(x = x, y = y)) +
  geom_moon(ratio = 1, size = 8, fill = "black") +
  geom_moon(aes(ratio = ratio, right = right),
            fill = "yellow", size = 8)

出力例

少しでも、あなたの解析が楽になりますように！！