Nothing
# Clustering classification
bacisClassification<- function(numGroup = 5,
tau1 = NA,
tau2 = .001,
phi1 = 0.1,
phi2 = 0.3,
clusterCutoff = NA,
MCNum = 50000,
nDat = c(25, 25, 25, 25, 25),
xDat = c(2, 3, 7, 6, 10),
seed = NA
)
{
if (is.na(seed))
{
set.seed( as.integer((as.double(Sys.time())*1000+Sys.getpid()) %% 2^31) )
}
else{
set.seed(seed)
}
if (is.na(tau1))
{
sd <- (logit(phi2) - logit(phi1)) / 6
tau1 <- 1 / sd / sd
cat("The value of tau1 is set at:", tau1, "\n")
}
#print(tau1)
xLim <- max(xDat / nDat) + 0.3
#print(xDat)
t <- ModelOne(
xDat = xDat,
nDat = nDat,
numGroup = numGroup,
pp1 = phi1,
pp2 = phi2,
tau1 = tau1,
tau2 = tau2,
clusterCutoff = clusterCutoff,
MCNum = MCNum
)
return(list(highResponseGroup=t$highResponseGroup, lowResponseGroup=t$lowResponseGroup))
}
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.