Nothing
BM_gaussian_covariates <- setRefClass("BM_gaussian_covariates",
contains = "scalar_model_with_covariates",
methods = list(
initialize = function(membership_type,adj,covariates,...)
{
if(class(covariates)=="matrix")
{
covariates_list <- list(covariates)
}
else
{
covariates_list <- covariates
}
.self$initFields(membership_name = membership_type,
model_name = "gaussian_covariates",
adj = adj,
covariates = covariates_list,
...)
.self$postinit()
},
plot_parameters = function(Q)
{
matrixplot(.self$plot_transform(model_parameters[[Q]]$mu))
},
prediction = function(Q)
{
B <- matrix(0,nrow(adj),ncol(adj))
for(k in 1:length(covariates))
{
B <- B + model_parameters[[Q]]$beta[k] * covariates[[k]]
}
if(membership_name=='LBM')
{
return(
(
memberships[[Q]]$Z1
%*%
model_parameters[[Q]]$mu
%*%
t(memberships[[Q]]$Z2)
) + B
)
}
else
{
return(
(
memberships[[Q]]$Z
%*%
model_parameters[[Q]]$mu
%*%
t(memberships[[Q]]$Z)
) + B
)
}
}
)
)
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.