similar to: progressive slowdown during script execution? [Broadcast ]

I'm an R novice, so I hope my question is a valid one. I'm trying to run the following script in the current version of R. for (i in 1:1640){for (j in (i+1):1641){ if (i == 1 && j == 2){x <- cor(sage[i,],sage[j,],method="spearman"); y <- cor(frie[i,],frie[j,],method="spearman")} if (i != 1 || j != 2){x <-

Hi all, I am having some strange problems on a customer site : I am pretty sure I disabled oplocks entirely - I don't like the idea of them, and would rather have stable data. The customer (5-10 users) runs Sage Line 50, and until recently was getting numerous errors within sage. We have temporarily replaced their switch and this has reduced the fatality rate and a few

Hi. I am currently trying to run some Spearman correlations, and have encountered two issues. 1) When using cor.test() with a variable that includes ties, I get the "Cannot compute exact p-values with ties" error. I have read that this function now uses an asymptotic formula that allows for ties, so do not understand why I am getting this error. (I am running version 2.4.0.) I

Full_Name: Vladimir Obolonkin Version: tested in 2.0 to 2.4.1 OS: linux, win, mac Submission from: (NULL) (202.14.96.194) {{ Subject shortened manually -- to pass anti-spam filters Original-subject: Incorrect matrix of spearman correlations working \ large (24000 by 425 and 78 by 425 data frames) in 64-bit Linux machines;\ the same code gives correct results in 32-bits Win and

vobolonkin at lic.co.nz wrote: > Full_Name: Vladimir Obolonkin > Version: tested in 2.0 to 2.4.1 > OS: linux, win, mac > Submission from: (NULL) (202.14.96.194) > > {{ Subject shortened manually -- to pass anti-spam filters > > Original-subject: Incorrect matrix of spearman correlations working \ > large (24000 by 425 and 78 by 425 data frames) in 64-bit Linux

Hi, I have just encountered a strange behaviour from 'cor' with regards to the treatment of NAs when calculating Spearman correlations. I guess it is a subtle bug. If I understand the help page correctly, the two modes 'complete.obs' and 'pairwise.complete.obs' specify how to deal with correlation coefficients when calculating a correlation _matrix_. When calculating

Dear R-users: I am runing R 1.6.2 with Windows XP. I try to calculate Pearson correlation and Spearman correlation of any pairwise columns of 8000 x 80 data matrix with missing values and randomize the matix 1000 times and calculate this two correlations again. The code bellow for Pearson is working fine but for Spearman got the error bellow for randomized data matrix and R console is stop

Hallo! Is there an error in cor to calculate Spearman correlation with cor if there are NA's? cor.test gives the correct result. At least there is a difference. Or am I doing something wrong??? Does anybody know something about this? a<-c(2,4,3,NA) b<-c(4,1,2,3) cor(a, b, method="spearman", use="complete.obs") # -0.9819805 cor.test(a, b,

I'm aware that S N Krishna asked the same question. However, I have failed to implement the posted solution for running rank order correlations on multiple subsets of data using the by() function. Here is my problem: Take a set of data from two subjects, who provided numerical infant mortality (IM) estimates for five countries: sub <- c(1, 1, 1, 1, 1, 2, 2, 2, 2, 2)

Hi, How can I accomplish this in R. Example: I have the following data.frame: data <- data.frame(x=c(1,2,3,4,5,6,5,3,7,1,0,4,8),y=c(1,2,1,2,2,2,1,1,1,2,2,2,2),z=c(5,8,4,3,4,1,6,3,3,6,3,5,7)) Supposing that data$y is a factor, I would like to find the Spearman correlation between data$x and data$z indexing it by data$y. To be more specific, I want to find two correlations: between x and z

There seems to be an issue in computing rank correlations with missing values present. I think this comes from the way rank() function works but I am not sure how to go about this. Rank function places missing values at the end by default thus skewing the rank relationship between two vectors: Example: R : Copyright 2003, The R Foundation for Statistical Computing Version 1.8.1 (2003-11-21),

Hi, I'm currently trying to find/define a relationship between one dependent and several independant variables. The problem is that i cannot use the normal multiple regression/correlation in Spss because the data is not normal distributed. i calculated the spearman roh and Kendalls tau Correlation and also some partial correlations in R. Now i wanna find out the the multiple correlation

Arrrrr, I think I've found a bug in the behavior of the stats::cor function when NAs are present, but in case I'm missing something, could you look over this example and let me know what you think: > a = c(1,3,NA,1,2) > b = c(1,2,1,1,4) > cor(a,b,method="spearman", use="complete.obs") [1] 0.8164966 > cor(a,b,method="spearman",

Hello, I calculate the correlation between two matrices cor(x,y, method="spearman") and I am wondering if it possible to see only the significant correlations. I can do that for single OTUs with cor.test() command but I would like to have an output for whole matrix. Besides,I would like to plot the relative abundance of each OTU, against the number of significant pair-wise correlations

I have a data set where I want the correlations between 2 variables conditional on a students grade level. This code works just fine. by(tmp[,c('mtsc07', 'DCBASmathscoreSPRING')], tmp$Grade, cor, use='complete', method='pearson') However, this generates an error by(tmp[,c('mtsc07', 'DCBASmathscoreSPRING')], tmp$Grade, cor, use='complete',

Hello, I would like to do a PCA with dudi.pca or PCA, but also with the use of Spearman or Kendall correlations Is it possible ? Otherwise, how can I do, according to you ? Thanking you in advance Eric Bourgade RTE France [[alternative HTML version deleted]]

Dear R maintainers, R is great. Now that I have that out of the way, I believe I have encountered a bug, or at least an inconsistency, in how Spearman and Kendall rank correlations are handled. Specifically, cor() and cor.test() do not produce the same answer when the data contain NAs. cor() treats the NAs as data, while cor.test() eliminates them. The option use="complete.obs" has

Hello to all of you, R-expeRts! I am trying to compute the cor.test for a matrix that i labelled mydata according to mydata=read.csv... then I converted my csv file into a matrix with the mydata=as.matrix(mydata) NOW, I need to get the p-values from the correlations... I can successfully get the spearman's correlation matrix with: cor(mydata, method="s",

Hi The cor function in the stats package calculates the correlation between columns of data, does anyone know if it is at all possible to calculate the correlation between rows instead ? Or is there an appropriate package or function that is more appropriate I'd like to calculate spearman & pearson correlations between rows. Many thanks Jason -- --------------------------------

I have a quick statistical question and hoped somebody has a tip for me without me having to go to the local statistician on Monday. I assess 4 statements from 90 subjects. Each of the 4 statements receives one of three responses (say -1, 0, or 1). I can use Cramer's V or Spearman correlations to assess the correlation between each pair of statements, but I am looking for a measure of