假设我有以下数据:
datapoints1 = data.frame(categ=c(rep(1, n), rep(2, n)), vals1=c(rt(n, 1, 2), rnorm(n, 3, 4)))
datapoints2 = data.frame(categ=c(rep(1, n), rep(2, n)), vals2=c(rt(n, 5, 6), rnorm(n, 7, 8)))
使用 ggplot2，我如何使用 facet 功能在单个命令中创建两个 QQplots，即一个具有两个 t 样本, 另一个是两个高斯样本？
最佳答案
首先，合并两个数据框:
dat <- cbind(datapoints1, vals2 = datapoints2[ , 2])
然后，对数据进行排序:
dat_sort <- do.call("rbind", lapply(unique(dat$categ), FUN = function(x) {data.frame(categ = x, vals1 = sort(dat$vals1[dat$categ == x]), vals2 = sort(dat$vals2[dat$categ == x]))}))
如果两个样本向量的长度相同则很简单:
ggplot() +
geom_point(data = dat_sort, aes(x = vals1, y = vals2)) +
facet_wrap( ~ categ, scales = "free")
n = 1000 的示例:
关于r - 使用 ggplot2 的多面 qqplots，我们在Stack Overflow上找到一个类似的问题： https://stackoverflow.com/questions/12369377/