170 likes | 303 Views
Parameter Estimation & Model Fitting. More SAAMII Applications & Intro to DIMSUM. Experiment. Experimental point z is the sum y and z Variable y is the output of your model; deterministic Variable e is the error at any given point in time; stochastic. Parameters.
E N D
Parameter Estimation & Model Fitting More SAAMII Applications & Intro to DIMSUM
Experiment • Experimental point z is the sum y and z • Variable y is the output of your model; deterministic • Variable e is the error at any given point in time; stochastic
Parameters • Microparameters: kij’s, Vi’s so far have been given • Macroparameters: Ai’s, lambdai’s derived from our kij’s, Vi’s and knowing the structure of the model
DIMSUM • W3dimsum available online at http://biocyb.cs.ucla.edu • If working from home, must first d/l Java Start up and turn off proxy servers • Using a set of data (which one: z, y, or e?) finds the best macroparams (which are?) via WRSS (which is?)
Curve peeling v. Marquardt • Used to fit the model to data • Curve peeling: page 353 in your reader; basically strips away the individual exp curves composing the composite curve • Marquardt: requires initial guesses, more expensive in computation time but more robust.
WRSS as a determination of goodness of fit • Weighted Residual Sum of Squares • Minimizes residuals from the expected curve(s) (which is z,y,or e??)
Information Criteria • Akaike Info Criterion • Schwartz Criterion • Also…F-test…RMS….see pps 357+ AIC and SC are closely related but weight residuals differently.
Ready, set, go! • After starting W3DIMSUM, click on data and enter to start entering data values. • Enter number of data points and the t & z values. • Enter error model to be used….
HUH? WHAT’S AN ERROR MODEL • Given z’s, the error model assigns how much “credit” should be given to the z’s, I.e. how closely they should match the y • Also, determines how points are weighted – why would more allowance, for example, be given to earlier data points??
CV v. Var, Constant or not • CV: sd/mean • Var: sd squared • When would we choose to have error same for all data points? When no? • CV<10% is ok, <5% is even better.
Now run the algorithms • Go to Models on the top bar and select automated modeling. • The program will automatically find the best macroparameters for 1, 2, 3, and 4 exponential models using the Curve-Peeling algorithm, followed by the Marquardt algorithm.
Expert System • Expert system is like a little statistician in the program – analyzes AIC, etc and determines best fit • Notice that when you ran the error model you were given the Ai’s and lambdai’s.
SAAMII and parameter estimation • Given a compartmental layout and ranges for your microparameters, you can determine microparameters that best fit your model • Again, uses CV, var, etc to determine how closely the model should fit the data
How to do it… • Lay out model in compartments and fluxes. • Assign fixed or adjustable ranges (with an intial value) for all microparameters, including volume. • Open data window…
Entering data • First line: DATA • Second line: error model (FSD f) for fractional standard deviation – 5% cv is equivalent to f=.05 • Third line: headers of columns: t, smple, etc • …data separated by spaces & enters • Last line: END • Close the window
Fitting • Associate: go to show, then choose associations. Choose to associate s1, etc with whichever column you want. • Fit: click the red bar, and then the blue “fit” icon on the taskbar. • Clicking the values icon on gives your best fit parameters. Clicking on the sigma looking icon gives statistical breakdown.
One last note about fit • Given a set of errors, the line best fit through them should ideally look like what? ?