負の二項分布（negative binomial distribution）は，（a）その母数を成功回数（ないしサイズ母数）rと，成功確率pとする場合（こっちに親しみ）と，（b）サイズ母数rと平均μとするときの二種類があるんだそうだ。パラメーター化の方法が違うってのは，ガン…

自分で関数書けばいいってのもごもっともだし，他にも色々あるんだけど，bayesmパッケージにnmatっていう関数がある。 vcovm<-matrix(c(3,2,2,3),2,2) cm<-nmat(vcovm) matrix(cm,2,2)

12/17の名古屋，外国語教育メディア学会中部支部外国語教育基礎研究部会第4回年次例会で討論型ワークショップやります！ このワークショップは，（私の得てしてくだらない話題提供のあと）基本的に参加者同士でグループを組んでもらい，想いの丈を共有し合う…

外国語教育研究では強く歪んだ分布をもつ連続型確率変数を扱うときがある。 学習時間，ある変数上でのテスト成績，心理学的な実験における反応時間，読解時間などがそう。場合によっては，これを正規分布にしたがう確率変数にすることを目的に対数変換すると…

美しいリレーの話。 まずはRのスクリプトファイル（.r）を用意する（a.rとする）。 それでこのスクリプトファイルをRscriptで実行するバッチファイル（b.bat）を用意する。 このバッチファイルを実行するVBScript（c.vbs）を用意する。 vbsは，コマンドプロ…

そういう論文ってどれだけの数あるかわからないのだけど，こんなメタ分析があるって考える。 （潜在変数としての）英語の力と（同じく潜在変数としての）国語（母語？）の力の関係をメタ分析でもとめる，という話にしよう。基本的に，英語の力（の推定値）は…

これ，どうやってご説明差し上げたらいいか迷ってたけど，説明用のグラフ用意した。 当たり前だけど相関係数と密接な関係がある（なのであまり操作としてよくない） 相関が十分に高いとき，それぞれ分割した変数の方は当然正規性が満たされなくなる（←これが…

ある質的な成果指標（ある試験の合格・不合格）がある。テストが終わったあとに，合格した人と不合格だった人に，どの教科書（A，B）を使っていたかそれぞれ聞く。すると以下のようなクロス集計が得られる。このとき，使用した教科書と合格・不合格の関係を…

何故か知らんかった。plot関数って思ったより賢い。（いや，plot関数が賢いのは昔から当たり前なんだけど）ひとつのグループ化変数とひとつの連続量を入れたデータフレームがあるとする。こんなふうにして作ろう。 group<-c(rep("a",25),rep("b",25)) score<…

オンライン教材へのログイン時間といった学習履歴は，現在の外国語教育現場では容易に得られるようになってきていて，（主に探索的な目論見で）分析に組み込まれることが多くなってきている。しかし外国語教育研究ではこういったデータの分析の歴史は浅く，…

あるテストを実施したとき，クラスごとに合計得点をもとめたいとか，そういうことって外国語教育ではよくある（そういうことしかない気もする）。 男女とか，3クラスとかだったらね，記述統計をエクセルでそれぞれ手動でもとめてもいいかもだけど，20クラス…

テトラコリック相関係数でメタ分析したいときに，先行研究がオッズ比で報告してたりするからどうするんだって話があるそうな。もちろんクロス集計があればいいのだけど，それがないときはどうするか。オッズ比からテトラコリック相関係数の近似値を求める方…

今更だけど，BESTっていうパッケージなんだって。Bayesian Estimation Supersedes the T-testということでBESTパッケージ。うむうむ。https://cran.r-project.org/web/packages/BEST/vignettes/BEST.pdf CRAN - Package BESTなるほどね。基本はまずJAGSでMCM…

…という矢の問い合わせ。基本図の中はサンセリフでいいと思うんだけど，わざわざこういう指定のジャーナルもあるんだそう。 なるほど。 windowsFonts(TNR = windowsFont("Times New Roman")) par(family="TNR") x<-data.frame("How about"=rnorm(100),"this"…

信号検出理論における被験者反応を4つに分けて（H, Miss, FA, CR）2×2表にしたときの対数オッズ比は，0.6倍すると信号検出理論における弁別力指標d'にまあまあ近似する。なので，大雑把にいえば互いに代用できなくもない。 #データを作る dat<-matrix(c(100,…

背景 日本の外国語教育研究では，先達のたゆまぬ努力によって1990年代から認知心理学的な研究手法が欧米より徐々に輸入されはじめ，2000年代後半から広く一般化し，今に続く研究の流れを形成した。認知心理学といっても，特に外国語の語彙処理に関わる研究が…

なぜか一部の人はMANOVAを嫌うのだけど，MANOVAしたほうがいい場合っていっぱいあると思う。 で，ある統計の相談案件で効果量はどう報告したらいいのっていう話になった。「そういやあまり聞かない」…というのは国内の外国語教育研究で概説が見当たらないっ…

私たちの友達グループ ある友達のグループで，これまでの人生で付き合ってきたひとの数を集計してみたら，大体真ん中（50%点）ぐらいが12人だったという。世では4人くらいを「普通」だといっているんだけど，私たちの友達グループでは12人を「普通」と呼ぼう…

ある指導法Aの効果をみたいとする。 NEGD（群間計画）を組めたとする。 処遇後における処置群と統制群の中心傾向の差をもって効果と考えたい。 このとき，結果変数として2変数もっているとする。 文法テストと語彙テスト。 このとき，効果量としてもちいるべ…

Rでマハラノビス距離をもとめるためには，mahalanobis関数を使う。 勉強用に書いてみた。マハラノビス距離自体は，やっぱり英語のwiki先生が詳しい。Mahalanobis distance - Wikipedia #2変量のマハラノビス距離をもとめる #数値例の作成 library(MASS) roun…

逆正規分布（逆ガウス分布，Inverse Gaussian Distribution）は，反応時間の分析で使ったりすることがあるらしい。歪んでいるので。2母数の連続型分布で，形状パラミタ（λ）と平均（μ）をもつ。 相変わらず英語のwiki先生は統計と数学に詳しい。Inverse Gaus…

これで外国語教育でやるような処理は大体できると思う。ご自由にどうぞ。 library(retimes) #ない場合はインストールすること library(MASS) library(ks) x<-rexgauss(1000,300,200,500) #数値例の生成 #基礎 x<-x[!is.na(x)] #欠損の除外（ここにはない） s…

前よくわからないから 積分でやったんだけど，調べたらwikipediaにあったから関数作ったった。wikipediaすげえな。でも700円のコーヒーを飲む層ではないんだ俺。 pexgauss<-function(x,mu,sigma,tau){ tau<-1/tau u<-tau*(x-mu) v<-tau*sigma phi1<-pnorm(u,…

Ex-Gaussian分布の累積分布確率をもとめる方程式が分からないので（わたしが馬鹿だというのと調べるのがメンドイ），retimesパッケージの確率密度関数dexgaussで出したのをintegrate関数で数値積分して出してしまえ的な。 でもretimesは使ってない。このまま…

背景 手元に判断課題のデータがあって，信号検出理論の指標を出したいとする。大概の信号検出理論指標は計算が簡単なのでExcelでもできるのだけど，大量データだと結構面倒くさい。 信号検出理論では，反応を4種類に分ける。 Hit, Miss, CR, FA。 これを数え…

背景 効果量のCohen's dは2変数の分布の重なりの程度だとも解釈できる。 こういうこと。 くわしくは浦野先生（2013）など。 LET2013 workshop | ひとりごと ただし，これは（想定上の）母集団についての話。 さらにCohen's dは，プールされた標準偏差を使う…

最近受けた統計の質問。 カイ二乗じゃなくてG検定をやれ？ これはこの前行われたCELESの，私のワークショップのときに，ある先生からご質問頂いたのだけれども，ちょっと事例が個別すぎたので，その場ではお答えできなかったもの。こちらから連絡差し上げら…

これはすごい便利。エンドユーザー向けでだれでも使えそう。 orddomというの。ordinal dominance statisticsの略かな。 共通言語効果量とか，Cliff's deltaとか一気に出してくれる。 対比較とかこれで結構いける的な。 Rogmann, J. J. (2013). Ordinal Domin…

背景 こんなデータを作ります。対応なしの2群のデータ（x, y）です。 xは，平均50，標準偏差10の正規分布に従う乱数39個です。 yは，平均55，標準偏差10の正規分布に従う乱数39個です。 なので，理論的には，こんな感じになり， 効果量（標準化平均差）は0.5…