avg_param <- as.double(params_hat)
std_param <- as.double(stdev)
matplot(cbind(params[o],avg_param[o],avg_param[o]+std_param[o],avg_param[o]-std_param[o]),
avg_param <- as.double(params_hat)
std_param <- as.double(stdev)
matplot(cbind(params[o],avg_param[o],avg_param[o]+std_param[o],avg_param[o]-std_param[o]),