generate_data = function(N){ lst( n = N , pop_mean = ( # 0 rnorm(1) ) , pop_sd = ( # 1 rweibull(1,2,1) ) , y = rnorm(n,pop_mean,pop_sd) , .join_data = lst(pop_mean,pop_sd) ) }