Discussion :

[Sansideevraiment]

Bonjour à tous,
Dans un premier temps je vous explique le contexte de mon problème. J'ai une structure "model", qui contient elle-même plusieurs structures. Celle qui m'intéresse se nomme "reaction". Dans cette structure, outre diverses informations, je souhaite manipuler une autre structure appelée "products". Cette dernière contient plein de champs mais celui que je veux se nomme "species" qui retourne les noms des espèces participant en tant que "reactants" à la réaction.
J'ai donc écrit ce bout de code :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
model = TranslateSBML('recon_2.2.xml');
 
appears_in_reaction = zeros(1000,1000);
 
reaction = model.reaction;
 
 
for i=1:7785
    listReactions{i} = reaction(i).id; 
end
 
% for i=1:7785
%     for j=1:4
%     listReact{i,j} = reaction(i).reactant(j).species;
%     end
% end
 
 
 
for i=1:1000
    for j=1:4
    listProd{i,j} = reaction(i).product(j).species;
    end
end

Mais j'obtient toujours le fameux "Index exceeds matrix dimensions", et là je ne sais pas pourquoi...
Quelqu'un pour m'éclairer svp ?

PS : "reaction" est une structure 1*7785, et "reactant" ne dépasse jamais 1*4.

Merci d'avance

[duf42]

Bonjour,

Tu pourrais remplacer les valeurs écrites "en dur" par le calcul de la taille de tes données à l'aide de la fonction LENGTH:

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
 
for i=1:length(reaction)
    listReactions{i} = reaction(i).id; 
end

Sinon, tu peux aussi simplifier cette ligne pour éviter d'utiliser une boucle:
[code]
listReactions = {reaction.id};
[code]

Duf