Displaying 1 result from an estimated 1 matches for "zredef".
Did you mean:
predef
2009 Sep 23
1
Numerical integration problem
...es: distance r
and angle th); a, b0 and bt are model parameters
point.hazard <- function(pos,a,b0,bt) a * exp(-(pos$r^2/(2*b0^2))) *
exp(-(pos$th^2/(2*bt^2)))
#point.hazard for a point related to input point pos by time t and speed
v
hazard <- function(t, v, gam, pos, a, b0, bt)
{ pos2 <- zredef(pos,-t*v,gam)
pos2$th[pos2$th>pi] <- 2*pi-pos2$th[pos2$th>pi]
point.hazard(pos2,a,b0,bt)
}
#Returns a list of polar co-ordinates for a point defined by distance m
in direction gam from starting point pos
zredef <- function(pos, m, gam)
{ x <- pos$r*sin(pos$th) + m*sin(gam)
y <-...