У меня есть функция gentimes (10000), которая генерирует массив значений (w1, w2) 10000x2, например:
gentimes = function(nreps)
{
x=c(rexp(nreps,1/5),rexp(nreps,1/5))
x1=min(x)
x2=max(x)
y1=rexp(nreps,1/4)
y2=rexp(nreps,1/4)
w1=y1
w2 = ifelse(x1+y1>x2,x1+y1-x2+y2,y2)
z=array(0,c(nreps,2))
z[,1] = w1
z[,2] = w2
z
}
Я хотел бы просмотреть выходные данные и выбрать случаи, для которых w1 <1, а затем нарисуйте гистограмму значений w2 для этих случаев. Я пробовал это: </p>
z = gentimes(10000)
for (i in 1:10000)
{
if (z[1,i]<1)
{
new_array=array(0,c(10000,2))
new_array[i]=c(z[1,i],z[2,i])
}
new_array
}
Я новичок в R и пытаюсь создать новый массив новых дел. Очень благодарен за любую помощь:)