我希望該圖的列標題和行標題分別為phi = 0.8, phi = 0.9, phi = 0.95和sd = 1, sd = 3, sd = 5, sd = 10。本phi應顯示為希臘字母符號,而sd仍然是英文字母。
## simulate ARIMA(1, 0, 0)
set.seed(289805)
x1 <- arima.sim(n = 10, model = list(ar = 0.8, order = c(1, 0, 0)), sd = 1)
set.seed(671086)
x2 <- arima.sim(n = 10, model = list(ar = 0.9, order = c(1, 0, 0)), sd = 1)
set.seed(799837)
x3 <- arima.sim(n = 10, model = list(ar = 0.95, order = c(1, 0, 0)), sd = 1)
set.seed(289805)
x4 <- arima.sim(n = 10, model = list(ar = 0.8, order = c(1, 0, 0)), sd = 3)
set.seed(671086)
x5 <- arima.sim(n = 10, model = list(ar = 0.9, order = c(1, 0, 0)), sd = 3)
set.seed(799837)
x6 <- arima.sim(n = 10, model = list(ar = 0.95, order = c(1, 0, 0)), sd = 3)
set.seed(289805)
x7 <- arima.sim(n = 10, model = list(ar = 0.8, order = c(1, 0, 0)), sd = 5)
set.seed(671086)
x8 <- arima.sim(n = 10, model = list(ar = 0.9, order = c(1, 0, 0)), sd = 5)
set.seed(799837)
x9 <- arima.sim(n = 10, model = list(ar = 0.95, order = c(1, 0, 0)), sd = 5)
set.seed(289805)
x10 <- arima.sim(n = 10, model = list(ar = 0.8, order = c(1, 0, 0)), sd = 10)
set.seed(671086)
x11 <- arima.sim(n = 10, model = list(ar = 0.9, order = c(1, 0, 0)), sd = 10)
set.seed(799837)
x12 <- arima.sim(n = 10, model = list(ar = 0.95, order = c(1, 0, 0)), sd = 10)
xx <- 1:10
df <- data.frame(xx, x1, x2, x3, x4, x5, x6, x7, x8, x9, x10, x11, x12)
reshapp <- reshape2::melt(df, id = "xx")
NEWDAT <- data.frame(y = reshapp$value, x = reshapp$xx, sd = rep(rep(c(sd=1, sd=3, sd=5, sd=10), each = 10), each = 3),phi = rep(rep(c(.8, .9, .95), each = 10), 4))
ggplot(NEWDAT, aes(x = x, y = y)) geom_line() geom_point() labs(x = 'lb', y = 'RMSE') facet_grid(sd ~ phi, scales = "free_y")
theme_bw() ggplot2::scale_y_continuous(expand = c(0.0, 0.00))

uj5u.com熱心網友回復:
[編輯以說明第三行代碼中的重構]
一個快速的解決方法是使用以 UTF8 編碼表示的 Phi 創建或重命名分面鍵,加上創建您想要的字串(即“Phi|SD = N”)在這種情況下,我創建一個新變數:
NEWDAT %>%
mutate(phi_label = paste0("\U03D5 = ", phi), #"\U03D5" represents the character "?"
sd_label = fct_reorder(.f = paste0("sd =", sd), .x = sd)) %>%
ggplot(aes(x = x, y = y)) geom_line() geom_point() labs(x = 'lb', y = 'RMSE')
facet_grid(sd_label ~ phi_label, scales = "free_y")
theme_bw() ggplot2::scale_y_continuous(expand = c(0.0, 0.00))

uj5u.com熱心網友回復:
這是一個plotmath僅使用標準字符,沒有特殊轉義序列的解決方案。
NEWDAT <- data.frame(y = reshapp$value, x = reshapp$xx, sd = rep(rep(c(sd=1, sd=3, sd=5, sd=10), each = 10), each = 3),phi = rep(rep(c(.8, .9, .95), each = 10), 4))
NEWDAT$sd <- factor(NEWDAT$sd, levels = NEWDAT$sd, labels = paste("sd ==", NEWDAT$sd))
NEWDAT$phi <- with(NEWDAT, factor(phi, levels = phi, labels = paste("phi ==", phi)))
ggplot(NEWDAT, aes(x = x, y = y))
geom_line()
geom_point()
scale_y_continuous(expand = c(0.0, 0.00))
labs(x = 'lb', y = 'RMSE')
facet_grid(
sd ~ phi,
scales = "free_y",
labeller = label_parsed
)
theme_bw()

uj5u.com熱心網友回復:
部分答案是這樣的:使用此facet_grid呼叫而不是您現在正在使用的呼叫。
facet_grid(sd ~ phi, scales = "free_y", labeller = . %>% label_both(sep = " = "))
使用 bquote
facet_grid(sd ~ phi, scales = "free_y",
labeller =
label_bquote(rows = sigma == .(sd), cols = phi == .(phi))
)
這將產生正確的標簽。
靈感來自這篇文章的答案
轉載請註明出處,本文鏈接:https://www.uj5u.com/yidong/352299.html
