search

markusheinonen / bamfa

Bayesian metabolic flux analysis
9 stars 0 forks source link

Bmfa function error #1

Closed jmayoral1 closed 8 months ago

jmayoral1 commented 9 months ago

Hello Markus,

Thank you and your team for the developing this approach! It seems like it's been a few years since this code was created, but I'm nonetheless interested in giving Bayesian flux sampling a try. However, I encounter an error when calling the bmfa function, shown below:

model = readCbModel('Ec_core_flux1.xml'); %% attempted with some of my own genome-scale models as well, but similar error [sol,~] = bmfa(model,[],[],[],500,10,'gibbs',100); %% no experimentally measured fluxes of my own, so I left y, g, and sdx input empty

Optimising flux variances... Error using Incorrect dimensions for matrix multiplication. Check that the number of columns in the first matrix matches the number of rows in the second matrix. To perform elementwise multiplication, use '.'.

Error in bmfa (line 121) vC = Sv - A S Sv + A Sdx A';

Now admittedly, I don't fully grasp the nuances of matrix multiplication, so I was hoping you could provide some guidance on the possible sources of the error. Is element-wise multiplication the intended step here?

Any help would be greatly appreciated. In case it is relevant, I am using Matlab version R2021b and Cobra toolbox v3.0 Josh

jmayoral1 commented 8 months ago

Closing this issue - the code works fine when providing at least one flux prior, as demonstrated in the demos.