Nothing
recurse <- function(fy,tpm,ispd,lns)
{
#
# Function recurse to calculate the ``recursive probabilities'',
# given the parameters theta, and the observations y.
#
# Set a bunch of constants:
K <- nrow(tpm)
K2 <- K*K
L <- ncol(fy)
M <- K*L
nreps <- length(lns)
nxi <- L - nreps
N <- K*K*nxi
epsilon <- sqrt(.Machine$double.eps)
if(inherits(ispd,"matrix")) {
if(ncol(ispd) != nreps)
stop("Number of columns of \"ispd\" must equal \"nreps\".\n")
cis <- 0
nis <- nreps
} else {
cis <- 1
nis <- 1
}
# Recursive probabilities:
rp <- .Fortran(
'recurse',
fy=as.double(fy),
xispd=as.double(ispd),
tpm=as.double(tpm),
nreps=as.integer(nreps),
epsilon=as.double(epsilon),
lns=as.integer(lns),
nstate=as.integer(K),
nis=as.integer(nis),
cis=as.integer(cis),
wrk=double(K2),
xlc=double(L),
ntot=as.integer(L),
nxi=as.integer(nxi),
alpha=double(M),
beta=double(M),
gamma=double(M),
xi=double(N),
xisum=double(K2),
PACKAGE="hmm.discnp"
)
list(gamma=matrix(rp$gamma,nrow=K),xi=matrix(rp$xisum,nrow=K),llc=rp$xlc)
}
Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.