simulates track locations from the joint precision matrix of a ssm model fit. Currently, the joint precision of the SSM movement parameters is not included (ie. a full posterior simulation).

sim_post(x, what = "predicted", reps = 1, sim_only = FALSE)

Arguments

x

a ssm fit object with class ssm_df

what

simulate fitted or predicted locations

reps

number of replicate tracks to simulate from the ssm model fit object

sim_only

logical, do not include ssm estimated locations in output (default is FALSE)

Value

a fG_sim_post object containing the paths simulated from a ssm fit object

Examples

fit <- fit_ssm(ellie, model = "crw", time.step = 24)
#> fitting crw SSM to 1 tracks...
#> 
 pars:   1 1 0 -2.91153      
 pars:   0.38512 0.37191 -0.01692 -3.38813      
 pars:   -1.4595 -1.51236 -0.06769 -4.81792      
 pars:   -2.20058 -2.24383 -0.11527 -4.95732      
 pars:   -0.38346 -0.4056 -0.04615 -3.85456      
 pars:   -1.23664 -1.20844 -0.10327 -3.65984      
 pars:   -2.25711 -2.16871 -0.17159 -3.42693      
 pars:   -3.07869 -3.34599 -0.26131 -1.24778      
 pars:   -4.46622 -1.62752 -0.49092 0.12577      
 pars:   -3.18442 -2.86246 -0.30197 -0.80152      
 pars:   -3.40949 -3.30584 -0.49915 -0.40195      
 pars:   -3.06904 -3.1641 -1.05456 -0.36592      
 pars:   -3.55185 -2.81433 -1.34577 -0.44032      
 pars:   -3.30075 -3.01959 -1.10393 -0.39264      
 pars:   -3.27829 -3.03843 -1.38115 -0.39634      
 pars:   -3.2318 -2.97181 -2.21587 -0.40477      
 pars:   -3.28781 -3.0824 -5.29383 -0.34231      
 pars:   -3.30223 -3.07095 -5.32168 -0.38301      
 pars:   -3.30223 -3.07095 -5.32168 -0.38301      
psim <- sim_post(fit, "p", reps = 10)
plot(psim, type = "lines")