stanKOC2017の分析はじめにキングオブコント,2017が先日開催されましたね。いや,今年も面白かった。にゃんこスター,すごかったなあ,一本目。あんなんずるいやん,と言いながらずっとゲラゲラ笑ってました。みたいなのが生まれてくる土壌があるところだと思っています。審査員長の松本人志がいう,「ある程度いやらしさが必要」というのはそうところだと思うの。そんなの卑怯じゃない,そんなのマナー違反じゃない,と思いながら笑わずに得られないというかね。あかんと思いながら笑ってしまっている自分がいるというかね。2本同じ系列だったのは残念なところ。かまいたち・山内に言われなくても「一回冷静に」なっちゃったもんね。違う切り口で笑わせてもらいたいと思うわけですよ。あのやり方,二回連続では笑わない。いう意味では,かつての2700はよかった。右肘左肘の次にやった「キリンスマッシュ・キリンレシーブ」のネタってみなさん覚えてます?あれすごいことですよ。色物すぎて賞は取れなかったけどねえ。採点」で調べると,やれ松本人志の基準がぶれているだの,アンガールズに甘かっただのといった意見があるわけですが,個人的には我々が共感できないレベルでの評価軸を持っててこその審査員だと思うので,私はそういう反論に組みするものではないです。見て楽しい人と,舞台に立って実際にやる人,やってる経験がある人から考える上手い下手ってあると思うもんね。いうことで,審査員の好みとか,にゃんこスターの特殊性みたいなのがデータから観れるかどうか,分析してみます。まとめだらだら読むと大変なので,先に結論を。相関係数をベイズ推定するときはサンプルサイズをたっぷりねウィッシャートブンプを使ったらその辺すこしはましだわねMDSとMCMCの組み合わせで幻想的な花火が描けますよここから本題。94),ncol=5,byrow=T)まずは単純な審査者の評定の類似性。これは相関係数でいいと思うのです。いうことで,しれっとstanで推定してみます。parametersブロックのところでquad_form_diagという謎関数を使っていますが,これはStanマニュアルを読んでいただければ「あー,行列とベクトルを組み合わせた計算するのね」とお分りいただけると思います。あとは素直にコード化しているので,読んでいただければわかるかなと。相関係数はちょっと歪みがあるかもしれないので,get_posterior_meanで取り出すのではなくmedian(MAP)推定値が欲しい。いうことで,MCMCサンプルから中央値を取り出して行列にするコードをR側で準備。1.0000000これでできましたよっと。1.0000000む,ちょっと違うな。いうことでちょっとシミュレーションとかして確認してたんですが,サンプルサイズが小さい(今回15回分しかない)のが問題だったかも。irisデータとか使うと,まあまあな数字になるんですけど。それにしてもずれてるなあ。ツイッターで呟いたら,どこからともなく(@simizu706$$)「偏差積和行列渡して,ウィッシャート分布使いなはれ。あと標準偏差の推定やめなはれ」というアドバイスが。1.0000000ふぬぅ,まあいいでしょう。これをみると,バナナマンの二人はやっぱり評価が似てるんですよね。松本さんと大竹さんも評価が似ている。似てる,似てないからマップを作るのは,そうだね,MDSだね。cmdscale((1-tmp),2)plot(result.mds,type="n")text(result.mds,colnames(tmp))こうしてみると,左右はっきり分かれます。左側がバナナマン的な評価次元。設楽さんはわりとロジックの通ったシュールな世界とかを作るので,やや論理的な方向,右側は勢いを評価する次元というところかな。上下で考えると,上はシュールで下はインパクトみたいな感じもしますが。ちなみにこれはMAP推定値をプロットしているんであって,実際はたくさんのMCMCサンプルがある=座標も分布しているわけです。この分布をプロットしてみましょう。雲のように重ねてみました。してみると,三村さんのポジションがちょっと広く散らばっているようですね。他の人と相対的にみると,ちょっと立ち位置の確信度が低い,という言い方ができるかもしれません。相関行列に対してMDSをする,というのはない話ではないのですが,MCMCサンプルについてそういうことをやるっていうのは珍しいかな,と思ってブログにしてみました。いう遊び方もありだと思うんです。おまけ最後におまけです。行と列を入れ替えて,すなわち演者の方についてMDSをやってみましょう。公平を期すために?1回目のネタだけを対象にやってみたいと思います。自由度の関係でMCMCもできないから,記述統計量で。geom_point()+geom_text(family="HiraKakuProN-W3")gこうしてみると,やっぱりにゃんこスターは独特のところに位置づくんだな。さんから,相関行列のMCMCサンプルから行列で推定値を取り出す方法,教えてもらいました。今年も塾合宿の季節がやってまいりました。史上最大規模の,33人が参加。キャンセル投稿を送信できませんでした。メールアドレスを確認してください。メール送信チェックに失敗しました。このブログではメールでの投稿共有はできません。