1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
240
241
242
243
244
245
246
|
<html><h1>errbar1</h1><img src="errbar1.gif"><p>
<a href="gall.using.html">How to download and try this example</a><br>
<pre><b>Usage: pl -gif errbar1.htm
or for postscript: pl -ps errbar1.htm | lp
In this example the means and standard deviations were already calculated,
so proc bars is used to render the bars. An alternative is to use proc
rangebar with meanmode option to calculate mean & standard deviation
and render bars.
Categories are used to position bars in X; proc catslide is used to
adjust the positioning slightly so that pairs of bars may be displayed
for each category.
<hr></b>
#set TAILLEN = 0.05
#proc Page:
pagesize: 8.5 11
#if @DEVICE in gif,png
scale: 0.4
#endif
Title: Mean FES Scores
For Families in Emergency Shelters vs.
Families in Transitional Housing
/////////////////////////////
////// Top Panel
/////////////////////////////
///////// Left bars
// specify data for top panel, left of each pair
// (must be done before areadef because areadef uses data to set up category names)
#Proc Getdata:
Data:
// E (1-5)
//categ n obs mean stdev pval
//==== ===== ===== ===== =====
coh 23 37.83 18.66 p<.05
exp 23 42.17 9.00 p<.05
conf 23 54.48 13.67 p>.05
ind 23 42.52 14.71 p>.05
achiev 23 54.48 8.44 p>.05
// set up plot area
#Proc Areadef:
areaname: 2hi
xscaletype: categories
xcategories: datafield=1
yrange: 0 100
Xaxis.stubs: text
Cohesion
Expressiveness
Conflict
Independence
Achievement\nOrientation
yaxis.stubs: inc 10
Yaxis.label: Standard Score
// do horizontal lines..
#Proc Line
linedetails: width=0.4 dashscale=5
notation: locval
points: min 50(s) max 50(s)
#Proc Line
linedetails: width=0.4 dashscale=5 style=1
notation: locval
Points: min 40(s) max 40(s)
min 60(s) max 60(s)
// do centered p values
#Proc Scatterplot
Xfield: 1
ylocation: 5(s)
labelfield: 5
// prepare to do left bars
#proc catslide
axis: x
amount: -0.1
// render points
#Proc Scatterplot:
Xfield: 1
Yfield: 3
symbol: shape=pixdiamond style=solid color=black radius=0.05
// render bars
#Proc bars
locfield: 1
lenfield: 3
errbarfields: 4
tails: @TAILLEN
////////// Right bars ////////////////
// prepare to do right bars
#proc catslide
axis: x
amount: 0.1
// specify data for right bars
#Proc Getdata:
// T (1-5)
// cat n obs mean stdev
//===== ===== ===== =====
Data:
coh 39 52.38 15.44
exp 39 48.90 8.94
conf 39 50.90 10.99
ind 39 44.62 9.59
achiev 39 54.92 8.59
// render points
#Proc Scatterplot:
Xfield: 1
Yfield: 3
symbol: shape=pixsquare style=solid color=black radius=0.05
// render bars
#proc bars
locfield: 1
lenfield: 3
errbarfields: 4
tails: @TAILLEN
/////////////////////////////////////////////////
//// Bottom row
/////////////////////////////////////////////////
// specify data for left bars
#proc getdata:
// E (6-10)
//cat n obs mean stdev pval
//===== ===== ===== ===== =====
data:
int 23 45.78 14.30 p<.05
act 23 45.00 12.30 p>.05
mor 23 53.26 9.11 p<.05
org 23 49.17 10.43 p<.05
con 23 59.48 10.40 p>.05
#proc areadef:
areaname: 2lo
xscaletype: categories
xcategories: datafield 1
yrange: 0 100
#proc xaxis
stubs text
Intellectual-\nCultural\nOrientation
Active-\nRecreational\nOrientation
Moral-\nReligious\nEmphasis
Organization
Control
#proc yaxis
stubs: inc 10
label: Standard Score
// do horiz. lines
#Proc line
notation: locvalues
linedetails: width=0.4
Points: min 50(s) max 50(s)
#proc line:
notation: locvalues
linedetails: style=1 dashscale=5 width=0.4
Points: min 40(s) max 40(s)
min 60(s) max 60(s)
// do centered P values
#Proc Scatterplot
Xfield: 1
ylocation: 5(s)
labelfield: 5
#proc catslide
axis: X
amount: -0.1
// render points
#Proc Scatterplot:
Xfield: 1
Yfield: 3
symbol: shape=pixdiamond style=solid color=black radius=0.05
// render bars
#Proc bars
locfield: 1
lenfield: 3
errbarfields: 4
tails: @TAILLEN
//////////// Right bars
#proc catslide
axis: x
amount: 0.1
#proc getdata:
// T (6-10)
// cat n obs mean stdev
//==== ===== ===== =====
Data:
int 39 53.51 9.78
act 39 47.54 11.50
mor 39 58.67 6.34
org 39 55.33 9.86
con 39 56.31 8.97
// render points
#Proc Scatterplot:
Xfield: 1
Yfield: 3
symbol: shape=pixsquare style=solid color=black radius=0.05
// render bars
#proc bars
locfield: 1
lenfield: 3
errbarfields: 4
tails: @TAILLEN
// Define legend entries explicitly
#proc legendentry
sampletype: symbol
label: Emergency Shelters (N = 23)
details: shape=diamond style=filled fillcolor=black radius=0.05
#proc legendentry
sampletype: symbol
label: Transitional Housing (N = 39)
details: shape=square style=filled fillcolor=black radius=0.05
#proc legend:
location: 5.5 9
|
