个人中心
文章发布
退出
作者热门文章
- html - 出于某种原因，IE8 对我的 Sass 文件中继承的 html5 CSS 不友好？
- JMeter 在响应断言中使用 span 标签的问题
- html - 在 :hover and :active? 上具有不同效果的 CSS 动画
- html - 相对于居中的 html 内容固定的 CSS 重复背景？
25
4
我正在尝试使用 ggplot 绘制引导样本(1000 列)，但结果图不正确。
这是数据:
dat<- structure(list(YY = c(1961L, 1961L, 1961L, 1961L, 1961L, 1961L,
1961L, 1961L, 1961L, 1961L, 1961L, 1961L, 1961L, 1961L, 1961L,
1961L, 1961L, 1961L, 1961L, 1961L, 1961L, 1961L, 1961L, 1961L,
1961L, 1961L, 1961L, 1961L, 1961L, 1961L, 1961L, 1961L, 1961L,
1961L, 1961L, 1961L, 1961L, 1961L, 1961L, 1961L, 1961L, 1961L,
1961L, 1961L, 1961L, 1961L, 1961L, 1961L, 1961L, 1961L, 1961L,
1961L, 1961L, 1961L, 1961L, 1961L, 1961L, 1961L, 1961L, 1961L,
1961L, 1961L, 1961L, 1961L, 1961L, 1961L, 1961L, 1961L, 1961L,
1961L, 1961L, 1961L, 1961L, 1961L, 1961L, 1961L, 1961L, 1961L,
1961L, 1961L, 1961L, 1961L, 1961L, 1961L, 1961L, 1961L, 1961L,
1961L, 1961L, 1961L, 1961L, 1961L, 1961L, 1961L, 1961L, 1961L,
1961L, 1961L, 1961L, 1961L, 1961L, 1961L, 1961L, 1961L, 1961L,
1961L, 1961L, 1961L, 1961L, 1961L, 1961L, 1961L, 1961L, 1961L,
1961L, 1961L, 1961L, 1961L, 1961L, 1961L, 1961L, 1961L, 1961L,
1961L, 1961L, 1961L, 1961L, 1961L, 1961L, 1961L, 1961L, 1961L,
1961L, 1961L, 1961L, 1961L, 1961L, 1961L, 1961L, 1961L, 1961L,
1961L, 1961L, 1961L, 1961L, 1961L, 1961L, 1961L, 1961L, 1961L,
1961L, 1961L, 1961L, 1961L, 1961L, 1961L, 1961L, 1961L, 1961L,
1961L, 1961L, 1961L, 1961L, 1961L, 1961L, 1961L, 1961L, 1961L,
1961L, 1961L, 1961L, 1961L, 1961L, 1961L, 1961L, 1961L, 1961L,
1961L, 1961L, 1961L, 1961L, 1961L, 1961L, 1961L, 1961L, 1961L,
1961L, 1961L, 1961L, 1961L, 1961L, 1961L, 1961L, 1961L, 1961L,
1961L, 1961L, 1961L, 1961L, 1961L, 1961L, 1961L, 1961L, 1961L,
1961L, 1961L, 1961L, 1961L, 1961L, 1961L, 1961L, 1961L, 1961L,
1961L, 1961L, 1961L, 1961L, 1961L, 1961L, 1961L, 1961L, 1961L,
1961L, 1961L, 1961L, 1961L, 1961L, 1961L, 1961L, 1961L, 1961L,
1961L, 1961L, 1961L, 1961L, 1961L, 1961L, 1961L, 1961L, 1961L,
1961L, 1961L, 1961L, 1961L, 1961L, 1961L, 1961L, 1961L, 1961L,
1961L, 1961L, 1961L, 1961L, 1961L, 1961L, 1961L, 1961L, 1961L,
1961L, 1961L, 1961L, 1961L, 1961L, 1961L, 1961L, 1961L, 1961L,
1961L, 1961L, 1961L, 1961L, 1961L, 1961L, 1961L, 1961L, 1961L,
1961L, 1961L, 1961L, 1961L, 1961L, 1961L, 1961L, 1961L, 1961L,
1961L, 1961L, 1961L, 1961L, 1961L, 1961L, 1961L, 1961L, 1961L,
1961L, 1961L, 1961L, 1961L, 1961L, 1961L, 1961L, 1961L, 1961L,
1961L, 1961L, 1961L, 1961L, 1961L, 1961L, 1961L, 1961L, 1961L,
1961L, 1961L, 1961L, 1961L, 1961L, 1961L, 1961L, 1961L, 1961L,
1961L, 1961L, 1961L, 1961L, 1961L, 1961L, 1961L, 1961L, 1961L,
1961L, 1961L, 1961L, 1961L, 1961L, 1961L, 1961L, 1961L, 1961L,
1961L, 1961L, 1961L, 1961L, 1961L, 1961L, 1961L, 1961L, 1961L,
1961L, 1961L, 1961L, 1961L, 1961L, 1961L, 1961L, 1961L, 1961L,
1961L, 1961L, 1961L, 1961L, 1961L, 1961L, 1961L, 1961L), MM = c(1L,
1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L,
1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 2L,
2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L,
2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 3L, 3L, 3L, 3L, 3L, 3L,
3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L,
3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 4L, 4L, 4L, 4L, 4L, 4L, 4L,
4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L,
4L, 4L, 4L, 4L, 4L, 4L, 4L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L,
5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L,
5L, 5L, 5L, 5L, 5L, 5L, 6L, 6L, 6L, 6L, 6L, 6L, 6L, 6L, 6L, 6L,
6L, 6L, 6L, 6L, 6L, 6L, 6L, 6L, 6L, 6L, 6L, 6L, 6L, 6L, 6L, 6L,
6L, 6L, 6L, 6L, 7L, 7L, 7L, 7L, 7L, 7L, 7L, 7L, 7L, 7L, 7L, 7L,
7L, 7L, 7L, 7L, 7L, 7L, 7L, 7L, 7L, 7L, 7L, 7L, 7L, 7L, 7L, 7L,
7L, 7L, 7L, 8L, 8L, 8L, 8L, 8L, 8L, 8L, 8L, 8L, 8L, 8L, 8L, 8L,
8L, 8L, 8L, 8L, 8L, 8L, 8L, 8L, 8L, 8L, 8L, 8L, 8L, 8L, 8L, 8L,
8L, 8L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L,
9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L,
10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L,
10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L,
10L, 10L, 10L, 10L, 10L, 11L, 11L, 11L, 11L, 11L, 11L, 11L, 11L,
11L, 11L, 11L, 11L, 11L, 11L, 11L, 11L, 11L, 11L, 11L, 11L, 11L,
11L, 11L, 11L, 11L, 11L, 11L, 11L, 11L, 11L, 12L, 12L, 12L, 12L,
12L, 12L, 12L, 12L, 12L, 12L, 12L, 12L, 12L, 12L, 12L, 12L, 12L,
12L, 12L, 12L, 12L, 12L, 12L, 12L, 12L, 12L, 12L, 12L, 12L, 12L,
12L), DD = 1:365, RR = c(0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 5.3, 15.7, 3.6, 11.7, 40.1, 0, 0, 34, 9.4,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 16.3, 15, 0,
0, 0, 4.6, 0, 0, 0, 5.3, 14.7, 7.1, 0, 4.1, 0.5, 0, 1.8, 1.3,
0, 0, 0, 0, 0, 39.4, 1, 25.9, 0, 4.3, 4.6, 2, 9.7, 0, 0, 2, 0,
0, 0.3, 57.9, 31, 2.3, 0, 1.8, 0, 0, 0, 0, 20.8, 11.7, 1, 3.1,
0, 0, 0, 0.3, 8.6, 0, 0.3, 5.1, 0, 0, 0, 10.9, 7.1, 0, 4.1, 1.5,
2.8, 15.8, 109.5, 140, 25.7, 57.2, 9.4, 140.2, 119.6, 31.2, 0.8,
43.9, 2.8, 52.3, 45.5, 0.8, 12.4, 29, 66.3, 116.6, 155.7, 22.6,
6.9, 21.3, 39.1, 21.6, 19.8, 0, 14, 0, 11.2, 0, 14.2, 11.7, 52.8,
74.4, 34, 5.1, 0, 0, 0, 0, 2.5, 0.3, 6.6, 6.1, 31.8, 6.9, 7.6,
1.8, 6.6, 7.4, 8.4, 0, 20.8, 3, 15.5, 8.9, 6.4, 11.9, 11.7, 127,
100.3, 56.4, 6.1, 28.4, 84.6, 33.3, 6.4, 5.1, 30.5, 4.6, 4.1,
19.3, 39.6, 27.2, 7.9, 3.6, 3.3, 9.7, 10.7, 1.3, 59.7, 51.6,
65.8, 28.7, 13.2, 2.8, 32.5, 10.7, 57.2, 2.5, 0, 0, 29.5, 0,
0, 6.6, 31.2, 11.7, 1.5, 0, 0.5, 1.5, 29.7, 13.2, 20.3, 34.3,
9.7, 1, 1.5, 4.1, 13.5, 33, 0, 0, 10.9, 1.5, 7.4, 13, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 5.3, 0, 0, 8.9, 13.7, 0, 0, 0, 0,
0, 0, 18.3, 0, 0, 0, 0, 0, 0, 0, 3.6, 14.2, 7.6, 0, 0.5, 0, 0,
0, 0, 0, 0, 0, 0, 0, 6.6, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1), CUMSUM = c(0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 5.3, 21, 24.6,
36.3, 76.4, 76.4, 76.4, 110.4, 119.8, 119.8, 119.8, 119.8, 119.8,
119.8, 119.8, 119.8, 119.8, 119.8, 119.8, 119.8, 119.8, 119.8,
119.8, 119.8, 119.8, 119.8, 119.8, 119.8, 119.8, 119.8, 119.8,
119.8, 119.8, 119.8, 119.8, 119.8, 119.8, 119.8, 119.8, 119.8,
119.8, 119.8, 119.8, 119.8, 119.8, 119.8, 136.1, 151.1, 151.1,
151.1, 151.1, 155.7, 155.7, 155.7, 155.7, 161, 175.7, 182.8,
182.8, 186.9, 187.4, 187.4, 189.2, 190.5, 190.5, 190.5, 190.5,
190.5, 190.5, 229.9, 230.9, 256.8, 256.8, 261.1, 265.7, 267.7,
277.4, 277.4, 277.4, 279.4, 279.4, 279.4, 279.7, 337.6, 368.6,
370.9, 370.9, 372.7, 372.7, 372.7, 372.7, 372.7, 393.5, 405.2,
406.2, 409.3, 409.3, 409.3, 409.3, 409.6, 418.2, 418.2, 418.5,
423.6, 423.6, 423.6, 423.6, 434.5, 441.6, 441.6, 445.7, 447.2,
450, 465.8, 575.3, 715.3, 741, 798.2, 807.6, 947.8, 1067.4, 1098.6,
1099.4, 1143.3, 1146.1, 1198.4, 1243.9, 1244.7, 1257.1, 1286.1,
1352.4, 1469, 1624.7, 1647.3, 1654.2, 1675.5, 1714.6, 1736.2,
1756, 1756, 1770, 1770, 1781.2, 1781.2, 1795.4, 1807.1, 1859.9,
1934.3, 1968.3, 1973.4, 1973.4, 1973.4, 1973.4, 1973.4, 1975.9,
1976.2, 1982.8, 1988.9, 2020.7, 2027.6, 2035.2, 2037, 2043.6,
2051, 2059.4, 2059.4, 2080.2, 2083.2, 2098.7, 2107.6, 2114, 2125.9,
2137.6, 2264.6, 2364.9, 2421.3, 2427.4, 2455.8, 2540.4, 2573.7,
2580.1, 2585.2, 2615.7, 2620.3, 2624.4, 2643.7, 2683.3, 2710.5,
2718.4, 2722, 2725.3, 2735, 2745.7, 2747, 2806.7, 2858.3, 2924.1,
2952.8, 2966, 2968.8, 3001.3, 3012, 3069.2, 3071.7, 3071.7, 3071.7,
3101.2, 3101.2, 3101.2, 3107.8, 3139, 3150.7, 3152.2, 3152.2,
3152.7, 3154.2, 3183.9, 3197.1, 3217.4, 3251.7, 3261.4, 3262.4,
3263.9, 3268, 3281.5, 3314.5, 3314.5, 3314.5, 3325.4, 3326.9,
3334.3, 3347.3, 3347.3, 3347.3, 3347.3, 3347.3, 3347.3, 3347.3,
3347.3, 3347.3, 3347.3, 3347.3, 3347.3, 3347.3, 3352.6, 3352.6,
3352.6, 3361.5, 3375.2, 3375.2, 3375.2, 3375.2, 3375.2, 3375.2,
3375.2, 3393.5, 3393.5, 3393.5, 3393.5, 3393.5, 3393.5, 3393.5,
3393.5, 3397.1, 3411.3, 3418.9, 3418.9, 3419.4, 3419.4, 3419.4,
3419.4, 3419.4, 3419.4, 3419.4, 3419.4, 3419.4, 3419.4, 3426,
3426, 3426, 3426, 3426, 3426, 3426, 3426, 3426, 3426, 3426, 3426,
3426, 3426, 3426, 3426, 3426, 3426, 3426, 3426, 3426, 3426, 3426,
3426, 3426, 3426, 3427)), .Names = c("YY", "MM", "DD", "RR",
"CUMSUM"), class = "data.frame", row.names = c(NA, -365L))
这是我的脚本:
library(reshape2)
library(ggplot2)
library(meboot)
set.seed(1234)
dat[is.na(dat)]<- " "
out <- meboot(x=dat$CUMSUM,reps=1000,trim=0.10,elaps=TRUE)
f<-out$ens
Z<-out$ens
colnames(Z)<-c(seq(1:1000))
ZZ<-cbind(dat[3],Z)
df.melted<-melt(ZZ,id="DD")
p<-ggplot()
p <- p + geom_line(data=df.melted,aes(x=DD,y=value),color="black",alpha=0.1,size=0.5)
p <-p + geom_line(data=dat,aes(x=DD,y=CUMSUM),color="red",size=1)
p<-p + theme(panel.background=element_rect(fill="white"),
plot.margin = unit(c(0.5,0.5,0.5,0.5),"cm"),
panel.border=element_rect(colour="black",fill=NA,size=1),
axis.line.x=element_line(colour="black"),
axis.line.y=element_line(colour="black"),
axis.text=element_text(size=15,colour="black",family="serif"),
axis.title=element_text(size=15,colour="black",family="serif"),
legend.position = "none")
p <- p + labs(x="Day of Year",y="Cummulative Rainfall (mm)")
outImg <- paste0("bootrapped",".png")
ggsave(outImg,p,width=8,height=5)
生成的图像是这样的:
问题/问题
melt。关于如何在 ggplot 中正确执行此操作的任何建议？非常感谢任何帮助。
最佳答案
您需要按 variable 分组.
p <- p + geom_line(data=df.melted,aes(x=DD,y=value, group = variable),color="black",alpha=0.1,size=0.5)
关于r - 使用 ggplot2 绘制自举样本(多列)，我们在Stack Overflow上找到一个类似的问题： https://stackoverflow.com/questions/51077729/
25
4
0
我正在处理不平衡的面板数据，我想从中抽取一个随机样本，该样本不受每个单位不同观察次数的影响。例如，在下面的代码中，IBM 被选中的可能性是 GOOG 的两倍，被选中的可能性是 MSFT 的五倍。有没有
有人可以指出实现 CollectionChanged 的示例。我正在使用 wpf mvvm 灯。我试图谷歌，没有找到任何足够好的东西。 最佳答案 public ObservableCollecti
我有 df我想对某些变量的分布进行一些抽样。比方说 df['type'].value_counts(normalize=True)返回: 0.3 A 0.5 B 0.2 C 我想做类似 sampled
我仍然无法理解样本/像素/片段之间有什么区别。 由于片段着色器按像素执行，我认为片段只是指一个像素，这是正确的吗？谁能给我一个例子和每个人的定义？ 最佳答案 片段着色器按片段执行并发出像素。它们非常相
我正在尝试理解这个名为“The Amazing Audio Engine”的 GitHub 项目，它简化了在 iOS 上处理音频的过程。 我从麦克风捕捉并使用这种方法: id receiver = [
如何在诺基亚 Qt SDK(用于手机)中使用 QTableView。我引用了一些文档，但我仍然不清楚 QTableView。请任何人建议如何使用 QTableView。 我想显示具有三列的 QTabl
我已经能够获取 JmDNS 附带的示例来编译和运行，但是我无法获取任何类来发现我的服务。 我正在运行一个 Windows 环境，多台 PC 运行 VNC、SSH 和 Apache，我一直在尝试让 Jm
问题的具体实例 我的整数范围是 1-100。我想生成此范围内的 n 个总数，这些数字尽可能均匀分布并包括第一个和最后一个值。 示例 start = 1, end = 100, n = 5 Outp
我在线程组中有几个带有脚本的 JSR 223 采样器，它们在执行在调用 HTTP 请求之前进行一些工作。 问题在于 JSR 233 采样器包含在最终摘要报告中 我的问题是如何从最终计算中排除那些 JS
我需要有关存储后端歌曲预览的好方法的建议(现在正在查看iTunes，也许还有spotify和soundcloud)。我的想法是，我需要预下载并可能缓存30秒及更少的音频文件，以方便召回。 然后，我需要
我刚刚从 Github 下载了 Atmosphere 样本。当我在聊天样本上运行 jetty:run goal 时，我遇到了一些问题。 我可以使用浏览器访问该页面(http://localhost:9
我有一个包含五个项目的向量。 my_vec 有更换，当我需要时没有更换。最有效的方法是什么？ 请注意，在我的向量中，我有两次值“a” - 因此，在返回的打乱向量集中，它们都应该有两次“a”。 最佳答案
我正在尝试学习如何阅读规范。让我们看看尝试压缩后会得到什么:1) 一个空缓冲区和 2) 一个感叹号: >>> zlib.compress(b'', 0) b'x\x01\x01\x00\x00\xff
我是 cuda 的新手，几周前才开始阅读有关并行编程和 cuda 的内容。在我安装了 cuda 工具包之后，我正在浏览 sdk 示例(安装工具包时附带的)并想尝试其中的一些。我从 0_Simple 文
我正在使用 FsCheck 生成自定义数据的 Gen . 假设你有一个函数返回 Gen : let chooseRectangle widthMax heightMax offset = gen
我有一个包含大约 800 万个观察值的数据框。我需要从中提取样本，但想从多个列中采样。 我尝试了以下方法，但不起作用: import pandas as pd state = ['mi', 'mi',
我学习安卓图形，我遇到了一个奇怪的问题: 我发现很多提到“FingerPaint”样本的地方，但我在样本文件夹中找不到它。 只有 47 个示例项目，没有一个是关于图形的。 如何下载此示例？我按照这
假设我有一个 pandas 数据框 rid category 0 0 c2 1 1 c3 2 2 c2 3 3 c3 4 4
我想用 Pandas sample功能，但具有不分组或过滤数据的标准。 import pandas as pd import numpy as np df = pd.DataFrame(np.rand
我正在尝试从 MNIST 数据集中绘制 10 个样本。每个数字之一。这是代码: import sklearn import pandas as pd import matplotlib.pyplot
我是一名优秀的程序员,十分优秀！