SOCR EduMaterials Activities ExpCarTimeExperiment
Contents
Exponential Time Car Activity
Description
TBD.
Goal
To provide an interactive simulation demonstrating exponential processes, exponential decay, simulation and model fitting. Additional aims are to:
- Draw analogy between Gamma, Exponential and Poisson (discrete!)
- Provide the means for (large) data generation/sampling
- Provide the infrastructure for model-fitting and assessment.
Experiment
Go to the SOCR Experiments and select the Exponential Car-Time Experiment from the drop-down list of experiments on the top left. The image below shows the initial view of this experiment with pointers and highlights of the most important components and controls:
Select the number of cars (k) and the rate of the exponential process (r). When pressing the play button, one trial will be executed that involves sampling k random exponential time variables (representing the between car times). The results will be displayed graphically in the two graph windows (timeline and histogram) for this process. Summary statistics and model vs. data distributions are shown in the distribution table below. The stop button ceases any activity and is helpful when the experimenter chooses a large number of cars (k is large). The reset button will clear and re-initialize the entire experiment, deleting all previous information and data collected.
As the number of cars (k) increases, the distribution of the between-car times will have exponential-decay shape. Note that as the number of cars (k) increases, the empirical histogram (graph) and moments (table) of the random exponential times begin to converge to their corresponding Exponential distribution counterparts. Arrows point to the match of the empirical histogram of the between-car times and theoretical Exponential distribution, the random sample of between-car times (Y column of left table) and the similarity between the first two moments of the sample data and model (Exponential) distributions.
Applications
Relations between Gamma, Exponential and Poisson distributions
Number of Arrivals in a fixed time-period follows Poisson distribution.
- \({\Gamma}(k=1, \theta=1/\lambda)\,\) reduces to an Exponential distribution with rate parameter \(\lambda\).
- For an integer \(k\) parameter, the Gamma distribution \({\Gamma}(k=1, \theta=1/\lambda)\,\) becomes the probability distribution of the waiting time until the \(k^{th}\) arrival of a one-dimensional Poisson process with expected number of arrivals in a fixed time interval equal to \(\lambda = {1\over \theta}\).
- If one observes a Poisson (Number of arrivals in a fixed time-interval), then the Exponential distribution describes the lengths of the inter-arrival times in the homogeneous Poisson processes.
- The Exponential distribution may also be considered as a continuous analog of the Geometric distribution. Geometric distribution describes the number of Bernoulli trials until the first change (success). Similarly, Exponential distribution describes the time for a continuous process to change state and is sometimes called the Mean-Time-to-Failure distribution.
Exponential data sampling
We can use this applet to geenrate a sample of an arbitrally size from an Exponential Distribution with specified rate parameter. Then these data may be used for simulations or exploratory data analsys. For example, suppose we ran the experiment with k=1,000, as shown onthe image above. We can copy the sample from the Y column and past it in any of the SOCR Charts. The image below shows the construction of a Statistical Bar Chart using the 1,000 Exponential (between-car) time measurements. Notice how the sumamry statistics on this Bar Chart coinside with the data-statistics reported on the distribution table of the Exponential Time Car Experiment applet. The last graph (bottom-right) illustrates a (Q-Q) plot of the Normal vs. Exponential-sample quantiles. The highly non-linear shape of this QQ plot indicates the significant departure of normality of these Exponential time measurements (between-car times).
Exponential model-fitting and assessment
The image below demonstrates how we can fit a Gamma-model to the sample of (k=1,000) Exponential measurements that we simulated as between-car times. Go to SOCR Modeler and click the Raw Data check-box. Then paste in all of your raw Exponential measurements recorder in the Y collumn of the Exponential Time Car Experiment. Choose Gamma-Fit, from the drop-down list of SOCR modelers on the top-left. Click on the Estimate parameters and Scale up check-boxes. Then go to the Graph tab and use the sliders on the top to select an approptiate view for the data and the model-fit. Under the Results tab you will find the Gamma model parameter estimates (these are MLE estimates). In this case, the best Gamma distribution model is estimated to have parameters Shape = 0.8942640916359458; Scale = 1.0443144529333228.
Translate this page:
