df$freq <- as.factor(df$freq)

df$time <- as.factor(df$time)

我學(xué)會(huì)了如何使用循環(huán)繪制一個(gè)因子級(jí)別的直方圖：

windows(width=19, height=10)

par(las=1, cex.lab=0.75, cex.axis=0.6, bty="n", mgp=c(1, 0.6, 0),

oma=c(2, 4, 2, 0) + 0.1, mar=c(4, 0, 3, 3) + 0.1)

a <- layout(matrix(c(1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17,

18, 19, 20, 21), nrow=3, ncol=7, byrow=T))

layout.show(a)

for (i in 1:length(unique(df$freq))) {

value <- subset(df, freq == unique (df$freq)[i])

hist(value$thr, main=paste0("freq: ", unique(df$freq)[i]))

}

我為兩個(gè)因素嘗試了這個(gè)循環(huán)的變體，但遺憾的是這不起作用：

for (i in 1:length(unique(df[c("freq", "time")]))) {

value <- subset(df, freq == unique (df$freq)[i] & time == unique(df$time)[i])

hist(value$thr, main=paste0("freq: ", unique(df$freq)[i]))

}

我還想學(xué)習(xí)如何根據(jù)兩個(gè)因素（不僅僅是一個(gè)）的水平標(biāo)記每個(gè)直方圖......