小编Geo*_*rge的帖子

您好我使用数据框中的ggplot2 geom_histogram函数绘制了一个图,请参阅下面的示例并链接到ggplot直方图需要使用嵌套的ddply函数和facet wrap来标记每个geom_vline的因子

现在我需要创建一个包含上图中使用的汇总数据的数据框.有人可以帮忙吗？提前致谢.G

Sector2 Family  Year    Length
BUN Acroporidae 2010    332.1300496
BUN Poritidae   2011    141.1467966
BUN Acroporidae 2012    127.479
BUN Acroporidae 2013    142.5940556
MUR Faviidae    2010    304.0405
MUR Faviidae    2011    423.152
MUR Pocilloporidae  2012    576.0295
MUR Poritidae   2013    123.8936667
NTH Faviidae    2010    60.494
NTH Faviidae    2011    27.427
NTH Pocilloporidae  2012    270.475
NTH Poritidae   2013    363.4635

可以帮我改变这些图的标题文字大小.即使它们变大？

脚本

ggplot(NMPSCMOR, aes(Length,  fill=Year)) + 
  geom_histogram(position="dodge", binwidth=60,  colour="black") + xlim(0, 600) +
  scale_fill_grey(start = 1, end = 0)+ 


  geom_vline(data=ddply(NMPSCMOR, Year~Morphology~Sector2, numcolwise(mean)), 
             mapping=aes(xintercept=Length,color=Year),  linetype=2, size=1) + 
  scale_color_grey(start=1,end=0)+

  xlab("Length Class") +
  ylab(expression(paste("Total Count"))) + #( ", m^2, ")", sep = 
  facet_wrap( ~ Morphology + Sector2, ncol=3, scales = "free") +
  theme(

    panel.grid.minor = element_blank(),               #removes minor grid lines
  panel.grid.major = element_blank())

我在这里发布了一个问题,并且能够重现 克劳斯的答案 ,使用虹膜数据上的tidyverse计算加性模型中每个物种的多个r平方值.但是,包发生了更新,现在没有计算R-sq值.不确定为什么......
这是子句响应和输出

library(tidyverse)
library(broom)
iris %>% nest(-Species) %>% 
  mutate(fit = map(data, ~mgcv::gam(Sepal.Width ~ s(Sepal.Length, bs = "cs"), data = .)),
         results = map(fit, glance),
         R.square = map(fit, ~ summary(.)$r.sq)) %>%
  unnest(results) %>%
  select(-data, -fit)

#      Species  R.square       df    logLik      AIC      BIC deviance df.residual
# 1     setosa 0.5363514 2.546009 -1.922197 10.93641 17.71646 3.161460    47.45399
# 2 versicolor 0.2680611 2.563623 -3.879391 14.88603 21.69976 3.418909    47.43638
# 3  virginica 0.1910916 2.278569 -7.895997 22.34913 28.61783 4.014793    47.72143

然而,我的代码和输出会产生这些 …