mean = c(0.3533024,1.5799308,0.8018807)

n =100

m = data.frame()
v=data.frame()
for(i in 1:3){
  method= rep(i,n)
  value = rnorm(mean[i],0.5,n)

  m = rbind(m,method)
  v = rbind(v,value)
}

df = cbind(m,v)
colnames(df) = c("method", "value")

op <- par(no.readonly=TRUE)
par(bty = "l",mar = rep(3,4))
boxplot(value ~ method, data = df,
        col = "lightgray", pch = 4)