similar to: scree plot

Hello, I am using function princomp from the package psych. I have my principle component object mypc: mypc <- princomp(covmat=mycor) plot(mypc) # shows me a screeplot Question: how could I actually see the values displayed in the screeplot. I don't mean on the graph - I just want to know the actual value for each component (e.g., 10, 3.2, 1.8, etc.) I need to know how much variance,

hello, I've two questions today. 1) I'm trying to do a scree diagram, I did a Google for a specific command I could used to do so. All I could find is a screeplot. Are they the same command? 2) what command can I used to present a PC scores, eigenvectors of the PC scores, and component correlations? thanks! -- View this message in context:

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:

I'm trying to make a screeplot including the Cumulative Proportion of the Variance, something that can easily be done in S-Plus with 'screeplot(pc.object,cumulative=T)'. How can I access the Proportion of Variance in an princomp object and how could I get the Cumulative Proportion of the Variance on the screeplot? Many thanks in advance, Jan:-) --

Dear all, I'm using R 2.8.1 under Vista. I programmed a Simulation with the code enclosed at the end of the eMail. After the simulation I want to analyse the columns of the single simulation-runs, i.e. e.g. Simulation[[1]][,1] sth. like that but I cannot address these columns... Can anybody please help? Best, Thomas ############################ CODE ############################

Hello, I am trying to complete a PCA on a set of standardized ring widths from 8 different sites (T10, T9, T8, T7, T6, T5, T3, and T2). The following is a small portion of my data: T10 T9 T8 T7 T6 T5 T3 T2 1.33738 0.92669 0.91146 0.98922 0.9308 0.88201 0.92287 0.91775 0.82181 1.05319 0.92908 0.97971 0.95165 0.98029 1.14048 0.77803 0.88294 0.96413 0.90893 0.87957 0.9961 0.74926 0.71394 0.70877

My previous post about prcomp and princomp was done in some haste as I had long ago indicated to Kurt that I would try to have this ready for the June release, and it appeared that I would miss yet another release. I also need to get it out before it becomes hopelessly buried by other work. Brian Ripley kindly pointed out some errors, and also pointed out that I was suggesting replacing some

Hi, I am using R to do a principal components analysis for a class which is generally using SPSS - so some of my question relates to SPSS output (and this might not be the right place). I have scoured the mailing list and the web but can't get a feel for this. It is annoying because they will be marking to the SPSS output. Basically I'm getting different values for the component

--=-hiYzUeWcRJ/+kx41aPIZ Content-Type: text/plain; charset="UTF-8" Content-Transfer-Encoding: 8bit Hi, on March 10 I filed a wishlist bug report asking for the inclusion of some changes to factanal() and the associated print method. The changes were originally proposed by John Fox in 2005; they make print.factanal() display factor correlations if factanal() is called with rotation =

********************************************************************** This email and any files transmitted with it are confidential and intended solely for the use of the individual or entity to whom they are addressed. If you have received this email in error please notify the system manager. ********************************************************************** Hi Just installing R and some

Sorry, sorry.... of course levels(testf)[c(2,1,3)] will do the job My excuses to all Anne PS I will meditate the following saying "la parole est d'argent et le silence est d'or" BONNES FETES A TOUS SEASONAL GREETINGS ---------------------------------------------------- Anne Piotet Tel: +41 79 359 83 32 (mobile) Email: anne.piotet@m-td.com

Hello R helpers! I looked but did not find a table-lookup R-utility. I could use a loop to do the job (old FORTRAN/C habits die hard) but if I have a big table in which I have to search for the values corresponding to a vector, I end up logically with a double loop. Is there already such a utility? Otherwise, is there a way without loops? Thanks as always Anne

I've read through many postings about principle component analysis in the R-help archives, but haven't been able to piece together the information I need. I'd like to recreate an SPSS-like experience of factor analysis using R. Here's what SPSS produces: 1. Scatterplots of all possible variable pairs, with regression lines. xyplot(my.dataframe) is perfect but for the lack of

Hello! I would like to know if there is a simple way to reorder levels of a given factor.Let's say that the vector testf<-factor(c("red","red","red","blue","blue","white")) levels(testf) : blue red white should have reordered levels such as levels(testf) : red blue white (this is for presentation purposes) I guess

Is there an upgrate of the vcd library (visualisation of categorical data) for the latest R version? Trying to download it from CRAN I get URL /data/WWW/ftp/pub/R/bin/windows/contrib/r-release/vcd_0.1-3.4.zip was not found on this server. googling it, I found it for instance on http://www.sourcekeg.co.uk/cran/bin/windows/contrib/1.9/ but trying to install it gave me the message >

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)

hello to all! I have a small problem wit predict() for lm Let's say I have predictors x1 and x2, response y I want to predict for a new ds say dn<-data.frame(x1= seq(min(x1),max(x1),length=10),x2=rep(median(x2),10)) predict(lm(y~x1+x2),dn,se.fit=T) Error message > Error: variables 'x1', 'x2' were specified differently from the fit (I looked in the help and found

Hi! I'm trying to do some ploytomous logistic regression using multinom() in the nnet package, but am a bit confused about interpretation of the results Is it possible to get the following quantities: I: maximum likelihood estimates to test for fit of model and significance of each predictor (I would like to produce a table of the following type) Analysis of Variance: MLE (values are

Hello R helpers! I'm using the avas function form package acepack (called from areg.boot package Hmisc) to estimate automatically transformations of predictors (in this case monotonous) and response. Well, it seems to work quite well, but I have 3 basic questions: - which set of basis functions is used in this procedure? - how do I back transform my estimate (y hat ) to the originasl scale?

Dear R-helpers, I have a problem which I suppose is trivila, but... I have included NA values as factors ( (to be able to make nice printed summaries with NAs % ba category ) with the following code dat$x.f<-factor(dat$x, exclude=NULL); levels(dat$x.f)<-c("A1","A2","A3","A4","NA"); length(dat$x.f) Now, I want to impute the missing values.