similar to: 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

Hi, I input this command on a data matrix: > apply( traffic, 2, function(z){z-mean(z)}) and got the following error messages: Error in sum(..., na.rm = na.rm) : invalid "mode" of argument Can anyone tell me what is wrong here? Running R1.2.2 on Windows ME. Cheers, Kevin ------------------------------------------------- Ko-Kang Kevin Wang Statistical Analysis Division

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

I have two problems where I've come up with some code that will do the analysis that I want, but it looks pretty clumsy. In the first case, I calculate the variance on five different columns for each of 14 clusters and get them into one matrix. I get the job done, but I would have thought that it could be done in one or two lines, not six, and be generalized so that it didn't matter how

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, I suggest addition of a

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,

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'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:

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

I'm wondering if I'm seeing a regression of 5616 with 3.4.8. I'm trying to set-up pptpd with winbind, which I'm doing on two machines on Debian lenny, and I'm trying on Debian Squeeze now. The Windows client is saying "Error 778: It was not possible to verify the identity of the server." The logs say that everything is ok, and that the client is hanging up the

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

Hi, I run the example for princomp for R211 I got the following error for biplot > ## The variances of the variables in the > ## USArrests data vary by orders of magnitude, so scaling is appropriate > (pc.cr <http://pc.cr> <- princomp(USArrests)) # inappropriate Erreur dans cov.wt(z) : 'x' must contain finite values only > princomp(USArrests, cor = TRUE) # =^=

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

At the moment I've no host with which to test ppp/lp things on. If you (or anyone) could provide complete regexes, (each beginning with ^ and ending with $) I will patch the rules accordingly. If not, sending the full log lines is a good start. Thanks! -- [ Todd J. Troxell ,''`. Student, Debian GNU/Linux Developer, SysAdmin, Geek :

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

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

Dear all, as the bug # 7991 is flagged not-reproducible, let me give you some pieces of code, as I have the same or similar problem. The problem always shows up with the first example (a small text file) and only sometimes (but without obvious pattern) with the second example, which is a binary file. > download.file("ftp://ftp.nhc.noaa.gov/pub/atcf/btk/bal012006.dat",

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 ?

Dear all, I'm trying to produce 2 png files, one consisting of an image plot and a color-table (also an image plot) and the other one consisting of 4 image plots and a color table. I'd like the color table to be exactly the same. The way I proceded is the following: for one plot and the color-table png(file = png.file, width = 650, height = 800, pointsize = 16) layout(matrix(c(1, 2),

Dear all, I have some C code using the zlib uncompress. It compiles fine under Win XP (with SP2) with R up to 2.3.1 and generates four files (Makedeps, *.d, *.dll, and *.o): C:\Data\R-packages\cmg>rcmd SHLIB cmg.c -o cmg.dll latex: not found making cmg.d from cmg.c gcc -IC:/PROGRA~1/R/R-2.3.1/include -Wall -O2 -c cmg.c -o cmg.o cmg.c: In function `readcmg': cmg.c:33: warning: