[R-sig-Geo] Mann-Kendall test on raster data

[Tiffany Smith]

Hello,

I'm trying to do a Mann-Kendall test on my RasterBrick data. So far, I've been able to complete the Mann-Kendall test (via the "Kendall" package) by extracting a series like so:

onelocation<-rev(extract(myraster,SpatialPoints(cbind(-86,33))))
onelocation.mk<-MannKendall(onelocation)

I need help either figuring out a way to do the MK test with raster data, or if I loop through all points and use the above commands, how can I return that data into a Raster so I can visualize the results by making maps?

Thanks,
Tiffany