gausskernxy:

In galuardi/hmmwoa: bones of a hidden markov model for psat data

Usage

gausskernxy(siz, xsigma, ysigma, muadv = 0)

Arguments

siz
xsigma
ysigma
muadv

Examples

##---- Should be DIRECTLY executable !! ----
##-- ==>  Define data, use random,
##--	or do  help(data=index)  for the standard data sets.

## The function is currently defined as
function (siz, xsigma, ysigma, muadv = 0) 
{
    x = 1:round(siz)
    mu = c(mean(x), mean(x)) + muadv
    fx = (matrix(exp((-0.5 * (x - mu[1])/xsigma)^2))/(sqrt(2 * 
        pi) * xsigma))
    options(digits = 5)
    fx = exp(-0.5 * ((x - mu[1])/xsigma)^2)/sqrt((2 * pi) * xsigma)
    fy = exp(-0.5 * ((x - mu[2])/ysigma)^2)/sqrt((2 * pi) * ysigma)
    fx[!is.finite(fx)] = 0
    fy[!is.finite(fy)] = 0
    kern = (fx %*% t(fy))
    kern = kern/(sum(sum(kern, na.rm = T), na.rm = T))
    kern[is.nan(kern)] = 0
    kern
  }

galuardi/hmmwoa documentation built on May 16, 2019, 5:37 p.m.