1要因の場合

  1. Excelを起動し、1列目に要因、2列目にデータを入力する。
Treatment Measurement
A 2
A 3
A 4
B 8
B 9
B 7
  1. Rに以下のスクリプトをペーストする。ただし、リターンを含めていはいけない。
a <- read.delim("clipboard", header=F)
  1. Excelに戻りデータを選択し、コピーする。ただし、ヘッダー（タイトル）を含んではいけない。
  2. Rに移動しリターンを押す。
  3. Rで以下のスクリプトをペーストする。
library(psych)
library(agricolae)
anova <- aov(V2 ~ V1, data=a)
summary(anova)
(describeBy(a$V2, list(a$V1), mat=T))
(HSD.test(anova, "V1", group=T))
 
結果をコピーしExcelのシートに張り付け、「テキストウィザードで編集」を選択し、スペース区切りを選ぶ。
 
すぐに傾向を見たいのなら
b <- describeBy(a$V2, list(a$V1), mat=T)
par(mar=c(5, 4, 4, 2)) #x軸のラベルが切れる場合は1要素目を変更
c <- barplot(b$mean, names.arg=b$group1, ylim=c(0,max(b$mean+b$se)), las=3)
arrows(c, b$mean-b$se, c, b$mean+b$se, angle=90, length=0.1)
arrows(c, b$mean+b$se, c, b$mean-b$se, angle=90, length=0.1)
 

2要因の場合

  1. Excelを起動し、1列目と2列目に要因、3列目にデータを入力する。
  2. Rを起動し、ファイル->作業スペースの保存をクリックし保存先フォルダを選択する。
  3. Rに以下のスクリプトをペーストする。ただし、リターンを含めていはいけない。

a <- read.delim("clipboard", header=F)

  1. Excelに戻りデータを選択し、コピーする。ただし、ヘッダー（タイトル）を含んではいけない。
  2. Rに移動しリターンを押す。
  3. Rで以下のスクリプトをペーストする。
library(psych)
library(agricolae)
anova <- aov(V3 ~ V1*V2, data=a)
summary(anova)
a2 <- cbind(paste(a$V1, a$V2, sep="_"), a)
colnames(a2) <- c("treatment", "treatment1", "treatment2", "measurement")
anova2 <- aov(measurement ~ treatment, data=a2)
(describeBy(a$V3, list(a$V1, a$V2), mat=T))
(HSD.test(anova2, "treatment", group=T))
 
結果をコピーしExcelのシートに張り付け、「テキストウィザードで編集」を選択し、スペース区切りを選ぶ。
 
すぐに傾向を見ないのなら

b <- describeBy(a$V3, list(a$V1, a$V2), mat=T)
par(mar=c(5, 4, 4, 2)) #x軸のラベルが切れる場合は1要素目を変更
c <- barplot(b$mean, names.arg=paste(b$group2, b$group1, sep=", "), ylim=c(0,max(b$mean+b$se)), las=3)
arrows(c, b$mean-b$se, c, b$mean+b$se, angle=90, length=0.1)
arrows(c, b$mean+b$se, c, b$mean-b$se, angle=90, length=0.1)

3要因の場合

  1. Excelを起動し、1列目と2列目、3列目に要因、4列目にデータを入力する。
  2. Rを起動し、ファイル->作業スペースの保存をクリックし保存先フォルダを選択する。
  3. Rに以下のスクリプトをペーストする。ただし、リターンを含めていはいけない。
a <- read.delim("clipboard", header=F)
  1. Excelに戻りデータを選択し、コピーする。ただし、ヘッダー（タイトル）を含んではいけない。
  2. Rに移動しリターンを押す。
  3. Rで以下のスクリプトをペーストする。
library(psych)
library(agricolae)
anova <- aov(V4 ~ V1*V2*V3, data=a)
summary(anova)
a2 <- cbind(paste(a$V1, a$V2, a$V3, sep="_"), a)
colnames(a2) <- c("treatment", "treatment1", "treatment2","treatment3", "measurement")
anova2 <- aov(measurement ~ treatment, data=a2)
(describeBy(a$V4, list(a$V1, a$V2, a$V3), mat=T))
(HSD.test(anova2, "treatment", group=T))
 
結果をコピーしExcelのシートに張り付け、「テキストウィザードで編集」を選択し、スペース区切りを選ぶ。

すぐに傾向を見たいのなら

b <- describeBy(a$V4, list(a$V1, a$V2, a$V3), mat=T)
par(mar=c(5, 4, 4, 2)) #x軸のラベルが切れる場合は1要素目を変更
c <- barplot(b$mean, names.arg=paste(b$group3, b$group2, b$group1, sep=", "), ylim=c(0,max(b$mean+b$se)), las=3)
arrows(c, b$mean-b$se, c, b$mean+b$se, angle=90, length=0.1)
arrows(c, b$mean+b$se, c, b$mean-b$se, angle=90, length=0.1)