Nothing
multi_tpfit_ml <-
function(data, coords, tolerance = pi/8, rotation = NULL, mle = "avg") {
# Estimation for matrixes of transition rates
#
# data vector of data
# coords coordinates matrix
# tolerance angle for tolerance (in radians)
# rotation vector of rotation angles (in radians)
# mle argument to pass to the function tpfit
if (!is.factor(data)) data <- as.factor(data)
if (!is.matrix(coords)) coords <- as.matrix(coords)
n <- dim(coords)[1]
nc <- dim(coords)[2]
nk <- nlevels(data)
if (n < (nk + nk^2 * nc)) stop("there are not enough data to estimate the parameters")
dire.mat <- diag(, nc)
if (!is.null(rotation)) {
if (length(rotation) != nc - 1) stop("wrong length of rotation vector; the right length is ", nc - 1)
storage.mode(rotation) <- "double"
dire.mat <- .C('rotaxes', nc = as.integer(nc), ang = as.double(rotation),
res = as.double(dire.mat), PACKAGE = "spMC")$res
dire.mat <- t(matrix(dire.mat, nc, nc))
}
res <- list()
res$coordsnames <- colnames(coords)
if (is.null(res$coordsnames)) {
res$coordsnames <- paste("X", 1:nc, sep = "")
}
res$coefficients <- apply(dire.mat, 1, function(d) {
md.coef <- list()
md.coef$coefficients <- tpfit_ml(data, coords, d, tolerance, mle)$coefficients
class(md.coef) <- "tpfit"
return(md.coef)
})
res$prop <- table(data)
res$prop <- res$prop / sum(res$prop)
res$tolerance <- as.double(tolerance)
res$rotation <- rotation
class(res) <- "multi_tpfit"
return(res)
}
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.