This package provides the ability to draw joyplot or ridgeline plots in Stata. It is based on the Joyplot Guide that I released in October 2020.
The package can be installed via SSC or GitHub. The GitHub version, might be more recent due to bug fixes, feature updates etc, and may contain syntax improvements and changes in default values. See version numbers below. Eventually the GitHub version is published on SSC.
The package (v1.2) is available on SSC and can be installed as follows:
ssc install joyplot, replace
Or it can be installed from GitHub (v1.4):
net install joyplot, from("https://raw.githubusercontent.com/asjadnaqvi/stata-joyplot/main/installation/") replace
The palettes
package is required to run this command:
ssc install palettes, replace
ssc install colrspace, replace
Even if you have the package installed, make sure that it is updated ado update, update
.
If you want to make a clean figure, then it is advisable to load a clean scheme. These are several available and I personally use the following:
ssc install schemepack, replace
set scheme white_tableau
You can also push the scheme directly into the graph using the scheme(schemename)
option. See the help file for details or the example below.
I also prefer narrow fonts in figures with long labels. You can change this as follows:
graph set window fontface "Arial Narrow"
The syntax for v1.4 is as follows:
joyplot y [x] [if] [in], over(variable)
[ overlap(num) bwidth(num) color(str) alpha(num) normglobal lines lwidth(num) lcolor(str)
xangle(str) xlabsize(num) xlabcolor(str) xticks(str)
ylabsize(num) ylabcolor(str) ylabposition(str) offset(num)
yline(str) ylinecolor(str) ylinewith(str) ylinepattern(str)
xtitle(str) ytitle(str) yreverse xreverse
title(str) subtitle(str) note(str) xsize(num) ysize(num) scheme(str) name(str) ]
See the help file help joyplot
for details.
The most basic use is as follows:
joyplot y, over(variable)
or
joyplot y x, over(variable)
where y
is the variable we want to plot, and x
is usually the time dimension. The over
variable splits the data into different groupings that also determines the colors. The color schemes can be modified using the color()
option. Here any scheme from the colorpalettes
package can be used.
Set up the data:
clear
set scheme white_tableau
graph set window fontface "Arial Narrow"
use "https://github.com/asjadnaqvi/The-Stata-Guide/blob/master/data/OWID_data.dta?raw=true", clear
keep if group10==1
keep country date new_cases
We can generate basic graphs as follows:
joyplot new_cases date if date > 22460, over(country)
joyplot new_cases date if date > 22460, over(country) yline
joyplot new_cases date if date > 22460, over(country) lc(black) color(white) alpha(100)
joyplot new_cases date if date > 22460, over(country) lc(white) color(black) alpha(50) lw(0.05)
joyplot new_cases date if date > 22460, over(country) lines lw(0.2)
joyplot new_cases date if date > 22460, over(country) lines lw(0.2) color(black)
joyplot new_cases date if date > 22460, over(country) lines lw(0.2) color(black)
joyplot new_cases date if date > 22460, over(country) lines lw(0.2) ylabpos(right)
joyplot new_cases date if date > 22460, over(country) lw(0.2) ylabpos(right) xsize(5) ysize(7)
joyplot new_cases date if date > 22460, over(country) yrev
joyplot new_cases date if date > 22460, over(country) yrev xrev offset(20)
joyplot new_cases date if date > 22460, over(country) normg
joyplot new_cases date if date > 22460, over(country) normg overlap(15) xangle(45)
joyplot new_cases date if date > 22460, over(country) normg overlap(15) xangle(45) lines
joyplot new_cases date if date > 22460, over(country) bwid(0.1) off(-20) overlap(10) lw(none)
We can also customize dates, increase the overlap of the layers, change the palette, and change the color intensity:
qui summ date if date > 22460
local xmin = r(min)
local xmax = r(max)
joyplot new_cases date if date > 22460, over(country) overlap(8) color(CET C1) alpha(100) ///
lc(white) lw(0.2) xticks(`xmin'(60)`xmax') off(-30) ///
xtitle("Date") ytitle("Countries") ///
title("{fontface Arial Bold:My joyplot}") subtitle("Some more text here") ///
note("Some text here", size(vsmall))
Next we modify the scheme and make sure the colors are passed correctly. We use neon
from schemepack which has a black background:
qui summ date if date > 22460
local xmin = r(min)
local xmax = r(max)
joyplot new_cases date if date > 22460, over(country) overlap(8) color(CET C1) alpha(90) ///
lc(black) lw(0.1) xticks(`xmin'(60)`xmax') off(-30) ///
ylabc(white) xlabc(white) ///
xtitle("Date") ytitle("Countries") ///
title("{fontface Arial Bold:My joyplot}") subtitle("a subtitle here", color(white)) ///
note("Some text here", size(vsmall)) scheme(neon)
qui summ date if date > 22425
local xmin = r(min)
local xmax = r(max)
joyplot new_cases date if date > 22425, over(country) overlap(8) color(black) alpha(100) bwid(0.1) ///
lc(white) lw(0.2) xticks(`xmin' `xmax') off(+20) ///
ylabc(none) xlabc(none) xangle(0) ///
xtitle("") ytitle("") ///
title("{fontface Arial Bold:The Joy Division look}") scheme(neon)
Load the data that contains average USA state-level monthly temperatures for the period 1991-2020:
use "https://github.com/asjadnaqvi/The-Stata-Guide/blob/master/data/us_meantemp.dta?raw=true", clear
lab de month 1 "Jan" 2 "Feb" 3 "Mar" 4 "Apr" 5 "May" 6 "Jun" 7 "Jul" 8 "Aug" 9 "Sep" 10 "Oct" 11 "Nov" 12 "Dec", replace
lab val month month
joyplot meantemp, over(month) yline
joyplot meantemp, over(month) yline yrev
joyplot meantemp, over(month) yline ylw(0.2) ylc(blue) ylp(dot) ylabpos(right)
joyplot meantemp, over(month) bwid(1.5) xlabs(3) ylabs(3) overlap(3) yline normg yrev ///
ytitle("Month") xtitle("degrees Centigrade") ///
title("Mean average temperature in the USA") subtitle("2009-2020 average") ///
note("Source: World Bank Climate Change Knowledge Portal (CCKP).", size(vsmall)) ///
xsize(3) ysize(5)
qui summ meantemp
local xmin = r(min)
local xmax = r(max)
joyplot meantemp, over(month) bwid(1.5) xlabs(3) ylabs(3) overlap(3) yline normg color(cividis) alpha(100) ///
ytitle("Month") xtitle("degrees Centigrade") xticks(`xmin'(10)`xmax') ///
title("Mean average temperature in the USA") subtitle("2009-2020 average") ///
note("Source: World Bank Climate Change Knowledge Portal (CCKP).", size(vsmall)) ///
xsize(3) ysize(5)
Please open an issue to report errors, feature enhancements, and/or other requests.
v1.4 (26 Apr 2022)
- option to reverse x and y axes added
- option save graph names added
- optimization to x-axis scales
- various bug fixes
v1.3 (24 Apr 2022)
- option to add y-axis lines.
- the syntax is now allows for variable level densities.
- label placement optimized.
v1.21 (15 Apr 2022)
- y label position option added. It can now be position on the left or right.
- xsize and ysize options added to change the graph dimensions.
v1.2 (14 Apr 2022)
- x-axis angle option added.
- Global normalization option added.
- Draw lines only option added.
- Major fixes to syntax, defaults values, and how lines and areas are drawn.
v1.1 (08 Apr 2022)
- Public release. Several options and features added.
v1.0 (13 Dec 2021)
- Beta version.