1 
  p1 <- c(.25,.025,.025,.1,.2,.4)

  T <- sample(1:6,size=N,replace=TRUE, prob=someprobabilityvector)
  Y <- rbinom(N,1,p1[c(T)])

こんにちは皆さん、私はRとプログラミング全般に不慣れで、基本を理解するのに助けが必要です。上記のベクトルYで何が起こっているのか誰かが私に説明してもらえますか？私は上記のことを理解してp1[c(T)]います。しかし、ベクトルYが何をしているのかわかりません。すべての助けは事前に感謝されます。

4

1 に答える 1

2

コードの最初の行は、6 つの確率のベクトルを作成します。

p1 <- c(.25,.025,.025,.1,.2,.4)

N2 行目では、1 から 6 の数字から値をランダムに選択します (置換あり)。各値の確率は に指定されていsomeprobabilityvectorます。Nしたがって、関数は1 ～ 6 の値を含む長さのベクトルを返します。

T <- sample(1:6,size=N,replace=TRUE, prob=someprobabilityvector)

3 行目ではN、1 回の試行と確率で指定された 2 項分布から乱数p1[c(T)]が生成されます。c(T)と同じTです: 1 から 6 までの値を含むベクトル。ベクトルは、ベクトルのインデックス付けに使用されp1ます。したがって、 vectorからの値をp1[c(T)]含むベクトルを返します。Np1

Y <- rbinom(N,1,p1[c(T)])

指定された二項分布には試行が 1 つしかないため、ベクトルYにはゼロと 1 が含まれます。

于 2013-02-03T09:16:15.170 に答える