反復測定二元配置分散分析をR（carパッケージ）でやってみる（S. Kubota）

例えば，業務ソフトAとBを３カ月ずつ使用し，どちらが業務効率が良かったかを検証したいとする．つまり，業務ソフト違いによる業務効率の影響と使い続けることでの時系列での業務効率の変化を検討する．このとき，業務ソフトの違い，時系列の変化をそれぞれ要因という．また，この検証を同一の社員で，AとBのソフトをそれぞれ３か月ずつ使い検証したこととする．この場合は，同一の社員で2つの要因について検証を行ったことになるので，反復測定二元配置分散分析というのを行う．今回はRとcarパッケージを使ってみる．データの並びは以下の通り．エクセルにあるデータを読みこんでいく前提で記述する．データは行に被験者，列に要...

混合計画の分散分析

混合計画（混合モデルではない）の二元配置分散分析。要するに被験者間、被験者内計画。car packageでMANOVAで実行してみる。詳細は下記の通り。> dat<-read.delim("clipboard") #データ読み込み> head(dat) #データの冒頭部だけ見せるとこんな感じ．データの中身はこんな感じ． gender trial1 trial2 trial31 male 10.0121 14.0931 15.87682 male 16.1474 15.9080 20.84173 male 57.0574 47.2407 75.45794 male ...

Rを使うために入れておくべきソフトとパッケージ

データサイエンスにRを使うのに便利なソフトとパッケージを挙げておく．Rはフリーで高度な統計学的手法を用いることができるのが特徴．なお以下は最小限の構成．・Rstudio　Rのための統合開発環境　一度使うとやめられなくなる．・dplyrとtidyr　データラングリングにはかかせないパッケージ．・ggplot2　作図のためのパッケージ．このパッケージを使って作図するときれいな図が作れる．Written by S. Kubota

相関行列の分析をしてデータの見通しを立てる-Hmisc，PerformanceAnalytics　パッケージを使用- (S. Kubota)

相関行列の見通しを立てるのに，corrplotを用いるのもいいが，他の方法も記述しておく．まずはデータ，以下のようなデータをエクセルから読み込んだ（corrplotのメモと同じデータ）．dat<-read .delim=("clipboard")dat EMG1 EMG2 EMG3 EMG4 time force fatigue comfort1 0.020 0.007 0.015 0.020 1.70 361.132 68.6 60.62 0.039 0.008 0.010 0.012 1.36 293.362 65.0 48.23 0.0...

ezパッケージで簡単に反復測定分散分析-一元配置と二元配置- (S. Kubota)

まずはezパッケージをインストールしておく．データの読み込み．Excelなんかでデータを成形しておき，読み込む．データは以下の通り．同一被験者で2つの要因について組み合わせて実験を行った，という設定のデータである． subjectsは被験者，factor1と2が要因，valueが値．twoway<-read.delim("clipboard")twoway subjects factor1 factor2 value1 a A N 4.152 b A N 4.083 c A...