我正在测量每天的剧集持续时间(ep.dur
以分钟为单位的向量),为期数T=364
天的观察期。该向量ep.dur
有一个length(ep.dur)
,T=364
在没有发生事件的天数中为零,并且range(ep.dur)
在 0 到 1440 之间
T 周期内的情节持续时间之和为a<-sum(ep.duration)
现在我有一个向量den
,带有length(den)=99
. 向量 den 显示每 1% (1%, 2%, 3%, ...) 的开发需要多少天a
现在给出 den
and a
,我想模拟多个ep.dur
这可能吗?
澄清 1: : (danas.zuokas 的第一条评论) 的元素den
代表持续时间而不是确切的天数。这意味着,例如 1,1a
天开发 1%(=1195.8),2 天开发 2%,3 天开发 3%,4 天开发 4%,5 天开发 5%,5天开发6% ......)。这些剧集可以在 T 中的任何地方发生
澄清 2:(danas.zuokas 的第二条评论)不幸的是,无法假设持续时间如何发展。这就是为什么我必须模拟大量 ep.dur 向量。但是,如果这有任何帮助,我可以将 den 向量扩展为更有限的分辨率(即:而不是 1% 的跳跃,0.1% 的跳跃)。
算法描述
算法应该满足 den 向量提供的所有信息。我想象算法如下(示例 3): a 的每 1% 跳跃是 335,46 分钟。den[1]
告诉我们 1% 的 a 是在 1 天内开发的。所以假设我们生成ep.dur[1]
=335,46。好的。我们去den[2]
: 2% 的 a 是在d[2]
=1 天内开发的。因此,ep.dur[1]
不能为 335,46 并被拒绝(一天内仍应出现 2% 的 a)。可以说产生了ep.dur[1]
= 1440。d[1]
满意,满意d[2]
(至少 2% 的总持续时间在dur[2]
=1 天内开发),dur[3]
=1 也满意。守门员?但是,dur[4]
如果 ep.dur[1]=1440 不满足 =2,因为它表明 4% 的 a (=1341) 应该在 2 天内发生。所以ep.dur[1]
被拒绝。现在让我们说ep.dur[1]
= 1200。dur[1:3]
被接受。然后我们生成ep.dur[2]
等等,确保生成的 ep.dur 都满足 den 提供的信息。
这在编程上可行吗?我真的不知道从哪里开始这个问题。一旦赏金开始期结束,我将提供慷慨的赏金
示例 1:
a<-119508
den<-c(1, 2, 3, 4, 5, 5, 6, 7, 8, 9, 10, 10, 11, 12, 13, 14, 15, 15,
16, 17, 18, 19, 20, 20, 21, 22, 23, 24, 25, 25, 26, 27, 28, 29,
30, 30, 31, 32, 33, 34, 35, 35, 36, 37, 38, 39, 40, 40, 41, 42,
43, 44, 45, 45, 46, 47, 48, 49, 50, 50, 51, 52, 53, 54, 55, 55,
56, 57, 58, 59, 60, 60, 61, 62, 63, 64, 65, 65, 66, 67, 68, 69,
70, 70, 71, 72, 73, 74, 75, 75, 76, 77, 78, 79, 80, 80, 81, 82,
83)
示例 2:
a<-78624
den<-c(1, 2, 2, 3, 3, 4, 4, 5, 5, 6, 7, 7, 8, 8, 9, 9, 10, 10, 11,
11, 12, 13, 13, 14, 14, 15, 15, 16, 16, 17, 18, 19, 21, 22, 23,
28, 32, 35, 36, 37, 38, 43, 52, 55, 59, 62, 67, 76, 82, 89, 96,
101, 104, 115, 120, 126, 131, 134, 139, 143, 146, 153, 160, 165,
180, 193, 205, 212, 214, 221, 223, 227, 230, 233, 234, 235, 237,
239, 250, 253, 263, 269, 274, 279, 286, 288, 296, 298, 302, 307,
309, 315, 320, 324, 333, 337, 342, 347, 352)
示例 3
a<-33546
den<-c(1, 1, 1, 2, 4, 6, 8, 9, 12, 15, 17, 21, 25, 29, 31, 34, 37,
42, 45, 46, 51, 52, 56, 57, 58, 59, 63, 69, 69, 71, 76, 80, 81,
87, 93, 95, 102, 107, 108, 108, 112, 112, 118, 123, 124, 127,
132, 132, 132, 135, 136, 137, 150, 152, 162, 166, 169, 171, 174,
176, 178, 184, 189, 190, 193, 197, 198, 198, 201, 202, 203, 214,
218, 219, 223, 225, 227, 238, 240, 246, 248, 251, 254, 255, 257,
259, 260, 277, 282, 284, 285, 287, 288, 290, 294, 297, 321, 322,
342)
示例 4
a<-198132
den<-c(2, 3, 5, 6, 7, 9, 10, 12, 13, 14, 16, 17, 18, 20, 21, 23, 24,
25, 27, 28, 29, 31, 32, 34, 35, 36, 38, 39, 40, 42, 43, 45, 46,
47, 49, 50, 51, 53, 54, 56, 57, 58, 60, 61, 62, 64, 65, 67, 68,
69, 71, 72, 74, 75, 76, 78, 79, 80, 82, 83, 85, 86, 87, 89, 90,
91, 93, 94, 96, 97, 98, 100, 101, 102, 104, 105, 107, 108, 109,
111, 112, 113, 115, 116, 120, 123, 130, 139, 155, 165, 172, 176,
178, 181, 185, 190, 192, 198, 218)