2016-05-30 51 views
0

我想計算p值並將其與相關性柵格圖 合併,但我不知道如何計算柵格p值並將其與相關性柵格圖進行比較。 ý13至24是NDVI光柵圖像如何從相關性柵格中找到p值

預先感謝您的任何幫助我的 - [R編程代碼如下所示

r.stack <- stack 
(y13,y14,y15,y16,y17,y18,y19,y20,y21,y22,y23,y24,y1,y2,y3,y4,y5,y6,y7,y8,y9,y10,y11,y12) 
z1 <- r.stack[[1:12]] 
z2 <- r.stack[[13:24]] 
z3 <- stack(z1,z2) 


r <- calc(z3, fun=function(x) cor(x[1:12], x[13:24], method='pearson')) 

plot(r,col=colorRampPalette(c("red", "yellow", "blue"))(255)) 

# I have no idea how to get P- value from the Correlation # 

而Y1至Y12是乾旱(SEPI)光柵圖像和 .. ...

+1

請分享代碼,而不是截圖,甚至不可讀。 – mtoto

+3

尋求調試幫助的問題(「爲什麼這個代碼不工作?」)必須包含所需的行爲,特定的問題或錯誤以及在問題本身中重現問題所需的最短代碼。沒有明確問題陳述的問題對其他讀者無益。請參閱:[如何創建最小,完整和可驗證示例](http://stackoverflow.com/help/mcve)。 –

+0

@mtoto ....我更正了....感謝您的建議^^ – faris

回答

2

我已經得到了答案,並希望它可以幫助其他人。 我corLocal打算corr。的代碼如下所示,

r.stack <- stack (y13,y14,y15,y16,y17,y18,y19,y20,y21,y22,y23,y24 ,y3, y4,y5,y6, y7,y8,y9,y10,y11,y12,y1,y2)

z1 <- r.stack[[1:12]] z2 <- r.stack[[13:24]] z3 <- stack(z1,z2)

r <- corLocal(r.stack[[1:12]], r.stack[[13:24]], test=TRUE) plot(r) # r and p-value map will be shown

# only correlation cells where the p-value < 0.05 will be ommited

r.mask <- mask(r[[1]], r[[2]] < 0.05, maskvalue=FALSE)

plot(r.mask)

相關問題