add check for log space operation

MunskyGroup · Jul 11, 2024 · 0ccf377 · 0ccf377
1 parent 92d2528
commit 0ccf377
Showing 1 changed file with 8 additions and 2 deletions.
diff --git a/WorkSpace/AlexP/SSIT_BurstingGene_MCMC.m b/WorkSpace/AlexP/SSIT_BurstingGene_MCMC.m
@@ -13,7 +13,7 @@
 
     %% Solve using FSP
     F1 = F1.formPropensitiesGeneral('BurstingGene');
-    F1.tSpan = [1:1:20];
+    F1.tSpan = [1:0.1:5];
     F1.initialTime = 0;
     F1.solutionScheme = 'FSP';    % Set solutions scheme to FSP.
     [FSPsoln, F1.fspOptions.bounds] = F1.solve;  % Solve the FSP analysis
@@ -152,9 +152,15 @@ function plotMHTraces(samples, paramNames)
 
 % Get computation time for MCMC
 tic
-MHResults = performMLEandMCMC(F1,true,0.1,3000);
+MHResults = performMLEandMCMC(F1,false,0.1,1500);
 toc
 
+tic
+logMHResults = performMLEandMCMC(F1,true,0.1,1500);
+toc
+
+transformed_logMHSamples = exp(logMHResults.mhSamples);
+
 % Assume MHResults.mhSamples is a 1000x4 matrix
 samples = MHResults.mhSamples;