iteratedMap_ioNoise: Time series with input and measurement noise

In PhilippVWC/myBayes: Functions for time series analysis based on Bayesian methods

Description

Time series with input and measurement noise

Usage

iteratedMap_ioNoise(x0, r, n, funct = function(x, r) myBayes::gsm_cpp(x =
  x0, r = r, alpha = 2), irfunct = identity, orfunct = identity)

Arguments

x0	double - starting value of timeseries
r	double - bifurcation parameter
n	integer - length of time series
funct	function of two scalar values - map
irfunct	function of scalar value - function for input noise
orfunct	function of scalar value - output noise function

Details

Iterated map with input and output noise

Value

the noisy time series

Examples

x0 = 0.7
r = 0.8
N = 100
y = iteratedMap_ioNoise(x0 = x0,
                        r = r,
                        n = N,
                        funct=function(x,r) myBayes::gsm_cpp(x,r,alpha = 1.7),
                        orfunct = function(u) min(1,max(0,rnorm(n = 1,
                                                                mean = u,
                                                                sd = 0.1))))
plot(x = 1:N,
     y = y,
     type = "l",
     ylim = c(0,1))

PhilippVWC/myBayes documentation built on Oct. 2, 2020, 8:25 a.m.