我正在R中绘制散点图,并想要添加箱线图。x轴是日期,y轴权重。我想随时间显示体重,因此每周都应该有一个箱线图。但是,箱形图的宽度太小。这是我目前拥有的:
箱图绘制如下:
data <- read.table("weight2.txt", header=TRUE)
data$datetime <- strptime(paste(data$Date, data$Time), "%d.%m.%Y %H:%M")
data$week <- strftime(data$datetime,format="%W")
data$timestamp <- as.numeric(data$datetime)
plot(data$datetime, data$Weight, xlab="Date",
ylab="Weight", ylim=c(61,68))
library(plyr)
dt <- data.table(data)
aggrWeek <- ddply(dt,~week,summarise,
min=min(Weight),
firstQ=quantile(Weight,0.25),
mean=mean(Weight),
thirdQ=quantile(Weight,0.75),
max=max(Weight),
timestamp=mean(timestamp))
aggrWeek$datetime <- as.POSIXct(aggrWeek$timestamp, origin="1970-01-01")
boxplotData <- t(
data.frame(aggrWeek$min, aggrWeek$firstQ, aggrWeek$mean, aggrWeek$thirdQ, aggrWeek$max))
bxp(list(
stats=data.matrix(boxplotData), n=rep(1,ncol(boxplotData))),
add=TRUE, at=aggrWeek$datetime, show.names=FALSE)
如何使箱形图的宽度更大?
就像我想的那样,您的x值是POSIXt日期,自1970年1月1日起以秒为单位存储。因此,如果宽度为1,则实际上是在绘制一个“ 1秒”宽的箱形图,该范围很小。如果您想要宽约“ 1天”的箱形图,请尝试
bxp(list(
stats=data.matrix(boxplotData), n=rep(1,ncol(boxplotData))),
add=TRUE, at=aggrWeek$datetime, show.names=FALSE, boxwex=60*60*24)
通过一个可复制的示例,看看这有多容易!
本文收集自互联网,转载请注明来源。
如有侵权,请联系[email protected] 删除。
我来说两句