1

私はRが初めてで、観測値の表を期待値の1つと比較して、chisqを計算しようとしています。課題の一環として、期待値テーブルを、観測値からランダムな順列を使用して作成した 999 個のテーブルのセットと比較する必要があります。各テーブルの chisq 値 (nsim=999) を計算し、観測データからの実際の chisq とともにすべての chisq 値のヒストグラムをプロットする必要があります。私が使用しているデータとコードは次のとおりです。

> survival=table(titanic[,c("CLASS","SURVIVED")])
> survival
      SURVIVED
CLASS   no yes
  1st  122 203
  2nd  167 118
  3rd  528 178
  crew 673 212

> expected=expected(survival) #library(epitools)
> expected
      SURVIVED
CLASS        no       yes
  1st  220.0136 104.98637
  2nd  192.9350  92.06497
  3rd  477.9373 228.06270
  crew 599.1140 285.88596

>nsim=999
>random= rep(survival,nsim)

そして今、私は立ち往生しています！

4

1 に答える 1

1

順列を生成する最も簡単な方法は、列で次のsampleコマンドを使用することです。"SURVIVED"

sample(titanic[,"SURVIVED"])

その列のラベルをシャッフルし、yes/noこれを 999 回繰り返すことができます。

replicate(999, {
  permSurvival <- sample(titanic[,"SURVIVED"])
  # Code to measure chi square test goes here
})
于 2013-09-25T01:57:57.607 に答える