-
Notifications
You must be signed in to change notification settings - Fork 0
Description
I put the code here and will submit multiple PR in the coming days but you can already use it
Part of code to introduce around line 260
detailScoring={};
for j=1:SampleNumber
incR=1;
for i=1:210
if ~isempty(essentialRxns{i})
rxns=essentialRxns{i};
rxnID=findRxnIDs(model,rxns);
rxnID(rxnID==0)=[];
if ~isempty(rxnID)
for k=1:length(rxnID)
%1st column = sample ID
detailScoring{incR,((j-1)*8)+1}=j;
%2nd column = task ID
detailScoring{incR,((j-1)*8)+2}=i;
%3th column = task score for this sample
detailScoring{incR,((j-1)*8)+3}=ScorebyTask(i,j);
%4th column = task score in binary version for this sample
detailScoring{incR,((j-1)*8)+4}=ScorebyTask_binary(i,j);
%5th column = essential reaction associated to this
%task
detailScoring{incR,((j-1)*8)+5}=rxns(k);
%6th column = expression score associated to the
%reaction listed in the 5th column
detailScoring{incR,((j-1)*8)+6}=expression.Rxns(rxnID(k),SampleNumber);
%7th column = gene used to determine the expression of the
%reaction listed in the 5th column
geneName=expression.gene_used(rxnID(k),SampleNumber);
detailScoring{incR,((j-1)*8)+7}=geneName;
%8th column = original expression value of the gene
%listed in the 7th column
detailScoring{incR,((j-1)*8)+8}=data.value(strcmp(data.gene,geneName),SampleNumber);
incR=incR+1;
end
end
end
end
and add the output detailScoring in the call function
[score, score_binary ,taskInfos, detailScoring]=CodeWebTool_new(data,ThreshType,SampleNumber,ref,percentile_or_value,percentile,value,LocalThresholdType,percentile_low,percentile_high,value_low,value_high)