日本一簡単な医療統計を目指すblog

筆者の覚書のためにできるだけ簡単に統計が使えるようになることを目指すブログです。そういう意味では細かい統計事項には触れませんし、統計家がかいているわけではありませんので間違いもあると思います。

Rを使って無料でとにかく最速で累積発生率(0から上に階段状に上がるKaplan-Meier曲線)を作成する方法

 

前回は上から下に階段状に落ちるKaplan-Meier曲線

1 カプランマイヤーの元データをCSVでデータ保存

デスクトップにエクセルを使ってCSVでデータを保存します。

データは英語で1行目にMonth Survival  Groupとして、Survivalは生存・打ち切りが0,死亡が1でデータを用意します。

Kaplan Meier 元データ

カプランマイヤーデータ

2 Rのディレクトリーの設定

Rを立ち上げて、ディレクトリーをデスクトップにします。

 

f:id:USMLE:20201101114544p:plain

Rでカプランマイヤー

3 以下をコピペしてRに打ち込めばとりあえずはKMカーブが完成。

library(survival)

x <- read.csv("test.csv")

survfit1 <- survfit(Surv(month, survival)~group, data=x)

plot(survfit1, fun="event", col=1:2)

 

Rで作成 累積発生率(下から上へのKaplan-Meier曲線)

Rで作成 累積発生率(下から上へのKaplan-Meier曲線)

 赤字が前回との違いです。FUNのところで”EVENT”を指定すると下から上への累積での積み上げ表示になります。COLは色を指定しています。

ちなみに前回の記事はコチラ

easystats.hatenablog.com

ちなみに最初のlibrary(survival)はsurivalというライブラリー(プログラムセット)を読み込むコマンドです。

 

無料ソフトEZRを使用するともっと楽にできます。

 EZRの詳しい説明はコチラ

 

EZRでやさしく学ぶ統計学~EBMの実践から臨床研究まで~ 2版

EZRでやさしく学ぶ統計学~EBMの実践から臨床研究まで~ 2版

  • 作者:神田善伸
  • 発売日: 2015/04/17
  • メディア: 単行本(ソフトカバー)
 

 

 

エクセルを使ってLog-rank testで二つのKaplan-Meier Curveを比較する。

 

Kaplan-Meier曲線を描いてLog-rank testで評価

前回までの内容はコチラからどうぞ

easystats.hatenablog.com

 

easystats.hatenablog.com

Log-rank testのためのデータをExcelで集める 

ここからは2群の

その場合にはGroupの列(column)を用意します。

そしてLog-rank testを行えば有意差検定が行えます。

この前提としてはHazard性が保たれていることが必要条件ですが、通常がんなどのsurvival dataではKaplan-MeierとLog-rank testを用いることが許されています。

 

Case No Group Months from surgery Event
1 0 5 0
2 1 12 1
3 0 14 0
4 0 14 1
5 1 24 1
6 1 30 0
7 0 60 0

エクセルでカプランマイヤーを実際に描く

このあとにエクセルでカプランマイヤーを描きログランクテストで有意差検定をするのであれば、エクセル統計がいいでしょう。エクセル統計はインストールするとエクセルにプラグインされますので、ほとんどエクセルの操作感覚で論文や学会に使えるKaplan Meierが完成します。

操作としては上の表のGroupからEventまでの3列を選んでボタンを押すだけです。

 

 

エクセル統計 通常版1年1台

エクセル統計 通常版1年1台

  • 発売日: 2018/06/06
  • メディア: CD-ROM
 

 これがあれば、あっという間にカプランマイヤーカーブを完成できます。

 

エクセル統計のカプランマイヤーカーブ、ログランクテストの実際のサンプル画面

実際にはエクセルでのデータの入力方法は2パターンあります。上述の方法ともう一つの方法が記載されていますので参照してみてください。カプランマイヤーカーブのサンプルイメージも見られます。

bellcurve.jp

 

Log-rankもsurvival packageを使えばできますが、簡単にはEZRでしょう。

もう一つの解析方法はWilcoxonです。

 

EZRでやさしく学ぶ統計学~EBMの実践から臨床研究まで~ 2版

EZRでやさしく学ぶ統計学~EBMの実践から臨床研究まで~ 2版

  • 作者:神田善伸
  • 発売日: 2015/04/17
  • メディア: 単行本(ソフトカバー)
 

 

ゼロから始めて一冊でわかる!みんなのEBMと臨床研究

ゼロから始めて一冊でわかる!みんなのEBMと臨床研究

  • 作者:神田 善伸
  • 発売日: 2016/10/17
  • メディア: 単行本
 

 

Rを使って無料でとにかく最速でKaplan-Meier曲線を作成する方法

 

1 カプランマイヤーの元データをCSVでデータ保存

デスクトップにエクセルを使ってCSVでデータを保存します。

データは英語で1行目にMonth Survival  Groupとして、Survivalは生存・打ち切りが0,死亡が1でデータを用意します。

Kaplan Meier 元データ

カプランマイヤーデータ

2 Rのディレクトリーの設定

Rを立ち上げて、ディレクトリーをデスクトップにします。

 

f:id:USMLE:20201101114544p:plain

Rでカプランマイヤー

3 以下をコピペしてRに打ち込めばとりあえずはKMカーブが完成。

library(survival)

x <- read.csv("test.csv")

survfit1 <- survfit(Surv(month, survival)~group, data=x)

plot(survfit1,lty=1:2,mark.t=T)

 

R カプランマイヤーの最も簡単な描き方

R カプランマイヤー

ちなみに最初のlibrary(survival)はsurivalというライブラリー(プログラムセット)を読み込むコマンドです。

完成 カプランマイヤー

完成 カプランマイヤー

 

無料ソフトEZRを使用するともっと楽にできます。

 EZRの詳しい説明はコチラ

 

EZRでやさしく学ぶ統計学~EBMの実践から臨床研究まで~ 2版

EZRでやさしく学ぶ統計学~EBMの実践から臨床研究まで~ 2版

  • 作者:神田善伸
  • 発売日: 2015/04/17
  • メディア: 単行本(ソフトカバー)
 

 

 

英文論文の書き方 文献ソフト Endnoteエンドノートの簡単な使い方2

 

はじめにー論文の文献に関して

英文論文を書くのはなかなか大変です。

そんな際にMUSTなのが、Endnoteです。今回はEndnoteX9を使用しますが、他のバージョンでもほとんど以下の方法が使えると思います。

 

 文献の付け方は雑誌によって様々ですが、このEndnoteがあれば、簡単に雑誌に合わせた体裁に変換することができ、大幅な時間削減に役立ちます。

 

前回までは文献リストの作り方などを説明しました。

easystats.hatenablog.com

リファレンスの形式の整え方

文献リストの形式は雑誌によって異なります。

Endnoteで有名な雑誌の文献リストを選ぶ

ある程度有名な雑誌の場合にはEndnote側ですでに形式を用意していますので、簡単です。

Edit>output styleからopen style managerを選んで雑誌を選びます。

Endnoteで新たにスタイルを作成する

Edit>output styleからnew styleを選んで作成します。

しかしかなり大変です。

そこで、

Endnoteで他のスタイルを選び、変更する

のが楽です。

以下の図のように

Edit>output styleからここではEdit"Vancouver"を選びます。

Endnote 文献形式の選び方

Endnote 文献形式の選び方

Endnoteエンドノートで本文の文献サイテーションの形式を整える

EndnoteでEditをおすと以下のような画面に移る。

Citation >Templatesを選んだ画面が以下の通り。

Vancouverでは論文の番号が()で囲まれた形になる。

もう一つ多くあるパターンは上付きであるが、この形式にしたければ

以下の画面のBibliography Numberを選択したまま図のまるで囲まれたボタンを押せば上付きになる。

項目を新たに追加するには Insert Fieldから選ぶ。

Endnote 本文中の文献サイテーションの形式変更

Endnote 本文中の文献サイテーションの形式変更

 

Endnoteエンドノートで文献リストの形式を変更する

  1. Edit"Vancouver"からBibliography>Templatesを選ぶ
  2.  一番多いと思われる雑誌の形式を例に下の図に示す。
  3. 右側のInsertから項目を選び、順番やつなぎの記号(カンマ、コロンなど)を選ぶ。

Endnote 文献リストの形式の変更方法

Endnote 文献リストの形式の変更方法


Endnoteエンドノートで文献リストの著者数(Author数)を変更する

  1. 下図の左側でCitation>Author Listsを選ぶ
  2. 選べる項目はわかりやすいので直感的に最適化できる。

 

  1. Endnote 著者数を変更する

    Endnote 著者数を変更する

    変更を新しいスタイルで保存

 最後に保存を忘れずに

FILE>SAVE ASで好きな名前に変更します。

 

私は読んだことはありませんが、日本語の解説本もあるようです。

 

最新EndNote活用ガイド デジタル文献整理術 第7版

最新EndNote活用ガイド デジタル文献整理術 第7版

  • 作者:讃岐 美智義
  • 発売日: 2018/02/15
  • メディア: 単行本(ソフトカバー)
 

 

 

英文論文の書き方 文献ソフト Endnoteエンドノートの簡単な使い方1

 

はじめにー論文の文献に関して

英文論文を書くのはなかなか大変です。

そんな際にMUSTなのが、Endnoteです。今回はEndnoteX9を使用しますが、他のバージョンでもほとんど以下の方法が使えると思います。

 

 文献の付け方は雑誌によって様々ですが、このEndnoteがあれば、簡単に雑誌に合わせた体裁に変換することができ、大幅な時間削減に役立ちます。

 

なぜ文献が重要か?

ちなみに文献は適当でよいと思っている人も多いかと思いますが、

Reviewerのチェック項目に大事な文献が入っているかという項目があります。

少しでもインパクトファクターの高い雑誌を目指すのであれば、重要文献を必ず入れる必要があります。

文献が必要なのは主にIntroductionとDiscussionです。

Endnoteエンドノートで文献リストを簡単に作成する方法

Endnoteでまずは文献リストを作成する。

下の図のように

  1. 地球のマークをクリックしてEndnoteエンドノートをOnline Search modeにする
  2. Pubmed(NLM)を選択
  3. 検索する。後述するようにPMIDから選ぶと便利
  4. 文献を右クリックしてCopy reference toから使用しているライブラリーにコピーする。

最も簡単なEndnoteの使い方1

最も簡単なEndnoteの使い方1



Endnoteエンドノートでワードに最も簡単に文献をひもづける方法

  1. EndnoteエンドノートをLocal library modeにするために図のフォルダーマークをクリックする
  2.  使用する文献を右クリックしてCopyする。
  3. ワードに移り、サイトする箇所にペーストする。

最も簡単なEndnoteの使い方2

最も簡単なEndnoteの使い方2

最も簡単なEndnoteの使い方3

最も簡単なEndnoteの使い方3

ワードで英文論文に文献リストを作成する。

  1. 下図の丸で囲った部分を投稿する雑誌(または同じ体裁の雑誌)を選ぶ。
  2. そのすぐしたのUpdate Citations and Bibliographyをおすと、サイト部分は数字に変わり、文献リストが作成される。

最も簡単なEndnoteの使い方4

最も簡単なEndnoteの使い方4

 

 

私は読んだことはありませんが、日本語の解説本もあるようです。

 

最新EndNote活用ガイド デジタル文献整理術 第7版

最新EndNote活用ガイド デジタル文献整理術 第7版

  • 作者:讃岐 美智義
  • 発売日: 2018/02/15
  • メディア: 単行本(ソフトカバー)
 

 

次回、第2回は、文献リストの体裁の整え方の予定です。

Cox-regressionで多変量解析をする方法

 

Cox regressionで多変量解析

Cox regressionは時間が関係するアウトカム(生存率、再発率など)に関する多変量解析の方法です。今回はそのデータの集め方を説明します。

Kaplan-Meier曲線を描くのに必要な用語(復習)

Kaplan-Meier曲線は例えば癌の生存率曲線などを描くときに必要です。

その際に必要なデータは何か?一人につき最低2つの情報が必要です。

 

生死、そして期間です。

イベントEvent (生死)

まずは、患者さんが生存しているかどうかというデータです。

通常は死亡しているとイベントEventといい、たとえばEZRでは1を割り当てます*1

生存しているか、情報が不明の場合は0です。

 期間

それと手術から死亡日までの日数、または生存が確認されている最終日までの日数が必要です。

手術日と最終生存確認日(または死亡日)の二つのデータを引き算してもよいでしょう。

打ち切り

例えば引越しなどである時期までは生存していて、その後不明というときには

最終生存確認日をいれて、それまでは生きていたというデータにします。これを打ち切りといいます。

Kaplan-Meierカーブのためのデータをエクセルで集める

Case No Months from surgery Event
1 5 0
2 12 1
3 14 0
4 14 1
5 24 1
6 30 0
7 60 0

 


 1人目は経過5ヶ月で生存中

2人目は経過12ヶ月で死亡

3人目は14ヶ月に外来にきてその後引っ越して状況不明と行った場合もEventは0です

 

これをRで読み込んでsurvival packageを使えばタダでKMカーブがかけます。

 

もっと簡単にはEZRでかけます。

英語の説明はhttp://www.jichi.ac.jp/saitama-sct/SaitamaHP.files/manual.html

EZRでやさしく学ぶ統計学~EBMの実践から臨床研究まで~ 2版

EZRでやさしく学ぶ統計学~EBMの実践から臨床研究まで~ 2版

  • 作者:神田善伸
  • 発売日: 2015/04/17
  • メディア: 単行本(ソフトカバー)
 

Cox regression testのためのデータをExcelで集める 

ここからはいろいろな予後に関係しそうなFactorを入力します。

今回の場合にはSex、Smokingの列(column)を用意します。

それぞれの単変量はLog-rank testを行えば有意差検定が行えます。

easystats.hatenablog.com

この前提としてはHazard性が保たれていることが必要条件ですが、通常がんなどのsurvival dataではKaplan-MeierとLog-rank testを用いることが許されています。

そして多変量を行うのがCox regressionです。

 

Case No Sex Smoking Months from surgery Event
1 0 1 5 0
2 1 0 12 1
3 0 1 14 0
4 1 1 14 1
5 1 0 24 1
6 1 0 30 0
7 0 1 60 0


Cox-regressionもRを使えばできますが、簡単にはEZRでしょう。 

 

EZRでやさしく学ぶ統計学~EBMの実践から臨床研究まで~ 2版

EZRでやさしく学ぶ統計学~EBMの実践から臨床研究まで~ 2版

  • 作者:神田善伸
  • 発売日: 2015/04/17
  • メディア: 単行本(ソフトカバー)
 

 

ゼロから始めて一冊でわかる!みんなのEBMと臨床研究

ゼロから始めて一冊でわかる!みんなのEBMと臨床研究

  • 作者:神田 善伸
  • 発売日: 2016/10/17
  • メディア: 単行本
 

 

Kaplan-Meier:エクセルの打ち切り・イベントの入力をすべて逆にしてしまったときの対処法

エクセルで打ち切り・イベントを逆に入力してしまった、、、

前回はKaplan-Meierのデータの集め方を説明しました。

easystats.hatenablog.com

Kaplan-Meier曲線を描いてみたら、、、

Kaplan-Meierのデータを描いてみたら、生存率0%に墜落したグラフが出ることがあります。全員死亡していないはずなのに、、、というときはイベントと打ち切りの数字が逆になっていることがあります。

ソフトによっては打ち切りを1にするまたは0にするといったチェックを入れるだけで解決するものもありますが、そうはいかないソフトもあります。

 

なぜあなたは論文が書けないのか?

なぜあなたは論文が書けないのか?

  • 作者:佐藤 雅昭
  • 発売日: 2016/07/08
  • メディア: 単行本
 

 

エクセルで簡単に0を1、1を0にする方法

イベントと打ち切りはそれぞれ片方が1であればもう一方は0という関係です。これをすべて入力し直すのは間違える可能性もありますし、時間もかかります。

今回は打ち切りCensorが必要だったのにイベント(死亡)を1として入力してしまった場合を想定します。C列イベントが死亡の際に1となってしまっていたときを想定します。

1からもう一方を引けば、0と1が入れ替わりますので、

f:id:USMLE:20200518081331p:plain

DにCensor列を作り計算式で1からC列を引き算します。

上のように入力して下の図の蛍光ペンがついたポッチをダブルクリックすると

f:id:USMLE:20200518081701p:plain

ダブルクリックすると左側のC列が埋まっていれば自動計算されてD列が自動入力される

f:id:USMLE:20200518081930p:plain

数値に変えたければコピーの後、形式を選択して数値のみペーストする。

鬼滅の刃、人気ですが連載終了らしいですね。

 

 

鬼滅の刃 20 (ジャンプコミックスDIGITAL)

鬼滅の刃 20 (ジャンプコミックスDIGITAL)

 

 

 

Log-rank testで二つのKaplan-Meier Curveを比較する。

 

Kaplan-Meier曲線を描いてLog-rank testで評価

Kaplan-Meier曲線を描くのに必要な用語

Kaplan-Meier曲線は例えば癌の生存率曲線などを描くときに必要です。

その際に必要なデータは何か?一人につき最低2つの情報が必要です。

 

生死、そして期間です。

イベントEvent (生死)

まずは、患者さんが生存しているかどうかというデータです。

通常は死亡しているとイベントEventといい、たとえばEZRでは1を割り当てます*1

生存しているか、情報が不明の場合は0です。

 期間

それと手術から死亡日までの日数、または生存が確認されている最終日までの日数が必要です。

手術日と最終生存確認日(または死亡日)の二つのデータを引き算してもよいでしょう。

打ち切り

例えば引越しなどである時期までは生存していて、その後不明というときには

最終生存確認日をいれて、それまでは生きていたというデータにします。これを打ち切りといいます。

Kaplan-Meierのためのデータをエクセルで集める

 

今回は実際に描いてみたいと思います。

仮にデータを以下のように用意します。

 

 

Case No Months from surgery Event
1 5 0
2 12 1
3 14 0
4 14 1
5 24 1
6 30 0
7 60 0

 1人目は経過5ヶ月で生存中

2人目は経過12ヶ月で死亡

3人目は14ヶ月に外来にきてその後引っ越して状況不明と行った場合もEventは0です

 

これをRで読み込んでsurvival packageを使えばタダでKMカーブがかけます。

 

もっと簡単にはEZRでかけます。

英語の説明はhttp://www.jichi.ac.jp/saitama-sct/SaitamaHP.files/manual.html

EZRでやさしく学ぶ統計学~EBMの実践から臨床研究まで~ 2版

EZRでやさしく学ぶ統計学~EBMの実践から臨床研究まで~ 2版

  • 作者:神田善伸
  • 発売日: 2015/04/17
  • メディア: 単行本(ソフトカバー)
 

Log-rank testのためのデータをExcelで集める 

ここからは2群の

その場合にはGroupの列(column)を用意します。

そしてLog-rank testを行えば有意差検定が行えます。

この前提としてはHazard性が保たれていることが必要条件ですが、通常がんなどのsurvival dataではKaplan-MeierとLog-rank testを用いることが許されています。

 

Case No Group Months from surgery Event
1 0 5 0
2 1 12 1
3 0 14 0
4 0 14 1
5 1 24 1
6 1 30 0
7 0 60 0

 

Log-rankもsurvival packageを使えばできますが、簡単にはEZRでしょう。

もう一つの解析方法はWilcoxonです。

 

EZRでやさしく学ぶ統計学~EBMの実践から臨床研究まで~ 2版

EZRでやさしく学ぶ統計学~EBMの実践から臨床研究まで~ 2版

  • 作者:神田善伸
  • 発売日: 2015/04/17
  • メディア: 単行本(ソフトカバー)
 

 

ゼロから始めて一冊でわかる!みんなのEBMと臨床研究

ゼロから始めて一冊でわかる!みんなのEBMと臨床研究

  • 作者:神田 善伸
  • 発売日: 2016/10/17
  • メディア: 単行本
 

 

Kaplan-Meier曲線を描きたい

 

Kaplan-Meire曲線を描いてみたいと思います。

本日はPart1データ収集です。

使用するのは無料ソフトEZRです。

 

EZRの詳しい説明はコチラ

 

EZRでやさしく学ぶ統計学~EBMの実践から臨床研究まで~ 2版

EZRでやさしく学ぶ統計学~EBMの実践から臨床研究まで~ 2版

  • 作者:神田善伸
  • 発売日: 2015/04/17
  • メディア: 単行本(ソフトカバー)
 

 

Kaplan-Meier カプランマイヤー曲線に必要なデータを集める

Kaplan-Meier曲線は例えば癌の生存率曲線などを描くときに必要です。

その際に必要なデータは何か?一人につき最低2つの情報が必要です。

 

生死、そして期間です。

イベント(カプランマイヤー曲線)

まずは、患者さんが生存しているかどうかというデータです。

通常は死亡しているとイベントEventといい、たとえばEZRでは1を割り当てます*1

生存しているか、情報が不明の場合は0です。

 

それと手術から死亡日までの日数、または生存が確認されている最終日までの日数が必要です。

手術日と最終生存確認日(または死亡日)の二つのデータを引き算してもよいでしょう。

 

例えば引越しなどである時期までは生存していて、その後不明というときには

最終生存確認日をいれて、それまでは生きていたというデータにします。これを打ち切りといいます。

 

このデータがあれば、EZRで簡単にKaplan-Meierカーブを描くことができます。

EZRの詳しい説明はコチラ

 

EZRでやさしく学ぶ統計学~EBMの実践から臨床研究まで~ 2版

EZRでやさしく学ぶ統計学~EBMの実践から臨床研究まで~ 2版

  • 作者:神田善伸
  • 発売日: 2015/04/17
  • メディア: 単行本(ソフトカバー)
 

 

*1:統計ソフトによってはイベントが0というソフトもあります。

PCRは何の略

新型コロナウィルスが猛威を振るい、多くの方がお困りの状況と察します。

早く、自体が収束することを願うばかりです。

 

さて、医学英語にはいろいろな略語があります。
今回はPCRです。



PCRはポリメラーゼ連鎖反応

Polymerase (ポリメレース)

Chain

Reaction

の略です。
Chain Reaction は連鎖反応ですが、

-aseは酵素につく接尾語です。
poly-多数の、mono-は単数です。

 

home.p08.itscom.net

 

First Aid for the USMLE Step 1 2020

First Aid for the USMLE Step 1 2020

 

 

目次

統計ソフトに関して

統計ソフトは何を使う?

Rをダウンロードしよう

EasyRをダウンロードする

エクセル統計を使う

2X2の統計

カイ二乗?Fisher正確テスト?

Rで実際にFisher正確テスト

Rで実際にカイ二乗検定

Kaplan Meierを描いてみた

 

 

 

 

エクセル統計

 

統計ソフトはいろいろあってどれがいいか迷うところです。

エクセル統計

わかりやすさという面ではエクセル統計も人気である。

エクセルにプラグインするので、エクセルに入力しボタンを押せば、グラフができあがる。ただし、入力方法はある程度決められているのでそのあたりがデメリットだ。

 値段的には1年で2万円強い、4年2台なら6万円ほどで△といったところか

エクセル統計 通常版1年1台

エクセル統計 通常版1年1台

  • 発売日: 2018/06/06
  • メディア: CD-ROM
 

 

4Stepエクセル統計 第4版

4Stepエクセル統計 第4版

  • 作者:柳井久江
  • 発売日: 2015/10/20
  • メディア: 単行本
 

 

 

 

 
easystats.hatenablog.com

Rで実際にカイ二乗検定

 

Rで実際にカイ二乗検定をしてみる

こんどはRでカイ二乗検定してみます。

  160cm以上 160cm未満
男性 30 8
女性 24 12

 

データはすべてのカラムで6以上なのでカイ二乗が適応できます(以下参照)。

easystats.hatenablog.com

Rを立ち上げると以下のようなマークが出てくると思います。

細かいことはともかくとして、ここに以下のように打ち込んでみます。

 

chisq.test(rbind(c(30,8),c(24,12)))

そうすると結果が出てきます。

 

前回のFisherと比べてみましょう。

fisher.test(rbind(c(30,8),c(24,12)))

 

前回と一緒で少し解説すると

c(30,4)やc(24,12)は行列の一列分で

rbind()で2X2の行列を作っています。

「このテーブルをFisher正確検定しなさい」というコマンドです。

 

簡単にカイ二乗検定を行いたい人はコチラ

難しくてむりという人はEasyRが良いかもしれません。

 

EZRでやさしく学ぶ統計学~EBMの実践から臨床研究まで~ 2版

EZRでやさしく学ぶ統計学~EBMの実践から臨床研究まで~ 2版

  • 作者:神田善伸
  • 発売日: 2015/04/17
  • メディア: 単行本(ソフトカバー)
 

 

easystats.hatenablog.com

Rで実際にFisher正確テストを試してみる。

さて、前回のデータを実際にRで統計処理してみます。

 

Rを立ち上げると以下のようなマークが出てくると思います。

細かいことはともかくとして、ここに以下のように打ち込んでみます。

fisher.test(rbind(c(30,4),c(24,12)))

 

そうすると結果が出てきます。

少し解説すると

c(30,4)やc(24,12)は行列の一列分で

rbind()で2X2の行列を作っています。

「このテーブルをFisher正確検定しなさい」というコマンドです。

 

難しくてむりという人はEasyRが良いかもしれません。

easystats.hatenablog.com

統計ソフトは何を使う?

医療用の統計ソフトはいろいろとあるがどれが良いのか?

思いつくだけでもいくつも種類がある。

  1. JMP
  2. R
  3. EasyR
  4. SPSS
  5. エクセル統計
  6. PRISM

もし、施設でJMPやSPSS、PRISMが利用できれば比較的使いやすいだろう。しかし自分で購入するとなるとかなり効果だ。

 

この中でフリーソフトはRとEasyRだ。

ダウンロードは以下を参照

easystats.hatenablog.com

easystats.hatenablog.com