- html - 出于某种原因,IE8 对我的 Sass 文件中继承的 html5 CSS 不友好?
- JMeter 在响应断言中使用 span 标签的问题
- html - 在 :hover and :active? 上具有不同效果的 CSS 动画
- html - 相对于居中的 html 内容固定的 CSS 重复背景?
这是我的数据示例
library(tidyr)
library(dplyr)
library(ggplot2)
resource <- c("good","good","bad","bad","good","good","bad","bad","good","good","bad","bad","good","good","bad","bad")
fertilizer <- c("none", "nitrogen","none","nitrogen","none", "nitrogen","none","nitrogen","none", "nitrogen","none","nitrogen","none", "nitrogen","none","nitrogen")
t0 <- sample(1:20, 16)
t1 <- sample(1:20, 16)
t2 <- sample(1:20, 16)
t3 <- sample(1:20, 16)
t4 <- sample(1:20, 16)
t5 <- sample(1:20, 16)
t6 <- sample(10:100, 16)
t7 <- sample(10:100, 16)
t8 <- sample(10:100, 16)
t9 <- sample(10:100, 16)
t10 <- sample(10:100, 16)
replicates <- c(1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16)
data <- data.frame(resource, fertilizer,replicates, t0,t1,t2,t3,t4,t5,t6,t7,t8,t9,t10)
data$resource <- as.factor(data$resource)
data$fertilizer <- as.factor(data$fertilizer)
data.melt <- data %>% ungroup %>% gather(time, value, -replicates, -resource, -fertilizer)
data.melt$predict <- sample(1:200, 176)
其中,资源和肥料有 2 个因素,因此有效有 4 个处理和 4 x 4 = 16 个重复。时间是一个有 10 个级别的因子。我运行了一个模型,并预测了预测
列中的值。
现在我想绘制一个时间序列,其中 x 轴为时间,y 轴为拟合值(预测)的平均值,原始值(值)为每种类型的资源和肥料( 4 个处理)[即 4 个地 block ]。我还想为每个时间点的藻类生长添加一个置信区间。这是我对代码的尝试。
ggplot(df, aes(x=time, y=predicted)) + geom_point(size=3)+ stat_summary(geom = "point", fun.y = "mean") + facet_grid(resource + fertilizer ~.)
通过这个简单的代码,我仍然只得到 2 个图表,而不是 4 个。此外,预测函数的均值没有绘制出来。我不知道如何将值
和预测
绘制在一起,以及相应的置信区间。
如果有人也可以展示所有四种处理方法如何在一个图上进行,以及如果我可以将其分面(如上所示),那将会很有帮助
最佳答案
我建议的解决方案是创建第二个 data.frame ,其中包含所有摘要统计信息,例如平均预测值。我展示了一种使用 dplyr
包中的 group_by
和 summarize
来实现此目的的方法。摘要数据需要具有与主要数据匹配的resource
、肥料
和time
列。摘要数据还包含带有附加 y
值的列。
然后,需要将主要数据和摘要数据分别提供给适当的 ggplot 函数,但不要在主 ggplot() 调用中提供。 facet_grid
可用于将数据分成四个图。
# Convert time to factor, specifying correct order of time points.
data.melt$time = factor(data.melt$time, levels=paste("t", seq(0, 10), sep=""))
# Create an auxilliary data.frame containing summary data.
# I've used standard deviation as place-holder for confidence intervals;
# I'll let you calculate those on your own.
summary_dat = data.melt %>%
group_by(resource, fertilizer, time) %>%
summarise(mean_predicted=mean(predict),
upper_ci=mean(predict) + sd(predict),
lower_ci=mean(predict) - sd(predict))
p = ggplot() +
theme_bw() +
geom_errorbar(data=summary_dat, aes(x=time, ymax=upper_ci, ymin=lower_ci),
width=0.3, size=0.7, colour="tomato") +
geom_point(data=data.melt, aes(x=time, y=value),
size=1.6, colour="grey20", alpha=0.5) +
geom_point(data=summary_dat, aes(x=time, y=mean_predicted),
size=3, shape=21, fill="tomato", colour="grey20") +
facet_grid(resource ~ fertilizer)
ggsave("plot.png", plot=p, height=4, width=6.5, units="in", dpi=150)
关于r - 绘制 2x2x2 时间序列的原始值和预测值,我们在Stack Overflow上找到一个类似的问题: https://stackoverflow.com/questions/47745478/
我的问题:非常具体。我正在尝试想出解析以下文本的最简单方法: ^^domain=domain_value^^version=version_value^^account_type=account_ty
好吧,这就是我的困境: 我正在为 Reddit 子版 block 开发常见问题解答机器人。我在 bool 逻辑方面遇到了麻烦,需要一双更有经验的眼睛(这是我在 Python 中的第一次冒险)。现在,该
它首先遍历所有 y 值,然后遍历所有 x 值。我需要 X 和 y 同时改变。 For x = 3 To lr + 1 For y = 2 To lr anyl.Cells(x, 1)
假设我有一个包含 2 列的 Excel 表格:单元格 A1 到 A10 中的日期和 B1 到 B10 中的值。 我想对五月日期的所有值求和。我有3种可能性: {=SUM((MONTH(A1:A10)=
如何转换 Z-score来自 Z-distribution (standard normal distribution, Gaussian distribution)到 p-value ?我还没有找到
我正在重写一些 Javascript 代码以在 Excel VBA 中工作。由于在这个网站上搜索,我已经设法翻译了几乎所有的 Javascript 代码!但是,有些代码我无法准确理解它在做什么。这是一
我遇到过包含日期格式的时间戳日期的情况。然后我想构建一个图表,显示“点击”项目的数量“每天”, //array declaration $array1 = array("Date" => 0); $a
我是scala的新手! 我的问题是,是否有包含成员的案例类 myItem:Option[String] 当我构造类时,我需要将字符串内容包装在: Option("some string") 要么 So
我正在用 PHP 创建一个登录系统。我需要用户使用他或她的用户名或电子邮件或电话号码登录然后使用密码。因为我知道在 Java 中我们会像 email==user^ username == user 这
我在 C++ 项目上使用 sqlite,但是当我在具有文本值的列上使用 WHERE 时出现问题 我创建了一个 sqlite 数据库: CREATE TABLE User( id INTEGER
当构造函数是显式时,它不用于隐式转换。在给定的代码片段中,构造函数被标记为 explicit。那为什么在 foo obj1(10.25); 情况下它可以工作,而在 foo obj2=10.25; 情况
我知道这是一个主观问题,所以如果需要关闭它,我深表歉意,但我觉得它经常出现,让我想知道是否普遍偏爱一种形式而不是另一种形式。 显然,最好的答案是“重构代码,这样你就不需要测试是否存在错误”,但有时没有
这两个 jQuery 选择器有什么区别? 以下是来自 w3schools.com 的定义: [attribute~=value] 选择器选择带有特定属性,其值包含特定字符串。 [attribute*=
为什么我们需要CSS [attribute|=value] Selector根本当 CSS3 [attribute*=value] Selector基本上完成相同的事情,浏览器兼容性几乎相似?是否存在
我正在解决 regx 问题。我已经有一个像这样的 regx [0-9]*([.][0-9]{2})。这是 amont 格式验证。现在,通过此验证,我想包括不应提供 0 金额。比如 10 是有效的,但
我正在研究计算机科学 A 考试的样题,但无法弄清楚为什么以下问题的正确答案是正确的。 考虑以下方法。 public static void mystery(List nums) { for (
好的,我正在编写一个 Perl 程序,它有一个我收集的值的哈希值(完全在一个完全独立的程序中)并提供给这个 Perl 脚本。这个散列是 (string,string) 的散列。 我想通过 3 种方式对
我有一个表数据如下,来自不同的表。仅当第三列具有值“债务”并且第一列(日期)具有最大值时,我才想从第四列中获取最大值。最终值基于 MAX(DATE) 而不是 MAX(PRICE)。所以用简单的语言来说
我有一个奇怪的情况,只有错误状态保存到数据库中。当“状态”应该为 true 时,我的查询仍然执行 false。 我有具有此功能的 Controller public function change_a
我有一个交易表(针对所需列进行了简化): id client_id value 1 1 200 2 2 150 3 1
我是一名优秀的程序员,十分优秀!