カイ二乗検定のサンプルサイズ計算を SPSS で行う方法 拡張機能のインストール 拡張機能 → 拡張ハブ を選択 検索窓に power と入力して、Power Analysis for Chi-Square Models を見つけて、拡張の取得にチェックを入れ、OK をクリック 規約に同意して、イ...

カイ二乗検定のためのサンプルサイズ計算の方法 パッケージのインストール R の pwr パッケージをインストールして使用する install.packages('pwr') # 1 回だけインストール library(pwr) サンプルサイズ計算をする分割表の例 サンプルサイズ計算をする分...

多重代入法で欠測値を補完し、Fine-Gray 回帰 を実施する方法 サンプルデータ 競合リスクを含むイベント変数、時間変数、着目したい変数、交絡因子を以下のとおりとする 競合リスクを含むイベント変数：CompRisk 時間変数：DaysFPS 着目したい変数：PS34 ...

Mac で EZR を使い始めるときに、EZR プラグインを当てて起動する方法 Mac で EZR を起動するために EZR プラグインの適用 EZR のインストールが終了したあと、ツール → Rcmdr プラグインのロード を選択 RcmdrPlugin.EZR が選ばれていることを確認して OK...

因子分析の推定法にはどんなものがあり、どう使い分けるか 因子分析の推定法 因子分析の推定法には、主に、最尤法、最小二乗法、重み付き最小二乗法の 3 つがある 因子分析の推定法 最尤法 最尤法は、精度よく推定できる サンプルサイズが十分大きいときに...

ロジスティック回帰の説明変数の数はいくつ入れてよいのか？ ロジスティック回帰の説明変数の数 ロジスティック回帰分析をする場合、サンプルサイズによっていくつの説明変数を投入できるか サンプルサイズ諸論 の記述がもっともクリアでわかりやすい 二項...

説明変数を対数変換した場合の重回帰分析の結果の書き方 重回帰分析の結果の書き方 重回帰分析をはじめ、様々な回帰分析の結果の書き方を見比べられる、リンク集はこちら 重回帰分析の結果の書き方 対数変換した説明変数がある場合どう書くか 例 1 ： ( lo...

IPTW で重み付けした競合リスク累積発生曲線を書く方法 IPTW 競合リスク累積発生率曲線の書き方の前に IPTW、競合リスク、累積発生曲線とは？ IPTW に関しては以下の過去記事を参照 競合リスクに関しては以下の過去記事を参照 累積発生率（累積イベント）...

SPSS で 3 群比較を二項検定を繰り返して行う方法 SPSS で 3 群比較する方法 多重比較 ここで言う、3 群比較＆多重比較とは、3 群が同じ割合かどうかを、2 群ずつ確認する方法 3 群のうち、2 群のみを取り出し、二項検定を繰り返す 4 群以上も同様に行える...

SPSS で 3 群が同じ割合かどうかを検定する方法 SPSS で 3 群比較する方法 検定の種類 ここで言う、3 群比較とは、3 群の割合が同じかどうかを検討することである すべての割合が同じと仮定した分布への適合度検定を行う ちなみに、4 群以上でも同様である...

SPSS で 2 群の比較をする方法の解説 SPSS で 2 群比較する方法 ここで言う 2 群比較というのは、2 群の割合の比較のことであり、半々かどうかの検定という意味である つまり、二項分布を使った二項検定のことである まず、分析 → ノンパラメトリック検定 ...

反復測定データを用いて、時点ごとの中央値を計算した後、折れ線グラフを書く方法 中央値の折れ線グラフの前に 平均値の折れ線グラフ 中央値の折れ線グラフの前に、平均値の折れ線グラフを示す ここでは、平均値の折れ線グラフとは、以下のようなグラフを...

EZR で反復測定データの時点ごとの平均値を結んだ折れ線グラフを書く方法の解説。 EZR で平均値の折れ線グラフを書くデータの準備 まず、下記のようなデータを用意する。 一人一行で、列方向に反復測定データを並べたデータである。 EZR で平均値の折れ線...

SPSS で連続データを区切ってカテゴリカルデータを作る方法を解説 SPSS で連続データをカテゴリカルデータにする方法 例 1 ：変数の計算 年齢を 10 歳の区切りでカテゴリカルデータにしたい場合 10 で割って小数点以下を切り捨てるという方法を使う まず、...

割合のメタアナリシスの方法の解説 割合のメタアナリシスデータ準備 R で割合のメタアナリシスを行う準備 metafor と meta というパッケージをインストールして準備する install.packages(c("metafor", "meta")) library(metafor) library(meta) サンプル...

SPSS でグループごとの分析を行う場合、どのようにしたらよいか？ SPSS でグループごとの分析を行う方法 1： ケースの選択 分析するデータを読み込んだ後、以下の手順で設定する データ → ケースの選択 をクリック 例えば性別 gender が 1 のグループだけ...

SPSS で ROC 曲線を群間比較する方法 2 つの検査値を比較する方法も紹介 ROC 曲線を比較する方法を解説するデータ アウトカムデータ outcome, 検査値データ s100b, 性別データ gender を含むデータ こちらからダウンロードして試してみることが可能 aSAH.s...

EZR で ROC 曲線を群間比較する方法 ROC 曲線を比較するためのサンプルデータ R の pROC パッケージに含まれる aSAH データセットを使う outcome データで s100b という検査値のカットオフ値を求める gender で female と male の群間比較を行う pROC パッ...

カイ二乗検定は、2つのカテゴリカル変数の間に統計的に有意な関連があるかどうかを判断するのに役立つ手法である。しかし、統計的有意性だけでは、その関連性の「強さ」についてはわからない。そこで重要になるのが効果量である。効果量は、群間の差や変数...

サンプルの連続データが正規分布していない場合に、平均値の 95 % 信頼区間を計算する方法 正規分布していない連続データが対数正規分布だった場合 正規分布していない連続データを対数変換すると正規分布に見える場合、その連続データは対数正規分布のデ...

ステップワイズ変数選択の際に ある変数だけは残っていてほしい場合の指定方法 ステップワイズ変数選択を EZR で実行する方法 ステップワイズ変数選択を EZR で実行するには、データを読み込んだ後、例えば重回帰のメニューからステップワイズを選択する ...

対数変換した後の数値の要約値の論文表への書き方 実際、どのように記載されているか 幾何平均を書いている論文 Does the association of the triglyceride to high-density lipoprotein cholesterol ratio with fasting serum insulin differ by race/eth...

分散分析のサンプルサイズ計算をSPSSで行う方法 分散分析のサンプルサイズをSPSSで計算する前に効果量を見積もる 分散分析のサンプルサイズをSPSSで計算する際に、必要なのは効果量を見積もること 先行研究の分散分析表を使うとよい 今回は tension の高さ...

重回帰分析で、独立変数にカテゴリ変数を使う方法を解説する。 SPSSでは、ダミー変数を作成しておく必要がある。 ダミー変数とは何か？ ダミー変数の作り方は？ 重回帰分析をSPSSで行う方法 重回帰分析をSPSSで行う場合、２つの方法がある。 ひとつは、「...

IPTW ログランク検定をEZRとRを使って行う方法。 IPTWとは？ IPTW とは、Inverse Probability of Treatment Weights の頭文字語。 日本語では、逆確率重みづけと言う。 交絡因子調整方法の一つ。 詳しくは他の記事も参照のこと。 IPTW ログランク検定とは...

Rでヒストグラムと箱ひげ図をggplot2 を使って重ねて描く方法の紹介。 Rでヒストグラムをggplot2 を使って描く方法 Rでヒストグラムをggplot2を使って描く方法を紹介する。 まずggplot2 パッケージをインストールして、呼び出しておく。 install.packages(...

カイ二乗検定をわかりやすく解説。 エクセルでの計算方法も紹介。 カイ二乗検定とは？ カイ二乗検定とは、サンプルのグループ間の「何か」の割合を比較することで、グループの母集団の割合が統計学的に違うかどうかを調べる方法である。 グループ間で何か...

主成分分析をカテゴリ変数が含まれるデータセットで実施する方法を紹介する。 主成分分析をカテゴリ変数で実施可能か？ カテゴリ変数でも主成分分析はできるか？ 答えはYesである。 SPSSで実施可能である。 ただし、追加のパッケージ (Complex Sample) が...

Brain Network Analysisでよく使われている手法について理解したところをまとめてみる。 Brain Network Analysisで登場する用語 Brain Network Analysisで登場する用語は、一般的なネットワーク分析で登場する用語と被るものが多い。 ネットワーク分析の基...

EZR で箱ひげ図を書くにはどうしたらよいか？ 色分けするにはどうしたらよいか？ EZR で箱ひげ図を書く方法 データを読み込んだ後、グラフと表から箱ひげ図メニューを選択する。 変数のほうで連続データを指定し、群別する変数の枠で群を表すカテゴリカル...

ネットワーク分析の基礎知識として用語の簡単な解説。 ネットワーク分析とは？ ネットワーク分析とは、社会学や通信ネットワークなどの分野で用いられる分析手法である。 脳画像を用いた脳内のネットワーク分析も行われている。 このネットワーク分析は、...

主成分分析で軸の回転をすることはできるか 主成分分析とは？ 主成分分析とは、多数の変数の情報を少数の合成変数に縮約する解析手法である。 たくさん変数があると何がどうなっているかわからないが、それら変数同士の相関関係を活用して、代表的な方法と...

MMRMは、Mixed Model with Repeated Measuresの略である。 反復測定の連続データを線形混合モデルで解析することを意味している。 この試験・研究デザインのとき、サンプルサイズ計算はどのようにしたらよいか？ MMRMのサンプルサイズ計算サイト 以下のサ...

正規性の確認は、検定は不要で、グラフで確認すればよい。 ヒストグラムが一番簡単な方法だが、もう一つグラフで確認する方法がある。 それがQQプロットだ。 QQプロットとは何か？ QQプロットの前に 正規性とは？ 統計解析の検定において正規性が条件であ...

R におけるネットワーク分析のパッケージ network の簡単な使い方。 network パッケージのインストール まず統計ソフトRにパッケージをインストールして呼び出す。 install.packages("network") library(network) 準備するデータ ネットワークのつながりあ...

対数変換した説明変数の単位数当たりのハザード比を、真数に戻した説明変数の単位数当たりのハザード比に計算できるだろうか？ 論文で発表された数値を使うメタアナリシスを行うことはできるか？ つまり、生データがない場合である。 どういう状況か？ 例...

正準相関分析とは、どんな分析で、どんなときに使うのか？ 正準相関分析とはどんな分析か 正準相関分析 canonical correlation analysis は、一言で言うと、２つの変数群を２つの合成変数にまとめ、その合成変数間の相関係数を最大にする方法である。 正準...

線形判別分析とは？ 数量化II類とは？ 違いは何か？ 線形判別分析とは？ 線形判別分析とは、いくつかの変数で、群分けできる「判別関数」という係数の組み合わせを求めることである。 いくつかの変数と最適な係数の掛け合わせ＋足し合わせで、合成変数を作...

前後2時点の対応のあるデータを群間比較するときは、共分散分析を使う EZR と SPSS の方法を解説 対応のあるデータの群間比較の考え方 対応のあるデータの群間比較の際に、考え方が２つある。 一つは、後のデータを目的変数に、群別変数を着目する説明変数...

重回帰分析に交互作用を含んだ解析をSPSSで計算する方法。 SPSSの重回帰分析では、交互作用項を自動で作ってくれない。 どのように作って、どのように解析すればよいか？ SPSSで交互作用項の作り方 重回帰モデル内で交互作用の検討をするためには、交互作...

ケンドールの順位相関係数はどのように計算するか紹介する ケンドールの順位相関係数 ケンドール（Kendall）の順位相関係数 $ \tau $ は、順位を使わない相関係数である。 何を使うかというと、値の大小を使う。 変数Xと変数Yの相関を計算する際に、2つのX...

IPTW カプランマイヤー曲線の書き方。 RのパッケージとEZRを使うと書ける。 IPTWとは？ IPTW とは、Inverse Probability of Treatment Weights の頭文字語。 逆確率重み付けの意味。 傾向スコア (Propensity Score) を使った、交絡要因調整の一種。 処方さ...

EZRでグラフが出力されたときに、余白が狭くて軸の文字が見えないことがある。 余白を調整するにはどうしたらよいか？ EZRでグラフの余白調整が必要な例 もっとも典型的なグラフは、分散分析でTukeyのオプションを入れたときに出力されるグラフである。 左...

サブグループ解析の結果をフォレストプロットで書く方法。 ロジスティック回帰の場合で、EZR で行う方法。 Version 1.62（2023 年 11 月 1 日公開）から、EZR 本体に組み込まれたので、以下の方法を使う必要なし 具体的な方法はこちらの動画参照。 フォレ...

重回帰分析は偏回帰係数を求めて予測式（回帰式）を作り、標準化偏回帰係数の計算、回帰式の有意性の分散分析、決定係数の算出、偏回帰係数の検定、など行うが、実際どんな計算をしているのか？ 実際どんな計算をしているかわかると、偏回帰係数や標準化偏...

Dunn検定と言えば、順位和を使った検定が有名である。 これはSPSSでクラスカルウォリス検定をしたあとのPost-hoc検定として有名。 どんな検定なのか？ Dunn検定の明確化 Olive Jean Dunn先生は、1961年と1964年に二つのDunn検定を発表している。 1961年のD...

クラスカルウォリス検定を、原著に立ち返って、わかりやすく解説。 各群のサンプル数が偏ってはいけないなどの条件があるのかどうか、また、タイデータ（同順位データ）の影響について。 クラスカルウォリス検定とは？ クラスカルウォリス検定とは、３群以...

主成分分析は、多変量情報の縮約と言われるが、実際にはどんな計算をしているのか？ 数学的に少し詳しくわかりたい人向け。 主成分分析の計算上の目標は合成変数の作成と分散の最大化 主成分分析の数学的な計算の目的は、合成変数の作成と、その合成変数の...

「共分散構造分析（SEM）は、なんだか難しそう…」と感じていないだろうか。そして、「どれくらいのデータがあれば、この分析ができるのだろうか？」と悩んだことはないだろうか。論文や研究でSEMを使いたいけれど、統計や数学は苦手だという方、安心してほ...

目的変数が正規分布している必要はない。 説明変数も正規分布している必要はない。 前提知識 実際に測定された目的変数を実測値と言う。 回帰式で計算された目的変数を予測値と言う。 実測値と予測値の差を残差と言う。 正式には母集団の誤差項の話になる...

カーネル密度推定とはどんな推定で、どんな時に使われるのか？ バンド幅の決め方も紹介 カーネル密度推定とは何か？ とても分かりやすい記事を見つけたので、そのまま引用する。 さて、カーネル密度推定法を大雑把に言うと、「確率密度が全然わからない確...

繰り返し計測された連続データを適切に扱うために必要な線形混合モデル。 ランダム切片モデルが一番単純で有名だが、その切片はどんなふうになっているか。 通常の線形モデルの切片はどんな計算になっているか？ 例えば、連続データ CDratio1 を 3 群のFCZ...

Jonckheere-Terpstra 検定（ヨンクヒール タプストラ検定、ヨンキー検定）とは？ Jonckheere-Terpstra検定とは？ ３群以上の群で、単調に増加するか、単調に減少するかなどの傾向を検定するノンパラメトリックな検定法。 Jonckheere-Terpstra 検定は、デー...

マンホイットニーのU検定を行う際に最低限何例あればよいか？ マンホイットニーのU検定はどんな計算をしているか？ マンホイットニーのU検定はそもそもどんな計算をしているか、という点から、各群最低何例必要かを考えたい。 数学的には同値になる、方法...

相関係数の比較はどうやるか？ 相関係数の差の検定とは？ 相関係数の比較　相関係数の差の検定 （この記事で言う相関係数は、すべてPearsonの積率相関係数を指している） 相関係数の差の検定は、以下の検定統計量 T を計算し、標準正規分布で有意確率を求...

EZRの出力をExcelにきれいに貼り付ける方法。 標準の機能を使う方法と、テキストファイルウィザードを使う方法を紹介する。 ベースラインの患者背景の表 「グラフと表」メニュー中の「サンプルの背景データのサマリー表の出力」を使うと、ベースラインの患...

重回帰分析を実施するのにサンプル数はどのくらい必要か？ ここでいうサンプル数とは、サンプルサイズのこと。 重回帰分析のサンプルサイズはどのように計算するのか？ 重回帰分析の効果量とは何か？ 重回帰分析のサンプル数計算に必要な値は？ サンプル数...

連続データの代表的な要約値、平均値と標準偏差を統合するにはどうやったらよいか？ 連続データを層別で要約値を求めた後に統合する方法。 連続データの層別解析の統合とは？ 目的変数としての連続データに対して、交絡を調整しながら分析するには、共分散...

層別ログランク検定を SPSS と EZR で行う方法を解説する 層別ログランク検定の使い分け 生存時間データで交絡を調整したいということなら、コックス比例ハザードモデルを使えばよいのだが、なんらかの理由で層別解析＆統合をしたい場合、どうしたらよいか...

マンテル・ヘンツェル検定は、層別の分割表を統合する方法だ。 SPSSではどのように行うか。 SPSSでマンテル・ヘンツェル検定を行う方法 分析→記述統計→クロス集計表 と進む。 行、列、層にそれぞれ、説明変数、目的変数、層別変数を投入する。 説明変数と...