{smcl} {* *! version 1.5.0 04oct2012}{...} {hline} {cmd:help paran} {hline} {title:Title} {p 4 21 2} {cmd:paran} {hline 2} HornŐs Test of Principal Components/Factors {title:Syntax} {pstd} Parallel analysis of data {p 8 17 2} {cmd:paran} [{it:varlist}] [{it:weight}] [{cmd:if} {it:exp}] [{cmd:in} {it:range}] [{cmd:,} {it:options}] {synoptset 20 tabbed}{...} {synopthdr} {synoptline} {syntab :Model} {synopt :{opt ite:rations(#)}}specify the number of iterations{p_end} {synopt :{opt cent:ile(#)}}specify using a centile value of instead of the mean{p_end} {synopt :{opt factor(factor_type)}}use {help factor:factor} instead of {help pca:pca} (defaults to {cmd:pca} if left blank){p_end} {synopt :{opt cit:erate(#)}}communality re-estimation iterations ({help ipf:ipf} only){p_end} {syntab :Reporting} {synopt :{opt q:uietly}}suppresses {help pca:pca} or {help factor:factor} output{p_end} {synopt :{opt nost:atus}}suppresses the status indicator{p_end} {synopt :{opt all}}report all eigenvalues (default reports only those retained){p_end} {syntab :Graphing} {synopt :{opt gr:aph}}graphs unadjusted, adjusted, and random eigenvalues{p_end} {synopt :{opt color}}renders graph in color (default is black and white){p_end} {synopt :{opt lcolors(3 x rgb)}}specifies colors using three {help colorstyle:rgb triples} for observed, random, and adjusted eigenvalues (overrides the {cmd:color} option){p_end} {synopt :{opt saving(filename)}}saves graph as a {cmd:.gph} file{p_end} {synopt :{opt replace}}replaces an existing file when {cmd:saving()}{p_end} {syntab :Miscellaneous} {synopt:{opt pr:otect(#)}}perform {it:#} optimizations and report the best solution ({cmd:ml} only with {help factor:factor}){p_end} {synopt :{opt seed(#)}}seed the random number generator with the supplied integer{p_end} {synopt :{opt mat(matrix name)}}option to provide a correlation matrix instead of the {it:varlist}{p_end} {synopt :{opt n(#)}}specifies the required sample size when using the {cmd:mat()} option{p_end} {synopt :{opt copyleft}}displays the GPL license for {cmd:paran}{p_end}{synoptline} {p2colreset}{...} {p 4 4 2} {cmd:fweights}, and {cmd:aweights} are allowed when using {it:varlist}. See help {help weights:weights}. {title:Description} {p 4 4 2} {cmd:paran} is an implementation of HornŐs technique for evaluating the components or common factors retained in a principle component analysis (PCA) or a common factor analysis (FA). According to Horn, a common interpretation of non-correlated data is that they are perfectly non-collinear, and one would expect therefore to see eigenvalues equal to 1 in a PCA of such data (or equal to 0 in the case of a common factor analysis, as with {cmd: pf}). However, Horn notes that multi-colinearity occurs due to sampling error and least-squares "bias," even in uncorrelated data, and therefore actual PCAs of such data will reveal eigenvalues of components greater than and less than 1. His strategy is to contrast eigenvalues produced through a PCA on a random dataset (uncorrelated variables) with the same number of variables and observations as the experimental or observational dataset to produce eigenvalues for components or factors that are adjusted for the sample error-induced inflation. Values greater than zero are retained in the adjustment given by: {pstd} For principal component analysis: {p 8 8} Observed Data Eigenvalue_p - (Random Data Eigenvalue_p - 1) {pstd} For common factor analysis: {p 8 8} Observed Data Eigenvalue_p - Random Data Eigenvalue_p {p 4 4 2} {cmd:paran} is used in place of a {cmd:pca} {it:varlist} command (or {cmd:factor}). The user may also specify how many times to make the contrast with a random dataset (default is 30 per variable). Values less than 1 will be ignored (or less than 0 for {cmd:factor}), and the default value assumed. Random datasets are generated using the {cmd:uniform()} function. The program returns the estimated mean eigenvalues of random data if the {cmd:centile} option is unspecified, otherwise it returns the specified {cmd:centile}. biases on each eigenvector are also returned. {cmd:paran} may be thus be used to conduct parallel analysis following Glorfeld's suggestions to reduce the likelihood of over-retention. (Glorfeld, 1995) {p 4 4 2} When the {cmd:all} option is not used, {it:unadjusted} eigenvalues greater than 1 (for prinicpal components) or 0 (for factors) are reported, with retained adjusted eigenvalues printed in yellow, and unretained adjusted eigenvalues printed in red. {title:Options} {phang} {cmdab:iter:ation(}{it:#}{cmd:)} sets the number of contrast datasets to evaluate. The default value is 30 * the number of variables, and values less than 1 are ignored. For large datasets with large numbers of variables many iterations may be time consuming. The greater the number of iterations the more accurate the estimates of sample bias will be. {phang} {cmd:centile(}{it:#}{cmd:)} specifies that supplied centile value is to be used instead of the mean (assumed median, since the distribution is symmetrical) in estimating bias. Values above the mean/median, such as the 95th percentile, are more conservative estimates of chance bias in the eigenvalues from a PCA of sample data. This option supercedes the older {cmd:pnf}, which was equivalent to {cmdab:cent:ile(}{inp:95}{cmd:)}. Values of {cmdab:cent:ile()} must be greater than 0 and less than 100. Non-integer values will be rounded to the nearest integer value. Running {cmd:paran} without this option uses the mean value (very close to {cmdab:cent:ile(}{inp:50}{cmd:)}). (see Glorfeld, 1995) {phang} {cmdab:q:uietly} suppresses output for the PCA or factor analysis. This option is only used if a {it:varlist} is specified in the {cmd:paran} command. {phang} {cmdab:nost:atus} By default {cmd:paran} indicates when every tenth percent of the computation has been completed. {cmd:nostatus} eliminates this behavior. {phang} {cmd:factor(}{it:factor_type}{cmd:)} selects one of the factor estimation types: {cmd:pf}, {cmd:pcf}, {cmd:ipf}, or {cmd:ml} (for principal factors, principal component factors, iterated principal factors, or maximum likelihood factors, respectively). If you specify anything but one of these four abbreviations, you will be warned and the program will halt. CAVEAT: Conducting parallel analysis using factor methods other than {cmd:pf} is unorthodox. Interpret such results at your own risk. If {cmd:factor} is not used, or if one of the factor estimation types is not used {cmd:paran} performs parallel analysis using {help pca:pca} by default. {phang} {cmdab:cit:erate(}{it:#}{cmd:)} sets how many iterations will be used to re-estimate communalities for the iterated principal factor type. (see {help factor:factor}) {phang} {cmdab:pr:otect(}{it:#}{cmd:)} sets the number of optimizations for starting values option for the maximum likelihood factor type. (see {help factor:factor}) {phang} {cmd:all} reports all components or factors, not just those with unadjusted eigenvalues greater than one (or greater than zero for factor). The default is not to report all components or factors. {phang} {cmdab:gr:aph} draws a graph of the observed eigenvalues, the random eigenvalues, and the adjusted eigenvalues much like the graphs presented by Horn in his 1965 paper. {phang} {cmd:color} renders the graph in color (only with {cmdab:gr:aph}) with unadjusted eigenvalues drawn in red, adjusted eigenvalues drawn in black, and random eigenvalues drawn in blue, and all lines drawn solid. Without the {cmd:color} option, the graph is rendered in black and white, and the line connecting the unadjusted eigenvalues is dashed, the line connecting the random eigenvalues is dotted, and the line connecting the adjusted eigenvalues is solid. {phang} {cmd:lcolors(}{it:# # # # # # # # #}{cmd:)} specifies the colors of each line on the graph using three {help colorstyle:rgb triples} (only with {cmdab:gr:aph}). The first triple indicates the {cmd:R}, {cmd:G} and {cmd:B} components of the observed eigenvalues, the second triple sets the values for the mean or centile random eigenvalues, and the third triple sets the values for the adjusted eigenvalues. These settings override the default (red, blue, and black) colors of the {cmd:color} option. {phang} {cmd:saving(}{it:filename}{cmd:)} outputs the graph to the specified {it:filename} as a {cmd:.gph} file (only with {cmdab:gr:aph}). {phang} {cmd:replace} overwrites an existing {it:filename} when the {cmd:saving()} option is used with {cmdab:gr:aph}. {phang} {cmd:seed(}{it:#}{cmd:)} specifies an integer seed for the random number generator (see {help generate:set seed}) so that results of {cmd:paran} upon a specific data set can be exactly reproduced. The default behavior of {cmd:paran} is not to specify a seed. {phang} {cmd:mat(}{it:matrix name}{cmd:)} specifies an optional correlation matrix to be used instead of the {it:varlist}; requires the {cmd:n(}{it:#}{cmd:)} option also be specified. This option is not compatible with {cmd:aweights} or {cmd:fweights}. {phang} {cmd:n(}{it:#}{cmd:)} specifies the sample size when using the {cmd:mat(}{it:matrix name}{cmd:)} option. {phang} {cmd:copyleft} displays the copying permission statement for {cmd:paran}. {cmd:paran} is free software, licensed under the GPL. The full license can be obtained by typing: {p 12 8 2} {inp: . net describe paran, from (http://www.doyenne.com/stata)} {phang} and clicking on the {net "describe paran, from (http://www.doyenne.com/stata)":click here to get} link for the ancillary file. {title:Remarks} {p 4 4 2} Hayton, et al. (2004) urge a parameterization of the random data to approximate the distribution of the observed data with respect to the middle ("mid-point") and the observed min and max. However, the PCA as I understand it is insensitive to standardizing transformations of each variable, and any linear transformation of all variables, and produces the same eigenvalues used in component or factor retention decisions. This is born by the notable lack of difference between analyses conducted using a variety of simulated distributional assumptions (Dinno, 2009). The central limit theorem would seem to make the selection of a distributional form for the random data moot with any sizable number of iterations. Former functionality implementing the recommendation by Hayton et al. (2004) has been removed, since parallel analysis is insensitive to it, and it only adds to the computation time required to conduct parallel analysis. {title:Examples} {p 4 8}{inp:. paran var1-var16} {p 4 8}{inp:. paran var1-var26, iter(5000) q centile(95)} {p 4 8}{inp:. paran var1-var10, iter(1) factor(ipf) cit(50)} {title:Author} {p 4 8} Alexis Dinno{p_end} {p 4 8} alexis dot dinno at pdx dot edu{p_end} {p 4 8} I am receptive to comments and requests.{p_end} {title:Reference} {phang} Dinno A. 2009 "Exploring the Sensitivity of HornŐs Parallel Analysis to the Distributional Form of Simulated Data" {it:Multivariate Behavioral Research}. 44: 362-388 {phang} Glorfeld LW. 1995. "An Improvement on HornŐs Parallel Analysis Methodology for Selecting the Correct Number of Factors to Retain. {it:Educational and Psychological Measurement}. 55: 377-393 {phang} Hayton JC, Allen DG, and Scarpello V. 2004. "Factor Retention Decisions in Exploratory Factor Analysis: A Tutorial on Parallel Analysis" {it:Organizational Research Methods}. 7: 191-205 {phang} Horn JL. 1965. "A rationale and a test for the number of factors in factor analysis." {it:Psychometrika}. 30: 179-185 {phang} Zwick WR, Velicer WF. 1986. "Comparison of Five Rules for Determining the Number of Components to Retain." {it:Psychological Bulletin}. 99: 432-442 {title:Saved results} {cmd:paran} saves the following 1 by P matrices in {cmd:e()}: {synoptset 20 tabbed}{...} {p2col 5 20 24 2: Matrices}{p_end} {synopt:{cmd:e(UnadjustedEv)}}Unadjusted eigenvalues from the {cmd:pca} or {cmd:factor} command{p_end} {synopt:{cmd:e(AdjustedEv)}}Eigenvalues from the analysis adjusted by subtracting the estimated bias{p_end} {synopt:{cmd:e(MeanRandomEv)}}The mean of the eigenvalues of random data sets of size N by P (only if the {cmd:centile()} option is unspecified){p_end} {synopt:{cmd:e(CentRandomEv)}}The centile of the eigenvalues of random data sets of size N by P as given by the {cmd:centile()} option (and only if that option is specified){p_end} {synopt:{cmd:e(Bias)}}The estimated bias (which {it:is} the mean of the eigenvalues of random data sets of size N by P when using the {cmd:factor} option{p_end} {p2colreset}{...} {title:Also See} {p 4 8} On-line:{space 2}help for: {help pca:pca}, {help factor:factor}