## Not run:
# # simulate mediational data
# a <- .5
# b <- .6
# t_prime <- .3
#
# X <- rnorm(50,0,1)
# M <- a*X + rnorm(50,0,1)
# Y <- t_prime*X + b*M + rnorm(50,0,1)
#
# # save jzs_med output
# res <- jzs_med(independent=X,dependent=Y,mediator=M)
#
# # plot results
# plot(res$main_result)
# ## End(Not run)
Run the code above in your browser using DataLab