MATLAB: Help with simple probability calculation (loop?)

loopMATLAB

I have a 2-column matrix that looks like this:

"1" indicates an occurrence and "0" does not. I need to compute a probability for each unique number in column 1. So the answers should be:

30 = 1
29 = .6667
28 = .50

Now, I have the "count_unique" function from the file exchange that counts the unique numbers and identifies how many occurrences there are for each number. Using this, I have tried to write a for loop that computes this problem, but I am stuck in the mud. I'm grateful for any assistance. Happy Thanksgiving.

Best Answer

I don’t know how robust this is, but it works here:

M = [30     1
    30     1
    30     1
    30     1
    30     1
    29     0
    29     1
    29     1
    28     1
    28     1
    28     0
    28     0];
raw_freq = accumarray(M(:,1), 1);
adj_freq = accumarray(M(M(:,2)>0), 1);
Prob = adj_freq./raw_freq;
Mu = unique(M(:,1));
Result = flipud([Mu Prob(~isnan(Prob))])
Result =
           30            1
           29      0.66667
           28          0.5

Related Solutions

MATLAB: Select data from a text file

Reading the file is easy with dlmread or readtable.

For calculating the means of identical epoch, you can use the Split-Apply-Combine workflow introduced in R2015b, or unique with accumarray:

 [epoch, ~, subs] = unique(m(:, [1 2]), 'rows'); %find unique epochs as rows of column 1 and 2
 meancol5 = accumarray(subs, m(:, 5), [], @mean); %mean of column 5 by epoch

MATLAB: Error with extract results from tomRun

This is the code in an embedeed block

function Wnlop  = loadsharing_Linearization_MPC(Wnlop3)%,e)
   Wnl = 1;
   dw = zeros (2,10);
   dw = [Wnl Wnl Wnl Wnl Wnl Wnl Wnl Wnl Wnl Wnl;
   Wnl Wnl Wnl Wnl Wnl Wnl Wnl Wnl Wnl Wnl];
FF = [ 945.8209  868.6367
       868.6367  797.7697];
S = [ -273.5275
 -251.1619];
Alpha = [ 0         0
    7.9634    7.3155
         0         0
    0.1303    0.1307
    7.9634    7.3155
   -7.9634   -7.3155
         0         0];
b_L = [29.1083   28.1583   28.1583   28.1583   28.1583   28.1583   28.1583   28.1583   28.1583   28.1583
   29.1083   28.1583   28.1583   28.1583   28.1583   28.1583   28.1583   28.1583   28.1583   28.1583
      -Inf   29.1540   29.1561   29.1561   29.1561   29.1561   29.1561   29.1561   29.1561   29.1561
      -Inf   29.1561   27.1562   -0.0366   14.0596   14.0597   14.0597  -14.0597         0         0
   57.3266   56.3166   56.3166   56.3166   56.3166   56.3166   56.3166   56.3166   56.3166   56.3166
      -Inf         0         0         0         0         0         0         0         0         0
   28.1583   28.1583   28.1583   28.1583   28.1583   28.1583   28.1583   28.1583   28.1583   28.1583];
b_U =[ 29.2083   28.1583   28.1583   28.1583   28.1583   28.1583   28.1583   28.1583   28.1583   28.1583
   29.2083   28.1583   28.1583   28.1583   28.1583   28.1583   28.1583   28.1583   28.1583   28.1583
   31.6583   29.1540   29.1561   29.1561   29.1561   29.1561   29.1561   29.1561   29.1561   29.1561
   31.4061   29.1561   27.1562   -0.0366   14.0596   14.0597   14.0597  -14.0597         0         0
       Inf   56.3166   56.3166   56.3166   56.3166   56.3166   56.3166   56.3166   56.3166   56.3166
    0.9700         0         0         0         0         0         0         0         0         0
       Inf   28.1583   28.1583   28.1583   28.1583   28.1583   28.1583   28.1583   28.1583   28.1583];
     eml.extrinsic('extractResult');
eml.extrinsic('qpAssign','tomRun');
    Prob = qpAssign(FF, S, Alpha, b_L(1:7), b_U(1:7), [], [], [], [],...
            [], [], [], [], [], [], []);
 Result = tomRun('qpSolve', Prob, 1)
 eml.extrinsic('extractResult');
 dwp = extractResult(Result)
 dwopt = [0 0];
 dwopt = dw
 wstart = [0.98 0.99]
  Wnlop = wstart + dwopt ;

Best Answer

Related Solutions

MATLAB: Select data from a text file

MATLAB: Error with extract results from tomRun

Related Question