簡體   English   中英

如何使用facet_grid(ggplot2)在每個構面上繪制多條線

[英]How to plot multiple lines per facet using facet_grid (in ggplot2)

我的問題

理想情況下,我想使用ggplot2創建一個具有6個小平面(角度和組的每個唯一組合為3 x 2,請參見下面的數據)和每個小平面有4條線(並排為2 x 2,請參見下面的數據)的圖),其中每行的x和y值為旋轉和predvalue。 但是,將兩者都包含在一個圖中會使數據消失。

抱歉,由於信譽欠佳,我無法發布圖片,因此我已將其編輯出去,以后將嘗試插入它們-當然,如果有人對其進行編輯也可以。

我的數據如下所示(請參閱本文末尾的完整數據說明):

> 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)

暫無
暫無

聲明:本站的技術帖子網頁,遵循CC BY-SA 4.0協議,如果您需要轉載,請注明本站網址或者原文地址。任何問題請咨詢:yoyou2525@163.com.

 
粵ICP備18138465號  © 2020-2024 STACKOOM.COM