Displaying 2 results from an estimated 2 matches for "272x976".
2012 Dec 20
4
comparison of large data set
...= as.matrix(read.ascii.grid(file1)$data)
modeldepth1[modeldepth1!=0]<-1
modeldepth2 = as.matrix(read.ascii.grid(file2)$data)
modeldepth2[modeldepth2!=0]<-1
...
modeldepth54 = as.matrix(read.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 indicat...
2012 Dec 21
0
variable names in numeric list and Bayesian inference
...lt;- sum(!obs & model)
s1/(s1 + s2 + s3)
}
modeldepths <- lapply(modeldepths, function(x) f2(observeddepth, x))
modeldepthnew <- modeldepth[modeldepth > 0.7] # Here the objective is to
pick the 31 of the 54 files (31 of the 54 elements have an f2 greater
than 0.7) with their complete 272x976 cells because that information
will then be used in the next step.
The next step will be to perform a Bayesian inference where the
modeldepthnew represents the prior likelihood.
Are there any suggestions on a specific R package to use in the Bayesian
inference steps?
Thank-you.
Irucka Embry...