Nothing
# DESP/R/DESP_RML.R by A. S. Dalalyan and S. Balmand Copyright (C) 2015-
#
# This program is free software; you can redistribute it and/or modify
# it under the terms of the GNU General Public License (version 3) as published by
# the Free Software Foundation.
#
# 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.
#
# A copy of the GNU General Public License is available at
# http://www.r-project.org/Licenses/
#
DESP_RML <-
function(X,B,Theta=NULL) {
# estimation of the diagonal of the precision matrix by likelihood maximization, when the true value of B is known or has already been estimated
# the observations of the data matrix X are assumed to have zero mean
# read the sample size and the number of variables
D = dim(X);
n = D[1]; # n is the sample size
p = D[2]; # p is the dimension
# compute the sample cov matrix
if(is.null(Theta))
{
S = crossprod(X)/n;
}
else
{
S = crossprod(X - Theta %*% MASS::ginv(B))/n;
}
Phi = diag(crossprod(S,B));
Phi = Phi*(Phi>0);
return(1/Phi)
}
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.