草薙の研究ログスマートフォン用の表示で見る草薙の研究ログ英語の先生をやってます。統計これで外国語教育でやるような処理は大体できると思う。library(retimes)#ない場合はインストールすることlibrary(MASS)library(ks)x<-rexgauss(1000,300,200,500)#数値例の生成#基礎x<-x[!is.na(x)]#欠損の除外(ここにはない)summary(x)#5数要約+平均skew(x)#歪度quantile(x,seq(0,1,.1))#10%点ずつquantile(x,seq(0,1,.05))#5%点ずつdata.frame("Counts"=hist(x)$counts,"Breaks"=hist(x)$breaks[-1])#度数分布表#外れ値y1<-x[x<mean(x)+2*sd(x)]#2SD以上を除去length(y1)/length(x)#除去率y2<-x[x<mean(x)+3*sd(x)]#3SD以上を除去y3<-ifelse(x>mean(x)+2*sd(x),mean(x),x)#2SD以上を平均値に置き換えy4<-ifelse(x>mean(x)+3*sd(x),mean(x),x)#3SD以上を平均値に置き換えy5<-ifelse(x>mean(x)+2*sd(x),mean(x)*2+sd(x),x)#2SD以上を2SDに置き換えy6<-ifelse(x>mean(x)+3*sd(x),mean(x)*3+sd(x),x)#3SD以上を3SDに置き換え#変換y6<-log(x,10)#対数変換(底を10)y7<-log(x,2)#対数変換(底を2)y8<-1/x#逆数変換#最尤法による分布へのフィットm1<-timefit(x)#Ex-Gaussian分布へのフィットm1@par#パラミターm1@logLik#対数尤度AIC(m1)#AICBIC(m1)#BICm1<-timefit(x,iter=100)#100個のブートストラップ標本m1@bootPar#ブートストラップ標本ごとの推定値apply(m1@bootPar,2,quantile,c(.025,.975))#パーセンタイル法による信頼区間apply(m1@bootPar,2,sd)#ブートストラップ標準誤差m2<-fitdistr(x,densfun="normal")#正規分布へのフィットm3<-fitdistr(x,densfun="lognormal")#対数正規分布へのフィットm4<-fitdistr(x,densfun="gamma")#ガンマ分布へのフィットm5<-fitdistr(x,densfun="weibull")#ワイブル分布へのフィットm2$estimate#推定値m2$sd#標準誤差m2$loglik#対数尤度AIC(m2)#AICBIC(m2)#BIC#モーメント法によるEx-Gaussian分布へのフィットmexgauss(x)#可視化boxplot(x,range=F,horizontal=F)#箱ひげ図k1<-kde(x)#カーネルplot(k1)#図示hist(x)#ヒストグラム#推定した母数をもつEx-Gaussian分布の確率密度曲線q<-seq(0,max(x)*1.2,1)plot(q,dexgauss(q,m1@par[[1]],m1@par[[2]],m1@par[[3]]),type="l",ylab="p",xlab="RT")#2017/4/6追記ある読者の方から,外れ値の除去を行列でやるにはどうしたらよいか?という質問をいただきました。ありがとうございます。やり方はありますが,以下のようにすれば最も単純ではないでしょうか。#このような実験では最も一般的なデータ形式,#つまり,行に項目,列に被験者の情報が入っている行列#またはデータフレームdatがあるとして,hazure<-function(x){x[x<mean(x)+2*sd(x)]}apply(dat,2,hazure)#これでOK#これで除去したもののみを返しているが,#これではデータの整理がつかないので,#この外れ値除去後に外れ値の平均を取る。#つまり,外れ値を除去した被験者平均。23:43Rで反応時間データの基礎的処理まとめTweet広告を非表示にする関連記事2017-02-23ゼロ過剰ポアソン分布をデータにフィットさせるゼロ過剰ポアソン分布では2つのプロセスが考えられている。...2017-01-11最尤推定した母数のもとでの確率密度曲線をヒストグラムに描き足すえっと,リクエストがあったのでここに書く。size)...これはすごい便利。エンドユーザー向けでだれでも使えそう。(2)はてなブログをはじめよう!kusanagikさんは、はてなブログを使っています。ブログを報告するスターをつけました引用をストックしましたストック一覧を見る閉じる引用するにはまずログインしてくださいログイン閉じる引用をストックできませんでした。再度お試しください閉じる限定公開記事のため引用できません。