Dear nmusers: Recently I am trying to build a cell-cycle based tumor growth model, i.e. tumor size = G1+ S + G2 + M. Meanwhile, I have quite a few simultaneous biomarkers to characterize the cycle behaviors. Those biomarkers could be calculated based on the population of G1, S, G2, and M. I have the following puzzles needed to be addressed:
1) In my data file, what is the CMT number I should give to those observations (such as tumor size and many biomarkers)? By the way, I need CMT to specify/initialize the corresponding compartment. Since I am not quite sure how to use L2 and PCMT, I tried many times, all failed except the following foolish way: I create one more compartment to hold the tumor size: DADT(iTUMROR) = DADT(iG1) + DADT(iS) + DADT(iG2) + DADT(iM) and in my data file, I put the CMT compartment created (iTUMOR) into the observation of tumor size. However, for all other biomarkers, there is no way to write this kind of DADT equations. How could I solve this in a smart way? 2) Since some biomarker data are collected based on certain time point, saying t0. In other words, they are relative fold-changes compared to the vehicle at t0. Therefore, I need to save the intermediate population of G1, S, G2 and M at t0, as constants, which will be used to scale the later-on populations so that the model predictions are comparable with the observations. However, in the both block $DES and $ERROR, there is no way to save these intermediate populations as global variables or constants. I guess, I have to use MSFO to separate the simulation to many sections? But it am not sure how to do it. Your thoughts and feedback are really appreciated! Thanks a lot! Feng
