gausskern:

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

Usage

gausskern(siz, sigma, muadv = 0)

Arguments

siz
sigma
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, sigma, muadv = 0) 
{
    x = 1:round(siz)
    mu = c(mean(x), mean(x)) + muadv
    fx = (matrix(exp((-0.5 * (x - mu[1])/sigma)^2))/(sqrt(2 * 
        pi) * sigma))
    options(digits = 5)
    fx = exp(-0.5 * ((x - mu[1])/sigma)^2)/sqrt((2 * pi) * sigma)
    fy = exp(-0.5 * ((x - mu[2])/sigma)^2)/sqrt((2 * pi) * sigma)
    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.