r - 为什么 Facet_Grid 输出与其数据不同
问题描述
我有这些时间序列数据,我绘制如下
## 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
# ggplot for x1
plot1 <- ggplot2::ggplot(NULL, aes(y = x1, x = xx)) + ggplot2::geom_line(color = "#F2AA4CFF") + ggplot2::geom_point(color = "#101820FF") + xlab('lb') + ylab('RMSE') + ggplot2::theme_bw()+ ggplot2::scale_y_continuous(expand = c(0.0, 0.00))
# ggplot for x2
plot2 <- ggplot2::ggplot(NULL, aes(y = x2, x = xx)) + ggplot2::geom_line(color = "#F2AA4CFF") + ggplot2::geom_point(color = "#101820FF") + xlab('lb') + ylab('RMSE') + ggplot2::theme_bw()+ ggplot2::scale_y_continuous(expand = c(0.0, 0.00))
# ggplot for x3
plot3 <- ggplot2::ggplot(NULL, aes(y = x3, x = xx)) + ggplot2::geom_line(color = "#F2AA4CFF") + ggplot2::geom_point(color = "#101820FF") + xlab('lb') + ylab('RMSE') + ggplot2::theme_bw()+ ggplot2::scale_y_continuous(expand = c(0.0, 0.00))
# ggplot for x4
plot4 <- ggplot2::ggplot(NULL, aes(y = x4, x = xx)) + ggplot2::geom_line(color = "#F2AA4CFF") + ggplot2::geom_point(color = "#101820FF") + xlab('lb') + ylab('RMSE') + ggplot2::theme_bw()+ ggplot2::scale_y_continuous(expand = c(0.0, 0.00))
# ggplot for x5
plot5 <- ggplot2::ggplot(NULL, aes(y = x5, x = xx)) + ggplot2::geom_line(color = "#F2AA4CFF") + ggplot2::geom_point(color = "#101820FF") + xlab('lb') + ylab('RMSE') + ggplot2::theme_bw()+ ggplot2::scale_y_continuous(expand = c(0.0, 0.00))
# ggplot for x6
plot6 <- ggplot2::ggplot(NULL, aes(y = x6, x = xx)) + ggplot2::geom_line(color = "#F2AA4CFF") + ggplot2::geom_point(color = "#101820FF") + xlab('lb') + ylab('RMSE') + ggplot2::theme_bw()+ ggplot2::scale_y_continuous(expand = c(0.0, 0.00))
# ggplot for x7
plot7 <- ggplot2::ggplot(NULL, aes(y = x7, x = xx)) + ggplot2::geom_line(color = "#F2AA4CFF") + ggplot2::geom_point(color = "#101820FF") + xlab('lb') + ylab('RMSE') + ggplot2::theme_bw()+ ggplot2::scale_y_continuous(expand = c(0.0, 0.00))
# ggplot for x8
plot8 <- ggplot2::ggplot(NULL, aes(y = x8, x = xx)) + ggplot2::geom_line(color = "#F2AA4CFF") + ggplot2::geom_point(color = "#101820FF") + xlab('lb') + ylab('RMSE') + ggplot2::theme_bw()+ ggplot2::scale_y_continuous(expand = c(0.0, 0.00))
# ggplot for x9
plot9 <- ggplot2::ggplot(NULL, aes(y = x9, x = xx)) + ggplot2::geom_line(color = "#F2AA4CFF") + ggplot2::geom_point(color = "#101820FF") + xlab('lb') + ylab('RMSE') + ggplot2::theme_bw()+ ggplot2::scale_y_continuous(expand = c(0.0, 0.00))
# ggplot for x10
plot10 <- ggplot2::ggplot(NULL, aes(y = x10, x = xx)) + ggplot2::geom_line(color = "#F2AA4CFF") + ggplot2::geom_point(color = "#101820FF") + xlab('lb') + ylab('RMSE') + ggplot2::theme_bw()+ ggplot2::scale_y_continuous(expand = c(0.0, 0.00))
# ggplot for x11
plot11 <- ggplot2::ggplot(NULL, aes(y = x11, x = xx)) + ggplot2::geom_line(color = "#F2AA4CFF") + ggplot2::geom_point(color = "#101820FF") + xlab('lb') + ylab('RMSE') + ggplot2::theme_bw()+ ggplot2::scale_y_continuous(expand = c(0.0, 0.00))
# ggplot for x12
plot12 <- ggplot2::ggplot(NULL, aes(y = x12, x = xx)) + ggplot2::geom_line(color = "#F2AA4CFF") + ggplot2::geom_point(color = "#101820FF") + xlab('lb') + ylab('RMSE') + ggplot2::theme_bw()+ ggplot2::scale_y_continuous(expand = c(0.0, 0.00))
# plot in a 3 by 5 grid by using plot_layout
plot1 + plot2 + plot3 + plot4 + plot5 + plot6 + plot7 + plot8 + plot9 + plot10 + plot11 + plot12 + patchwork::plot_layout(ncol = 3, byrow = TRUE)
我现在想给它一个整容facet_grid
,我尝试如下
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(c(1,3,5,10), each = 3), phi = rep(c(.8, .9, .95), 4))
ggplot(NEWDAT, aes(x = x, y = y)) + geom_line() + geom_point() + facet_grid(sd ~ phi, scales = "free_y") +
theme_bw() + ggplot2::scale_y_continuous(expand = c(0.0, 0.00))
我想要的是
我如何使之Facet_Grid Plot
成为现实Grid_Arranged Plot
?
指针
可以看到lowest
第一个系列中的数据点是-0.2505125
while the highest
is 4.6956397
。在和Facet-Grid Plot
中。lowest is -20+
highest is 30+
解决方案
当我阅读评论时,我可以看到您将删除此问题。但是,我想请您花点时间看看我的解决方案。我认为它更简单,更透明。
我们将从一个简单tibble
的计算参数开始。
library(tidyverse)
df = tibble(
seed = rep(c(289805, 671086, 799837), 4),
sd = rep(c(1, 3, 5, 10), each=3),
ar = rep(c(.8, .9, .95), 4)
)
输出
# A tibble: 12 x 3
seed sd ar
<dbl> <dbl> <dbl>
1 289805 1 0.8
2 671086 1 0.9
3 799837 1 0.95
4 289805 3 0.8
5 671086 3 0.9
6 799837 3 0.95
7 289805 5 0.8
8 671086 5 0.9
9 799837 5 0.95
10 289805 10 0.8
11 671086 10 0.9
12 799837 10 0.95
现在我们正在准备一个简单的farima
函数。
farima = function(seed, n, ar, order, sd){
set.seed(seed)
model = arima.sim(n=n, model=list(ar = ar, order = order), sd = sd)
tibble(
name = paste("sd -", sd, "ar -", ar),
x = model %>% as.numeric(),
lb = 1:length(x)
)
}
之后你所要做的就是一个简单的突变!
df %>% nest_by(seed, ar, sd) %>%
mutate(data = list(farima(seed, 10, ar, c(1, 0, 0), sd))) %>%
unnest(data) %>%
ggplot(aes(lb, x))+
geom_line(color = "#F2AA4CFF")+
geom_point(color = "#101820FF")+
theme_bw()+
scale_y_continuous(expand = c(0, 0))+
facet_wrap(.~name, ncol = 3, scales = "free")
但你也可以做一些不同的事情。
df %>% nest_by(seed, ar, sd) %>%
mutate(data = list(farima(seed, 10, ar, c(1, 0, 0), sd))) %>%
unnest(data) %>%
ggplot(aes(lb, x))+
geom_line(color = "#F2AA4CFF")+
geom_point(color = "#101820FF")+
theme_bw()+
scale_y_continuous(expand = c(0, 0))+
facet_grid(vars(sd), vars(ar), scales = "free")
我希望你喜欢我的解决方案,而不是删除这个问题,你会将它标记为有用。
@Daniel James 的小更新
您必须很好地了解使用该nest_by
功能时会发生什么。我会一步一步地告诉你。
df %>% nest_by(seed, ar, sd) %>% #Step 1
mutate(data = list(farima(seed, 10, ar, c(1, 0, 0), sd))) %>% #Step 2
unnest(data) %>% #Step 3
ggplot(aes(lb, x))+ #Step 4
geom_line(color = "#F2AA4CFF")+
geom_point(color = "#101820FF")+
theme_bw()+
scale_y_continuous(expand = c(0, 0))+
facet_wrap(.~name, ncol = 3, scales = "free")
在第一步之后,您会得到一个 shattered tibble
,但是,不是按变量,而是按行。变量出现在date
这里,在这种情况下基本上是空的(它包含零值)。
# A tibble: 12 x 4
# Rowwise: seed, ar, sd
seed ar sd data
<dbl> <dbl> <dbl> <list<tibble[,0]>>
1 289805 0.8 1 [1 x 0]
2 289805 0.8 3 [1 x 0]
3 289805 0.8 5 [1 x 0]
4 289805 0.8 10 [1 x 0]
5 671086 0.9 1 [1 x 0]
6 671086 0.9 3 [1 x 0]
7 671086 0.9 5 [1 x 0]
8 671086 0.9 10 [1 x 0]
9 799837 0.95 1 [1 x 0]
10 799837 0.95 3 [1 x 0]
11 799837 0.95 5 [1 x 0]
12 799837 0.95 10 [1 x 0]
如果您现在进行突变,它将应用于每一行中的值。
所以让我们看看我们在第二步之后有什么。
# A tibble: 12 x 4
# Rowwise: seed, ar, sd
seed ar sd data
<dbl> <dbl> <dbl> <list>
1 289805 0.8 1 <tibble [10 x 3]>
2 289805 0.8 3 <tibble [10 x 3]>
3 289805 0.8 5 <tibble [10 x 3]>
4 289805 0.8 10 <tibble [10 x 3]>
5 671086 0.9 1 <tibble [10 x 3]>
6 671086 0.9 3 <tibble [10 x 3]>
7 671086 0.9 5 <tibble [10 x 3]>
8 671086 0.9 10 <tibble [10 x 3]>
9 799837 0.95 1 <tibble [10 x 3]>
10 799837 0.95 3 <tibble [10 x 3]>
11 799837 0.95 5 <tibble [10 x 3]>
12 799837 0.95 10 <tibble [10 x 3]>
请注意,我们的farima
函数在mutate
函数中使用 arguments seed
,ar
并sd
在每一行上接收单个值!还要注意使用list
您将包装我们的函数的farima
函数。
现在我们要做的就是解决我们的data
变量。我们将在第 3 步中执行此操作。
# A tibble: 120 x 6
# Groups: seed, ar, sd [12]
seed ar sd name x lb
<dbl> <dbl> <dbl> <chr> <dbl> <int>
1 289805 0.8 1 sd - 1 ar - 0.8 -0.251 1
2 289805 0.8 1 sd - 1 ar - 0.8 0.886 2
3 289805 0.8 1 sd - 1 ar - 0.8 0.198 3
4 289805 0.8 1 sd - 1 ar - 0.8 1.54 4
5 289805 0.8 1 sd - 1 ar - 0.8 4.70 5
6 289805 0.8 1 sd - 1 ar - 0.8 4.44 6
7 289805 0.8 1 sd - 1 ar - 0.8 3.02 7
8 289805 0.8 1 sd - 1 ar - 0.8 1.81 8
9 289805 0.8 1 sd - 1 ar - 0.8 1.17 9
10 289805 0.8 1 sd - 1 ar - 0.8 1.16 10
# ... with 110 more rows
第 4 步及以后发生的一切只是创建适当的图表。
推荐阅读
- python - Tensorflow 对象检测:使用 .h5 (hdf5) 文件从头开始训练
- javascript - 显示来自 mustache 的 javascript 变量内容
- amazon-web-services - 负载均衡器/目标组关系查询
- python - Python 脚本无权更改文件
- python - Google Colaboratory 上带有 GPU 的 XGBoost
- java - install4j:如何使用在安装期间添加的 jdbc 驱动程序检查 JDBC 连接
- regex - 正则表达式 - 如何在保留原始字符串的同时替换字符串中的特定模式?
- repository - jcaptcha 依赖项不可用
- scala - 在 Gatling 测试中传递配置
- plot - Gnuplot 填充曲线未正确显示