我正在使用软件包VGAM的帮助页面中的以下示例
library(VGAM)
fit4 <- vgam(BMI ~ s(age, df = c(4, 2)), lms.bcn(zero = 1), data = bmi.nz, trace = TRUE)
qtplot(fit4, percentiles = c(5,50,90,99), main = "Quantiles", las = 1, xlim = c(15, 90), ylab = "BMI", lwd = 2, lcol = 4)
我正在用它获得一个合适的图表:
如何避免从图表中绘制点?此外,我需要在每个年龄段20,30,40 ... 80打印出这些百分位数的值(单独作为表格)。如何才能做到这一点?我可以在这里使用ggplot()格式命令而不是qtplot()命令吗?谢谢你的帮助。
答案 0 :(得分:1)
这样的事情怎么样:
# required packages
library(VGAM)
require(reshape2)
require(ggplot2)
# fitted values from vgam
fit4 <- vgam(BMI ~ s(age, df = c(4, 2)), lms.bcn(zero = 1), data = bmi.nz, trace = TRUE)
fitted.values <- data.frame(qtplot.lmscreg(fit4, percentiles = c(5,50,90,99))$fitted.values)
fitted.values[, 'age'] <- bmi.nz[, 'age']
# melt data.frame
dmelt <- melt(fitted.values, id.vars='age')
# ploting
ggplot(dmelt, aes(age, value, group=variable)) +
geom_line(color='blue') +
annotate(geom='text',
x = max(bmi.nz[, 'age']) + 3,
y = unlist(fitted.values[which.max(fitted.values[, 'age']), -ncol(fitted.values)]),
label=c(' 5%', '50%', '90%', '99%')) +
lapply(2:8*10, function(i) {
annotate(geom='text',
x = i,
y = 1+unlist(fitted.values[which.min(abs(fitted.values[, 'age'] - i)), -ncol(fitted.values)]),
label=paste0(round(unlist(fitted.values[which.min(abs(fitted.values[, 'age'] - i)), -ncol(fitted.values)]),1), '%'))
}) +
scale_y_continuous('BMI') +
theme_bw(base_size=16)