相関係数は標準化された共分散。 共分散は偏差積の平均値。 偏差は平均との距離。なので、片方の変数が平均から離れている度合いと、もう片方の変数の平均から離れている度合いをかけたものの平均から、スケールの大きさを取り除いたもの。数学的な含意は、…

識別力が低い項目ばかりのとき、被験者の（古典的テスト理論的な意味での）正答率と反応時間の中心の相関（普通は負の値）は必然的に小さくなる。スピードに関わる課題操作が正答率に及ぼす影響も小さくなる。刺激の条件間で2つの相関に違いがあるとき、まず…

正答率は0から1までの範囲しか取らない。 正答率は平均が0か1に近づくと分散は小さくなる。反応時間は非負で最小値がある。 反応時間は平均が高くなると分散も大きくなる。当たり前だけど大事なことなのでメモ。

一般に被験者間でみたときに、定数であるところの能力が高い人は、低い人に比べて早く正確に問題を解く。この関係（被験者の回答時間と正答率の共分散）は、問題が難しくなるにつれ大きくなる。一般に項目間でみたときに、定数であるところの難しい問題は簡…

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

時代とその時々のインフラによって物事のスタイルはよく変わる。物事自体はあまり変わらないのだけど。私が学部生のときには既に、プレゼンといえばスライド投影だった。Microsoftのパワポを使うのがかなり一般化していたので、当時からスライド投影のことを…

直接ではないが師匠筋にあたる（むしろ私の指導教官とみなされることがかなり多い）O先生が、私の後輩たちに増井先生の富豪的プログラミングについて語ったという。私はプログラミングに関して素人同然なので必然的に富豪的なんだけど、富豪的ってのは時と場…

RはCUIだから絶対に嫌という層は一定数いる。ならUIつけたらいいじゃない。 HSPは書くのがハイパー簡単だし簡単なUIを作るのに適している。ただ高度な計算がめんどい。なのでHSPでUIを作って、そのUIからRを走らせたい。基本的な方針としてはHSPで被せた実行…

ご尊名をお書きしてもよいものかと悩んでいたけども，結局伏せさせていただくことにする。 K先生とおっしゃる先生がいらして，近くの大学にお勤めなのだけど，指導を頂いた師匠筋というわけではなく，あくまでも私淑するに留まる関係性。でも幸いにも，先生…

ある種の動機づけの強さについて，性別に由来する差を知りたいと思った。それを測る質問紙（10項目）があるとする。これら10の観測値の平均をもって，この動機づけの強さの値だとみなそう。それで，男女100人ずつくらい集めて，その値の2群の平均値を比べた…

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

前回のおさらい さて，続き続き。 「観測変数の相関が高いってだけで，同じものを測っているとはいえない」という話でした。 構成概念妥当性観のことを言い出すと私は延々と止まらくなるので自重するけど，これは，観測変数の相関係数をもとめることが実質科…

外国語教育の研究者の中には，「相関係数の値が高ければそれらふたつの変数は同じものを測っている（またはその度合いが高い）」っていう信念を強く持っているひとが多い。 そんなにひどい間違いではないと思うの。 でも，それは微妙にヒューリスティックが…

バッティングセンターにいくとして，100kmから160kmまですこしずつ球速を上げていくとする。100kmだと60%は打てるんだんけど，110kmになればけっこうきついわ，20%は切るな。160kmなんてほぼ無理だ。 ――って考えるとする。 理想的にはこんな感じ。 でもこの…

検定力のシミュレーションに使えるかも。 はっきり言って意味はない。 #2標本のt検定（対応なし，Welch） mcpower2<-function(b,n1,n2,m1,m2,sd1,sd2){ p<-numeric(b) for(i in 1:b){ s1<-rnorm(n1,m1,sd1) s2<-rnorm(n2,m2,sd2) p[i]<-t.test(s1,s2)[[3]] }…

自分でもたいそう長いことそういうふうに書いてたのだけど、例えば「ある実験的操作が反応時間に影響を及ぼすか」みたいな研究課題は、論文として許されるのだろうけど途方もないくらい馬鹿らしい。そんなことはある筈がない。いや、あるといえばあるかもだ…

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

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

そういえば去年のFrontiers in Psychologyに反応時間データの処理についての論文があった。 いろんな変換とか外れ値除去したりする方法をシミュレーションで試したけど，Box-Cox変換をλ = -1でやるのがいいぞ！という変な推しだった。まあ，いいけど。反応時…

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

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

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

背景 卒業論文や修士論文で，指導教官や先輩，または投稿論文で査読者から「標準偏差」を報告しなさいと言われたことがある方も多いと思います。 ただ，「標準偏差とはなにか」を理解することは簡単じゃありません（と考えるひともいるようです）。 ここでは…

心情としてはLinux派（のエンドユーザー）なのだが，これまでどうしてもMS Officeでやらなきゃいけないことが多いので，Windows 7も使っていた。 ThinkPad Tシリーズなのだけど，RAMが4GB余っていたので，RAMDISKとしてTMPファイルなどいろいろ割り当ててい…

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

外国語教育の研究者で，Rとりあえずインストールしたけど，何ができるかわからない，という方。とりあえず以下の様なパッケージを入れておいてはどうでしょうか。 psych:汎用心理学系分析パッケージ，主に信頼性係数や因子分析に使う lavaan:潜在変数分析用…

自分でもきちんとできていないのだけど，もちろん，OSはLinux系だと思う。それで基本的には，（GNU）emacsの中にひきこもっちゃう。やはりemacs。 emacsでメールを送受信して，emacsでスケジュール管理して，emacsでpdf開いちゃう。論文を書くときは野鳥（ya…

まだモックアップ版です。RでS-P表を返します。 ヘッダーつきのデータフレームを入れます。 S-P表を返します。 オプションのplot=TでS CurveとP Curveも描画します。 このデフォルトはTです。 sp<-function(dat,plot=T){ dat[order(apply(dat,1,sum),decreas…

lubuntuかわいいよ好きよ 背景 文系院生に必要な研究環境は，はっきりいって以下の要件に尽きる。 pdfが閲覧できること ネット・ブラウジングができること メールの送受信ができること オフィススイートを備えていること （人によって）ターミナル上で文字処…

第二言語習得研究（一応はそういう名前の授業の単位は大学院で取ってあるレベル）が一般的な教育従事者にとってさほど魅力的でないということはさまざまな先生がいっておられる。そして2015年の今日では，それはさほど真新しい見方でもない。でも，「なぜ」…