統計ER

統計ソフトRの使い方を中心に、統計解析方法の解説をするブログ。ありそうでなかなか見つからないサンプルサイズ計算などニッチな方法について紹介しています。

EZRで比例ハザード性を確認するには?

生存時間解析でCox回帰を行う前提として比例ハザード性が確認できていないといけない。 比例ハザード性はどのように確認すればいいのか?

EZRで連続量を正規分布にする方法 Box-Cox 変換

連続量を何らかの方法で正規分布に近づける方法はいくつかある。 よく使う方法は自然対数をとる方法だ。 ここではBox-Cox変換の方法をまとめてみた。

EZRで描いたカプランマイヤー曲線の横軸を変更したい

EZRでカプランマイヤー曲線は描けたのだが、横軸(時間軸、X軸)を調整したい。 どうやれば、変更できるか?

EZRで描いたグラフを修正する方法 ― Macでは?Windowsでは?

EZRでグラフを描いたあと、一部文字を消したり、文字を加えたり、図形を加えたり、ちょっとした修正を行いたいことがある。 そういうときに、Mac版であれば、どうしたらよいか? Windows版であれば、どうすればよいか?

EZRでグループ別の散布図を描く方法

EZRでグループ別の散布図はどのようにして描くのか? 色はどのようにして変更するのか? 点の形は変更できるのか?

偏相関係数のサンプルサイズ計算

お互いにある変数と相関がある変数同士の相関係数を算出する場合、お互いに相関がある変数の影響を排除したうえでの相関を検討したいことがある。 そのようなときに使う方法が偏相関係数だ。 偏相関係数を算出する試験デザインを考える際に、サンプルサイズ…

EZRでフィッシャーの正確検定が有意なのに95%信頼区間が1をまたいでいる場合どうしたらよいか?

オッズ比の95%信頼区間が1をまたいでいないことと、有意水準5%で統計学的有意であることがイコールであることは知っていると思う。 しかし、EZRでフィッシャーの正確検定を実施すると、有意確率が5%未満であるにもかかわらず、オッズ比の95%信頼区間が…

EZRのFisherの正確検定でオッズ比を計算するとオッズ比が違う件

EZRのFisherの正確検定の関数でオッズ比を計算すると、通常のオッズ比と異なる数値が計算される。 これはなぜなのか? どう考えたらよいか? EZRのFisherの正確検定を使わないオッズ比の計算 EZRのFisherの正確検定を使った場合のオッズ比の計算 EZRのFisher…

EZRのフィッシャー正確検定でオッズ比の信頼区間のパーセントを変えるには?

統計ソフトEZRは、無料で使えて、しかもメニューから選んで分析できる、簡単便利な統計ソフト。 フィッシャーの正確検定でオッズ比の信頼区間のパーセントを変えることはできるか?

EZRでグラフの凡例の位置を調整するには?

統計ソフトEZRは、無料で使えて、しかもメニューから選択するタイプの使いやすい統計ソフトだ。 グラフを描く機能も優れている。 たまに、グラフの凡例がグラフと重なることがある。 それを調整するにはどうしたらよいか?

EZRの起動時のフォントやフォントサイズを変更するには?

EZR

統計ソフトEZRは無料で使える、医療統計に最適なソフトウェア。 起動時のフォントやフォントサイズを変更するにはどうやったらいいか。

尤度(ゆうど)とは?

読み方も意味しているところも分かりにくい尤度(ゆうど)。 理解したところを共有したい。

EZRのグラフの色を変えるには?

統計ソフトEZRは無料で使えるのに、メニューから統計解析を選んで操作できる、初心者に優しいソフトウェアだ。 ただし、グラフの色を変えるには少しコツがいる。 少しだけRスクリプトを知る必要がある。 この記事を読めば、どこをどのように変えれば、色を変…

KH Coderでテキストマイニングするためにどのようなデータをどのように準備したらよいか?

無料テキストマイニングソフト KH Coder でテキスト分析をするにあたり、どのようなデータをどのように準備すればいいのか? そんな悩みを抱えているあなたに。

相関係数に必要なデータ数は?

相関係数計算に必要なデータ数はいくつか? 相関係数の検定との関係から紹介。

時系列データの自己相関を計算してみる。時系列データへの変換。統計ソフトRでやってみる

時系列データの自己相関を計算する方法を調べたので紹介したい。 時系列データ解析法ARIMAの基礎的な部分になる。

反復測定分散分析 群分けあり Split Plot Design のサンプルサイズ計算

反復測定分散分析 群分けあり Split Plot Designと呼ばれるデザインの分析におけるサンプルサイズ計算は、いろいろ探してもなかなか出てこない。 ようやく方法を見つけて、やっとのことで読み解いたので、共有したい。 統計ソフトRの関数の形になっているの…

分散分析のサンプルサイズ計算―pwr.anova.test()を使う方法

分散分析のサンプルサイズ計算を統計ソフトRで行う方法を解説。分散分析は三群以上の平均値を比較する分析方法。少なくともどれか一つの群がほかの群とは異なることを証明するための分析方法だ。

反復測定分散分析でグループ比較をしたいときに生じる疑問

複数回測定しているデータがあり、複数グループに分けられるとする。そうすると、反復測定分散分析を用いてグループ比較を行いたくなる。 反復測定分散分析でグループ比較を行いたいと思ったときに生じるいくつかの疑問に答える記事。 反復測定分散分析でグ…

ウェルチの検定 Welch's t-test

平均値の差の検定、いわゆるt検定は、等分散かどうかなんて気にしないでいい。 どんな時もウェルチの検定を行えばいい。

おすすめ統計本―医学統計を中心に

統計ER的おすすめ統計本のリストアップ。 このブログでたびたび参考書籍として登場している本を厳選してチョイス。

統計ソフトRでサポートベクターマシンを使ってみよう!

統計ソフトRでサポートベクター分類器・サポートベクターマシンを使って機械学習させるにはどうすればよいか?

割合の差の検定のサンプルサイズ計算

二群の割合を比較したいときは割合の差の検定を使う。割合の差の検定を使う場合のサンプルサイズ計算はどのようにすればいいか? 統計ソフトRでやってみるとどうなるか? エクセルではできないか?

Clopper-Pearson法で割合の信頼区間を計算するには?

割合の信頼区間の推定法としてClopper-Pearson法というものを目にしたので、簡易的な割合の信頼区間推定法とどう違うのか統計ソフトRで計算してみた。

統計ソフトRのパッケージをsourceからインストールする方法とBioconducorのパッケージの場合

統計ソフトRはパッケージを追加すると新しい機能が追加できる。その方法もとても簡単だ。 統計ソフトRのパッケージの追加方法として例外的な方法があることを知ったので、共有したい。

時系列データ「ドル円レート」を分析してみた!

時系列データとは何か?時系列データを分析する方法は? 日時と値が得られていて、その後の予測をしたい場合、時系列データ分析を検討されたし。

成人喫煙率はいつゼロパーセントになるか?

成人の喫煙率は年々減少している。 このままだと近い将来ゼロパーセントになりそう。 実際のところ、いつごろゼロパーセントになるか。 時系列データをもちいた未来予測を行ってみたい。

Bonferroni 型 多重比較のサンプルサイズ計算 - 統計ソフトRのスクリプト

三群以上の平均値や割合の差を比較したいデザインの場合、サンプルサイズはどのようにすればよいか? Bonferroni (ボンフェローニ)型の比較は、保守的過ぎて有意水準が厳しすぎるとは言われるが、逆にBonferroniの条件を満たすサンプルサイズなら、HolmやH…

生存期間中央値の求め方 - 5年生存率から逆算で求める方法

がん患者さんの生存期間中央値を5年生存率から逆算で求めるにはどうしたらよいか?統計ソフトRで実施してみた。

ロジスティック回帰&トレンド検定&多重比較 - 統計ソフトRでタイタニック号データの分析

統計ソフトRに組み込みのタイタニック号の生存・死亡データで、生存・死亡に関するロジスティック回帰分析と、独立変数の多重比較を実行してみる。 今から100年以上前の未曽有の大惨事の鎮魂を祈りつつ、解析してみる。