[R] Calculating Kendall's tau
Desta Yoseph
desta_yo at yahoo.com
Tue Mar 31 18:08:56 CEST 2015
I am analyzing trend using Mann-kendall test for 31 independent sample, each sample have 34 years dataset. I supposed to find Kendall “tau” for each sample. The data is arranged in column wise (I attached the data).To find Kendall tau, I wrote R script as:
desta<-read.csv("rainfall.csv", header=T, sep=",") require(Kendall) MK<-function(y) { nc<-ncol(y) MannKendalltau<- numeric(nc) for(i in 2:nc){ MannKendalltau[i]<-MannKendall(y[,i]) } MannKendalltau } MK(desta)
The displayed result showed both “tau” and “2-sided p-value”in unorganized way. But, I want only “tau” value that is presented in organized manner. Anyone can tell me how can I get orderly displayed “tau” value? here is my sample result: [[1]][1] 0
[[2]][1] 0.4352941attr(,"Csingle")[1] TRUE
[[3]][1] 0.5462185attr(,"Csingle")[1] TRUE
[[4]][1] 0.4218487attr(,"Csingle")[1] TRUE....Thank you for your guidance
[[alternative HTML version deleted]]
More information about the R-help
mailing list