This sounds like homework. Homework is discouraged on this list (but you might get lucky).
Cheers, Bert Bert Gunter Genentech Nonclinical Biostatistics (650) 467-7374 "Data is not information. Information is not knowledge. And knowledge is certainly not wisdom." Clifford Stoll On Tue, Mar 31, 2015 at 9:08 AM, Desta Yoseph via R-help <r-help@r-project.org> wrote: > 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. ______________________________________________ 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.