So I replaced the replic
function with following, but didn’t work,
# with parallel
function replic_par(R, N1, N0, B, N1_b, N0_b, seedn)
Random.seed!(seedn)
R_means = SharedArray{Float64}(R)
R_bootmeans = SharedArray{Float64}(R)
R_bootvars = SharedArray{Float64}(R)
@sync @distributed for i in 1:R
treated, control = create_data(N1, N0) # create the data
matched = match_tc(treated, control) # match and create matched data
ATET = mean(matched[:, 1] - matched[:, 2]) # estimate
R_means[i] = ATET
# Bootstrap replication is 100
ATET_b = boot_ATET(B, treated, control, N1_b, N0_b)
boot_mean = mean(ATET_b)
R_bootmeans[i] = boot_mean
boot_var = sum((ATET_b .- ATET).^2)/B
R_bootvars[i] = boot_var
end
R_means, R_bootmeans, R_bootvars
end
thanks for all the help anyways,