similar to: Note on PCA (not directly with R)

Hi all, I am using the vegan package to run a prcincipal components analysis on forest structural variables (tree density, basal area, average height, regeneration density) in R. However, I could not find out how to extract factor loadings (correlations of each variable with each pca axis), as is straightforwar in princomp. Do anyone know how to do that? Moreover, do anyone knows

Dear R-list, Among the (R-)tools, I've seen on the net, for (bivariate) Principal Component scatter plots (+histograms), "plotpc" [1] is the one I like most. By default it performs PCA on a bivariate dataset based on R's "princomp()" (which is the eigenvector-based algebraic solution to PCA). I would like to modify "plotpc()" in order be able, as an

Dear all, I have a question about the 'sdev' value returned by the princomp function (which does principal components analysis). On the help page for princomp it says 'sdev' is 'the standard deviations of the principal components'. However, when I calculate the principal components for the USArrests data set, I don't find this to be the case: Here is how I

After doing a PCA using princomp, how do you view how much each component contributes to variance in the dataset. I'm still quite new to the theory of PCA - I have a little idea about eigenvectors and eigenvalues (these determine the variance explained?). Are the eigenvalues related to loadings in R? Thanks, Paul -- View this message in context:

Dear expe-R-ts, I try some test data for a factorAnalysis (resp. pca) in the sense of Prof. Ripley's MASS ? 11.1, p. 330 ff., just to prepare myself for an analysis of my own empirical data using R (instead of SPSS). 1. the data. ## The test data is (from the book of Backhaus et al.: Multivariate ## Analysemethoden. Springer 2000 [9th ed.], p. 300 ff):

Dear All, I want to know if there is some easy and reliable way to estimate the number of dominant eigenvalues when applying PCA on sample covariance matrix. Assume x-axis is the number of eigenvalues (1, 2, ....,n), and y-axis is the corresponding eigenvalues (a1,a2,..., an) arranged in desceding order. So this x-y plot will be a decreasing curve. Someone mentioned using the elbow (knee)

Hey, All In principal component analysis (PCA), we want to know how many percentage the first principal component explain the total variances among the data. Assume the data matrix X is zero-meaned, and I used the following procedures: C = covriance(X) %% calculate the covariance matrix; [EVector,EValues]=eig(C) %% L = diag(EValues) %%L is a column vector with eigenvalues as the elements percent

Hello All, The basic premise of what I want to do is the following: I have 20 "entities" for which I have ~500 measurements each. So, I have a matrix of 20 rows by ~500 columns. The 20 entities fall into two classes: "good" and "bad." I eventually would like to derive a model that would then be able to classify new entities as being in "good

Hi, I want to include a vector representing the sites - environmental data correlation in a PCA. I currently use prcomp (no scaling) to perform the PCA, and envfit to retrieve the coordinates of the environmental data vector. However, the vector length is different from the one obtained in CAnoco when performing a species - environmental biplot (scaling -2). How can I scale the vector in order to

Hi All, would appreciate an answer on this if you have a moment; Is there a function (before I try and write it !) that allows the input of a covariance or correlation matrix to calculate PCA, rather than the actual data as in princomp() Regards Glenn [[alternative HTML version deleted]]

Dear all: When I have more variables than units, say a 195*10896 matrix which has 10896 variables and 195 samples. prcomp will give only 195 principal components. I checked in the help, but there is no explanation that why this happen. Can we get more than 195 PCs for this case? Thank you very much. Best! Alan Aug-12-2005

Hi all, I am attempting to run PCA on a matrix (nrow=66, ncol=84) using 'prcomp' (stats package). My data (referred to as 'Q' in the code below) are separate river streamflow gaging stations (columns) and peak instantaneous discharge (rows). I am attempting to use PCA to identify regions of that vary together. I am entering the following command:

Hi all, I'm trying to do model reduction for logistic regression. I have 13 predictor (4 continuous variables and 9 binary variables). Using subject matter knowledge, I selected 4 important variables. Regarding the rest 9 variables, I tried to perform data reduction by principal component analysis (PCA). However, 8 of 9 variables were binary and only one continuous. I transformed the data by

Dear all, I have read various descriptions of employing resampling techniques, such as the bootstrap, to estimate the uncertainties of the eigenvectors computed by PCA. When I try

Hi All, I have been trying to get to plot my PCA legend outside of the PCA plot, but success still alludes me. Can you guys please advise how I can achieve this. I used locater() to obtain coordinates for below the Comp.1 axis. Using these coordinates the legend disappears. Below is the code for the PCA and legend. Thanks in advance for the help. Regards Tinus r.cols <-

Hi R-community, I am doing a PCA and I need plots for different combinations of axes (e.g., PC1 vs PC3, and PC2 vs PC3) with the arrows indicating the loadings of each variables. What I need is exactly what I get using biplot (pca.object) but for other axes. I have plotted PC2 and 3 using the scores of the cases, but I don't get the arrows proportional to the loadings of each variables on

Dear R users, i would like to apply a PCA on image data for data reduction. The image data is available as three matrices for the RGB values. At the moment i use x <- data.frame(R,G,B)#convert image data to data frame pca<-princomp(x,retx = TRUE) This is working so far. >From this results then i want to create a new matrix from the first (second..) principal component. Here i stuck.

Hi all, I am doing bootstrap on a distance matrix, in which samples have been drawn with replacement. After that I do PCA on a resulted matrix, and these 2 steps are repeated 1000 times. pca(x) is a vector where I wanted to store all 1000 PCAs; and x is from 1 to 1000 SampleD is a new matrix after resampling; I am getting the following error message, which I don't understand: ....

Dear R-helpers, Using the package Lattice, I performed a PCA. For example pca.summary <- summary(pc.cr <- princomp(USArrests, cor = TRUE)) The Output of "pca.summary" looks as follows: Importance of components: Comp.1 Comp.2 Comp.3 Comp.4 Standard deviation 1.5748783 0.9948694 0.5971291 0.41644938 Proportion of Variance 0.6200604

Hi, I have to construct an ellipse interval region on a PCAbiplot, I have my parameters "a" and "b" and I would apply the formula: draw.ellipse(x, y, a = , b = ) I have done a PCA on my data so I have my scores and loading for the first and second component, but my answer is what I have to choose as X and Y into the formula? if "a" and "b" are scalars or