Replicating graphs to show means, standard deviations and distributions.

[Aisha Khan]

I recently saw the below graph posted on twitter and I am wondering if the same can be replicated in Stata for my data.
I have data on farm hsoueholds for three-time periods, a total of around 7000 observations overall. Some of these households access credit and some don’t. I have credit as a dummy variable that takes the value 1 if a houhold access credit and 0 otherwise. I want to see crop yields of each households based on whether or not they access credit.
Can I replicate all three graphs for households that access credit and those that don’t??

[Nick Cox]

You can get close with twoway bar using ci results; graph box; and dotplot.

[Aisha Khan]

Nick Cox Thank you for your reply.
Following your suggestion I use my data to generate the graphs as follows:

* Example generated by -dataex-. For more info, type help    dataex
clear
input double Household_No float Survey_Round str12 Credit float Yield
1 1 "Non-adopters" 4118.3335
1 2 "Non-adopters"    6177.5
1 3 "Adopters"      8086.909
2 3 "Adopters"        9142.7
3 1 "Non-adopters"     741.3
3 2 "Non-adopters"  6816.552
4 1 "Non-adopters"  2671.351
4 3 "Adopters"      6678.378
5 1 "Non-adopters" 1317.8667
5 2 "Non-adopters"  4612.533
5 2 "Non-adopters"  1372.778
5 3 "Adopters"          7413
9 1 "Non-adopters"  4561.846
10 1 "Non-adopters" 1560.6316
12 2 "Non-adopters"   3162.88
12 3 "Adopters"      5814.118
14 1 "Non-adopters" 1900.7693
14 1 "Non-adopters" 1900.7693
14 2 "Non-adopters"      3530
14 2 "Non-adopters" 4706.6665
14 2 "Non-adopters"      4942
14 2 "Non-adopters"      4942
14 2 "Non-adopters" 4667.4443
14 2 "Non-adopters"      3530
14 3 "Adopters"        5683.3
14 3 "Adopters"        5683.3
15 1 "Non-adopters" 2635.7334
15 2 "Non-adopters"  2534.359
15 2 "Non-adopters"    5930.4
15 2 "Non-adopters" 2196.4443
15 3 "Adopters"      8086.909
17 3 "Adopters"      5068.718
20 1 "Non-adopters"     741.3
20 2 "Adopters"        4800.8
25 1 "Non-adopters"    6177.5
25 2 "Non-adopters"    3953.6
25 2 "Non-adopters"  5702.308
25 2 "Non-adopters"  5765.667
27 1 "Non-adopters"      3530
27 2 "Non-adopters"    4447.8
27 2 "Non-adopters"  3154.468
27 2 "Non-adopters" 2208.1277
27 2 "Non-adopters"  2280.923
27 2 "Non-adopters"    8401.4
27 2 "Non-adopters" 2353.3333
27 2 "Non-adopters" 3801.5386
27 2 "Non-adopters"  2280.923
27 2 "Non-adopters"    4447.8
27 2 "Non-adopters"    3706.5
27 3 "Non-adopters"    1482.6
27 3 "Non-adopters"  946.3405
27 3 "Non-adopters"  950.3846
28 1 "Adopters"      4007.027
28 1 "Non-adopters"  3294.667
28 1 "Non-adopters"  5491.111
28 1 "Non-adopters" 2635.7334
28 2 "Non-adopters" 2196.4443
28 2 "Non-adopters" 3184.8445
28 2 "Non-adopters" 2196.4443
28 2 "Non-adopters"   4324.25
28 2 "Non-adopters" 3205.6216
28 2 "Non-adopters"  2671.351
28 2 "Non-adopters"      2471
28 2 "Non-adopters"  3075.022
28 2 "Non-adopters" 2270.6487
28 2 "Non-adopters"  2137.081
28 3 "Non-adopters"    6036.3
28 3 "Non-adopters"    3953.6
28 3 "Non-adopters" 4173.2446
31 1 "Non-adopters"  2280.923
33 1 "Non-adopters"     988.4
33 1 "Non-adopters"    3953.6
33 1 "Non-adopters"      2471
33 1 "Non-adopters"    1235.5
33 2 "Non-adopters" 3801.5386
33 2 "Non-adopters"    3706.5
33 2 "Non-adopters"  4633.125
33 2 "Non-adopters"      4942
33 2 "Non-adopters"      4942
33 2 "Non-adopters"    3706.5
33 3 "Non-adopters"      4942
33 3 "Adopters"        5930.4
33 3 "Non-adopters"    3953.6
33 3 "Non-adopters"    3953.6
33 3 "Non-adopters"   6325.76
36 1 "Non-adopters" 4392.8887
36 1 "Non-adopters"  5271.467
36 1 "Non-adopters" 2196.4443
36 1 "Non-adopters"      2471
36 2 "Non-adopters" 1098.2222
36 2 "Non-adopters" 4392.8887
36 2 "Non-adopters"  5491.111
36 3 "Non-adopters"  6589.333
36 3 "Non-adopters"    5930.4
37 1 "Non-adopters"  7282.947
37 2 "Non-adopters"  6415.096
37 2 "Non-adopters"      7413
37 3 "Adopters"       3688.06
37 3 "Adopters"      4492.727
37 3 "Adopters"          4942
end

[CODE
]***generate Yields for adopters and non-adopters
gen Yield_0=Yield if Credit=="Non-adopters"
gen Yield_1= Yield if Credit=="Adopters"
mvencode Yield_0 Yield_1, mv(0) override

***** Graph bar
graph bar Yield_0 Yield_1, percent scheme(s2mono) ytitle("Percent") graphr(c(white) lc(white)) plotregion(lcolor(black)) scale(1) legend(label(1 "Non-adopters") label(2 "Adopters ")) graphregion(color(white)) bgcolor(white) name(f1, replace)blabel(total, format(%9.2f))
***Graph box
graph box Yield , over(Credit) name(f2, replace)
***dotplot
dotplot Yield, over(Credit) center scheme(s2mono) ytitle("Yield (Kg/ha)") graphr(c(white) lc(white)) plotregion(lcolor(black)) scale(1) graphregion(color(white)) bgcolor(white) name(f3, replace)
graph combine f1 f2 f3, row(1)
*combine graphs
graph combine f1 f2 f3, row(1)
[/CODE]

I get the below graph. How can I include the error bars in the first graph?