Nothing
setMethod("Roblda",
signature(x = "IData"),
function(x, grouping, prior="proportions", CVtol=1.0e-5, egvtol=1.0e-10, subset=1:nrow(x),
# CovCase=1:4, SelCrit=c("BIC","AIC"), silent=FALSE, CovEstMet=c("Pooled","Globdev"), SngDMet=c("fasttle","fulltle"),
CovCase=1:4, SelCrit=c("BIC","AIC"), silent=FALSE, CovEstMet=c("Pooled","Globdev"), SngDMet=c("fasttle","fulltle"), k2max=1e6,
Robcontrol=RobEstControl(), ...)
{
CovEstMet <- match.arg(CovEstMet)
SelCrit <- match.arg(SelCrit)
SngDMet <- match.arg(SngDMet)
if (SngDMet=="fulltle") { stop("Roblda with fulltle not implemented yet.\n") }
limlnk2 <- log(k2max)
if (length(subset) < x@NObs)
{
x <- x[subset,]
grouping <- grouping[subset]
}
grouping <- factor(grouping,exclude=NULL)
grplvls <- levels(grouping)
n <- x@NObs
p <- 2*x@NIVar
if (p==2) CovCase <- q1CovCase(CovCase)
k <- length(grplvls)
if (length(grplvls)==1)
{
errmsg <- "The data belongs to one single group. A partition into at least two different groups is required\n"
if (silent==FALSE) {
stop(errmsg)
} else {
warning(errmsg)
return(NULL)
}
}
MxtDEst <- RobMxtDEst(x,grouping,Mxt="Hom",CovCase=CovCase,SelCrit=SelCrit,CovEstMet=CovEstMet,...)
nk <- as.numeric(table(grouping))
if (sum(nk)!=n) { stop("Dimensions of the x and grouping arguments do not agree with each other\n") }
glbmeans <- c(colMeans(x@MidP),colMeans(x@LogR))
mugdev <- scale(MxtDEst@RobNmuE,center=glbmeans,scale=FALSE)
vnames <- unlist(dimnames(MxtDEst@RobNmuE)[2])
B <- matrix(0.,nrow=p,ncol=p,dimnames=list(vnames,vnames))
for (g in 1:k) B <- B + (nk[g]/n) * outer(mugdev[g,],mugdev[g,])
selmodel <- BestModel(MxtDEst)
# Ilda(Conf=selmodel,p=p,nk=nk,prior=prior,means=MxtDEst@RobNmuE,W=MxtDEst@CovConfCases[[selmodel]]$RobSigE,B=B,egvtol=egvtol,limlnk2=limlnk2)
Slda(Conf=selmodel,p=p,nk=nk,prior=prior,means=MxtDEst@RobNmuE,W=MxtDEst@CovConfCases[[selmodel]]$RobSigE,B=B,egvtol=egvtol,limlnk2=limlnk2)
}
)
setMethod("Robqda",
signature(x = "IData"),
function(x, grouping, prior="proportions", CVtol=1.0e-5, subset=1:nrow(x),
CovCase=1:4, SelCrit=c("BIC","AIC"), silent=FALSE, SngDMet=c("fasttle","fulltle"), k2max=1e6,
Robcontrol=RobEstControl(), ...)
{
SelCrit <- match.arg(SelCrit)
SngDMet <- match.arg(SngDMet)
if (x@NIVar==1) CovCase <- q1CovCase(CovCase)
limlnk2 <- log(k2max)
if (length(subset) < x@NObs)
{
x <- x[subset,]
grouping <- grouping[subset]
}
grouping <- factor(grouping,exclude=NULL)
grplvls <- levels(grouping)
if (length(grplvls)==1)
{
errmsg <- "The data belongs to one single group. A partition into at least two different groups is required\n"
if (silent==FALSE) {
stop(errmsg)
} else {
warning(errmsg)
return(NULL)
}
}
MxtDEst <- RobMxtDEst(x,grouping,Mxt="Het",CovCase=CovCase,SelCrit=SelCrit,...)
selmodel <- BestModel(MxtDEst)
# Iqda(Conf=selmodel,p=2*x@NIVar,nk=as.numeric(table(grouping)),lev=grplvls,
Sqda(Conf=selmodel,p=2*x@NIVar,nk=as.numeric(table(grouping)),lev=grplvls,limlnk2=limlnk2,
prior=prior,means=MxtDEst@RobNmuE,Wg=coef(MxtDEst,selmodel)$Sigma)
}
)
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.