統計計算を行うのに何かと乱数が必要になります。Fortran には組み込み関数として random_number() がありますが、 コンパイラによって実装が異なります*1*2。環境によって品質が変わると困るので、自分で実装したものを使うのが安全だと思います。私はメル…

何事も実際にやってみるのは大切です．経験することで直観が養われます．本記事は分散の大きさを感覚的に理解することが目的です．分散は平均とともに統計学で最初に習うものですから，定義を知っている人は多いでしょう，しかし，散布図を見て分散がどのく…

本稿ではベイズ統計の基本要素である真の分布，統計モデル，事前分布を説明し，真の分布のベイズ推測とは何かを説明します．参考書にしているのは渡辺（2012）*1，Watanabe（2018）*2および，渡辺澄夫先生の解説記事や講義資料です． 野外観測や実験等で我々…

統計学では真の分布（一般的に未知）からサンプルが生成されると考えます．分析者は統計モデルを構築し，その統計モデルとサンプルから真の分布を推測します．これが統計学における問題設定です*1．同じ統計モデルでも推測方法が違えば予測分布（推測された…

Fortranの配列用組み込み関数の中にpack関数というものがあります．直観的に配列を扱うことができ，データをまとめたりする際に大変便利な関数ですpack関数は配列xの要素の内，条件を満たす要素（正確にはmaskで指定した論理型配列の成分が真である部分の要…

配列にデータ（生物の個体数や水位など）を格納して平均を求める場合，doループを使ってもよいですが，Fortranを使うなら配列用の組み込み手続きを活用したいものです．大きさ10の整数値配列xの成分を x(:) = (/1, 2, 3, 4, 5, 6, 7, 8, 9, 10/) で与えます…