# intial example
data("USstates")
head(USstates@data )
statePolys <- create_map_table(USstates, 'ST')
head(statePolys)
## Not run:
# data("edPov")
#
# # basic figure 1
# lmplot(stat.data=edPov,
# map.data=statePolys,
# panel.types=c('labels', 'dot', 'dot','map'),
# panel.data=list('state','pov','ed', NA),
# ord.by='pov',
# grouping=5, median.row=TRUE,
# map.link=c('StateAb','ID'))
#
# # publication figure 1a
# lmplot(stat.data=edPov, map.data=statePolys ,
# panel.types=c('labels', 'dot', 'dot','map'),
# panel.data=list('state','pov','ed', NA),
# ord.by='pov',
# grouping=5,
# median.row=TRUE,
# map.link=c('StateAb','ID'),
#
# plot.height=9,
# colors=c('red','orange','green','blue','purple'),
# map.color2='lightgray',
#
# panel.att=list(list(1, header='States', panel.width=.8, align='left', text.size=.9),
# list(2, header='Percent Living Below \n Poverty Level',
# graph.bgcolor='lightgray', point.size=1.5,
# xaxis.ticks=list(10,15,20), xaxis.labels=list(10,15,20),
# xaxis.title='Percent'),
# list(3, header='Percent Adults With\n4+ Years of College',
# graph.bgcolor='lightgray', point.size=1.5,
# xaxis.ticks=list(0,20,30,40), xaxis.labels=list(0,20,30,40),
# xaxis.title='Percent'),
# list(4, header='Light Gray Means\nHighlighted Above',
# inactive.border.color=gray(.7), inactive.border.size=2,
# panel.width=.8)))
#
# edPov$points <- 0
#
# # publication figure 1b
# lmplot (stat.data=edPov, map.data=statePolys,
# panel.types=c('dot', 'labels', 'dot', 'dot', 'map'),
# panel.data=list('points', 'state', 'pov', 'ed', NA),
# map.link=c('StateAb','ID'),
# ord.by='pov',
# grouping=5,
# median.row=TRUE,
#
# plot.height=9,
#
# colors=c('red','orange','green','blue','purple'),
# map.color2='lightgray',
#
# panel.att=list(list(1, panel.width=.15, point.type=20,
# graph.border.color='white',
# xaxis.text.display=FALSE, xaxis.line.display=FALSE,
# graph.grid.major=FALSE),
#
# list(2, header='States', panel.width=.8,
# align='left', text.size=.9),
#
# list(3, header='Percent Living Below\nPoverty Level',
# graph.bgcolor='lightgray', point.size=1.5,
# xaxis.ticks=list(10,15,20), xaxis.labels=list(10,15,20),
# xaxis.title='Percent'),
#
# list(4, header='Percent Adults With\n4+ Years of College',
# graph.bgcolor='lightgray', point.size=1.5,
# xaxis.ticks=list(20,30,40), xaxis.labels=list(20,30,40),
# xaxis.title='Percent'),
#
# list(5, header='Light Gray Means\nHighlighted Above',
# inactive.border.color=gray(.7), inactive.border.size=2,
# panel.width=.8)))
#
# # publication figure 1c
# myPlot <- lmplot(stat.data=edPov, map.data=statePolys,
# panel.types=c('map', 'dot', 'labels', 'dot', 'dot'),
# panel.data=list(NA, 'points', 'state', 'pov', 'ed'),
# map.link=c('StateAb','ID'),
# ord.by='pov',
# grouping=5,
# median.row=TRUE,
#
# plot.height=9,
#
# colors=c('red','orange','green','blue','purple'),
# map.color2='lightgray',
#
# panel.att=list(list(2, panel.width=.15, point.type=20,
# graph.border.color='white',
# xaxis.text.display=FALSE, xaxis.line.display=FALSE,
# graph.grid.major=FALSE),
#
# list(3, header='States', panel.width=.8,
# align='left', text.size=.9),
#
# list(4, header='Percent Living Below\nPoverty Level',
# graph.bgcolor='lightgray', point.size=1.5,
# xaxis.ticks=list(10,15,20), xaxis.labels=list(10,15,20),
# xaxis.title='Percent'),
#
# list(5, header='Percent Adults With\n4+ Years of College',
# graph.bgcolor='lightgray', point.size=1.5,
# xaxis.ticks=list(20,30,40), xaxis.labels=list(20,30,40),
# xaxis.title='Percent'),
#
# list(1, header='Light Gray Means\nHighlighted Above',
# inactive.border.color=gray(.7), inactive.border.size=2,
# panel.width=.8)))
#
# print(myPlot, name='myExhibit.tiff', res=300)
# ## End(Not run)
Run the code above in your browser using DataLab