x <- rbind(cbind(rnorm(10, mean=0), rnorm(10, mean=0,), rnorm(10, mean=0)),
cbind(rnorm(10, mean=10), rnorm(10, mean=10), rnorm(10, mean=0)),
cbind(rnorm(10, mean=10), rnorm(10, mean=0), rnorm(10, mean=10)),
cbind(rnorm(10, mean=10), rnorm(10, mean=10), rnorm(10, mean=10))
)
res <- ockc(x, k=4, nboot=4, order=c(1:10, 21:40, 11:20))
res
Run the code above in your browser using DataLab