r - 如何使用 facet_grid 为每个构面绘制多条线(在 ggplot2 中)
问题描述
我的问题
理想情况下,我想使用 ggplot2 创建一个具有 6 个面(3 x 2,对于角度和组的每个独特组合 - 请参阅下面的数据)和每个面四条线(2 x 2 用于并排视图,请参阅下面的数据)的图),其中每条线的 x 和 y 值是旋转和预置值。但是,将两者都包含在一个图中会使数据消失。
抱歉,由于缺乏声誉,我无法发布图像,所以我现在已经将它们编辑出来,稍后会尝试插入它们 - 当然,如果其他人编辑它们也可以。
我的数据看起来像这样(请参阅本文末尾的完整数据的 dput):
> head(sub.y)
side view rotation angle age predval
1706 l back 120 0 old 1.322787
1694 l back 120 300 old 1.847914
1739 l back 120 60 old 1.332836
1744 l back 240 0 old 1.157399
1725 l back 240 300 old 1.540411
1713 l back 240 60 old 1.165085
我可以在一个方面产生(荒谬的,因为它们包含重复的旋转)线:
tmp.fig = ggplot(sub.y, aes(x = rotation, y = predval))
tmp.fig + geom_line(aes(group = view:side, colour = side, linetype = view))
我也可以创建构面,但如果我将它们添加到绘图中,所有线条都会消失:
tmp.fig = ggplot(sub.y, aes(x = rotation, y = predval))
tmp.fig + geom_line(aes(group = view:side, colour = side, linetype = view))
tmp.fig + facet_grid(rows = angle ~ age)
如果我通过将每个交互减少到一个因素(不包括图像)来简化绘图,也会发生同样的情况:
sub.y2 = droplevels(subset(sub.y,age == "young" & side == "r"))
tmp.fig = ggplot(sub.y2, aes(x = rotation, y = predval))
tmp.fig + geom_line(aes(group = view))
tmp.fig + facet_grid(rows = vars(angle))
我最接近期望的结果是按角度和年龄进行子集化,并产生一个单一的方面:
sub.single = subset(sub.y, angle =="60" & age == "young")
tmp.fig = ggplot(sub.single, aes(x = rotation, y = predval, group = comb, colour = s))
tmp.fig + geom_line(aes(group = comb, colour = side, linetype = view))
我尝试过的其他事情
- 旋转作为数值变量或无序因子也是如此。
- 使用 facet_grid 的新符号,即:
facet_grid(row = vars(angle), cols = vars(age))
- 将 side 和 view 组合成一个变量并将其用于分组,例如:
sub.y$comb = factor(paste0(sub.y$side, sub.y$view))
tmp.fig = ggplot(sub.y, aes(x = rotation, y = predval))
tmp.fig + geom_line(aes(group = view))
tmp.fig + facet_grid(rows = angle ~ age)
- 子集以消除角度和年龄的每个独特组合内的旋转重复:
sub.lpalm = subset(sub.y,comb == "lpalm")
- 还结合指定颜色和线型的固定值:
tmp.fig = ggplot(sub.lpalm, aes(x = rotation, y = predval))
tmp.fig + geom_line(aes(x = rotation, y = predval), colour = "red", linetype = 2)
tmp.fig + facet_grid(rows = angle ~ age)
我的问题
- 是否可以直接在 ggplot2 中执行此操作?
- 如果是:如何?
- 如果不是:最好的解决方法是什么?ggplot2中有一个吗?还是我需要切换到 par(mfrow)/layout/split.screen?
这是完整的数据集:
sub.y = structure(list(side = structure(c(1L, 1L, 1L, 1L, 1L, 1L, 1L,
1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L,
1L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L,
2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 1L, 1L, 1L, 1L, 1L, 1L, 1L,
1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L,
1L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L,
2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L), .Label = c("l", "r"), class = "factor"),
view = structure(c(1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L,
1L, 1L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 1L,
1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L,
2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 1L, 1L, 1L, 1L, 1L, 1L, 1L,
1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L,
2L, 2L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 2L,
2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L), .Label = c("back",
"palm"), class = "factor"), rotation = structure(c(2L, 2L,
2L, 3L, 3L, 3L, 4L, 4L, 4L, 1L, 1L, 1L, 2L, 2L, 2L, 3L, 3L,
3L, 4L, 4L, 4L, 1L, 1L, 1L, 2L, 2L, 2L, 3L, 3L, 3L, 4L, 4L,
4L, 1L, 1L, 1L, 2L, 2L, 2L, 3L, 3L, 3L, 4L, 4L, 4L, 1L, 1L,
1L, 2L, 2L, 2L, 3L, 3L, 3L, 4L, 4L, 4L, 1L, 1L, 1L, 2L, 2L,
2L, 3L, 3L, 3L, 4L, 4L, 4L, 1L, 1L, 1L, 2L, 2L, 2L, 3L, 3L,
3L, 4L, 4L, 4L, 1L, 1L, 1L, 2L, 2L, 2L, 3L, 3L, 3L, 4L, 4L,
4L, 1L, 1L, 1L), .Label = c("60", "120", "240", "300"), class = c("ordered",
"factor")), angle = structure(c(1L, 2L, 3L, 1L, 2L, 3L, 1L,
2L, 3L, 1L, 2L, 3L, 1L, 2L, 3L, 1L, 2L, 3L, 1L, 2L, 3L, 1L,
2L, 3L, 1L, 2L, 3L, 1L, 2L, 3L, 1L, 2L, 3L, 1L, 2L, 3L, 1L,
2L, 3L, 1L, 2L, 3L, 1L, 2L, 3L, 1L, 2L, 3L, 1L, 2L, 3L, 1L,
2L, 3L, 1L, 2L, 3L, 1L, 2L, 3L, 1L, 2L, 3L, 1L, 2L, 3L, 1L,
2L, 3L, 1L, 2L, 3L, 1L, 2L, 3L, 1L, 2L, 3L, 1L, 2L, 3L, 1L,
2L, 3L, 1L, 2L, 3L, 1L, 2L, 3L, 1L, 2L, 3L, 1L, 2L, 3L), .Label = c("0",
"300", "60"), class = "factor"), age = structure(c(1L, 1L,
1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L,
1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L,
1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L,
1L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L,
2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L,
2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L,
2L, 2L, 2L, 2L), .Label = c("old", "young"), class = "factor"),
predval = c(1.32278735306685, 1.84791387068062, 1.33283630791921,
1.15739906981041, 1.54041056424184, 1.16508497061474, 1.07521745696964,
1.39817944551224, 1.08184750079453, 1.12479365463523, 1.48318811113428,
1.1320512540569, 1.39380248667788, 1.26377942966172, 1.30745255320685,
1.69059385166861, 1.50302783162468, 1.56520862954539, 1.54867304027266,
1.38979681736988, 1.44279652008074, 1.43351747463353, 1.29634372314341,
1.34233748386223, 1.12608718760191, 1.48543811583713, 1.13336154334054,
1.21240995016786, 1.63941177663727, 1.22084649460943, 1.07662100847547,
1.40055372523695, 1.08326842636178, 0.994888143137882, 1.26532710222997,
1.00056192025627, 1.77947866922719, 1.57287633657913, 1.64110190502596,
1.30305741118218, 1.18871937048205, 1.2272796256123, 1.32008919689553,
1.20287711129925, 1.24237660794742, 1.6083792173394, 1.43769159033688,
1.4944817898592, 1.21942340929136, 1.65226153894706, 1.22795818578818,
1.0774857498948, 1.40201746905525, 1.08414388349372, 1.00591000333728,
1.28320933888409, 1.01171055708072, 1.04917234385345, 1.35445633709304,
1.05548413456299, 1.27952165685907, 1.16910164994662, 1.20637971253239,
1.52534725655376, 1.37098239147552, 1.42253022462099, 1.40885870946257,
1.27614512965268, 1.32069215403839, 1.31291299288743, 1.19691581894265,
1.23601841416658, 1.05029770489711, 1.35633247083872, 1.05662308396455,
1.12500641475433, 1.48355807827798, 1.13226676891647, 1.00713834011054,
1.28520892651887, 1.01295310979221, 0.935262602075179, 1.17042599521691,
0.94027496475633, 1.59733529275344, 1.4288609015962, 1.48494199557954,
1.20263686648971, 1.10457968208415, 1.13779812715232, 1.21713008591452,
1.11679383676805, 1.15076228140728, 1.4580995773513, 1.31641344398338,
1.36386842665239)), row.names = c(1706L, 1694L, 1739L, 1744L,
1725L, 1713L, 1733L, 1720L, 1729L, 1728L, 1730L, 1708L, 1761L,
1763L, 1742L, 1804L, 1787L, 1712L, 1778L, 1722L, 1709L, 1714L,
1747L, 1699L, 1751L, 1737L, 1779L, 1726L, 1935L, 1781L, 1719L,
1905L, 1723L, 1735L, 1760L, 1762L, 1780L, 1836L, 1773L, 1767L,
1768L, 1785L, 1753L, 1701L, 1695L, 1774L, 1783L, 1765L, 5621L,
5618L, 5592L, 5598L, 5596L, 5557L, 5612L, 5565L, 5532L, 5540L,
5562L, 5578L, 5539L, 5567L, 5544L, 5560L, 5704L, 5611L, 5607L,
5599L, 5603L, 5555L, 5561L, 5576L, 5552L, 5604L, 5615L, 5564L,
5693L, 5593L, 5569L, 5556L, 5614L, 5580L, 5610L, 5591L, 5714L,
5589L, 5602L, 5597L, 5590L, 5729L, 5546L, 5617L, 5716L, 5726L,
5720L, 5642L), class = "data.frame")
解决方案
中间代码中有一个小错字,您打算使用tmp.fig = tmp.fig +...
但缺少tmp.fig =
. 这意味着+ geom_line()
第二行中的不包含在第三行中,因此图中没有行。
sub.y2 = droplevels(subset(sub.y,age == "young" & side == "r"))
tmp.fig = ggplot(sub.y2, aes(x = rotation, y = predval))
tmp.fig + geom_line(aes(group = view))
tmp.fig + facet_grid(rows = vars(angle))
选项 1:添加tmp.fig = ...
以增量修改绘图对象
sub.y2 = droplevels(subset(sub.y,age == "young" & side == "r"))
tmp.fig = ggplot(sub.y2, aes(x = rotation, y = predval))
tmp.fig = tmp.fig + geom_line(aes(group = view))
tmp.fig = tmp.fig + facet_grid(rows = vars(angle))
选项 2:在一个连接的流中构建绘图
ggplot(sub.y, aes(x = rotation, y = predval)) +
geom_line(aes(group = view:side, colour = side, linetype = view)) +
facet_grid(rows = angle ~ age)
推荐阅读
- jquery - 从页面上所有表格中的特定单元格获取值
- node.js - 一旦自适应卡出现在微软聊天机器人中,如何禁用输入字段?
- vue.js - 无论数据中的值如何,都看不到带有 v-if 的 div
- angular - 如何访问拦截器内部的错误响应主体?
- android - 如何在 BottomNavigationView 上方显示 Snackbar 并移动 FAB?
- ios - 如何在 Swift 的 collectionView 中显示每个类别的多个类别?
- python - 如何比较 2 个数据框并生成新的数据框
- java - 如何停止下载管理器下载我的下载中已经存在的文件?
- php - 自动完成引导 php mysql
- firebase - firebase_messaging 如何清除通知?