Monte Carlo simulation to compare two vectors, one of which has no variance

zerocor.R

library(MonteCarlo)

df <- data.frame(x=c(2.2, 2.2, 2.2, 2.2, 2.2), y=c(2.2, 1.2, 1.3, 2.4, 2.6))
df$trend <- seq(1,length(df$x))

pv <- function(noparams) {
    return(list("p.value"=cor(jitter(df$x), df$y)))
}

param_list <- list("noparams"=c(0))

MC_result<-MonteCarlo(func=pv, nrep=100000, param_list=param_list)

mean(MC_result$results$p.value)