Nothing
options( width=150, max.print=1000000 )
# some seeds to make more seeds
Seeds1 <- c( 826847827, 679365679, 779140483, 146694967, 539272129 )
# --------------------------------------
SGETID <- Sys.getenv(c("SGE_TASK_ID"))
SGENID <- as.numeric(SGETID)
samp <- SGENID
n <- floor(log(samp,10))+1
oname <- paste("RDS/Sample2/results_", paste(rep(0,5-n),sep="",collapse=""), samp, ".rds",sep="")
sd <- Seeds1[2]
set.seed(sd)
seedList <- ceiling( 1000000 * runif( 1000 ) )
seed <- seedList[SGENID]
# only do the jackknife on the main data once
MJK <- ( SGENID == 1 )
# --------------------------------------
library(samon, lib.loc="../../samlib")
data("samonPANSS2")
PANSS2 <- samonPANSS2
alphaList <- -20:20
Results <- samon(
mat = PANSS2, # imput matrix
Npart = 10, # number of partitions
InitialSigmaH = 10.0, # initial value
HighSigmaH = 50.0, # high value for p
InitialSigmaF = 8.0, # initial value
HighSigmaF = 50.0, # high value for q
lb = 30, # parameters for
ub = 210, # cumulative
zeta1 = 4.0, # beta distribution
zeta2 = 7.0,
NSamples = 50, # bootstraps
seed0 = seed,
MJackknife = MJK,
SJackknife = TRUE,
alphaList = alphaList
)
# --------------------------------------
saveRDS(Results, file=oname)
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.