old.metropolis function

Older versions of solar location estimation

Some deprecated functions, originally used purely for light level estimation before the sampling algorithm was generalized for satellite models as well.

mkElevationSeg(segments, day)

mkNLPosterior(segments, day, light, calib)

old.dist.gc(x1, x2 = NULL)

old.find.init(mask, nseg, nlpost, pars = c("Lon", "Lat", "k"))

old.metropolis(nlpost, lookup, p0, cov0, start, end, iter = 1000, step = 100)

old.mkLookup(x, binArray = TRUE)

k.prior(seg, ps)

Arguments

• segments: vector identifying the segment of each time and light value

• day: date-time values in POSIXct

• light: vector of light data

• calib: calibration function for light levels

• x1: matrix of track locations

• x2: matrix of track locations (optional second part)

• mask: image object of masked areas

• nseg: number of (twilight) segments

• nlpost: negative log posterior function

• pars: names of parameters

• lookup: lookup function for masked areas

• p0: initial locations for sampler

• cov0: covariance matrix for sampler

• start: known start parameters

• end: known end parameters

• iter: number of iterations

• step: number of thinning iterations per iter

• x: image-like object of matrix or array of binary masks

• binArray: logical: are the masks compressed into bits?

• seg: segment

• ps: light attenuation value

Details

These functions are included for legacy purposes, this was the original implementation.

Returns

If it is a LIST, use

Author(s)

Michael D. Sumner

See Also

Please use the more up to date function metropolis, with the models such as solar.model or satellite.model.