Simulation with cpm package

Hi,
I am running the following code based on the cpm vignette's code. I believe
the code is syntactically correct but it just seems to hang R. I can get
this to run if I set the sims to 100 but with 2000 it just hangs. Any ideas
why?

Thanks,
Chris

library(cpm)
cpmTypes <- c("Kolmogorov-Smirnov","Mann-Whitney","Cramer-von-Mises")
changeMagnitudes <- c(1, 2, 4, 5)
changeLocations <- c(50,100,300)
sims <- 2000
ARL0 <- 500
startup <- 20
results <- list()
for (cpmType in cpmTypes) {
   results[[cpmType]] <- matrix(numeric(length(changeMagnitudes) *
                                          length(changeLocations)), nrow =
length(changeMagnitudes))
   for (cm in 1:length(changeMagnitudes)) {
     for (cl in 1:length(changeLocations)) {
       print(sprintf("cpm:%s magnitude::%s location:%s",
                     cpmType, changeMagnitudes[cm], changeLocations[cl]))
       temp <- numeric(sims)
       for (s in 1:sims) {
         x <-c(rchisq(changeLocations[cl], df=3), rchisq(2000,

df=changeMagnitudes[cm]))
         temp[s] <-detectChangePoint(x, cpmType,
                                     ARL0=ARL0,
startup=startup)$detectionTime
       }
       results[[cpmType]][cm,cl] <- mean(temp[temp > changeLocations[cl]]) -
         changeLocations[cl]
     }
   } }

	[[alternative HTML version deleted]]