Nothing
esagda.pred <- function(ynew, y, ina) {
## xnew is the new observation(s)
## x is the data set
## ina is the group indicator variable
ynew <- as.matrix(ynew)
if ( ncol(ynew) == 1 ) ynew <- t(ynew)
ina <- as.numeric(ina)
g <- max(ina)
mat <- matrix(0, dim(ynew)[1], g)
for (j in 1:g) {
mod <- Directional::esag.mle( y[ina == j, ] )
mat[, j] <- Directional::desag(ynew, mod$mu, mod$gam, logden = TRUE )
}
Rfast::rowMaxs(mat)
}
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.