正規分布近似を用いた信頼区間とClopper-Pearson 信頼区間を計算する方法。 割合の区間推定（正規分布近似を用いた簡易的な方法） 割合の分散は、割合をp、サンプルサイズをnとすると、 $$ \frac{p (1 - p)}{n} $$ で計算できる。 標準誤差はその平方根で...

R はパッケージを追加すると新しい機能が追加できる。 その方法もとても簡単だ。 R のパッケージの追加方法として source からインストールする方法があるので、その解説。 R のパッケージを source からインストールする方法 R のパッケージは、Windows用...

R で時系列データを分析するための基礎的な方法の解説。 時系列データとは何か？ 時系列データとは、日時のデータと値のデータのペアのこと。 例としては、 四半期ごとの売上高。 毎月の経費。 週当たりの納入量。 毎日の株価・ドル円レート。 など。 こう...

多重比較のサンプルサイズ計算を R で行う方法 ボンフェローニ型 多重比較のサンプルサイズ計算 Bonferroni（ボンフェローニ）型多重比較とは、比較する数で有意水準を割って、割った有意水準より小さい有意確率の場合、統計学的有意と考える方法。 三群あ...

がん患者さんの生存期間中央値を5年生存率から逆算で求めるにはどうしたらよいか？ R で計算する方法。 5 年生存率から逆算で生存期間中央値を求める計算式 生存確率は、一般論として指数関数で近似できる。 実際には、がんの治療後すぐになくなったり、数...

R に組み込みのタイタニック号の生存・死亡データで、生存・死亡に関するロジスティック回帰分析と、独立変数の多重比較を実行してみる。 解析するタイタニック号の生存・死亡データの確認 解析するデータは、carData パッケージの TitanicSurvival という...

タイタニック号は、1912年4月14日の夜、氷山に激突し、北大西洋の底に1,500名以上の命と一緒に沈んだ。 乗客乗員の生存・死亡のデータを用いて、ロジスティック回帰分析を実行してみる。 タイタニック号の生存・死亡データはどこにあるのか？ R はインスト...

偏相関係数について。 順位相関係数は偏相関係数が計算できるか。 偏相関係数が計算したい場面 新規の検査値Aと既存の検査値Bの関連性を検討したい。 nは40例程度。 解析するに当たり、患者背景が問題になると思っている。 既存の検査値 Bと身長が関連する...

相関係数は、相関関係の強さを示す指標。 一方が大きいときにもう一方が大きければ、正の相関関係で、相関係数は1に近い。 一方が大きいときにもう一方が小さい場合は、負の相関関係で相関係数は－1に近い。 では、偏相関係数とは何か？ 違いは何か？ 偏相...

回帰分析をする際に、説明変数や目的変数が正規分布をしていないことで悩んでいる人は多い。 悩むところはそこじゃない。 重回帰分析では、残差が正規分布している必要がある。 重回帰分析の前提は何か？ 重回帰分析の前提は4つある。 独立性（データそれ...