snippet of data.frame snippet of list of data.frames
我的列表如下所示:
> dput(lapply(list1, head))
list(structure(c(NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA,
NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA), .Dim = c(6L,
4L), .Dimnames = list(NULL, c("FPAR", "gpp", "LAI", "npp"))),
structure(c(63L, 83L, 66L, 84L, 92L, 85L, 9976L, 3318L, 9456L,
9435L, 9002L, 9395L, 21L, 32L, 18L, 34L, 50L, 36L, 6742L,
5228L, 5405L, 5136L, 5387L, 5339L), .Dim = c(6L, 4L), .Dimnames = list(
NULL, c("FPAR", "gpp", "LAI", "npp"))), structure(c(NA,
NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA,
NA, NA, NA, NA, NA, NA, NA, NA), .Dim = c(6L, 4L), .Dimnames = list(
NULL, c("FPAR", "gpp", "LAI", "npp"))), structure(c(0L,
95L, 0L, 82L, 0L, 82L, 10306L, 10205L, 10306L, 10627L, 10306L,
10627L, 0L, 64L, 0L, 31L, 0L, 31L, 6396L, 6340L, 6396L, 6396L,
6396L, 6396L), .Dim = c(6L, 4L), .Dimnames = list(NULL, c("FPAR",
"gpp", "LAI", "npp"))), structure(c(NA, NA, NA, NA, NA, NA,
NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA,
7331L, NA, NA), .Dim = c(6L, 4L), .Dimnames = list(NULL,
c("FPAR", "gpp", "LAI", "npp"))), structure(c(96L, 96L,
100L, 88L, 93L, 100L, 12734L, 12215L, 11383L, 11446L, 12672L,
11531L, 61L, 61L, 70L, 40L, 58L, 70L, 7807L, 7357L, 7695L,
6400L, 6009L, 7735L), .Dim = c(6L, 4L), .Dimnames = list(
NULL, c("FPAR", "gpp", "LAI", "npp"))), structure(c(95L,
95L, 96L, 96L, 96L, 95L, 10829L, 10829L, 9652L, 10321L, 9652L,
10829L, 61L, 61L, 62L, 65L, 62L, 61L, 5144L, 5144L, 5188L,
5422L, 5188L, 5144L), .Dim = c(6L, 4L), .Dimnames = list(
NULL, c("FPAR", "gpp", "LAI", "npp"))), structure(c(NA,
NA, 87L, 63L, 87L, 87L, NA, NA, 7891L, 7891L, 0L, 7891L,
NA, NA, 46L, 19L, 46L, 46L, NA, NA, NA, NA, NA, NA), .Dim = c(6L,
4L), .Dimnames = list(NULL, c("FPAR", "gpp", "LAI", "npp"
))), structure(c(100L, 100L, 100L, 78L, 100L, 100L, 4011L,
4011L, 4112L, 4306L, 3664L, 4112L, 70L, 70L, 70L, 21L, 70L,
70L, 2253L, 2425L, 2479L, NA, 2253L, 2253L), .Dim = c(6L,
4L), .Dimnames = list(NULL, c("FPAR", "gpp", "LAI", "npp"
))), structure(c(83L, 87L, 71L, NA, 0L, 82L, 5627L, 6626L,
5862L, NA, 6963L, 5541L, 29L, 37L, 16L, NA, 0L, 28L, 3877L,
3695L, 3635L, 1642L, 3692L, 3644L), .Dim = c(6L, 4L), .Dimnames = list(
NULL, c("FPAR", "gpp", "LAI", "npp"))), structure(c(88L,
89L, 61L, NA, 0L, 88L, 6132L, 5636L, 6166L, NA, 6510L, 6502L,
49L, 51L, 18L, NA, 0L, 50L, 3413L, 3788L, 3519L, NA, 3463L,
3755L), .Dim = c(6L, 4L), .Dimnames = list(NULL, c("FPAR",
"gpp", "LAI", "npp"))))
因此，我尝试为列表中的 11 个数据帧有效地绘制变量 1 (FPAR) 和变量 2-4 (gpp、LAI、npp) 之间的散点图。我附上了我的列表的片段。
这是我到目前为止想出的办法，但它不起作用。
library(ggplot2)
scatterplot <- function(yvar){
ggplot(list1, aes_(x=~age, y=as.name(yvar))) +
geom_point() +
facet_wrap(~dataframe)
}
lapply(list1, scatterplot)
这是我第一次使用列表和 ggplot，所以我有点卡住了..
最佳答案
这是一种可以做到这一点的方法。
library(tidyverse)
list1 %>%
map(., as_data_frame) %>%
bind_rows(., .id = "dataframe") %>%
gather(., key, value, gpp:npp) %>%
mutate(dataframe = factor(dataframe, levels = unique(dataframe))) %>%
ggplot(., aes(FPAR, value, col = key)) +
geom_point() +
facet_wrap(. ~ dataframe, ncol = 3)
warning: Removed 71 rows containing missing values (geom_point).
您可能想阅读:
关于r - 绘制数据框列表，每个数据框有 4 个变量，我们在Stack Overflow上找到一个类似的问题： https://stackoverflow.com/questions/53099525/