plotTraceback: Plot deformation traceback as a function of iterations

In guiludwig/bsplinedef: Spatial Deformation via Tensor Product of B-splines

Description

This function is a debugging tool.

Usage

plotTraceback(model, sleep = 2, ...)

Examples

# Example using artificially generated data
# Example using artificially generated data
set.seed(1)
m <- 10
x1 <- (0:m)/m
x2 <- (0:m)/m
x <- as.matrix(expand.grid(x1,x2))
n <- nrow(x)
F1 <- function(x1,x2, a = 2.5, b = 1.0) {
x <- x1 - 0.5; y <- x2 - 0.5
angle <- a*exp(-(x*x+y*y)/(b*b)) + 3*pi/2
return(cos(angle)*x + sin(angle)*y + 0.5)
}
F2 <- function(x1,x2, a = 2.5, b = 1.0) {
x <- x1 - 0.5; y <- x2 - 0.5
angle <- a*exp(-(x*x+y*y)/(b*b)) + 3*pi/2
return(-sin(angle)*x + cos(angle)*y + 0.5)
}
TIME <- 20
covModel <- RMexp(var = 1, scale = .25) + RMnugget(var = 1) # Independent in time
data <- RFsimulate(covModel, x = F1(x[,1],x[,2]), y = F2(x[,1],x[,2]), T = 1:TIME)
y <- as.numeric(unlist(data@data))
covModelM <- RMexp(var = NA, scale = NA) + RMnugget(var = NA)
test.def <- bdef(x, y, tim = 1:TIME, cov.model = covModelM, maxit = 10)
# Estimated deformation
plotGrid(test.def)
plotTraceback(test.def, margins = TRUE)

guiludwig/bsplinedef documentation built on May 16, 2020, 10:24 p.m.