powered by
summary(output, y, eps.diff = 1e-5)
y <- rbind(matrix(rnorm(120,0,1),ncol=4), matrix(rnorm(120,4,1), ncol=4), matrix(rnorm(120,0,1),ncol=4)) output <- parse(K = c(1:2), lambda = c(0,1), y=y) output$mu.hat.best summary(output, y)
Run the code above in your browser using DataLab