ep.dur数日間の観察期間で、1 日あたりのエピソードの持続時間 (分単位のベクトル) を測定していT=364ます。ベクトルep.durのlength(ep.dur)_ T=364_ range(ep.dur)_
T 期間にわたるエピソード期間の合計はa<-sum(ep.duration)
今、私は ベクトル を持っていdenますlength(den)=99。ベクトル den は、各 1% (1%、2%、3%、...) の開発に必要な日数を示しています。a
と が与えられた ので、複数をシミュレートしたいとden思いますaep.dur
これは可能ですか？
明確化 1: : (danas.zuokas の最初のコメント) の要素は、正確な日数ではなく期間denを表します。つまり、たとえば 1、1%(=1195.8) が 1 日で、2% が 2 日で、3% が 3 日で、4% が 4 日で、5% が5 日で、6% が5 日で開発されることを意味します。 .....）。エピソードは T のどこでも発生する可能性がありますa
明確化 2: (danas.zuokas の 2 番目のコメント) 残念ながら、期間がどのように発展するかについての仮定はありません。そのため、多数の ep.dur ベクトルをシミュレートする必要があります。ただし、これが役立つ場合は、デンベクトルをより有限の解像度に拡張できます (つまり、1% のジャンプではなく、0.1% のジャンプ)。
アルゴリズムの説明
アルゴリズムは、den ベクトルが提供するすべての情報を満たす必要があります。アルゴリズムが次のようになると想像しました (例 3): a の 1% の各ジャンプは 335.46 分です。den[1]は、a の 1% が 1 日で開発されることを示しています。ep.dur[1]=335,46を生成するとします。わかった。den[2]: a の 2% は =1 日で開発されますd[2]。そのため、ep.dur[1]335,46 にすることはできず、拒否されます (a の 2% は 1 日で発生するはずです)。ep.dur[1]=1440を生成したとしましょう。d[1]満足している、満足してd[2]いる (合計期間の少なくとも 2% がdur[2]=1 日で開発されている)、dur[3]=1 も満足しています。キーパー？ただし、dur[4]ep.dur[1]=1440 の場合、=2 は成立しません。これは、a (=1341) の 4% が 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)
