генерация наборов значений из нескольких распределений в цикле for в R - PullRequest
0 голосов
/ 25 апреля 2018

Я хочу составить список векторов / данных из нового набора значений, полученных из 3 разных распределений, два - lnorm, а 1 - смешанное равномерное распределение.Мне нужно сложить сгенерированные значения каждого из этих распределений в одно значение (v1 + v2 + v3), чтобы согласовать распределение с суммой из 3 распределений.

Мой код не суммирует сумму распределений в одном наборе данных (с функцией: list), и я не уверен, работает ли часть кода функции UnivarMixingDistribution?Я хочу, чтобы он сгенерировал одно значение, сохранил его в xv3 и добавил к двум другим значениям.

Отображается ошибка:

In e[i] <- list(v1, v2, xv3) :
  number of items to replace is not a multiple of replacement length

код:

e <- list(mode="vector",length=n)
for(i in 1:1000){
    v1=rlnorm(n, mean, sd) #simulate values for lnorm distribution
    v2=rlnorm(n, mean, sd)
    v3=myMixv4<-UnivarMixingDistribution(Unif(Min=, Max=), #1st distribution
                                         Unif(Min=, Max=), #2nd distribution
                                         mixCoeff=c((ND/n),
                                                    (Detected/n)))
    rmyMixv3<-r(myMixv3)
    xv3 <- rmyMixv3(1)
    e[i] <- list(v1,v2,xv3)
}
Добро пожаловать на сайт PullRequest, где вы можете задавать вопросы и получать ответы от других членов сообщества.
...