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]]
______________________________________________
R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see
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.
