MIT OpenCourseWare http://ocw.mit.edu 12.740 PaleoceanographySpring 2008For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms.12.740 Paleoceanography Spring 2006 Problem 2 Factor analysis primer. The things you are expected to include in your answer are indicated numbered and in boldface, e.g. III.(2). Eigenvalue/eigenvector analysis primer. You are given email attached files with 100 x,y,z data points and certain files with which to examine them. The data will be explored in part using Viz!on within Microsoft Excel 98; other aspects can be achieved using any recent version of Excel or other mathematical packages. It is assumed that you are minimally proficient at basic Macintosh operations (pointing, clicking, double-clicking, dragging, and quitting programs) or can find someone to tutor you in these. You can do the exercise on any computer that has these applications; however, Viz!on will not work on PCs. The files are: Double-clicking on any of these files will open them up. 1. On the desktop of the Macintosh 7500, double click on the Problem Set 2 icon on the desktop. Double click on the “Excel98” icon to launch Excel. After Excel opens, close the blank spreadsheet and double click on the Vis!on” icon to add Vis!on into Excel. 2. Open the file Data.txt. Highlight the three columns (ABC). Click on Viz!on in the menu bar and select Vizualize 3 Columns.3. Click on “Click Me” and select 3Quant. Select the x,y, and z axes by dragging them in turn from their listing on the right to their appropriate box on the left: (5) Click on "View the data in a rotating plot”. Click on “Add Color by: x” (this makes the plot a little easier to follow when it is spinning) (6) You can expand the screen either by clicking on the expand bar or using the resizing handle on the lower right. When you expand at first, the plot remains the same size. You can expand or contract the size of the plot by using the tool shown below. This tool works by pointing and clicking near the outer edge if you want to expand, and near the center if you want to contract.also click and drag without letting go to move it as you please. II.(1)Examine the data set and describe the characteristics you notice. 2. Open the file Eigen-analysis (XL). Near the top of the screen you will see the correlation matrix: A x y z x 1 0.93 0.86 y 0.93 1 0.88 z 0.86 0.88 1 Immediately below that matrix is a section where guesses for the eigenvalues can be entered (bold font), where the determinant of the matrix A-λI is calculated ( font), and where the pre-computed eigenvectors are shown (sshhaaddooww font). You can change the guesses for the eigenvalues until you arrive at a solution (). There are three solutions (roots) to the characteristic cubic polynomial that arises from the A-λI determinant. You can arrive at your guesses to the roots by examining the graph that is located below the eigenvector section: Characteristic Polynomial Equationλ-0.05-0.04-0.03-0.02-0.0100.010.020.030.040.050 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 2.2 2.4 2.6 2.8 3 If you would like to see this equation plotted on a different scale, you may change the scale axes of this graph by double-clicking on an axis and selecting the scale option, and entering the new values you wish for the scale. By interacting with this graph and the calculation section above, II.(2) find the eigenvalues of the data set (to three decimal places) and describe how you found them(*). Notice that when you find an eigenvalue, the numbers to the right of the pre-computed eigenvector x (which are the elements of the column vector (A-λI)xx ) approach zero [note: this assumes that you put the smallest eigenvalue in the top group and the largest in the bottom group].. What is the sum of the eigenvalues? This result is not a coincidence - the sum of the eigenvalues corresponds to the rank of the matrix (the minimum number of orthogonal vectors required to span the sample space). II.(3) Divide the eigenvalues (7) You start the plot spinning by using the (arrow) hand tool to click and drag (letting go) to give the plot a spin. You canby the sum of the eigenvalues. These numbers are equal to the fraction of the total variance "accounted" for by the eigenvector associated with that eigenvalue. II.(4) Look at the eigenvectors and describe the spatial location of each (i.e., where in the four x-y quadrants and what angle from the x-y plane (1° accuracy); e.g.: "x positive, y negative, 30° above the horizontal"). (* )[Note: if you are familiar with the "Solver" in Microsoft Excel, you can use it to find the roots - but note that this Newton-Rapson solver requires that your initial guess be close to the correct answer). Close this file. 3. Double click-on DataVec.txt. Open up the columns using Viz!on. This file illustrates the raw data set and normalized eigenvectors from step 2. Notice that the eigenvectors are orthogonal and follow the principle axes of elongation of the data set. Which eigenvalues are associated with each of these eigenvectors?III. Factor analysis example. You are given a data set consisting of a subset of the Atlantic core top foram census (50 samples selected at random with data for G. ruber, G. sacculifer, G. bulloides, N. pachyderma (R), N. pachyderma (L). You will follow the Q-mode factor analysis of this data set into three factors, and verify its correctness. 1. The raw data (# of specimens of each species in a sample is given at the top of the spreadsheet. The sum of squares (for unit-length normalization) is calculated. You are also given estimates for sea surface temperature. Immediately below the raw data data is the row-normalized data, W: The analysis then proceeds by calculation of the 5 x 5 minor Product Moment (mPM). Actually, you could also use the Major Product Moment (50 x 50), but the calculations are obviously easier for the mPM. Note: If you move the cursor to any cell of a matrix, it will show the matrix formula used to generate the matrix. We then calculate the eigenvalues and eigenvectors. Excel will not do this calculation, so the mPM was transferred from Excel into Mathematica, which gives eigenvalues and eigenvectors as a result of the commands:etc.. The results were then copied back into Excel: rub III.(1) Verify that these are the eigenvalues and eigenvectors of the mPM (i.e., demonstrate orthogonality of (A-λ) and I