I am estimating the probabilities over time.
for example, my data matrix is 6 by 5 (where 6 is the number of simulations and 5 is the number of time periods).
data = [2, 2 , 1, 1, 1, 1; 2, 2 , 1, 1, 1, 2; 2, 2 , 1, 2, 1, 2; 2, 2 , 1, 2, 1, 2; 2, 3 , 2, 2, 1, 2]States = unique([data(:,1); data(:,2)]) % Find the unique rows of A based on the data in the first two columns.
[TFinitial, fromstateAge27] = ismember(sequence(:,3), States) % logical arrays
[TF25, tostateAge27] = ismember(sequence(:,4), States) % logical array
went_from_to_countAge27 = accumarray( [fromstateAge27(:), tostateAge27(:)], 1, []) % frequency of transitions from one state to another
went_from_to_probAge27 = went_from_to_countAge27 ./ sum(went_from_to_countAge27) % probability of transition from one state to another
Now I want to make a for loop that runs the above code for each column of matrix "data".
For the above code (column 1 to 2), my output is "went_from_to_probAge27 "
for column 2 to 3, my output will be "went_from_to_probAge28"
from column 3 to 4, my output will be "went_from_to_probAge29" and so on.
Best Answer