- html - 出于某种原因,IE8 对我的 Sass 文件中继承的 html5 CSS 不友好?
- JMeter 在响应断言中使用 span 标签的问题
- html - 在 :hover and :active? 上具有不同效果的 CSS 动画
- html - 相对于居中的 html 内容固定的 CSS 重复背景?
我有一个函数可以使用移动窗口计算中位数和 90% 置信区间。因此,对于每个 x = seq(xmin, xmax, by = wStep)
,我返回所有 y
的中位数以及 5% 和 95% 分位数,其 x
值小于 wSize/2
。我想通过创建自定义平滑函数 stat_movingwindow()
使用 ggplot2 将其显示为线条和色带。我可以使用 geom_smooth(data = ..., stat = "identity")
创建我想要的结果:
moveWin <- function(d, wSize = 0.5, wStep = 0.1,
f = function(x) quantile(x, prob = c(0.05,0.50,0.95), na.rm = TRUE)
){
x <- seq(min(d$x), max(d$x), by = wStep)
y <- matrix(NA, ncol = 3, nrow = length(x))
for(i in seq_along(x)){
y[i, ] <- f(d[abs(d$x - x[i]) < wSize/2, ]$y)
}
y <- as.tibble(y)
colnames(y) <- c("ymin","y","ymax")
y$x <- x
return(as.tibble(y))
}
set.seed(123)
d <- tibble(
x= sqrt(seq(0,1,length.out = 50))*10,
y= rnorm(50)
)
ggplot(data = d) + aes(x = x, y = y) +
geom_smooth(
data = function(d) moveWin(d, wSize = 1, wStep = 0.1),
mapping = aes(ymin = ymin, ymax= ymax),
stat = "identity") +
geom_point() + scale_x_continuous(breaks = 1:10)
按照小插图 Extending ggplot2 ,这是我到目前为止提出的代码。然而,问题是这并没有显示功能区。也许我需要某种方式来声明此自定义统计信息提供了美观的 ymin 和 ymax。如何让下面的代码输出与上面类似的结果?
StatMovingWindow <- ggproto("StatMovingWindow", Stat,
compute_group = function(data, scales, wSize, wStep, fun) {
moveWin(data, wSize = wSize, wStep = wStep, f = fun)
},
required_aes = c("x", "y")
)
stat_movingwindow <- function(mapping = NULL, data = NULL,
fun = function(d) quantile(d, probs = c(0.05, 0.50, 0.95), na.rm = TRUE),
wStep = 0.1, wSize = 1,
geom = "smooth", position = "identity", show.legend = NA, inherit.aes = TRUE,
...
){
layer(
stat = StatMovingWindow, data = data, mapping = mapping, geom = geom,
position = position, show.legend = show.legend, inherit.aes = inherit.aes,
params = list(wStep = wStep, wSize = wSize, fun = fun, ...)
)
}
ggplot(data = d) + aes(x = x, y = y) +
stat_movingwindow(wStep = 0.1, wSize = 1) +
geom_point() + scale_x_continuous(breaks = 1:10)
最佳答案
在您的 stat_movingwindow
代码中,对应 geom 的行是 geom = "smooth"
:
stat_movingwindow <- function(mapping = NULL, data = NULL,
fun = function(d) quantile(d, probs = c(0.05, 0.50, 0.95), na.rm = TRUE),
wStep = 0.1, wSize = 1,
geom = "smooth", # <- look here
position = "identity", show.legend = NA, inherit.aes = TRUE,
...
){
layer(
stat = StatMovingWindow, data = data, mapping = mapping, geom = geom,
position = position, show.legend = show.legend, inherit.aes = inherit.aes,
params = list(wStep = wStep, wSize = wSize, fun = fun, ...)
)
}
检查 geom_smooth
的代码,我们看到它包含参数 se = TRUE
,并使用 GeomSmooth
作为它的 geom:
> geom_smooth
function (mapping = NULL, data = NULL, stat = "smooth", position = "identity",
..., method = "auto", formula = y ~ x, se = TRUE, # <- look here
na.rm = FALSE,
show.legend = NA, inherit.aes = TRUE)
{
params <- list(na.rm = na.rm, se = se, ...)
if (identical(stat, "smooth")) {
params$method <- method
params$formula <- formula
}
layer(data = data, mapping = mapping, stat = stat, geom = GeomSmooth, # <- and here
position = position, show.legend = show.legend, inherit.aes = inherit.aes,
params = params)
}
深入研究 GeomSmooth,我们看到它的 draw_group
函数(负责绘制平滑线)将 se = FALSE
作为其默认参数。
从代码来看,如果 se == FALSE
,has_ribbon
也将是 FALSE
,即使 ymax
> & ymin
存在于您的数据中,感谢 StatMovingWindow$compute_group
函数。这反过来意味着 GeomLine$draw_panel(path, panel_params, coord)
的唯一结果将单独返回,没有 GeomRibbon$draw_group(ribbon, panel_params, coord)
.
> GeomSmooth$draw_group
<ggproto method>
<Wrapper function>
function (...)
f(...)
<Inner function (f)>
function (data, panel_params, coord, se = FALSE) # <- look here
{
ribbon <- transform(data, colour = NA)
path <- transform(data, alpha = NA)
has_ribbon <- se && !is.null(data$ymax) && !is.null(data$ymin) # <- and here
gList(if (has_ribbon) GeomRibbon$draw_group(ribbon, panel_params, coord),
GeomLine$draw_panel(path, panel_params, coord))
}
简而言之,geom_smooth
的默认参数 se = TRUE
覆盖了 GeomSmooth$draw_group
中的默认行为,(同样适用于stat_smooth
也是如此),如果我们想要获得相同的结果,我们应该在 stat_movingwindow
中执行相同的操作。
如果您希望经常绘制色带,则可以将 se = TRUE
作为参数包含在 stat_movingwindow
的定义中。如果它是临时的,您可以在需要时将其包含在您的代码中。
关于r - 如何创建自定义 ggplot2 平滑统计(不仅仅是自定义 lm 或 glm 模型),我们在Stack Overflow上找到一个类似的问题: https://stackoverflow.com/questions/54825289/
好的,所以我编辑了以下... 只需将以下内容放入我的 custom.css #rt-utility .rt-block {CODE HERE} 但是当我尝试改变... 与 #rt-sideslid
在表格 View 中,我有一个自定义单元格(在界面生成器中高度为 500)。在该单元格中,我有一个 Collection View ,我按 (10,10,10,10) 固定到边缘。但是在 tablev
对于我的无能,我很抱歉,但总的来说,我对 Cocoa、Swift 和面向对象编程还很陌生。我的主要来源是《Cocoa Programming for OS X》(第 5 版),以及 Apple 的充满
我正在使用 meta-tegra 为我的 NVIDIA Jetson Nano 构建自定义图像。我需要 PyTorch,但没有它的配方。我在设备上构建了 PyTorch,并将其打包到设备上的轮子中。现
在 jquery 中使用 $.POST 和 $.GET 时,有没有办法将自定义变量添加到 URL 并发送它们?我尝试了以下方法: $.ajax({type:"POST", url:"file.php?
Traefik 已经默认实现了很多中间件,可以满足大部分我们日常的需求,但是在实际工作中,用户仍然还是有自定义中间件的需求,为解决这个问题,官方推出了一个 Traefik Pilot[1] 的功
我想让我的 CustomTextInputLayout 将 Widget.MaterialComponents.TextInputLayout.OutlinedBox 作为默认样式,无需在 XML 中
我在 ~/.emacs 中有以下自定义函数: (defun xi-rgrep (term) (grep-compute-defaults) (interactive "sSearch Te
我有下表: 考虑到每个月的权重,我的目标是在 5 个月内分散 10,000 个单位。与 10,000 相邻的行是我最好的尝试(我在这上面花了几个小时)。黄色是我所追求的。 我试图用来计算的逻辑如下:计
我的表单中有一个字段,它是文件类型。当用户点击保存图标时,我想自然地将文件上传到服务器并将文件名保存在数据库中。我尝试通过回显文件名来测试它,但它似乎不起作用。另外,如何将文件名添加到数据库中?是在模
我有一个 python 脚本来发送电子邮件,它工作得很好,但问题是当我检查我的电子邮件收件箱时。 我希望该用户名是自定义用户名,而不是整个电子邮件地址。 最佳答案 发件人地址应该使用的格式是: You
我想减小 ggcorrplot 中标记的大小,并减少文本和绘图之间的空间。 library(ggcorrplot) data(mtcars) corr <- round(cor(mtcars), 1)
GTK+ noob 问题在这里: 是否可以自定义 GtkFileChooserButton 或 GtkFileChooserDialog 以删除“位置”部分(左侧)和顶部的“位置”输入框? 我实际上要
我正在尝试在主页上使用 ajax 在 magento 中使用 ajax 显示流行的产品列表,我可以为 5 或“N”个产品执行此操作,但我想要的是将分页工具栏与结果集一起添加. 这是我添加的以显示流行产
我正在尝试使用 PasswordResetForm 内置函数。 由于我想要自定义表单字段,因此我编写了自己的表单: class FpasswordForm(PasswordResetForm):
据我了解,新的 Angular 7 提供了拖放功能。我搜索了有关 DnD 的 Tree 组件,但没有找到与树相关的内容。 我在 Stackblitz 上找到的一个工作示例.对比drag'ndrop功能
我必须开发一个自定义选项卡控件并决定使用 WPF/XAML 创建它,因为我无论如何都打算学习它。完成后应该是这样的: 到目前为止,我取得了很好的进展,但还有两个问题: 只有第一个/最后一个标签项应该有
我要定制xtable用于导出到 LaTeX。我知道有些问题是关于 xtable在这里,但我找不到我要找的具体东西。 以下是我的表的外观示例: my.table <- data.frame(Specif
用ejs在这里显示日期 它给我结果 Tue Feb 02 2016 16:02:24 GMT+0530 (IST) 但是我需要表现为 19th January, 2016 如何在ejs中执行此操作?
我想问在 JavaFX 中使用自定义对象制作 ListView 的最佳方法,我想要一个每个项目如下所示的列表: 我搜了一下,发现大部分人都是用细胞工厂的方法来做的。有没有其他办法?例如使用客户 fxm
我是一名优秀的程序员,十分优秀!