Hi all,
I'm looking to create some summary statistic as seen below
Is it possible to combine the results of 2 proc means? Using the data below, I want to calculate proc means with 2 different classes. Like so:
proc means data= have n min max median mean stddev maxdec=1 ; class con parameter_n parameter ; var score; output out=overall; run;
proc means data= have n min max median mean stddev maxdec=1 ; class con parameter_n parameter site ; var score; output out=bysite; run;
MEMBER
SITE
PARAMETER
PARAMETER_N
SCORE
CON
1
A
OSUT
1
12
ARI
1
A
OCDT
2
12
ARI
2
A
OSUT
1
32
ARI
2
A
OCDT
2
37
ARI
3
A
OSUT
1
31
ARI
3
A
OCDT
2
24
ARI
4
A
OSUT
1
34
ARI
4
A
OCDT
2
24
ARI
5
A
OSUT
1
37
ARI
5
A
OCDT
2
20
ARI
22
A
OSUT
1
21
ASI
22
A
OCDT
2
48
ASI
25
A
OSUT
1
21
ASI
25
A
OCDT
2
41
ASI
26
A
OSUT
1
23
ASI
26
A
OCDT
2
42
ASI
27
A
OSUT
1
32
ASI
27
A
OCDT
2
52
ASI
28
A
OSUT
1
30
ASI
28
A
OCDT
2
44
ASI
9
B
OSUT
1
39
ARI
9
B
OCDT
2
25
ARI
10
B
OSUT
1
29
ARI
10
B
OCDT
2
27
ARI
11
B
OSUT
1
29
ARI
11
B
OCDT
2
23
ARI
15
B
OSUT
1
33
ARI
15
B
OCDT
2
29
ARI
16
B
OSUT
1
27
ARI
16
B
OCDT
2
36
ARI
17
B
OSUT
1
29
ASI
17
B
OCDT
2
32
ASI
20
B
OSUT
1
13
ASI
20
B
OCDT
2
36
ASI
21
B
OSUT
1
39
ASI
21
B
OCDT
2
34
ASI
22
B
OSUT
1
6
ASI
22
B
OCDT
2
32
ASI
23
B
OSUT
1
18
ASI
23
B
OCDT
2
33
ASI
My results should like like this:
Parameter Site Statistics ARI ASI Diff ASI - ARI
OSUT Overall n xx xx -xx
mean (SD) xx.x(x.x) xx.x(x.x)
median xx.x xx.x
Min, Max xx, xx xx, xx
95% CI xx.x, xx.x xx.x, xx.x -xx.x, -x.x
P-value <0.0001
A n xx xx -xx
mean (SD) xx.x(x.x) xx.x(x.x)
median xx.x xx.x
Min, Max xx, xx xx, xx
95% CI xx.x, xx.x xx.x, xx.x -xx.x, -x.x
P-value <0.0001
B n xx xx -xx
mean (SD) xx.x(x.x) xx.x(x.x)
median xx.x xx.x
Min, Max xx, xx xx, xx
95% CI xx.x, xx.x xx.x, xx.x -xx.x, -x.x
P-value <0.0001
Repeat for parameter OCDT
Not sure if this is possible.
I get close to what I want using proc univariate but the output just looks like raw data. I prefer the cleaner look in proc means.
... View more