320 likes | 476 Views
SIGMA Workshop Part 3: Statistical Screening. Gönenç Yücel SESDYN Research Group Boğaziçi University, Istanbul. A brief introduction. What is stat screening? What is it good for? What does it rely on? Altering uncertain (exogenous ) model parameters
E N D
SIGMA WorkshopPart 3: Statistical Screening Gönenç Yücel SESDYN Research Group Boğaziçi University, Istanbul
A brief introduction • What is stat screening? • What is it good for? • What does it rely on? • Altering uncertain (exogenous) model parameters • Relating the value of an outcome of interest to changes in parameter values • Uses correlation coefficient to quantify the degree and direction of relationship between a parameter and the value of the outcome at a certain time point
Background literature • Ford, A., & Flynn, H. (2005). Statistical screening of system dynamics models. System Dynamics Review, 21(4), 273–303. • Taylor, T., Ford, D., & Ford, A. (2007). Model Analysis Using Statistical Screening: Extensions and Example Applications . 25th International Conference of the System Dynamics Society. Boston: System Dynamics Society. • Taylor, T. R. B., David N. Ford, & Ford, A. (2010). Improving model understanding using statistical screening. System Dynamics Review, 26(1), 73–87.
Key Concept: Correlation Coefficients • A measure of the linear correlation (dependence) between two variables X and Y, giving a value between +1 and −1 inclusive, where • 1 is total positive correlation, • 0 is no correlation, and • −1 is total negative correlation.
Demo Model • Bass diffusion model • See Business Dynamics by Sterman (2000) for specifications of the model • Vensim version of the model is available in the Stat Screening folder on your computers
Our demo task • Perform a Statistical Screening Analysis on the demo model (i.e. Bass diffusion model) to evaluate the relative influence of • Two of the exogenous variables (i.e. contact rate and adoption fraction) and • The one exogenous initial condition (i.e. initial potential adopters). • Using the adoption rate (sales) as the main outcome of interest (performance variable).
Procedure • Perform statistical screening to calculate correlation coefficients and to plot these over time • Select a time period for analysis • Identify high-leverage parameters. • High-leverage parameters are the parameters with the highest absolute correlation coefficient values during the selected time • Create a list of high leverage parameters and their related model structures • Use additional structure-behavior analysis methods (e.g. verbal reasoning, scenario analysis, behavioral analysis) to explain how each parameter the structures they influence drive the behavior of the system.
A. Calculating Correlation Coefficients • Select uncertain model input parameters and a single performance variable for analysis • Specify a distribution for each uncertain model parameter • Simulate the model using a combination of values from the specified distributions (e.g. Using Vensim’s Sensitivity Analysis feature) • Export the results of the simulation set • Pick up the Excel template that best fits the simulation set • Import the results from the simulation set to the Excel template, and observe the plot of the correlation coefficients
Steps A.1 & A.2 • Selected parameters, and distributions • Model parameters to be analyzed • Model output to be analyzed • Adoption rate
A. Calculating Correlation Coefficients • Select uncertain model input parameters and a single performance variable for analysis • Specify a distribution for each uncertain model parameter • Simulate the model using a combination of values from the specified distributions (e.g. Using Vensim’s Sensitivity Analysis feature) • Export the results of the simulation set • Pick up the Excel template that best fits the simulation set • Import the results from the simulation set to the Excel template, and observe the plot of the correlation coefficients
Conducting a set of simulations on Vensim • Monte Carlo option in Vensim • Lets us to specify ranges for the parameters as well as their distribution • We will need to specify 2 things • An input control file (.vsc file) • An output savelist file (.lst file)
A. Calculating Correlation Coefficients • Select uncertain model input parameters and a single performance variable for analysis • Specify a distribution for each uncertain model parameter • Simulate the model using a combination of values from the specified distributions (e.g. Using Vensim’s Sensitivity Analysis feature) • Export the results of the simulation set • Pick up the Excel template that best fits the simulation set • Import the results from the simulation set to the Excel template, and observe the plot of the correlation coefficients
A. Calculating Correlation Coefficients • Select uncertain model input parameters and a single performance variable for analysis • Specify a distribution for each uncertain model parameter • Simulate the model using a combination of values from the specified distributions (e.g. Using Vensim’s Sensitivity Analysis feature) • Export the results of the simulation set • Pick up the Excel template that best fits the simulation set • Import the results from the simulation set to the Excel template, and observe the plot of the correlation coefficients
A.5 Choosing an Excel template • Choosing the right template! • Number of parameters • Number of simulations • Number of data points in a single run • In our example, we have • 3 parameters • 200 simulations • 100 data points for each simulation • The right template would be StatScreenTemplate-3inputs200runs100saveperiods.xls
A. Calculating Correlation Coefficients • Select uncertain model input parameters and a single performance variable for analysis • Specify a distribution for each uncertain model parameter • Simulate the model using a combination of values from the specified distributions (e.g. Using Vensim’s Sensitivity Analysis feature) • Export the results of the simulation set • Pick up the Excel template that best fits the simulation set • Import the results from the simulation set to the Excel template, and observe the plot of the correlation coefficients
Self-study Practice • Repeat the statistical screening on a modified version of the simple Bass diffusion model • Modification: • Add a quitting flow that flows from the adopters to the potential adopters stock • The amount of the flow is defined as • Adopters * Quitting Fraction • Reference value of the quitting fraction is set to be 0.1