2017-02-22 20 views
-1

毎週、毎月、ワイブル分布で風力をシミュレートしています。毎週風力/風力の違いがあるので、それらをすべて組み合わせて、その結果を複製して特定の風の場所が特定の時間にどのくらいのエネルギーを作り出すかを確認します。最初のシミュレーションでは、1週間目の結果が生成され、2週間目の2回目のシミュレーションなどが行われ、1か月間のマージ結果が得られ、目的の月数の​​結果が複製されます。2つ以上のシミュレーションをマージする

first<-rweibull(7,3.28,2.14) #time series for days in first week 
second<-rweibull(7,3.33,2.112)# time series for days in second week 
combined<-replicate(2,c(first,second))#combined 

    [,1]  [,2] 
[1,] 2.3314627 2.3314627 
[2,] 2.4219863 2.4219863 
[3,] 2.2993156 2.2993156 
[4,] 0.9576821 0.9576821 
[5,] 0.8065639 0.8065639 
[6,] 2.2820711 2.2820711 
[7,] 0.6315118 0.6315118 
[8,] 1.8152144 1.8152144 
[9,] 0.9296449 0.9296449 
[10,] 1.6707807 1.6707807 
[11,] 1.7307920 1.7307920 
[12,] 2.2075004 2.2075004 
[13,] 2.0368642 2.0368642 
[14,] 2.5320678 2.5320678  

ここで問題は、私が同じ結果を得るが、結果が異なるはずだということです。そのための簡単で速い解決策がありますか?

+0

あなたはもう少し説明できますか?期待される成果はどうですか?私はあなたが複製する関数を渡しているとは思わないが、結合して複製を使用すると混乱する... 'cbind'も結合する必要がある – user5249203

+0

毎週と毎月のワイブル分布をシミュレートする年。毎週風力/風力の違いがあるので、それらをすべて組み合わせたいと思って、その結果を複製して特定の時間に特定の風の場所が作り出すエネルギーを確認してください。 – kelamahim

+0

また、このようなものを 'replicate(2、c (rweibull(7,3.28,2.14)、rweibull(7,3.33,2.112))) '。変数に格納するのではなく、後でそれを渡す – user5249203

答えて

1

replicateに変数を渡す代わりに、これを行うことができます。

コード

fun <- function(){ 
    first <- rweibull(7,3.28,2.14) 
    second <- rweibull(7,3.33,2.112) 
} 

replicate(2, fun()) 
関連する問題