我有一个包含两列的矩阵,我想要分别绘制这两列的箱线图,但是它们的条目数量不同。例如,第一列有10个条目,而第二列只有7个条目,剩下的3个条目为零。为了比较,我想将它们并排绘制。是否有办法告诉R只对第一列绘制整个箱线图,而对第二列仅绘制前7个条目的箱线图?
## dummy version of your data
mat <- matrix(c(1:17, rep(0, 3)), ncol = 2)
## create object suitable for plotting with boxplot
## I.e. convert to melted or long format
df <- data.frame(values = mat[1:17],
vars = rep(c("Col1","Col2"), times = c(10,7)))
## draw the boxplot
boxplot(values ~ vars, data = df)
在上面,我按照你的说法假设你有一个矩阵。如果你实际上拥有的是数据框,则需要
df <- data.frame(values = c(mat[,1], mat[1:7, 2]),
vars = rep(c("Col1","Col2"), times = c(10,7)))
我假设两列数据是可比较的,因为这两列数据的存在意味着它们是一个分类变量,可以通过该变量将值分割(例如男女身高,以性别为分类变量)。
生成的箱线图如下所示:
## Load data from CSV; first row contains column headers
dat <- read.csv( 'your-filename.csv', header = T )
## Set plot region (when set 'ylim' skip first row with headers)
plot(
1, 1,
xlim=c(1,ncol(dat)), ylim=range(dat[-1,], na.rm=T),
xaxt='n', xlab='', ylab=''
)
axis(1, labels=colnames(dat), at=1:ncol(dat))
for(i in 1:ncol(dat)) {
## Get i-th column
p <- dat[,i]
## Remove 0 values from column
p <- p[! p %in% 0]
## Instead of 0 you can use any values
## For example, you can remove 1, 2, 3
## p <- p[! p %in% c(1,2,3)]
## Draw boxplot
boxplot(p, add=T, at=i)
}
mat
的前17个值,但你的其余部分是正确的。 - Gavin Simpsonboxplot(values ~ var, data = df)
应该是boxplot(values ~ vars, data = df)
。否则,这个函数运行得很好。 - bdetweiler