ログランク検定、重み付きログランク検定たちについてのメモ。 【目次】 計算式等 ノート サンプルデータ まずは手計算してみる SASでの実行 Rでの実行 プログラムコード 参考 前回のKaplan-meier 法に続く内容。 cochineal19.hatenablog.com 計算式等 先に…

カプランマイヤー法（Kaplan-meier method）についてのメモ。 【目次】 計算式等 ノート 打ち切りとは サンプルデータ 手計算してみる Rでの実行 SASでの実行 プログラムコード 参考 累積生存率とその信頼区間について記載する。 信頼区間はRで扱えるPlain（…

最近SQL Serverの記事を挙げられていなかったので、 題目通り、任意の分類ごとに１つだけフラグを立てるやり方。 例えば、次のようなデータにID毎に１つだけフラグを立てたいとか、ID・日付毎に１つだけフラグを立てたいというような場面。 用途として、立て…

ARIMA、SARIMA、SARIMAXモデルをRで試してみる。 【目次】 用意したデータ データ加工 観測値の時系列の描画 単位根の検定 ARIMAモデルの実行 SARIMAモデルの実行 SARIMAXモデルの実行 モデルの比較 おわりに 参考 用意したデータ 今回、横浜市の月別火災発…

本記事は、AR、MA、ARMA、ARIMA、SARIMA、ARIMAXモデル等の特徴について。 【目次】 計算式等 ノート 参考 計算式等 ＜ 前提＞ ■ 弱定常性（Weakly stationary） 時系列の期待値、分散および自己共分散が常に一定である状態（時点 t および k に依存しない状…

自己相関、自己共分散について、視覚的に分かるようにメモ。 【目次】 計算式 ノート 参考 自己相関、「自己」とつくので、ちょっと分かりづらい。 計算式 ■ 自己共分散 ■ 自己相関 ※t=時期、s=ラグ ノート サンプルデータで動きを見てみる。 Rで準備したデ…

Rでのデータ分割の方法のメモ。 データ準備 rsample で分割 dplyr で分割 参考 データ準備 1500行の何の意味もないデータ。 ads <- data.frame(id=seq(1:1500)) rsample で分割 library(rsample) ads.sp1 <- rsample::initial_split(ads, prop=3/4) adtrain1 …

共分散分析についてめも。 【目次】 帰無仮説、対立仮設 計算式等 計算例 プログラムコード 交互作用項での平行の評価 参考 共分散分析は、分散分析と重回帰分析を組み合わせた手法と言われる。分散分析では説明変数に集団などカテゴリー変数（要因）のみを…

二元配置分散分析（Two-way ANOVA）について。 【目次】 帰無仮説、対立仮設 計算式等 計算例 プログラムコード 交互作用のイメージ 参考 一元配置分散分析 は１要因に対する分析手法。２要因の場合は、二元配置分散分析を用いる。※本記事の二元配置分散分析…

備忘録がてら、ggplot で svg形式ファイルを作る方法とリンクの埋め込み方法を書きます。 svg形式ファイルとは？ ggplotでの出力方法 svg形式ファイルへのリンク埋め込み お試し画像 svg形式ファイルとは？ svg形式ファイルは xmlをベースにしたベクター画像…

ロジスティック回帰の評価について。 【目次】 計算式等 計算例 プログラムコード 参考 前回で理論部分を取り扱った。その続きでモデルの評価方法。 cochineal19.hatenablog.com 計算式等 モデル評価には、大きく分けて「Discrimination（判別能力）」と「Ca…

ロジスティック回帰分析について。 【目次】 計算式等 ノート プログラムコード 交差エントロピー誤差の偏微分 参考 今回はモデル作成まで。モデルの評価は別記事で。 計算式等 ロジスティック回帰は二値分類（）で用いられる分析手法。モデル式は回帰分析と…

重回帰分析について。 【目次】 計算式等 ノート プログラムコード 偏回帰係数の詳細 参考 単回帰分析については次の記事。 cochineal19.hatenablog.com 計算式等 重回帰分析は、多変量 { } の関係を直線で求める分析方法。モデルは次のとおり。 ■ 重回帰モ…