差分
このページの2つのバージョン間の差分を表示します。
両方とも前のリビジョン前のリビジョン次のリビジョン | 前のリビジョン | ||
r:markovchain [2018/12/12 23:55] – [マルコフ連鎖のための統計的推測] watalu | r:markovchain [2018/12/17 11:47] (現在) – watalu | ||
---|---|---|---|
行 7: | 行 7: | ||
* [[https:// | * [[https:// | ||
* [[https:// | * [[https:// | ||
+ | |||
+ | このパッケージはR 3.5.0以上を必要とする。各自でダウンロードしてインストールする必要があるかもしれない。 | ||
=== マルコフ解析 === | === マルコフ解析 === | ||
行 41: | 行 43: | ||
</ | </ | ||
- | < | ||
- | mcList <- new(" | ||
- | name = "A list of Markov chains" | ||
- | </ | ||
=== markovchainのメソッド === | === markovchainのメソッド === | ||
行 136: | 行 134: | ||
</ | </ | ||
+ | 関数statesはマルコフ連鎖のオブジェクトmcWeatherの中の状態の一覧を表示させる。 | ||
< | < | ||
states(mcWeather) | states(mcWeather) | ||
</ | </ | ||
+ | 関数namesもstatesと同じ動作をする。 | ||
< | < | ||
names(mcWeather) | names(mcWeather) | ||
</ | </ | ||
+ | 関数dimは状態数を返す。 | ||
< | < | ||
dim(mcWeather) | dim(mcWeather) | ||
</ | </ | ||
+ | 関数nameはマルコフ連鎖の名前を返す。これはnewで作成するときにつけている。 | ||
< | < | ||
- | name(mcWeather) | + | name(mcWeather) |
</ | </ | ||
行 279: | 行 278: | ||
< | < | ||
table(data.frame(before=rain$rain[-1], | table(data.frame(before=rain$rain[-1], | ||
+ | </ | ||
+ | |||
+ | 出力は次の通り。 | ||
+ | |||
+ | < | ||
+ | after | ||
+ | before | ||
+ | | ||
+ | 1-5 126 90 79 | ||
+ | | ||
</ | </ | ||
行 299: | 行 308: | ||
{{: | {{: | ||
+ | |||
+ | < | ||
+ | library(dplyr) | ||
+ | library(stringr) | ||
+ | library(DiagrammeR) | ||
+ | library(networkD3) | ||
+ | </ | ||
+ | |||
+ | この準備で頑張ると、[[https:// | ||
この推移から、遷移行列を推定する。 | この推移から、遷移行列を推定する。 | ||
行 372: | 行 390: | ||
< | < | ||
mcrain = rain.mcfit$estimate | mcrain = rain.mcfit$estimate | ||
+ | </ | ||
+ | |||
+ | |||
+ | < | ||
+ | rain.mcfit <- markovchainFit(rain$rain, | ||
</ | </ | ||