Skip to content

Commit

Permalink
Update workspace data.
Browse files Browse the repository at this point in the history
  • Loading branch information
Munsky committed Jul 24, 2024
1 parent bbfc984 commit cd18373
Show file tree
Hide file tree
Showing 2 changed files with 13 additions and 6 deletions.
19 changes: 13 additions & 6 deletions WorkSpace/EricModel/RonModelPipelineJuly11_2024.m
Original file line number Diff line number Diff line change
Expand Up @@ -11,7 +11,7 @@
addpath(genpath('../../src'));

loadPrevious = true;
savedWorkspace = 'workspaceJuly22';
savedWorkspace = 'workspaceJuly24';
addpath('tmpPropensityFunctions');

%% STEP 0 -- Preliminaries.
Expand Down Expand Up @@ -836,8 +836,14 @@
logPriorAll = @(x)-sum((log10(x)-log10PriorMean([1:12,14:17])).^2./(2*log10PriorStd([1:12,14:17]).^2));

objTPL = @(x)-getTotalFitErr(OrganizationTPL,exp(x),false,extraObjs)-logPriorAll(exp(x));
parsTPL=parsAllandTS;
parsTPL(16) = log(2)/5;

if loadPrevious
varNamesTPL = {'parsTPL'};
load(savedWorkspace,varNamesTPL{:})
else
parsTPL=parsAllandTS;
parsTPL(16) = log(2)/5;
end

% Check that the function works.
tic
Expand All @@ -861,8 +867,8 @@
%%
%% Save Results for Easier Use in subsequent runs.
parsAll_GR_Dusp1_TS = [extendedMod.parameters{:,2}];
parsAll_GR_Dusp1_TS(indsODEmod) = parsAllandTS;
varNames = unique({'ModelGR',
parsAll_GR_Dusp1_TS(16) = parsAllandTS(end);
varNames = unique({'ModelGR'
'GRfitCases'
'log10PriorMean'
'log10PriorStd'
Expand All @@ -882,9 +888,10 @@
'ModelGRDusp100nM_ext_red'
'fullPars'
'parsAll_GR_Dusp1_TS'
'parsTPL'
});

save('workspaceJuly22',varNames{:})
save('workspaceJuly24',varNames{:})

%%
%% Additonal Codes for FIM and PDO analyses (ALEX Paper)
Expand Down
Binary file added WorkSpace/EricModel/workspaceJuly24.mat
Binary file not shown.

0 comments on commit cd18373

Please sign in to comment.