我有一个三维数据结构,反映特定经度,纬度和深度的数据。我想将函数应用于此数据。通常,假设我想找到深度平均值,我会执行以下操作:
apply(MyData, MAR = c(1, 2), mean)
这对我有意义。我正在努力的是,我想要应用一个取决于经度和纬度的功能。有没有办法让应用将元素索引传递给函数?
答案 0 :(得分:5)
我认为你想使用outer()
并利用词汇范围
这样你就不必将myData传递给函数了
以经度和纬度呼唤:
myData <- read.table(...) # or whatever
outer(seq.int(dim(mydata)[1]),
seq.int(dim(mydata)[2]),
function(longitude,lattitude){
do things that depend on
myData[longitude,lattitude,]
})