example2: This is just an example on how to use functions from...

In MobilePhoneESSnetBigData/destim: R package for mobile devices position estimation

Description

This is just an example on how to use functions from destim.

Usage

example2()

Examples

#' #Towers position
data(towers)
# Function S for Tennekes's model
S <- function(x) if (x > 5) return(0) else return(20*log(5/x))
# Complete events matrix
E <- createEM(c(20,20), towers, S)
# Load detection events
data(events)
# Probabilities of transition from maximum
#likelihood
lhood <- function(delta) logLik(createTM(c(20,20),
           mask = matrix(c((pi/4)*delta^2, delta,
           (pi/4)*delta^2, delta, 1 - 4*delta - pi * delta^2,
           delta, (pi/4)*delta^2, delta, (pi/4)*delta^2),
           ncol = 3)), E, events)
delta <- optimize(lhood, interval = c(0, 0.1))$minimum
#Complete transition matrix
P <- createTM(c(20,20),
           mask = matrix(c((pi/4)*delta^2, delta,
           (pi/4)*delta^2, delta, 1 - 4*delta - pi * delta^2,
           delta, (pi/4)*delta^2, delta, (pi/4)*delta^2),
           ncol = 3))
# Estimate observed states (no need)
OS <- ostates(P, E, events)
# Estimate filtered states
FS <- fstates(P, E, events)
# Estimate smooth states
SS <- sstates(P, E, events, FS)
# Load some required packages
library(raster)
library(RColorBrewer)
library(animation)
library(ggplot2)
# Create the animations
# Observed states
saveGIF({
 pal <- colorRampPalette(c("#00000000","#000000FF"), alpha = TRUE)
 ani.options(interval = 0.02)
 for (i in 1:198) {
   plot(raster(cbind(matrix(0,ncol = 13, nrow = 20), matrix(1,ncol=1,nrow=20), matrix(0,ncol=6,nrow=20))), breaks = c(0,0.5,1), col = c("white","red"), legend = FALSE)
   plot(raster(matrix(OS[,i], ncol = 20)), zlim = c(0,1), col = pal(100), add = TRUE)
   ani.pause()
 }},  movie.name = 'obsest.gif')
# Filtered states
saveGIF({
 pal <- colorRampPalette(c("#00000000","#000000FF"), alpha = TRUE)
 ani.options(interval = 0.02)
 for (i in 1:198) {
   plot(raster(cbind(matrix(0,ncol = 13, nrow = 20), matrix(1,ncol=1,nrow=20), matrix(0,ncol=6,nrow=20))), breaks = c(0,0.5,1), col = c("white","red"), legend = FALSE)
   plot(raster(matrix(FS[,i], ncol = 20)), zlim = c(0,1), col = pal(100), add = TRUE)
   ani.pause()
 }},  movie.name = 'filteredest.gif')
# Smooth states
saveGIF({
 pal <- colorRampPalette(c("#00000000","#000000FF"), alpha = TRUE)
 ani.options(interval = 0.02)
 for (i in 1:198) {
   plot(raster(cbind(matrix(0,ncol = 13, nrow = 20), matrix(1,ncol=1,nrow=20), matrix(0,ncol=6,nrow=20))), breaks = c(0,0.5,1), col = c("white","red"), legend = FALSE)
   plot(raster(matrix(SS[,i], ncol = 20)), zlim = c(0,1), col = pal(100), add = TRUE)
   ani.pause()
 }},  movie.name = 'smoothest.gif')
# The matrix GRID relates the states with coordinates
GRID <- matrix(1, nrow = 2, ncol = 400)
GRID[1,] <- rep(1:20,20)
GRID[2,] <- rep(1:20, each = 20)
# Calculate square distance mean
fdist <- sapply (1:198, function (x)
          sum(apply(GRID - matrix(c(x %/% 10 + 1,14), nrow = 2, ncol = 400), 2,norm,type = "2") * FS[,x]))
sdist <- sapply (1:198, function (x)
          sum(apply(GRID - matrix(c(x %/% 10 + 1,14), nrow = 2, ncol = 400),2,norm,type = "2") * SS[,x]))
dists <- data.frame(T = 1:198, fdist = fdist, sdist = sdist)
ggplot() +
geom_line(data = dists, aes(x = T, y = sdist), colour = "green") +
geom_line(data = dists, aes(x = T, y = fdist), colour = "red")

MobilePhoneESSnetBigData/destim documentation built on Dec. 7, 2020, 7:35 p.m.