Nothing
store.dht.results <- function(results, dht.results, i, clusters, data, obs.tab){
# Calculate the number of observations with missing distances
strata.names <- dimnames(results$individuals$summary)[[1]]
strata.names.ND <- as.character(dht.results$individual$N$Label)
for(strat in seq(along = strata.names)){
# Get the number of missed dists in strata
if(strata.names[strat] == "Total"){
object.ids <- obs.tab$object
}else{
object.ids <- obs.tab$object[obs.tab$Region.Label == strata.names[strat]]
}
strata.data <- data[data$object %in% object.ids,]
data.miss.dists <- strata.data[is.na(strata.data),]
if(clusters){
# Need to sum cluster sizes
results$individuals$summary[strat,c("n.miss.dist"),i] <- sum(data.miss.dists$size)
results$clusters$summary[strat,c("n.miss.dist"),i] <- nrow(data.miss.dists)
}else{
results$individuals$summary[strat,c("n.miss.dist"),i] <- nrow(data.miss.dists)
}
results$individuals$summary[strat,c("Area", "CoveredArea", "Effort", "n", "ER", "se.ER", "cv.ER"),i] <- as.matrix(dht.results$individuals$summary[dht.results$individual$summary$Region == strata.names[strat],c("Area", "CoveredArea", "Effort", "n", "ER", "se.ER", "cv.ER")])
results$individuals$N[strat,c("Estimate", "se", "cv", "lcl", "ucl", "df"),i] <- as.matrix(dht.results$individuals$N[dht.results$individual$N$Label == strata.names.ND[strat],c("Estimate", "se", "cv", "lcl", "ucl", "df")])
results$individuals$D[strat,c("Estimate", "se", "cv", "lcl", "ucl", "df"),i] <- as.matrix(dht.results$individuals$D[dht.results$individual$D$Label == strata.names.ND[strat],c("Estimate", "se", "cv", "lcl", "ucl", "df")])
}
if(clusters){
for(strat in seq(along = strata.names)){
results$clusters$summary[strat,c("Area", "CoveredArea", "Effort", "n", "k", "ER", "se.ER", "cv.ER"),i] <- as.matrix(dht.results$clusters$summary[dht.results$clusters$summary$Region == strata.names[strat],c("Area", "CoveredArea", "Effort", "n", "k", "ER", "se.ER", "cv.ER")])
results$clusters$N[strat,c("Estimate", "se", "cv", "lcl", "ucl", "df"),i] <- as.matrix(dht.results$clusters$N[dht.results$clusters$N$Label == strata.names.ND[strat],c("Estimate", "se", "cv", "lcl", "ucl", "df")])
results$clusters$D[strat,c("Estimate", "se", "cv", "lcl", "ucl", "df"),i] <- as.matrix(dht.results$clusters$D[dht.results$clusters$D$Label == strata.names.ND[strat],c("Estimate", "se", "cv", "lcl", "ucl", "df")])
results$expected.size[strat, c("Expected.S","se.Expected.S"), i] <- as.matrix(dht.results$Expected.S[dht.results$Expected.S$Region == strata.names.ND[strat], c("Expected.S","se.Expected.S")])[1:length(strata.names.ND[strat]),]
}
}
return(results)
}
Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.