RMPとIRPの差分ベクトルの積和平均による項目選択ルール
期間内に終わりそうもないのでかわいそうだからまとめる.
これに懲りたら今後は計画的に作業を進めること.
追記:書き終えたあたりでなんか式に矛盾を感じてきた.
とりあえず雰囲気だけ掴んで細かい式は先輩方の論文の方を参考にしてください.
やりたいこと
コンピュータ適応型テスト (略称CAT)
項目応答理論 (略称IRT)によって推定されたパラメータの付与された問題項目(アイテム)をアイテムプールに置き、受験者の回答パターンに応じて、困難度の異なるアイテムを選んで出題することにより、動的、短時間、高精度で受験者特性の測定を判定を行う方法.
これをIRTではなくLRTでやりたいですね.
LRTだと離散的なので例えば
A君はクラス1,B君はクラス3
みたいな感じで受験者をレベルに応じてクラス分けできるぞ.
やりかた概要
- 初期暫定RMPを算出
- 暫定 RMPとIRPの差分ベクトルの積和平均が最小となるものを選択し出題する(ここを何度も繰り返す)
- 条件を満たしたらテスト終了
以下に詳しいやり方を書く.
1番
まずは初期暫定RMPを出す.
RMPを出すために受験者に何問か問題を解かせる.
今回はランク(困難度)が3つなので3つ解かせる.(各ランクからランダムに1つずつ)
3つ解かせ終わったら実際にRMPを算出する.
ランク・メンバーシップ・プロファイル(rank membership profile, RMP)は,各被験者のそれぞれの潜在ランクに対する所属確率を見るのに便利
とりあえず先に説明しておくと,RMPはベクトルだよ.
次元数はランクの数と同じだよ
今回は3つだから最終的にRMPをとすると
になるね
つまり,はその人がランク1に属する確率だよ
それでは,初期暫定RMPを求めていくよ
RMPを求めるために必要なのは以下の2つ
- 受験者の回答
今回は最初に3問出題することにしたのでとする.
はランク1の問題に対する回答
はランク2の問題に対する回答
なお,回答は当たったかどうか(1 or 0)の2値とする
1 : 正答
0 : 誤答
- 出題した問題のIRP
IRPもベクトルだね.簡単に言うとある問題に対してのそれぞれのランク毎の正答率だね.
ランク i の問題のIRPをとすると今回はランクが3つなので
ってかんじになるね.
はに対するランク1の人の正答率.
また,も今回は最初に3問出題することにしたので
の3つができるね.
↓下図はランク数10のIRP
RMPの要素を求める式は
だよ
も分子を変えるだけで求められるね.
一応だけど分子の p は確率(probability)の p だぞ.
で,これが何を意味しているかというと
正直資料が全然みつからなくてわからないんだぜ.
だから,俺の予想を書くぜ
まず, っていうのは事前確率を表しているぞ.
= 事前にわかっている回答者がクラス1に分類される確率ってことだね.
さて,今事前にわかっている事前確率はいくつかな?
↓
そうだね,まだわからないね.
というわけでまだ分からない
↓
どのクラスに分類される確率も等しい
とうことでクラス数が3つなので
にしましょう.
同じようにもだね.
それじゃあ
の分子 について考えてみよう.
っていうのは
のとき が起こる確率だね.
これは文字に起こすと
IRPが のとき回答が の確率といえる.
つまり,
とあるクラス1の問題は各ランクの人の正答率が なんですけども,その時に正答(or 誤答)する確率は?
って意味だと思う.
で,それに事前確率 が掛けられているので,
確率で掛け算は○○かつ××って意味なので
クラス1に属する確率 の人が,とあるクラス1の問題(各ランクの人の正答率が )に,その時に正答(or 誤答)する確率
ってことです.
(余計わかりづらくなった気がする…)
話が逸れたから実際の計算に戻るよ
RMPベクトル :
←これを求めるんだったね.
追記:多分ここが間違ってる↑なんか多分総和とかとらないといけないハズ
で求めるよ
実は,ここがマジで資料見つからなかったから曖昧だよ.
なんかややこしいけど,これは結局
正答だったら
誤答だったら
になるから,指数の左側か右側の片方は消えるよ
で,結局
正答だったら → IRPから得られる正答確率
誤答だったら → IRPから得られる誤答確率
っぽいのが出てくるっぽい.
そんでもって, の分母
これはたぶん の要素の総和が1になるように正規化するためのものだと思うよ.
総和が1だったらそのまま確率(百分率)で表せるからね.便利.
というわけで,これで初期暫定RMPは求められるはず.
RMPの求め方(初期暫定RMP以外のRMP)
ここら辺あいまい
初期暫定以外の普通のRMPはどうやって求めるかというと,
初期暫定RMPとほとんど同じですよ.
ただし,一つだけ違うのが事前確率
こいつは最初はわからなかったから 1/3 にしたけど,
2回目からはこの値にRMPを使えばいいんじゃない?
だって,RMPって被験者の各クラスに対する所属確率だったもんね.
そんでもって,ここが一番わからないんだけど
RMPの更新についてのお話
例えば,「初期暫定RMPが求まり,その値からクラス2の問題が選ばれた」としよう
↑つまり4問目
その時その問題の回答結果からRMPを更新することとなるが,
俺の中でどっちの更新式が2パターンある.
どっちが正しいかは知らんし,どっちも正しくないのかもしれない.
①どんどん値が増えてくパターン
追記:絶対間違ってる↑
②出題した問題のクラスに対してのみ更新するパターン
追記:絶対間違ってる↑
ニューラルネットワークに対する理解
ニューラルネットワークが何をしているのか? 今の自分の考えを記す.
とりあえずニューラルネットワークの図
図の丸はそれぞれパーセプトロン. 左が入力層,中央が隠れ層,右が出力層となっている.
何がしたい?
→ 何かの式を解きたい
とある真のモデルからサンプルを取ってくる ↓ それをNNに食わせる ↓ モデルを予測する
例えば
こんな感じのサンプルがあります.
なんとなくではありますが,
こんな感じに線を引けそうですね. この線を予測するのがニューラルネットワークの働きです. これにより真のモデルを予測します. ちなみに,今回はサンプルの次元が身長とバストの2次元だけだったので線を予測しましたが,3次元のデータなら面を,4次元以上なら超平面のようになります.
モデルが予測出来たら何がうれしいか?
→ サンプルにないデータも予測できる
例えば,今回のデータを使ってサンプルにないアイドルの身長さえ分かればバストが予測できる という寸法です.
じゃあ,予測するにはどうすればいいの?
→ 誤差逆伝播法
これについては後日.
パーセプトロンってなに?
ニューラルネットワークを作ろうと思ったので、とりあえず基本となるパーセプトロンから。
図にするとこんな感じ。
xが入力、zが出力となっていて
- u = w * x
- z = f(u)
関数fは活性化関数と呼ばれるもの。
単調増加する非線形関数が用いられる。
よく使われるのがロジスティック関数とかシグモイド関数とか呼ばれるもの
- f(u) = 1 / (1 + exp(-u))
他には双曲線正接関数やランプ関数などが用いられる
結局のところ、「ある値付近で急に立ち上がるインパルスみたいな関数であればなんでもいい」と理解している。(たぶん)
パーセプトロンの話に戻る。
入力は何個でも構わない
入力が増えたらその分だけ計算が増える
- u = w1 * x1 + w2 * x2 + w3 * x3
- z = f(u)
これをたくさん並べたものがニューラルネットワークとなる。
続きはまた今度。