Nothing
################################################################################
#
# iprior: Linear Regression using I-priors
# Copyright (C) 2018 Haziq Jamil
#
# This program is free software: you can redistribute it and/or modify
# it under the terms of the GNU General Public License as published by
# the Free Software Foundation, either version 3 of the License, or
# (at your option) any later version.
#
# This program is distributed in the hope that it will be useful,
# but WITHOUT ANY WARRANTY; without even the implied warranty of
# MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
# GNU General Public License for more details.
#
# You should have received a copy of the GNU General Public License
# along with this program. If not, see <http://www.gnu.org/licenses/>.
#
################################################################################
iprior_fixed <- function(mod) {
# When all that is required is to obtain the posterior estimate of the
# regression function (i.e. no estimation of any hyperparameters), then this
# estimatation method does just that.
#
# Args: An ipriorKernel object (mod)
#
# Returns: A list containing the optimised theta and parameters, y.hat and w,
# loglik values, standard errors, number of iterations, time taken, and
# convergence information.
w <- loglik <- NULL
start.time <- Sys.time()
if (is.ipriorKernel_nys(mod)) {
loglik_nystrom(mod$thetal$theta, mod, trace = TRUE, env = environment())
} else {
loglik_iprior(mod$thetal$theta, mod, trace = TRUE, env = environment())
}
end.time <- Sys.time()
time.taken <- as.time(end.time - start.time)
# Calculate fitted values ---------------------------------------------------
w <- get_w(u, V, Vy.inv.y, psi) # NOTE: psi, V, u and Vy.inv.y have been
y.hat <- get_y.hat(u, V, w) # assigned to environment by loglik
param.full <- theta_to_collapsed_param(mod$thetal$theta, mod)
list(theta = mod$thetal$theta, param.full = param.full, se = NA,
loglik = loglik, w = w, y.hat = y.hat, niter = NA,
start.time = start.time, end.time = end.time, time = time.taken,
convergence = NA, message = NULL)
}
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.