Random draws from conditional multinomial distribution

3 次查看（过去 30 天）
Laura Freitas
Laura Freitas 2021-10-26
评论： Laura Freitas 2021-10-26
I have to simulate data for 100 agents over 10,000 time periods for a variable, x. The initial value of x at time 0 for all agents is 0. The distribution of is a multinomial distribution, conditional on the past period's value, . Basically, can remain unchanged with probability 0.2, increase by 2571 with probability 0.5, increase by 2571*2 with probability 0.2, and increase by 2571*3 with probability 0.1, at any period. can never be lower than . How do I generate random values for this distribution?
I have tried:
x = mnrnd()
command, but I can't fit it to these specifications. Any help would be greatly appreciated!

请先登录，再回答此问题。

采纳的回答

Alan Stevens
Alan Stevens 2021-10-26
Something like
u = rand(1)
if u<0.2
d = 0;
elseif u<0.7
d = 2571;
elseif u<0.9
d = 2571^2;
else
d = 2571^3;
end
x(t) = x(t-1) + d;
(the last line assumes t is an integer).
  5 个评论
显示 3更早的评论
Alan Stevens
Alan Stevens 2021-10-26
function xnew = fn(x,t)
u = rand;
if u<0.2
d = 0;
elseif u<0.7
d = 2571;
elseif u<0.9
d = 2571^2;
else
d = 2571^3;
end
x(t) = x(t-1) + d;
end
Note that functions go at the end of a script.

请先登录，再进行评论。

更多回答（0 个）

类别

在 Help Center 和 File Exchange 中查找有关 Multinomial Distribution 的更多信息

标签

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by