円周上の分布 直線上の分布の特性関数・フーリエ変換は容易に円周上に拡張できる(周期性の制約を入れる) 円周上に定義されている分布 指数関数モデルと、直線上分布の円周上への巻きつけ トーラス上の分布、円柱上の分布 2つの独立な円周上のvon Mises分布を使って、トーラス上の"von Mises 分布"を描いてみる ソースは任意次元の2つの球面の掛け合わせになっている d1 <- 2 d2 <- 2 r1 <- 1 r2 <- 0.3 mu1 <- c(1,rep(0,d1-1)) mu2 <- c(1,rep(0,d2-1)) k1 <- 1.2 k2 <- 1.8 n <- 1000 out1 <- RvonMisesFisher(n,mu1,k1) out2 <- RvonMisesFisher(n,mu2,k2) torus.coords <- function(x1,x2){ r
von Mises-Fisher分布 (von Mises-Fisher distribution)† von Mises-Fisher分布は,\(\|\mathbf{\mu}\|=1\) なる \(d\)次元ベクトル \(\mathbf{\mu}\) と \(\sigma\gt0\) をパラメータとし,次の確率密度分布をもつ \[f(\mathbf{x};\mathbf{\mu},\sigma)=\frac{\sigma^{d/2-1}}{(2\pi)^{d/2}I_{d/2-1}(\sigma)}\exp(\sigma\mathbf{\mu}^\top\mathbf{x})\] ただし,\(I_j(x)\) は第1種変形Bessel関数. 超球上の正規分布ともいえる分布.von Mises分布の多次元版.\(\mathbf{\mu}\) が最頻値になる. -- しましま ↑
周期性のある変数・循環する変数を含むモデリングを実践しましたので紹介します。 スライドは埋め込んで、ソースコードのコピペ&解説をメインにします。 とある病んだ院生の体内時計(サーカディアンリズム) from . . 使用したデータは以下。自由に使ってください。 元データ: data.txt 起床時刻だけ抜き出したもの: data_sleepout.txt 起床時刻から起床時刻の時間: data_out2out.txt 起床時刻の累積時刻: data_cum.txt まずはじめに起床時間の分布のパラメータを求めたStanコードは以下になります(model1.stan)。 data { int<lower=0> N; real Sleep_out[N]; } parameters { real<lower=0, upper=2*pi()> mu; real<lower=0, upper=100
神奈川県水産技術センター リアルタイム海況データ城ケ島沖観測ブイ 海況図等のウェブサイト移行について 移行作業のご案内 海況図をはじめとした海況関係のコンテンツは、平成29年8月1日から新しいサイトに逐次移行します。移行の状況を逐次ご案内しております。 農林水産情報センターでの情報提供終了について 海況図の情報提供を行っていた農林水産情報センターのウェブサイトは、平成29年7月31日をもって終了しました。 平成9年8月以来、20年の長きにわたりご利用頂きましてありがとうございました。
リリース、障害情報などのサービスのお知らせ
最新の人気エントリーの配信
処理を実行中です
j次のブックマーク
k前のブックマーク
lあとで読む
eコメント一覧を開く
oページを開く