在R的單個圖中使用ggplot2繪制多盒圖以重復測量238個受試者

Question

我有238名患者。 每位患者每月進行一次就診，並進行排尿前尿液估計和排尿后尿液估計。 現在，在近兩年的時間里，每位患者都進行了23次檢查。我想為整個人群繪制尿液前后尿液的箱線圖，以便在一張圖中查看隨時間變化的趨勢（每次訪問）。 前空隙體積通常為300+，而后空隙體積小於50。 有沒有辦法在ggplot2中做到這一點 ？ 我只能使用boxplot制作2個單獨的圖形（一個圖形隨着時間的流逝在一個圖形中顯示多個prevoid，中位數通過線連接，另一個圖形顯示與后空隙相似的東西）.....實際上，我什至無法用ggplot繪制一個方框圖因為我的數據沒有xy軸

我嘗試了plotrix，但沒有...什么都沒有。

這種代碼與boxplot()有點配合

boxplot((mydata)[,c(13:35)], outline=FALSE, xlab= "Time in months", ylab="Prevoid urine",las=2)       ##plots box plots of columns 13-35 in one graph(23 boxplots)with outliers suppressed

p<-lines(1:23,(boxplot((mydata)[,c(13:35)],outline=FALSE, col="yellow", xlab= "Time in months", ylab="prevoid urine", las=2,ylim = c(0, 600)))$stats[3,], col="blue", lwd=3)        ##plots blue line connecting medians of 23 boxplots shaded yellow

Answer 1

ggplot中的facet選項對您有用嗎？

Answer 2

考慮調用sapply計算median ，然后傳遞載體導入lines ，以避免調用boxplot再次進行第二不必要的情節：

## plots box plots of columns 13-35 in one graph(23 boxplots)with outliers suppressed
boxplot(mydata[,c(13:35)], outline=FALSE, xlab= "Time in months", ylab="Prevoid urine",las=2)   

## plots blue line connecting medians of 23 boxplots shaded yellow
lines(1:23, sapply(mydata[,c(13:35)], median), col="blue", lwd=3) 

為了演示隨機數據：

set.seed(8132019)
mydata <- data.frame(replicate(35, runif(500, 300, 500)))

boxplot(...)
lines(...)

Answer 3

您的數據非常混亂，我們沒有大量有關變量順序的信息。 對於您提供的刪除答案，這是我所能做的。

數據

mydata <-data.frame(   serialnumber = 1:238,
                       urinarysummarypreop = c(68.75, 87.5, 81.25, 81.25, 81.25, 81.25, 87.5, 68.75, 81.25, 81.25, 81.25, 87.5, 81.25, 87.5, 81.25, 87.5, 68.75, 87.5, 87.5, 81.25, 81.25, 68.75, 68.75, 81.25, 87.5, 87.5, 87.5, 87.5, 81.25, 81.25, 81.25, 68.75, 81.25, 87.5, 81.25, 68.75, 77.08333333, 81.25, 77.08333333, 68.75, 81.25, 81.25, 68.75, 81.25, 87.5, 81.25, 68.75, 81.25, 68.75, 68.75, 81.25, 81.25, 81.25, 81.25, 68.75, 87.5, 87.5, 87.5, 87.5, 87.5, 81.25, 81.25, 87.5, 87.5, 87.5, 87.5, 81.25, 77.08333333, 87.5, 77.08333333, 87.5, 87.5, 68.75, 87.5, 87.5, 81.25, 81.25, 87.5, 87.5, 87.5, 87.5, 87.5, 87.5, 87.5, 81.25, 87.5, 81.25, 81.25, 87.5, 87.5, 87.5, 81.25, 81.25, 81.25, 87.5, 81.25, 87.5, 87.5, 77.08333333, 77.08333333, 87.5, 87.5, 87.5, 81.25, 81.25, 87.5, 77.08333333, 87.5, 87.5, 87.5, 81.25, 87.5, 81.25, 68.75, 81.25, 87.5, 87.5, 68.75, 87.5, 81.25, 87.5, 87.5, 81.25, 87.5, 87.5, 87.5, 87.5, 81.25, 77.08333333, 81.25, 87.5, 77.08333333, 87.5, 81.25, 68.75, 81.25, 81.25, 81.25, 81.25, 77.08333333, 87.5, 87.5, 87.5, 87.5, 87.5, 77.08333333, 87.5, 87.5, 77.08333333, 81.25, 77.08333333, 87.5, 77.08333333, 81.25, 87.5, 81.25, 87.5, 81.25, 87.5, 87.5, 87.5, 87.5, 81.25, 81.25, 87.5, 77.08333333, 87.5, 81.25, 81.25, 77.08333333, 81.25, 87.5, 68.75, 81.25, 87.5, 68.75, 68.75, 87.5, 77.08333333, 87.5, 87.5, 81.25, 87.5, 77.08333333, 81.25, 87.5, 87.5, 77.08333333, 77.08333333, 81.25, 87.5, 87.5, 81.25, 87.5, 87.5, 81.25, 87.5, 81.25, 87.5, 87.5, 87.5, 87.5, 81.25, 81.25, 87.5, 87.5, 87.5, 81.25, 81.25, 81.25, 81.25, 87.5, 77.08333333, 87.5, 87.5, 77.08333333, 87.5, 81.25, 87.5, 87.5, 77.08333333, 81.25, 81.25, 81.25, 87.5, 81.25, 81.25, 68.75, 87.5, 87.5, 81.25, 68.75, 87.5, 81.25, 81.25, 81.25, 87.5, 81.25),
                       urinarysummary6wks = c(26.33333333, 26.33333333, 26.33333333, 22.16666667, 26.33333333, 26.33333333, 22.16666667, 26.33333333, 26.33333333, 26.33333333, 26.33333333, 26.33333333, 26.33333333, 26.33333333, 22.16666667, 26.33333333, 26.33333333, 26.33333333, 22.16666667, 26.33333333, 26.33333333, 22.16666667, 22.16666667, 26.33333333, 22.16666667, 26.33333333, 22.16666667, 26.33333333, 26.33333333, 26.33333333, 26.33333333, 26.33333333, 26.33333333, 26.33333333, 26.33333333, 26.33333333, 26.33333333, 26.33333333, 16.66666667, 26.33333333, 26.33333333, 22.16666667, 22.16666667, 26.33333333, 26.33333333, 26.33333333, 22.16666667, 26.33333333, 26.33333333, 26.33333333, 26.33333333, 26.33333333, 22.16666667, 26.33333333, 22.16666667, 26.33333333, 22.16666667, 26.33333333, 26.33333333, 26.33333333, 26.33333333, 26.33333333, 26.33333333, 26.33333333, 26.33333333, 26.33333333, 26.33333333, 22.16666667, 26.33333333, 26.33333333, 26.33333333, 26.33333333, 22.16666667, 26.33333333, 26.33333333, 26.33333333, 26.33333333, 26.33333333, 22.16666667, 26.33333333, 26.33333333, 26.33333333, 26.33333333, 26.33333333, 26.33333333, 26.33333333, 26.33333333, 26.33333333, 26.33333333, 26.33333333, 26.33333333, 26.33333333, 26.33333333, 26.33333333, 22.16666667, 22.16666667, 26.33333333, 26.33333333, 22.16666667, 22.16666667, 26.33333333, 26.33333333, 22.16666667, 26.33333333, 26.33333333, 26.33333333, 26.33333333, 26.33333333, 26.33333333, 26.33333333, 26.33333333, 26.33333333, 26.33333333, 22.16666667, 26.33333333, 26.33333333, 26.33333333, 22.16666667, 26.33333333, 26.33333333, 26.33333333, 22.16666667, 16.66666667, 22.16666667, 26.33333333, 26.33333333, 26.33333333, 26.33333333, 22.16666667, 26.33333333, 26.33333333, 26.33333333, 26.33333333, 16.66666667, 22.16666667, 26.33333333, 22.16666667, 26.33333333, 26.33333333, 26.33333333, 26.33333333, 26.33333333, 26.33333333, 26.33333333, 26.33333333, 26.33333333, 26.33333333, 26.33333333, 26.33333333, 26.33333333, 22.16666667, 26.33333333, 22.16666667, 26.33333333, 26.33333333, 26.33333333, 26.33333333, 22.16666667, 26.33333333, 26.33333333, 26.33333333, 26.33333333, 26.33333333, 26.33333333, 26.33333333, 22.16666667, 26.33333333, 26.33333333, 26.33333333, 22.16666667, 22.16666667, 26.33333333, 22.16666667, 26.33333333, 22.16666667, 22.16666667, 26.33333333, 22.16666667, 26.33333333, 26.33333333, 26.33333333, 26.33333333, 22.16666667, 22.16666667, 26.33333333, 26.33333333, 26.33333333, 16.66666667, 22.16666667, 26.33333333, 26.33333333, 26.33333333, 26.33333333, 26.33333333, 26.33333333, 26.33333333, 26.33333333, 26.33333333, 26.33333333, 26.33333333, 26.33333333, 26.33333333, 26.33333333, 26.33333333, 26.33333333, 26.33333333, 26.33333333, 26.33333333, 26.33333333, 26.33333333, 26.33333333, 26.33333333, 26.33333333, 26.33333333, 16.66666667, 26.33333333, 22.16666667, 26.33333333, 26.33333333, 26.33333333, 22.16666667, 26.33333333, 26.33333333, 16.66666667, 26.33333333, 26.33333333, 22.16666667, 22.16666667, 26.33333333, 26.33333333, 26.33333333, 22.16666667, 26.33333333, 26.33333333, 26.33333333, 26.33333333, 26.33333333, 26.33333333),
                       urinarysummary3months = c(36.16666667, 57.66666667, 48.66666667, 32, 48.66666667, 26.33333333, 44.5, 36.16666667, 48.66666667, 46.58333333, 36.16666667, 57.66666667, 57.66666667, 48.66666667, 32, 36.16666667, 36.16666667, 36.16666667, 48.66666667, 46.58333333, 36.16666667, 48.66666667, 30.5, 32, 44.5, 36.16666667, 32, 57.66666667, 48.66666667, 36.16666667, 48.66666667, 36.16666667, 57.66666667, 57.66666667, 36.16666667, 36.16666667, 48.66666667, 46.58333333, 22.16666667, 36.16666667, 57.66666667, 44.5, 30.5, 57.66666667, 48.66666667, 36.16666667, 36.16666667, 48.66666667, 36.16666667, 36.16666667, 36.16666667, 36.16666667, 32, 36.16666667, 32, 36.16666667, 32, 36.16666667, 48.66666667, 57.66666667, 36.16666667, 48.66666667, 36.16666667, 36.16666667, 36.16666667, 36.16666667, 46.58333333, 26.33333333, 36.16666667, 36.16666667, 36.16666667, 57.66666667, 48.66666667, 36.16666667, 48.66666667, 36.16666667, 48.66666667, 57.66666667, 32, 36.16666667, 46.58333333, 36.16666667, 36.16666667, 36.16666667, 48.66666667, 48.66666667, 36.16666667, 36.16666667, 46.58333333, 36.16666667, 57.66666667, 36.16666667, 36.16666667, 36.16666667, 32, 57.66666667, 48.66666667, 36.16666667, 26.33333333, 36.16666667, 48.66666667, 48.66666667, 57.66666667, 48.66666667, 57.66666667, 36.16666667, 36.16666667, 36.16666667, 36.16666667, 36.16666667, 46.58333333, 36.16666667, 48.66666667, 41.66666667, 36.16666667, 36.16666667, 36.16666667, 32, 36.16666667, 36.16666667, 48.66666667, 44.5, 36.16666667, 36.16666667, 46.58333333, 36.16666667, 36.16666667, 36.16666667, 26.33333333, 48.66666667, 36.16666667, 48.66666667, 36.16666667, 22.16666667, 32, 57.66666667, 32, 32, 48.66666667, 36.16666667, 36.16666667, 48.66666667, 36.16666667, 36.16666667, 36.16666667, 48.66666667, 48.66666667, 36.16666667, 36.16666667, 36.16666667, 36.16666667, 46.58333333, 46.58333333, 46.58333333, 36.16666667, 48.66666667, 36.16666667, 32, 36.16666667, 48.66666667, 30.5, 48.66666667, 36.16666667, 36.16666667, 48.66666667, 46.58333333, 36.16666667, 36.16666667, 46.58333333, 26.33333333, 36.16666667, 36.16666667, 32, 57.66666667, 32, 30.5, 48.66666667, 32, 57.66666667, 36.16666667, 48.66666667, 36.16666667, 32, 46.58333333, 36.16666667, 46.58333333, 48.66666667, 36.16666667, 46.58333333, 36.16666667, 48.66666667, 36.16666667, 36.16666667, 46.58333333, 36.16666667, 36.16666667, 36.16666667, 46.58333333, 57.66666667, 36.16666667, 46.58333333, 46.58333333, 48.66666667, 36.16666667, 48.66666667, 57.66666667, 36.16666667, 36.16666667, 46.58333333, 48.66666667, 32, 36.16666667, 36.16666667, 32, 16.66666667, 36.16666667, 48.66666667, 46.58333333, 36.16666667, 48.66666667, 41.66666667, 46.58333333, 36.16666667, 16.66666667, 36.16666667, 48.66666667, 32, 36.16666667, 32, 36.16666667, 57.66666667, 30.5, 36.16666667, 36.16666667, 36.16666667, 36.16666667, 36.16666667, 57.66666667 ),
                       urinaryfunctionpreop = c(100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L),
                       urinaryfunction6wks = c(16.5, 16.5, 16.5, 16.5, 16.5, 16.5, 16.5, 16.5, 16.5, 16.5, 16.5, 16.5, 16.5, 16.5, 16.5, 16.5, 16.5, 16.5, 16.5, 16.5, 16.5, 16.5, 16.5, 16.5, 16.5, 16.5, 16.5, 16.5, 16.5, 16.5, 16.5, 16.5, 16.5, 16.5, 16.5, 16.5, 16.5, 16.5, 0, 16.5, 16.5, 16.5, 16.5, 16.5, 16.5, 16.5, 16.5, 16.5, 16.5, 16.5, 16.5, 16.5, 16.5, 16.5, 16.5, 16.5, 16.5, 16.5, 16.5, 16.5, 16.5, 16.5, 16.5, 16.5, 16.5, 16.5, 16.5, 16.5, 16.5, 16.5, 16.5, 16.5, 16.5, 16.5, 16.5, 16.5, 16.5, 16.5, 16.5, 16.5, 16.5, 16.5, 16.5, 16.5, 16.5, 16.5, 16.5, 16.5, 16.5, 16.5, 16.5, 16.5, 16.5, 16.5, 16.5, 16.5, 16.5, 16.5, 16.5, 16.5, 16.5, 16.5, 16.5, 16.5, 16.5, 16.5, 16.5, 16.5, 16.5, 16.5, 16.5, 16.5, 16.5, 16.5, 16.5, 16.5, 16.5, 16.5, 16.5, 16.5, 16.5, 16.5, 0, 16.5, 16.5, 16.5, 16.5, 16.5, 16.5, 16.5, 16.5, 16.5, 16.5, 0, 16.5, 16.5, 16.5, 16.5, 16.5, 16.5, 16.5, 16.5, 16.5, 16.5, 16.5, 16.5, 16.5, 16.5, 16.5, 16.5, 16.5, 16.5, 16.5, 16.5, 16.5, 16.5, 16.5, 16.5, 16.5, 16.5, 16.5, 16.5, 16.5, 16.5, 16.5, 16.5, 16.5, 16.5, 16.5, 16.5, 16.5, 16.5, 16.5, 16.5, 16.5, 16.5, 16.5, 16.5, 16.5, 16.5, 16.5, 16.5, 16.5, 16.5, 16.5, 16.5, 16.5, 0, 16.5, 16.5, 16.5, 16.5, 16.5, 16.5, 16.5, 16.5, 16.5, 16.5, 16.5, 16.5, 16.5, 16.5, 16.5, 16.5, 16.5, 16.5, 16.5, 16.5, 16.5, 16.5, 16.5, 16.5, 16.5, 16.5, 0, 16.5, 16.5, 16.5, 16.5, 16.5, 16.5, 16.5, 16.5, 0, 16.5, 16.5, 16.5, 16.5, 16.5, 16.5, 16.5, 16.5, 16.5, 16.5, 16.5, 16.5, 16.5, 16.5),
                       urinaryfunction3months = c(33.5, 66.75, 58.5, 33.5, 58.5, 16.5, 58.5, 33.5, 58.5, 33.5, 33.5, 66.75, 66.75, 58.5, 33.5, 33.5, 33.5, 33.5, 58.5, 33.5, 33.5, 58.5, 16.5, 33.5, 58.5, 33.5, 33.5, 66.75, 58.5, 33.5, 58.5, 33.5, 66.75, 66.75, 33.5, 33.5, 58.5, 33.5, 16.5, 33.5, 66.75, 58.5, 16.5, 66.75, 58.5, 33.5, 33.5, 58.5, 33.5, 33.5, 33.5, 33.5, 33.5, 33.5, 33.5, 33.5, 33.5, 33.5, 58.5, 66.75, 33.5, 58.5, 33.5, 33.5, 33.5, 33.5, 33.5, 16.5, 33.5, 33.5, 33.5, 66.75, 58.5, 33.5, 58.5, 33.5, 58.5, 66.75, 33.5, 33.5, 33.5, 33.5, 33.5, 33.5, 58.5, 58.5, 33.5, 33.5, 33.5, 33.5, 66.75, 33.5, 33.5, 33.5, 33.5, 66.75, 58.5, 33.5, 16.5, 33.5, 58.5, 58.5, 66.75, 58.5, 66.75, 33.5, 33.5, 33.5, 33.5, 33.5, 33.5, 33.5, 58.5, 50, 33.5, 33.5, 33.5, 33.5, 33.5, 33.5, 58.5, 58.5, 33.5, 33.5, 33.5, 33.5, 33.5, 33.5, 16.5, 58.5, 33.5, 58.5, 33.5, 16.5, 33.5, 66.75, 33.5, 33.5, 58.5, 33.5, 33.5, 58.5, 33.5, 33.5, 33.5, 58.5, 58.5, 33.5, 33.5, 33.5, 33.5, 33.5, 33.5, 33.5, 33.5, 58.5, 33.5, 33.5, 33.5, 58.5, 16.5, 58.5, 33.5, 33.5, 58.5, 33.5, 33.5, 33.5, 33.5, 16.5, 33.5, 33.5, 33.5, 66.75, 33.5, 16.5, 58.5, 33.5, 66.75, 33.5, 58.5, 33.5, 33.5, 33.5, 33.5, 33.5, 58.5, 33.5, 33.5, 33.5, 58.5, 33.5, 33.5, 33.5, 33.5, 33.5, 33.5, 33.5, 66.75, 33.5, 33.5, 33.5, 58.5, 33.5, 58.5, 66.75, 33.5, 33.5, 33.5, 58.5, 33.5, 33.5, 33.5, 33.5, 0, 33.5, 58.5, 33.5, 33.5, 58.5, 50, 33.5, 33.5, 0, 33.5, 58.5, 33.5, 33.5, 33.5, 33.5, 66.75, 16.5, 33.5, 33.5, 33.5, 33.5, 33.5, 66.75),
                       urinarybotherpreop = c(53.125, 81.25, 71.875, 71.875, 71.875, 71.875, 81.25, 53.125, 71.875, 71.875, 71.875, 81.25, 71.875, 81.25, 71.875, 81.25, 53.125, 81.25, 81.25, 71.875, 71.875, 53.125, 53.125, 71.875, 81.25, 81.25, 81.25, 81.25, 71.875, 71.875, 71.875, 53.125, 71.875, 81.25, 71.875, 53.125, 65.625, 71.875, 65.625, 53.125, 71.875, 71.875, 53.125, 71.875, 81.25, 71.875, 53.125, 71.875, 53.125, 53.125, 71.875, 71.875, 71.875, 71.875, 53.125, 81.25, 81.25, 81.25, 81.25, 81.25, 71.875, 71.875, 81.25, 81.25, 81.25, 81.25, 71.875, 65.625, 81.25, 65.625, 81.25, 81.25, 53.125, 81.25, 81.25, 71.875, 71.875, 81.25, 81.25, 81.25, 81.25, 81.25, 81.25, 81.25, 71.875, 81.25, 71.875, 71.875, 81.25, 81.25, 81.25, 71.875, 71.875, 71.875, 81.25, 71.875, 81.25, 81.25, 65.625, 65.625, 81.25, 81.25, 81.25, 71.875, 71.875, 81.25, 65.625, 81.25, 81.25, 81.25, 71.875, 81.25, 71.875, 53.125, 71.875, 81.25, 81.25, 53.125, 81.25, 71.875, 81.25, 81.25, 71.875, 81.25, 81.25, 81.25, 81.25, 71.875, 65.625, 71.875, 81.25, 65.625, 81.25, 71.875, 53.125, 71.875, 71.875, 71.875, 71.875, 65.625, 81.25, 81.25, 81.25, 81.25, 81.25, 65.625, 81.25, 81.25, 65.625, 71.875, 65.625, 81.25, 65.625, 71.875, 81.25, 71.875, 81.25, 71.875, 81.25, 81.25, 81.25, 81.25, 71.875, 71.875, 81.25, 65.625, 81.25, 71.875, 71.875, 65.625, 71.875, 81.25, 53.125, 71.875, 81.25, 53.125, 53.125, 81.25, 65.625, 81.25, 81.25, 71.875, 81.25, 65.625, 71.875, 81.25, 81.25, 65.625, 65.625, 71.875, 81.25, 81.25, 71.875, 81.25, 81.25, 71.875, 81.25, 71.875, 81.25, 81.25, 81.25, 81.25, 71.875, 71.875, 81.25, 81.25, 81.25, 71.875, 71.875, 71.875, 71.875, 81.25, 65.625, 81.25, 81.25, 65.625, 81.25, 71.875, 81.25, 81.25, 65.625, 71.875, 71.875, 71.875, 81.25, 71.875, 71.875, 53.125, 81.25, 81.25, 71.875, 53.125, 81.25, 71.875, 71.875, 71.875, 81.25, 71.875),
                       urinarybother6wks = c(31.25, 31.25, 31.25, 25, 31.25, 31.25, 25, 31.25, 31.25, 31.25, 31.25, 31.25, 31.25, 31.25, 25, 31.25, 31.25, 31.25, 25, 31.25, 31.25, 25, 25, 31.25, 25, 31.25, 25, 31.25, 31.25, 31.25, 31.25, 31.25, 31.25, 31.25, 31.25, 31.25, 31.25, 31.25, 25, 31.25, 31.25, 25, 25, 31.25, 31.25, 31.25, 25, 31.25, 31.25, 31.25, 31.25, 31.25, 25, 31.25, 25, 31.25, 25, 31.25, 31.25, 31.25, 31.25, 31.25, 31.25, 31.25, 31.25, 31.25, 31.25, 25, 31.25, 31.25, 31.25, 31.25, 25, 31.25, 31.25, 31.25, 31.25, 31.25, 25, 31.25, 31.25, 31.25, 31.25, 31.25, 31.25, 31.25, 31.25, 31.25, 31.25, 31.25, 31.25, 31.25, 31.25, 31.25, 25, 25, 31.25, 31.25, 25, 25, 31.25, 31.25, 25, 31.25, 31.25, 31.25, 31.25, 31.25, 31.25, 31.25, 31.25, 31.25, 31.25, 25, 31.25, 31.25, 31.25, 25, 31.25, 31.25, 31.25, 25, 25, 25, 31.25, 31.25, 31.25, 31.25, 25, 31.25, 31.25, 31.25, 31.25, 25, 25, 31.25, 25, 31.25, 31.25, 31.25, 31.25, 31.25, 31.25, 31.25, 31.25, 31.25, 31.25, 31.25, 31.25, 31.25, 25, 31.25, 25, 31.25, 31.25, 31.25, 31.25, 25, 31.25, 31.25, 31.25, 31.25, 31.25, 31.25, 31.25, 25, 31.25, 31.25, 31.25, 25, 25, 31.25, 25, 31.25, 25, 25, 31.25, 25, 31.25, 31.25, 31.25, 31.25, 25, 25, 31.25, 31.25, 31.25, 25, 25, 31.25, 31.25, 31.25, 31.25, 31.25, 31.25, 31.25, 31.25, 31.25, 31.25, 31.25, 31.25, 31.25, 31.25, 31.25, 31.25, 31.25, 31.25, 31.25, 31.25, 31.25, 31.25, 31.25, 31.25, 31.25, 25, 31.25, 25, 31.25, 31.25, 31.25, 25, 31.25, 31.25, 25, 31.25, 31.25, 25, 25, 31.25, 31.25, 31.25, 25, 31.25, 31.25, 31.25, 31.25, 31.25, 31.25), urinarybother3months = c(37.5, 53.125, 43.75, 31.25, 43.75, 31.25, 37.5, 37.5, 43.75, 53.125, 37.5, 53.125, 53.125, 43.75, 31.25, 37.5, 37.5, 37.5, 43.75, 53.125, 37.5, 43.75, 37.5, 31.25, 37.5, 37.5, 31.25, 53.125, 43.75, 37.5, 43.75, 37.5, 53.125, 53.125, 37.5, 37.5, 43.75, 53.125, 25, 37.5, 53.125, 37.5, 37.5, 53.125, 43.75, 37.5, 37.5, 43.75, 37.5, 37.5, 37.5, 37.5, 31.25, 37.5, 31.25, 37.5, 31.25, 37.5, 43.75, 53.125, 37.5, 43.75, 37.5, 37.5, 37.5, 37.5, 53.125, 31.25, 37.5, 37.5, 37.5, 53.125, 43.75, 37.5, 43.75, 37.5, 43.75, 53.125, 31.25, 37.5, 53.125, 37.5, 37.5, 37.5, 43.75, 43.75, 37.5, 37.5, 53.125, 37.5, 53.125, 37.5, 37.5, 37.5, 31.25, 53.125, 43.75, 37.5, 31.25, 37.5, 43.75, 43.75, 53.125, 43.75, 53.125, 37.5, 37.5, 37.5, 37.5, 37.5, 53.125, 37.5, 43.75, 37.5, 37.5, 37.5, 37.5, 31.25, 37.5, 37.5, 43.75, 37.5, 37.5, 37.5, 53.125, 37.5, 37.5, 37.5, 31.25, 43.75, 37.5, 43.75, 37.5, 25, 31.25, 53.125, 31.25, 31.25, 43.75, 37.5, 37.5, 43.75, 37.5, 37.5, 37.5, 43.75, 43.75, 37.5, 37.5, 37.5, 37.5, 53.125, 53.125, 53.125, 37.5, 43.75, 37.5, 31.25, 37.5, 43.75, 37.5, 43.75, 37.5, 37.5, 43.75, 53.125, 37.5, 37.5, 53.125, 31.25, 37.5, 37.5, 31.25, 53.125, 31.25, 37.5, 43.75, 31.25, 53.125, 37.5, 43.75, 37.5, 31.25, 53.125, 37.5, 53.125, 43.75, 37.5, 53.125, 37.5, 43.75, 37.5, 37.5, 53.125, 37.5, 37.5, 37.5, 53.125, 53.125, 37.5, 53.125, 53.125, 43.75, 37.5, 43.75, 53.125, 37.5, 37.5, 53.125, 43.75, 31.25, 37.5, 37.5, 31.25, 25, 37.5, 43.75, 53.125, 37.5, 43.75, 37.5, 53.125, 37.5, 25, 37.5, 43.75, 31.25, 37.5, 31.25, 37.5, 53.125, 37.5, 37.5, 37.5, 37.5, 37.5, 37.5, 53.125)

重塑

我們需要重塑數據並清理變量名稱。 列名隱藏2到3個不同的變量。 這不是ggplot2喜歡事物的方式（這可能就是您遇到麻煩的原因）。

m_data <- reshape2::melt(mydata, id.vars=c("serialnumber")) %>%
  separate(variable, into = c("urinary", "timepoint"), "bother|summary|function", remove=FALSE) %>%
  mutate(variable = str_extract(variable, "bother|summary|function"))%>%
  dplyr::select(-urinary)

充分了解數據的外觀

 serialnumber variable timepoint value
1            1  summary     preop 68.75
2            2  summary     preop 87.50
3            3  summary     preop 81.25
4            4  summary     preop 81.25
5            5  summary     preop 81.25
6            6  summary     preop 81.25

創建有序的x軸

現在讓我們按順序整理訪問（不確定是否是這種方式，但是我沒有更多信息了……）。

警告：這樣做是要混合“ summary”，“ bother”和“ function”值

m_data <- m_data %>%
  group_by(serialnumber) %>%
  mutate(visit_num = 1:length(serialnumber))

情節

現在您可以進行繪圖了

m_data  %>%
  ggplot(aes(factor(visit_num), value, group=serialnumber)) +
  geom_boxplot(aes(factor(visit_num), value, group=visit_num))+
  geom_line(color="red", alpha=0.01)+
  theme_bw()+
  NULL

我很驚訝時間點4的可變性為0 ...您的數據如何處理？

另外，每個患者只有9個時間點，難道我們不應該有23個時間點嗎？