Confidence Region (Confidence ellipse) plot for Bivariate Data

Part 1: Calculate Rotation Angle.
  1. Calculate Chi-square critical value based on required confidence (eg. alpha of 0.05, is 95% confidence ellipse). 
    • Using CHISQ.INV() function in Excel. 2 argument in function is degree of freedom, which is number of variables measured as part of sample data, for bivariate data dof=2.
  2. Calculate Variance Covariance matrix of sample data. 
    • Using COV() function of real stats.
  3. Calculate eigen value of covariance matrix. 
    • Using eVALUES() function of real stats or using formula considering VarCov matrix.
    • sqrt(λ1) is the radius of the major axis (the longer radius) 
    • sqrt(λ1) is the radius of the major axis (the longer radius) 
  4. Calculate angle from positive X-axis to the ellipse's major axis, in counterclockwise direction
    • Using ATAN2() excel function.
    • atan2()'s first parameter is y and second is x. atan2(covar, λ1-Var.X). 
    • θ is the angle in radian.
Part 2: Get the points on Confidence Ellipse
  1. Calculate rotation matrix, Q, using angle in radian.
    • Diagonal elements are, Cos(θ)
    • Off diagonal elements are, Sin(θ) and -Sin(θ)
  2. Create an input column for angle with values 0 to 2π in increments of π/10.
    • First cell, angle=0.
    • 2nd cell, angel =PI()/10
    • 3rd onward, cumulative sum, get values from 0 to 6.28 (21 values)
  3. Calculate the lengths of the two axes of the ellipse.
    • Length of major axis, a, a=sqrt(λ1)*sqrt(Chi-square Critical value).
    • Length of minor axis, b, b=sqrt(λ2)*sqrt(Chi-square Critical value).
  4. Calculate the points on confidence ellipse using formula.
    • X-value is, Xmean+(Cos(angle)*dig.of Q*mj_axis)+sin(angle)*offdig.neg.value*mn_axis)
    • Y-value is, Ymean+(Cos(angle)*offdigQ*mj_axis)+sin(angle)*dig.value*mn_axis)
    • Using Matrix operations, to be demonstrated.
  5. Get pair of points on confidence ellipse.
Part 3: Plot the graphs of confidence ellipse and include individual data and its mean in plot.
Part 4: Determine if point is inside or outside ellipse using inverse of covariance matrix.

Reference

Comments

Post a Comment

Popular posts from this blog

Covariance Error Ellipse - sqrt of Eigen value as ellipse axis length

Image
  Axis-aligned confidence ellipses In general, the equation of an axis-aligned ellipse with a major axis of length   and a minor axis of length  , centered at the origin, is defined by the following equation: In our case, the length of the axes are defined by the standard deviations   and   of the data such that the equation of the error ellipse becomes: where s defines the scale of the ellipse and could be any arbitrary number (e.g. s=1).  The question is now how to choose  , such that the scale of the resulting ellipse represents a chosen confidence level (e.g. a 95% confidence level corresponds to s=5.991). The left hand side of equation ( 2 ) actually represents the sum of squares of independent normally distributed data samples, cov is zero.  The sum of squared Gaussian data points is known to be distributed according to a so called  Chi-Square distribution . A Chi-Square distribution is defined in terms of ‘degrees of fre...
Read more

Visualizing Ellipse rotation and associated parametric equations

Image
  Diagonal matrix Ellipse axes are parallel to coordinate axes (no rotation). Zeros in off-diagonals means zero correlation. Same values in diagonals Ellipse is rotated 45 degrees if correlation (off-diagonal) is positive regardless of its magnitude.  Similarly -45 degrees if correlation is negative. Correlation matrix Correlation matrix has 1 in diagonals and values between -1 and 1 inclusive in off-diagonals. Ellipse size remains the same (always touches square of side 2 units).  Whether the correlation is positive or negative can be observed by the orientation of the ellipse. The amount of correlation can be interpreted by how thin the ellipse is. How It Works Below is a list of parametric equations starting from that of a general ellipse and modifying it step by step into a prediction ellipse, showing how different parts contribute at each step. Reference How to Draw Ellipse of Covariance Matrix (cookierobotics.com)
Read more