CIS 2033 Base on text book: A Modern Introduction to Probability and Statistics Understanding Why and How By: F.M. Dekking, C. Kraaikamp, H.P.Lopulaa, L.E.Meester Temple University Spring 2012 Slides by: Wanwisa Smith Modified by: Dr. Longin Jan Latecki
Chapter 15 Exploratory data analysis: graphical summaries The set of observations is called a dataset. By exploring the dataset we can gain insight into what probability model suits the phenomenon. To graphically represent univariate datasets, consisting of repeated measurements of one particular quantity, we discuss the classical histogram, the more recently introduced kernel density estimates and the empirical distribution function. To represent a bivariate dataset, which consists of repeated measurements of two quantities, we use the scatterplot.
15.2 Histograms: The term histogram appears to have been used first by Karl Pearson.
How to construct the histogram? Denote a generic (univariate) dataset of size n by First we divide the range of the data into intervals. These intervals are called bins and denoted by The length of an interval Bi is denoted by ǀBiǀ and is called the bin width. We want the area under the histogram on each bin Bi to reflect the number of elements in Bi. Since the total area 1 under the histogram then corresponds to the total number of elements n in the dataset, the area under the histogram on a bin Bi is equal to the proportion of elements in Bi: The height of the histogram on bin Bi must be equal to
Choice of the bin width Consider a histogram with bins of equal width. In that case the bins are of the from where r is some reference point smaller than the minimum of the dataset and b denotes the bin width. Mathematical research, however, has provided some guide- line for a data-based choice for b or m.
A kernel K is a function K:RR and a kernel K typically satisfies the following conditions.
Scaling the kernel K Then put a scaled kernel around each element xi in the dataset Scale the kernel K into the function
The bandwidth is too big The bandwidth is too small
The function g in blue is a mixture of two Gaussians. We draw 200 samples from it,which are shown as blue dots. We use the samples to generate the histogram (yellow) and its kernel density estimate f (red). The Matlab script is twoGaussKernelDensity1.m
15.4 The empirical distribution function Another way to graphically represent a dataset is to plot the data in a cumulative manner. This can be done by using the empirical cumulative distribution function .
15.5 Scatterplot In some situation we might wants to investigate the relationship between two or more variable. In the case of two variables x and y, the dataset consists of pairs of observations: We call such a dataset a bivariate dataset in contrast to the univariate. The plot the points (Xi, Yi) for i = 1, 2, …,n is called a scatterplot.