Nothing
### efront.R (2013-09-01)
###
### Efron t Statistic (2001)
###
### Copyright 2006-2013 Rainer Opgen-Rhein and Korbinian Strimmer
###
###
### This file is part of the `st' library for R and related languages.
### It is made available under the terms of the GNU General Public
### License, version 3, or at your option, any later version,
### incorporated herein by reference.
###
### 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, write to the Free
### Software Foundation, Inc., 59 Temple Place - Suite 330, Boston,
### MA 02111-1307, USA
efront.stat = function (X, L, verbose=TRUE)
{
FUN = efront.fun(L=L, verbose=verbose)
score = FUN(X)
return( score )
}
efront.fun <- function (L, verbose=TRUE)
{
if (missing(L)) stop("Class labels are missing!")
function(X)
{
tmp = centroids(X, L, lambda.var=0, lambda.freqs=0, var.groups=FALSE, verbose=verbose)
# differences between the two groups
diff = tmp$means[,1]-tmp$means[,2]
# standard error of diff
n1 = tmp$samples[1]
n2 = tmp$samples[2]
v = tmp$variances[,1] # pooled variance
sd = sqrt( (1/n1 + 1/n2)*v )
# tuning parameter
a0 <- quantile(sd, probs=c(0.9))
if (verbose) cat("Fudge factor a0 =", a0, "\n")
# t statistic
t = diff/(sd+a0)
return(t)
}
}
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.