search for: screeplot

Full_Name: Gavin Simpson Version: 2.5.0 OS: Linux (FC5) Submission from: (NULL) (128.40.33.76) Screeplots are a common plot-type used to interpret the results of various ordination methods and other techniques. A number of packages include ordination techniques not included in a standard R installation. screeplot() works for princomp and prcomp objects, but not for these other techniques as it was not...

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...

Hi, Suppose I've got a data set that I found the eigenvalues and eigenvectors. Then I want to draw a screeplot for the eigenvalues. However it returns: Error in matrix(w.m, nc = NC) : negative extents to matrix Then I tried a plot() function on the eigenvalues of the data set, and I can successfully draw it! And the output looks like what the screeplot would show to me. Is there any difference betwee...

Hi Here is a suggestion/request for stats::screeplot * Please allow to change xlab and ylab. The attached file does that, and maintains current behavior. Reason: strictly speaking, ylab should be "Variance" (singular). It would also be better to have an explicit label for the X axis: "Component" In the would-be-nice-if category...

Hi! Is there an easy way to add to the scree-plot labels to each value pertaining to the cumulative proportion of explained variance? Thanks and a happy new year Anne ---------------------------------------------------- Anne Piotet Tel: +41 79 359 83 32 (mobile) Email: anne.piotet@m-td.com --------------------------------------------------- M-TD Modelling and Technology Development PSE-C

...is can occasionally give a different result. - prcomp.Rd is changed to reflect the changes to prcomp - summary.prcomp is defined to give a result more like summary.princomp. In Splus summary.prcomp does not exist and the default method is applied, but this difference seems reasonable? - screeplot is made generic - function(...) is changed to function(x, ...) for some plot methods (to be consistent with plot). - print.prcomp is changed to be more like the result from Splus (Splus just uses print.default, so this is slightly different) - plot.prcomp is renamed screeplot.prcomp t...

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, in total, a certain number of comonents explain. Thanks a lot! -- Dim...

...read.csv2("C:\\Users\\Thomas\\Desktop\\PCA_Kohle_2007.csv",sep=";") rownames(analyse)<-analyse[,1] analyse<-analyse[,-1] require(tseries) require(fSeries) st_rendite<-apply(apply(analyse,2,log),2,diff)*100 pca<-princomp(st_rendite) summary(pca) summary(pca)$loadings screeplot(pca,type="l",main="Screeplot der Prinzipalkomponenten") loadings_pca1<--summary(pca)$loadings[1:6] loadings_pca2<--summary(pca)$loadings[7:12] loadings_pca3<--summary(pca)$loadings[13:18] loadings_pca4<--summary(pca)$loadings[19:24] loadings_pca5<--summary(pca)$...

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: http://r.789695.n4.nabble.com/Scree-diagram-tp2122453p2122453.html Sent from the R help mailing list archive at Nabb...

Hi, I?ve written a script to run several multivariate statistical analysis automatically. As one result a biplot and screeplot is produced. Now I?d like to display the name of the inputdatset as part of the title of these graphics and I do not want to enter it each time I run the script. How can I extract the name of a dataset? An (shortened) extraction from the script: >Test1 <- function(pri, sec, Pca, Ca, Dca){ &g...

Dear all, I am having trouble plotting a PCA result. The plot doesn't appear!!! R goes through without any errors but doesn't make a plot appear!! Could it be wrong window parameters? In this case how do I change them? I am under red hat 9 with the latest version of R! Thanks. Vous manquez d??espace pour stocker vos mails ?

Hi! I am getting a lot of numbers in the background of the pca screeplots if i use call("plot") and eval(somecall). Til now, creating the calls and plotting later on this way worked fine. Example: pcaI<-prcomp(iris[,1:4]) plot(pcaI) x<-call("plot",pcaI) eval(x) Anyone got an idea how i can avoid that? (also it might take a second or so f...

...3096 1.14002 1.09911 0.77523 1.31543 This is what I have done in R: rm(list=ls(all=TRUE)) PCA<-read.csv("PCAsites.csv", header=TRUE) PCA attach(PCA) names(PCA) ###using correlation matrix##### p1 <- princomp(PCA, cor = TRUE) summary(p1) loadings(p1) plot(p1) biplot(p1) p1$scores screeplot(p1, npcs=4, type="lines") The problem is that the purpose of this analysis is to derive a new data set, using the first component, that I can work with. In other words, I am doing this to compress the data into one column of ring widths. How do I derive a new data set? Thanks in advance...

...ies are small but not zero Loadings: Comp.1 Comp.2 Comp.3 Murder 0.663 0.314 0.679 Assault 0.694 -0.715 UrbanPop 0.279 -0.946 0.165 Comp.1 Comp.2 Comp.3 SS loadings 1.000 1.000 1.000 Proportion Var 0.333 0.333 0.333 Cumulative Var 0.333 0.667 1.000 > plot(pc.cr <http://pc.cr>) # shows a screeplot. > biplot(pc.cr <http://pc.cr>) Erreur dans plot.window(xlim, ylim, log, asp, ...) : 'xlim' nécessite des valeurs finies > [[alternative HTML version deleted]]

...for each eigenvector. This is the best part of the SPSS output. I feel like I'm close to finding the right function in R , but I don't know how to look at the eigenvalues of each component in R. princomp() seems a step in the right direction. 8. Scree plot. No problem, princomp() and screeplot() seem to produce about the right result. 9. Component matrix (lists the variable loading on each factor) factanal() seems to do this, but again the results don't jive with SPSS and I'm unsure why. 10. Factor rotation No problem, factanal(rotation="varimax") does this. If...

Dear R users, psy version 0.65 is on CRAN now. psy provides several methods used in psychometry (kappa, icc, cronbach, screeplot (with simulations), non linear mapping, etc.) A bug has been fixed in function wkappa (weighted kappa): in particular circumstances, the 2*2 table presented levels in an order different to what was suggested in the help file. The fpca function (PCA plot with pvalues and compatible with dependant/...

Hi: How i can make principal component plots with R. thanks Jorge -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send "info", "help", or "[un]subscribe" (in the "body", not the subject !) To: r-help-request at stat.math.ethz.ch

Hi, is there any other routine for factor analysis in R then factanal? Basically I'am interested in another extraction method then the maximum likelihood method and looking for unweighted least squares. Thanks in advance Sigbert Klinke

Dear R users, psy version 0.65 is on CRAN now. psy provides several methods used in psychometry (kappa, icc, cronbach, screeplot (with simulations), non linear mapping, etc.) A bug has been fixed in function wkappa (weighted kappa): in particular circumstances, the 2*2 table presented levels in an order different to what was suggested in the help file. The fpca function (PCA plot with pvalues and compatible with dependant/...

Hello again, I was hoping one of you could help me with this problem. Consider the sample data from R: > summary(prcomp(USArrests)) Importance of components: PC1 PC2 PC3 PC4 Standard deviation 83.732 14.2124 6.4894 2.48279 Proportion of Variance 0.966 0.0278 0.0058 0.00085 Cumulative Proportion 0.966 0.9933 0.9991 1.00000 How do I access the