-1

イベントの事後確率を計算するのに役立つ R のパッケージを探しています。ありますか？

わかりました、私はそのようなデータセットに取り組んでいます

学年 学歴 合格
群1 初等 50 なし
群2 高等 20 不可
群1 中等 70 可
群2 中等 67 可
群1 中等 55 可
群1 中等 49 不可
群1 中等 76 可

生徒が試験に合格する事前確率は 0.6 です。次に、生徒の年齢、教育レベル、学年
を考慮して、生徒が合格する事後確率を取得する必要があります。最初に取得する必要があることがわかっています P(age=group1| pass=yes )* P(education=primary| pass=yes)* P(grade>50 |pass=yes)
しかし、これはケース (行) ごとに行う必要があり、1000 行で設定された日付がある
ので、取得できると思いました関数がこれに役立ちます！

4

2 に答える 2

-1

これは、1 つの変数 (教育) と変数 (合格) に対する答えです。 

# get the prior probabilities  
prior<- c(prior_no, prior_yes)

# get contingency table of values for mydata
edu_table<- with(mydata,table(mydata$pass, mydata$education))

# get the sum across (pass)
tots<- apply(edu_table,1,sum)

# create matrix of 0's
ppn<- edu_table*0  
post<- edu_table*0  

# use a loop to get the prior probabilities& posterior probabilities
for(i in 1:length(tots)){
  for( j in 1: 4){

    ppn[i,j]=edu_table[i,j]/tots[i]
    post[i,j]=prior[i]*ppn[i,j]/(prior[1]*ppn[1,j]+prior[2]*ppn[2,j])

  }
} 

ppn # probability of education=j given y=i
post # posterior probability of y=i given education=j
于 2013-06-25T20:43:39.757 に答える