R– category –

R で二項分布のグラフを書いてみた 二項分布のグラフの前に 二項とは？ 二項とは、0か1か、YesかNoか、表か裏か、など二つに一つが選ばれること。 二項分布は、二項の場合の分布という意味である。 二項分布のグラフを書いてみる たとえば、10例で、確率5...

t分布とは何か？t検定とどう関係するのか？ t分布とはいったい何なのか？ t検定とはよく聞くが、それとt分布は関係あるのか？ t分布は要するに何なのか？ 実は、t分布は標準正規分布の代わりだ。 どんなときに代わりとして使うか。 2つ条件がある。 母分散...

平均70点、標準偏差15点のテストの場合、90点以上の学生は上位何パーセントに当たるか？ という問題に使うパーセンタイル percentile とクォンタイル quantile。 それぞれ百分位数と分位数ともいわれる。 パーセンタイルを求めるには？ pnorm()を使う。 母...

正規分布を R で書く方法の解説 正規分布を R で書く方法 curve()を使う。 dnormが、正規分布にしたがう確率分布を計算する。 横軸の範囲は-4から4としている。 3つの分布とも平均0とする。 R のスクリプトは以下の通り。 curve(dnorm(x,mean=0,sd=1.0), f...

R を使ったポアソン分布の書き方の解説 ポアソン分布とは？ ポアソン分布は、数学者ポアソンが発表した、離散確率分布。 離散とは、離れているということ。 0回、1回、2回、3回・・・が生じる確率だ。 ラムダλ が1の場合の 0 から 10 までの確率を計算して...

重回帰分析で、変数選択において、考慮しないといけないのが、多重共線性。 多重共線性のチェックは統計ソフトRでどうやるのか？ 多重共線性とは何か？ 生まれた西暦年と年齢は、どのくらい相関するだろうか？ 全く一致するといっても言い。 生まれた月に...

負の二項回帰モデルを R で行う方法 負の二項モデルを R で行う方法 まず、データを読み込む X <- c(rep(1, 913), rep(2, 1345), rep(3, 519), rep(9, 863)) Y <- c(rep(0, 910), rep(1, 2), rep(2, 1), rep(4, 0), rep(0, 1333), rep(1, 11), rep(2...

傾向検定は、サンプルデータで観察された、平均値や割合が、だんだんに大きくなる、だんだんに小さくなるという傾向が、母集団でもその通りか検定するもの。 Rでどのようにやるかまとめてみた。 単変量の傾向検定：コクラン・アーミテージ検定 回答の選択...

コレスポンデンス分析（対応分析）とは、人の属性の類似性(もしくは対応)を図に表す分析。 結果を見るといろいろと解釈したくなる、面白い分析。 R では、caというパッケージで分析可能。 コレスポンデンス分析のパッケージ ca をインストールして準備する...

ギブスサンプリング（ギブス法）を R で実施する方法。 初期値を与えて、何度もサンプリングしていると、だんだんに形になってくる。 ギブスサンプリング用のデータの用意 heikinという名前のデータを入力する。 heikin <- c(6.0,10.0,7.6,3.5,1.4,2.5,...