1
イム、私のコードは次のようになります。TukeyHSD特定の条件
#----------------------------------------------------------------------------------------#
# RING data:
#----------------------------------------------------------------------------------------#
library(doBy)
# Set working directory
setwd("")
#### Read data & Converting factors ####
dat <- read.table("afstand.txt", header =TRUE)
str(dat)
dat$Vial <- as.factor(dat$Vial)
dat$Line <- as.factor(dat$Line)
dat$Fly <- as.factor(dat$Fly)
dat$Temp <- as.factor(dat$Temp)
str(dat)
datSUM <- summaryBy(X0.5_sec+X1_sec+X1.5_sec+X2_sec+X2.5_sec+X3_sec~Vial_nr+Concentration+Sex+Line+Vial+Temp,data=dat, FUN=sum)
fl<-levels(datSUM$Line)
aov1 <- aov(X0.5_sec.sum ~ Concentration*Sex*Line*Temp, data=datSUM)
summary(aov1) #Overview of model
TukeyHSD(aov1, 'Line',ordered = TRUE, conf.level = 0.95)
私がやりたいどのようなインスタンスのラインと温度の間の相互作用を見て、私場合であります[ reached getOption("max.print") -- omitted 3716 rows ] TukeyHSD（aov1）を実行すると、このエラーが発生します。[ reached getOption("max.print") -- omitted 3716 rows ] LineとTempの間だけでテストしたいのですが、すべての組み合わせではなく、重要な結果を表示する方法TukeyHSD（avo1）を実行するだけですか？ 私はTukeyHSD(aov1, 'Line,Temp',ordered = TRUE, conf.level = 0.95),TukeyHSD(aov1, 'Line':'Temp',ordered = TRUE, conf.level = 0.95)とTukeyHSD(aov1, 'Line'&'Temp',ordered = TRUE, conf.level = 0.95)を使用しようとしましたが、運はありません。変数LineとTemp間の唯一のショーの相互作用に
structure(list(Concentration = structure(c(2L, 7L, 7L, 1L, 7L,
1L, 2L, 1L, 7L, 1L, 4L, 2L, 2L, 1L, 2L, 4L, 7L, 2L, 2L, 1L), .Label = c("a",
"b", "c", "d", "e", "x", "y"), class = "factor"), Sex = structure(c(1L,
2L, 2L, 2L, 1L, 2L, 2L, 1L, 1L, 2L, 1L, 1L, 1L, 1L, 2L, 1L, 1L,
1L, 2L, 1L), .Label = c("f", "m"), class = "factor"), Line = structure(c(3L,
3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 2L, 3L, 3L, 3L, 3L, 2L, 3L,
3L, 3L, 3L), .Label = c("20", "23", "40", "73"), class = "factor"),
Temp = structure(c(1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L,
2L, 1L, 2L, 1L, 1L, 2L, 1L, 2L, 1L, 1L), .Label = c("23",
"29"), class = "factor"), X0.5_sec.sum = c(107.19, 46.17,
58.05, 75.87, 87.75, 71.55, 46.17, 47.25, 22.41, 31.05, 27.36,
79.11, 87.58, 21.33, 34.29, 60.4, 85.05, 72.47, 114.21, 67.77
)), .Names = c("Concentration", "Sex", "Line", "Temp", "X0.5_sec.sum"
), row.names = c(NA, 20L), class = "data.frame")
再現可能なコードなしで助けがたいですが、 'which = c（" Line "、" Temp "）'をあなたの呼び出しの中に入れるべきでしょう。 – mtoto
コードのdput20を追加しました。私も'which = c（ "Line"、 "Temp"） 'を試してみましたが、それらの組み合わせでそれらを並べて実行します。 –
：あなたの提案で、問題を解決しました。 –