similar to: cor.test() -> p-values may be incorrect due to tie

Full_Name: Ian Wilson Version: R-1.8.0 OS: Windows (but own compilation) Submission from: (NULL) (139.133.7.38) the p-value is incorrect for cor.test using method "spearman" in R-1.8.0. This was not the case in R-1.7.1. Version R-1.8.0 on Windows > cor.test(rnorm(50),rnorm(50),method="spearman") Spearman's rank correlation rho data: rnorm(50) and rnorm(50) S

Full_Name: David Simcha Version: 2.10 OS: Windows XP Home Submission from: (NULL) (173.3.208.5) > a <- c(1:10) > b <- c(1:10) > cor.test(a, b, method = "spearman", alternative = "greater", exact = TRUE) Spearman's rank correlation rho data: a and b S = 0, p-value < 2.2e-16 alternative hypothesis: true rho is greater than 0 sample estimates:

Hi weekend R helpers, *Is it possible to get a ylabel on the right hand side axis ==> axis(4)? Thankx, Jan opar <- par(mfrow = c(2,1)) plot(ts.Origi[,1],ylab='NDII', main=name) # Add the legend text in the right order legend.txt <- c("NDII", "Inverse KBDI") legend(2001.8,0.29, legend.txt,col=c(1,2), lty=1)

Hi, all. Yet another beginner to R : ) I wonder, how it's possible to get the value of a coefficient from the object produced by cor.test() ? > cor.test(a, b, method="spearman") Spearman's rank correlation rho data: a and b S = 21554.28, p-value = 2.496e-11 alternative hypothesis: true rho is not equal to 0 sample estimates: rho 0.6807955 Warning message: In

Greetings. There seems to be a problem with the P-value computation in the cor.test with method="spearman". In R1.8.0 (MS Windows) I seem to be getting intermittently nonsense P-values, but the rho's are OK. I can get this reproducibly with the toy example attached where the first use is OK and subsequent calls with the same data give nonsense. (I have also seen the problem

Hello, Is there a way to extract output objects from a stats test without viewing the entire output? I am trying to do so in the following: define a vector of length j for( i in 1: length (vector)) { vector[i] = cor.test (datavector1, datavector2[i], method=("spearman")) } I would like the reported Spearman's rho to be saved in a vector. I have tried a few different ways of

I would like to turn off the warnings from cor.test while retaining exact=NULL. Is that possible ? > cor.test(c(1,2,3,3,4,5), c(1,2,3,3,4,5), method = "spearman") Spearman's rank correlation rho data: c(1, 2, 3, 3, 4, 5) and c(1, 2, 3, 3, 4, 5) S = 0, p-value < 2.2e-16 alternative hypothesis: true rho is not equal to 0 sample estimates: rho 1 Warning message: In

Dear HelpeRs, I have some data: "ice" <- structure(c(0.386, 0.374, 0.393, 0.425, 0.406, 0.344, 0.327, 0.288, 0.269, 0.256, 0.286, 0.298, 0.329, 0.318, 0.381, 0.381, 0.47, 0.443, 0.386, 0.342, 0.319, 0.307, 0.284, 0.326, 0.309, 0.359, 0.376, 0.416, 0.437, 0.548, 41, 56, 63, 68, 69, 65, 61, 47, 32, 24, 28, 26, 32, 40, 55, 63, 72, 72, 67, 60, 44, 40, 32, 27, 28, 33,

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

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

Hello all, I've run a Spearman's Rank test to discern relationships between landscape characteristics and a specific aspect of river behaviour. I've executed a correlation matrix between the one dependent variable and all of the predictors, which gives me a nice output of Spearman's Rho values. However, I also need to somehow find the "p" value, to assess the strength

Hi, I am not sure which kind of test is applied to the data if you use cor.test(x, y) ? Is it an unpaired t-Test? Regards -- View this message in context: http://www.nabble.com/cor.test%28x%2Cy%29-tp22492993p22492993.html Sent from the R help mailing list archive at Nabble.com.

I want to give the "rho" value below to another variable.How ? > Spearman's rank correlation rho > > > > data: a[, 3] and a[, 2] > > S = 22, p-value = 0.001174 > > alternative hypothesis: true rho is not equal to 0 > > sample estimates: > > rho > > 0.8666667 -- View this message in context:

R: I got a warning message when running the cor.test function using both Spearman and Kendall rank correlations saying that the p-value may be incorrect due to ties in the data. My data has 35 obs and one series has 6 pairs of ties. Does anyone know if this would likely have a great effect on the p-values calculated.. The values look good; tau = -0.68 with p-value = 8e-9 and rho = =0.84

cor.test and SuppDists give me different P-values for the same Spearman's rho. Which is correct, or am I doing something wrong? > x <- c(44.4, 45.9, 41.9, 53.3, 44.7, 44.1, 50.7, 45.2, 60.1) > y <- c( 2.6, 3.1, 2.5, 5.0, 3.6, 4.0, 5.2, 2.8, 3.8) > cor.test(x,y,method="spearman") Spearman's rank correlation rho data: x and y S = 48, p-value =

Dear R users, Does the results below make any sense? Can the the interval of the correlation coefficient be between *-1.0185* and -0.8265 at 95% confidence level? Liviu > library(boot) > data(mtcars) > with(mtcars, cor.test(mpg, wt, met="spearman")) Spearman's rank correlation rho data: mpg and wt S = 10292, p-value = 1.488e-11 alternative hypothesis: true rho is not

Hello, I wuold like to propose my modifications of the original cor.test to you : I tried to calcolate the correct p-value for Spearman and Kendall's test with ties. Let me know what you think. Thanks you for your time. Antonietta di Salvatore test <- function(x, ...) UseMethod("test") test.default <- function(x, y, alternative = c("two.sided",

Hi, I have a dataset with four categories of data, the number of samples are not the same in each category. I want to find the Spearaman's Rho. Let me give an example. x=(14.22770439,26.49420624,46.7277932,19.02550707,23.37379361,16.97789862,19.77100085,23.11270162,13.72929843,33.54430621,14.4756979,70.15811106,11.22789833,NA,NA,NA)

Dear R, I'm looking for exact p-values for Spearman's rank correlation in the presence of ties. This is available in StatXact and SPSS, but I haven't yet found it in R. Has anyone implemented this? -- Paul H Artes Assist Prof, Ophth Vis Sci Dalhousie University, QEII Eye Care Centre Halifax, NS, CANADA

Hi, I am not really familiar with Correlation foundations, although I read a lot. So maybe if someone kindly help me to interpret the following results. I had the following R commands: correlation <-cor( vector_CitationProximity , vector_Impact, method = "spearman", use="na.or.complete") cor_test<-cor.test(vector_CitationProximity, vector_Impact,