Displaying 1 result from an estimated 1 matches for "f2all".
Did you mean:
t2all
2012 Dec 20
4
comparison of large data set
...ad.ascii.grid(file54)$data)
modeldepth54[modeldepth54!=0]<-1
Each modeldepth* and observeddepth is 272x976
f2 <- function(obs, dat){
obs <- as.logical(obs)
dat <- as.logical(dat)
s1 <- sum(obs & dat)
s2 <- sum(obs & !dat)
s3 <- sum(!obs & dat)
s1/(s1 + s2 + s3)
}
f2all <- f2(observeddepth, modeldepth[1:54]) # I am including the
"[1:54]" in the code to indicate that I want to compare the
observeddepth data to each of the modeldepth data files
modeldepth <- as.matrix(c(modeldepth1, modeldepth2, ... modeldepth54))
modeldepth <- modeldepth_index[...