Re: [R] kendall's w

From: David Barron <mothsailor_at_googlemail.com>
Date: Wed 13 Sep 2006 - 17:47:36 GMT

I assume you are using the function in the concord package. This seems to happen because there is a different way of calculating the p-value depending on whether or not there are more than 7 cases. If there are 7 or less, the function doesn't work unless there are more than two rows. Here is the relevant code fragment:

   if (datadim[2] > 7) {

        p.table <- NA
        x2df <- datadim[2] - 1
        p.chisq <- pchisq(datadim[1] * (datadim[2] - 1) * W,
            x2df, lower.tail = FALSE)

    }
    else {
        p.table <- ifelse(W > Wcrit01[datadim[2] - 2, datadim[1] -
            2], "<0.01", ifelse(W > Wcrit05[datadim[2] - 2, datadim[1] -
            2], "<0.05", "> 0.05"))
        x2df <- NA
        p.chisq <- NA

    }

On 13/09/06, Bianca Vieru <bianca.vieru@free.fr > wrote:
>
> Hi,
>
> I try to calculate Kendall's W coefficient and I have a bizarre error.
>
>
> little.app.mat
> <-matrix(c(1,3,4,2,6,5,2,4,3,1,5,6,3,2,5,1,5,4),nrow=3,byrow=TRUE)
> print(kendall.w(little.app.mat[-1,]))
> >>> Kendall's W for ordinal data
> >>> W = 0.7753623Error in if ( is.na(x$p.table)) { : argument is of
> length zero
>
> big.app.mat
> <-matrix(c(1,3,4,2,6,5,2,4,3,1,5,6,3,2,5,1,5,42,3,5,3,6,7,9,9,8,7),nrow=3,byrow=TRUE)
> print(kendall.w(big.app.mat[-1,]))
> >>>Kendall's W for ordinal data
> >>>W = 0.4568966 p(X2[8]) = 0.5035488
>
> Why is that working for the big matrix and not for the little one?
>
> Thanks,
> Bianca
>
> ______________________________________________
> R-help@stat.math.ethz.ch mailing list
> https://stat.ethz.ch/mailman/listinfo/r-help
> PLEASE do read the posting guide
> http://www.R-project.org/posting-guide.html
> and provide commented, minimal, self-contained, reproducible code.
>

-- 
=================================
David Barron
Said Business School
University of Oxford
Park End Street
Oxford OX1 1HP


-- 
=================================
David Barron
Said Business School
University of Oxford
Park End Street
Oxford OX1 1HP

	[[alternative HTML version deleted]]

______________________________________________
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
Received on Thu Sep 14 03:59:12 2006

Archive maintained by Robert King, hosted by the discipline of statistics at the University of Newcastle, Australia.
Archive generated by hypermail 2.1.8, at Wed 13 Sep 2006 - 18:30:05 GMT.

Mailing list information is available at https://stat.ethz.ch/mailman/listinfo/r-help. Please read the posting guide before posting to the list.