2017-02-22 13 views
-1

Я имитирую энергию ветра с распределением weibull на каждую неделю и каждый месяц в году. Поскольку есть недельные различия в скорости и мощности ветра, я хочу объединить их все, а затем воспроизвести результат, чтобы узнать, сколько энергии может определить конкретное местоположение ветра в определенное время. Первая симуляция дает результаты для первой недели второго моделирования для второй недели и т. Д., При этом цель объединяет результаты в течение одного месяца, а затем реплицирует результаты для желаемого количества месяцев.слияние двух или более симуляций

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

Я имитирую энергию ветра с распределением weibull на каждую неделю и каждый месяц в год. Поскольку есть недельные различия в скорости и мощности ветра, я хочу объединить их все, а затем воспроизвести результат, чтобы узнать, сколько энергии может определить конкретное местоположение ветра в определенное время – kelamahim

+0

, или вы должны передать ему что-то вроде этого 'replicate (2, c (rweibull (7,3,28,2,14), rweibull (7,3,33,2,121))). вместо того, чтобы хранить в переменной adn передать ее позже – user5249203

ответ

1

Вместо передачи переменных replicate вы можете это сделать.

код

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

replicate(2, fun()) 
Смежные вопросы