Displaying 20 results from an estimated 6000 matches similar to: "Bug in Kendall for n<4?"
2009 Jul 06
0
Error exit, tauk2. IFAULT = 12
Dear McLeod,
I am a PhD student working with the Kendall package in R. I used
this package to make the so-called TauKr matrix correlation analysis
(Hemelrijk, 1990). However, lately I've been getting this warning:
Error exit, tauk2. IFAULT = 12. From the forums I've been looking
into, I noticed that when this warning pops up, it means that the
p-value of the correlation may be wrong
2012 Aug 20
1
Kendall package tau-a, b, and c
Hi all,
I would like to ask a question related to Kendall package. I ran Kendall
(x,y) and
saw the results. But I am not sure which tau values R reported. I have
ties in
my data set, so I want tau-b. Can anybody tell how Kendall package is
calculating tau values? I have looked at the package PDF, but I could
not find
any useful information. As long as I see from the following
link, there
2013 Apr 30
0
How to handle errors generated by the package Kendall?
Hello everyone,
I'm trying to perform some trend analysis using the package "Kendall" and for some observations I get the following error "WARNING Error exit, tauk2, IFAULT = 12"
Could anyone tell me how to create a logical vector of observations having its element equal to TRUE if this error happened and FALSE if not?
I'm thinking about something similar as
2005 Aug 18
2
kendall tau correlation test for ties: Potential error (PR#8076)
Full_Name: Dirk Koschuetzki
Version: 2.1.1
OS: source code
Submission from: (NULL) (194.94.136.34)
Hello,
>From the source code (R-2.1.1, file: .../R-2.1.1/src/library/stats/R/)
******************************
cor.test.default <-
function(x, y, alternative = c("two.sided", "less", "greater"),
method = c("pearson", "kendall",
2011 Mar 21
1
Kendall v MannKendall Functions
Hi,
I am running a correlation analysis on a temporal dataset. I was wondering
if you would receive the same tau and p values running the function:
MannKendall(x), where x is the dependant variable that changes with time
as you would running:
Kendall(d,x), where x is the exact same dataset as the x entered into
MannKendall and d is the date on which the observation was made (assuming
that
2009 Dec 13
2
O(N log N) Kendall Tau
I've noticed that the implementation of Kendall's Tau in R is O(N^2).
The following reference describes how it can be done in O(N log N):
A Computer Method for Calculating Kendall's Tau with Ungrouped Data
William R. Knight
Journal of the American Statistical Association, Vol. 61, No. 314, Part
1 (Jun., 1966), pp. 436-439
http://www.jstor.org/pss/2282833
I'm interested in
2005 Aug 13
1
R/S-Plus/SAS yield different results for Kendall-tau and Spearman nonparametric regression
Colleagues,
I ran some nonparametric regressions in R (run in RedHat Linux), then
a colleague repeated the analyses in SAS. When we obtained different
results, I tested S-Plus (same Linux box). And, got yet different
results. I replicated the results with a small dataset:
DATA:
37.5
23
37.5
13
25
16
25
12
100
15
12.5
19
50
20
100
13
100
10
100
10
100
16
50
10
87.5
2007 Sep 06
3
kendall test
Hello,
I thougth that there is a function which does the kendall test in R,
I writed on the console apropos("kendall") and I didn't found anything can you tell me how could I do to use the kendall test?
Thanks.
_____________________________________________________________________________
l
[[alternative HTML version deleted]]
2006 Sep 12
1
Kendall's tau-c
Hello,
I can't find a package which calculates Kendall's tau-c. There is the
package Kendall, but it only calcuates Kendall's tau-b.
Here is the example from
ttp://www2.chass.ncsu.edu/garson/pa765/assocordinal.htm.
cityriots <- data.frame(citysize=c(1,1,2,2,3,3),
riotsize=c(1,2,1,2,1,2), weight=c(4,2,2,3,0,4))
cityriots <- data.frame(lapply(cityriots,function(x)
2005 Jun 28
1
faster algorithm for Kendall's tau
Hi,
I need to calculate Kendall's tau for large data
vectors (length > 100'000).
Is somebody aware of a faster algorithm or package
function than "cor(, method="kendall")"?
There are ties in the data to be considered (Kendall's
tau-b).
Any suggestions?
Regards
Ferdinand
2012 Jun 25
2
Fast Kendall's Tau
Hello.
Has any further action been taken regarding implementing David Simcha's fast Kendall tau code (now found in the package pcaPP as cor.fk) into R-base? It is literally hundreds of times faster, although I am uncertain as to whether he wrote code for testing the significance of the parameter. The last mention I have seen of this was in 2010
2010 Feb 08
2
Incorrect Kendall's tau for ordered variables (PR#14207)
Full_Name: Marek Ancukiewicz
Version: 2.10.1
OS: Linux
Submission from: (NULL) (74.0.49.2)
Both cor() and cor.test() incorrectly handle ordered variables with
method="kendall", cor() incorrectly handles ordered variables for
method="spearman" (method="person" always works correctly, while
method="spearman" works for cor.test, but not for cor()).
In
2014 Sep 26
2
summary
Grandiosa Comunidad
Saludos
Quiero pedirles ayuda en los siguientes puntos;
.- Pregunto si puedo sacar, aumentarle al summary también los siguientes
puntos
Intervalo de confianza, desviación estándar?
.- como puedo obtener la correlación en variables dicotómicas ej. tau de
kendall?
.- Como puedo cambiar los títulos del ingles al español de los encabezados
de un surfit(Surv(tiempo, estado))?
--
2009 Aug 16
1
Printing the null hypothesis
Dear R developers,
Currently many (all?) test functions in R describe the alternative
hypothesis, but not the the null hypothesis being tested. For example,
cor.test:
> require(boot)
> data(mtcars)
> with(mtcars, cor.test(mpg, wt, met="kendall"))
Kendall's rank correlation tau
data: mpg and wt
z = -5.7981, p-value = 0.000000006706
alternative hypothesis: true tau is not
2006 Jun 22
2
programming advice
Dear R users
I want to compute Kendall's Tau between two vectors x and y.
But x and y may have zeros in the same position(s) and I wrote the
following function to be sure to drop out those "double zeros"
"cor.kendall" <- function(x,y) {
nox <- c()
noy <- c()
#
for (i in 1:length(x)) if (x[i]!= 0 | y[i] != 0)
nox[length(nox)+1]<- x[i]
for (i in
2004 Mar 03
1
cor(..., method="spearman") or cor(..., method="kendall") (PR#6641)
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
2008 May 05
0
kendall tau a,b,c
hi,
i have 2 lists of ranks for which i'd like to compute kendall tau.
there are ties in the ranks which (to the best of my knowledge) means i
cant use tau a but rather b or c.
how does R handle that? are ties automatically detected (using
corr.test()) and is tau b/c computed instead of tau a?
also kendall does not work when values in list 1 do not occur in list 2
(and vice versa) - how does
2011 Apr 18
4
splom, plotmath: how to add three lines of information with alignment?
Dear expeRts,
I would like to create a scatter plot matrix with splom(). The lower panel should
contain some additional information about the samples shown in the upper panel
plot, see the splom() call below. Now two questions came up:
(1) The lower panels show "tau" and "alpha" on top of each other. How can I plot
*three* expressions on top of each other? I tried several
2009 Oct 20
1
kendall.global
Hi every body:
I need some help with kendall.global. The example in the manual seems not working well, and cannot used with my data, always the same error.
data(mite)
> mite.hel <- decostand(mite, "hel")
>
> # Reproduce the results shown in Table 2 of Legendre (2005), a single group
> mite.small <- mite.hel[c(4,9,14,22,31,34,45,53,61,69),c(13:15,23)]
>
2011 Sep 02
1
Mann Kendall Test for Trend
Hi there,
I'm trying to apply the Mann Kendall test for trend analysis of a time
series. I have downloaded and installed the package Kendall and
subsequently loaded it into the software.
My time series is a .txt file with 2 columns - column 1 is the year (1985 -
2009) and column 2 is the corresponding entry variable.
According to the R guidelines, the call should be:
MannKendall(x)