similar to: values for the scree plot (package psych)

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

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:

Hello! I am running principle components analysis using princomp function in pacakge psych. mypc <- princomp(mydataforpc, cor=TRUE) Question: I'd like to use pairwise deletion of missing cases when correlations are calculated. I.e., I'd like to have a correlation between any 2 variables to be based on all cases that have valid values on both variables. What should my na.action be in

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

Hi Rui, Many thanks. The R code works BUT the results I get are quite weird I guess ! MIC = 0.2650 Normal 95% CI = (0.9614, 1.0398) The MIC is not inside the confidence intervals ! Is there something wrong in the R code ? Here is the reproducible example : ########## C=c(2,4,5,6,3,4,5,7,8,7,6,5,6,7,7,8,5,4,3,2) D=c(3,5,4,6,7,2,3,1,2,4,5,4,6,4,5,4,3,2,8,9) library(minerva) mine(C,D)$MIC

You need: myCor <- function(data, index){ mine(data[index, ])$MIC[1, 2] } results=boot(data = cbind(C,D), statistic = myCor, R = 2000) boot.ci(results,type="all") Look at the differences between: mine(C, D) and mine(cbind(C, D)) The first returns a value, the second returns a symmetric matrix. Just like cor() David L. Carlson Department of Anthropology Texas A&M

Hello, First of all, when I tried to use function mic I got an error. mic(cbind(C, D)) Error in mic(cbind(C, D)) : could not find function "mic" So I've changed your function myCor and all went well, with a warning relative to BCa intervals. myCor <- function(data, index){ mine(data[index, ])$MIC } results=boot(data = cbind(C,D), statistic = myCor, R = 2000)

Dear R-Experts, Here below is my R code (reproducible example) to calculate the confidence intervals around the spearman coefficient. ########## C=c(2,4,5,6,3,4,5,7,8,7,6,5,6,7,7,8,5,4,3,2) D=c(3,5,4,6,7,2,3,1,2,4,5,4,6,4,5,4,3,2,8,9) cor(C,D,method= "spearman") library(boot) myCor=function(data,index){ cor(data[index, ])[1,2] } results=boot(data=cbind(C,D),statistic=myCor, R=2000)

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

I'd like to paste a superstring with a number in an object. Thanks for any help. Murray mycor <- cor(1:10,1:10) plot(1:10,1:10) text(8,2,paste(expression(R^2)," = ",mycor)) [[alternative HTML version deleted]]

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

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

Hi all I'm using Samba version 3.0623 on my rhel3 with kernel 2.4.21 I mounted a windows partition with the command "smbmount". Now I have the problem, that I cannot copying files which are larger than 2GB. The cp command starts and everything works fine until the 2GB are reached. Then the bash gives me the error: "The maximum file size is exceeded" the

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

Hello! I am doing Principle Componenets Analysis using "psych" package: mypc<-principal(mydata,5,scores=TRUE) However, I was asked to run a case-weighted PCA - using an individual weight for each case. I could use "corr" from "boot" package to calculate the case-weighed intercorrelation matrix. But if I use the intercorrelation matrix as input (instead of the

Dear R users,I have a list of equally structured lists, how can I access e.g. all 2nd compontents in those sub-lists?An example:lst <- list(rep(list(1:3),3), rep(list(4:6),3))> lst[[1]][[1]][[1]][1] 1 2 3[[1]][[2]][1] 1 2 3[[1]][[3]][1] 1 2 3[[2]][[2]][[1]][1] 4 5 6[[2]][[2]][1] 4 5 6[[2]][[3]][1] 4 5 6What I want to get are all second sub-lists, in this case:[[1]][[2]][1] 1 2

Hello list, I am configuring NUT on a Dell workstation that is connected to an APC Back UPS Pro 1300, and I am running into an issue where starting everything manually seems to be successful, but fails when using the init.d scripts. The short story is that I get the following error when trying to check the status of the UPS: [root at mypc:/etc/nut]# upsc apc_pro1300 at localhost Error:

Hi, I'm a PhD student and I want to use the function 'daisy' from the package 'cluster' to compute dissimilarities. My variables are of mixed types so I use the argument 'stand' in daisy to define the type of my variables. I have the following error message : Warning message: binary variable(s) 13, 16, 17, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32, 33, 34, 35,

Dear Tomas, Thanks for your email. This can be easily reproduced if there is access to a shared folder. I've reproduced this behaviour on both Windows 7 and Windows 10. Steps to reproduce: 1. Ensure Offline Files is turned on within Windows 2. Using Windows Explorer, browse to a folder shared on a network using a UNC path, e.g. \\mypc\myshare\ 3. Create a test file, e.g. test.txt 4. Within