Rizzo: Statistical Computing with R, p.123より
x <- seq(.1, 2.5, length = 10) m <- 10000 z <- rnorm(m) dim(x) <- length(x) p <- apply(x, MARGIN = 1, FUN = function(x, z) {mean(z < x)}, z = z) Phi <- pnorm(x) print(round(rbind(x, p, Phi), 3))
ここでベクトル x に apply 関数を適用するためには,このベクトルに次元(MARGIN=1,すなわち行属性)を持たせることが必要.