title 'Profile Analysis Example: PROBE2 (Timm, Ex. 3.16)'; data probe2; *-------------------------------------------------------------* TIMM - Example 3.16 (p.244) Response speed in probed recall of sentences. The task is identical to that in dataset PROBE1. Group 1: Subjects with low STM capacity. Group 2: Subjects with high STM capacity *-------------------------------------------------------------*; input subjno 2-3 p1-p5; label P1='Position 1' P2='Position 2' P3='Position 3' P4='Position 4' P5='Position 5' SUBJNO='Subject number'; if subjno<11 then group='1Lo_STM'; else group='2Hi_STM'; datalines; S01 20 21 42 32 32 S02 67 29 56 39 41 S03 37 25 28 31 34 S04 42 38 36 19 35 S05 57 32 21 30 29 S06 39 38 54 31 28 S07 43 20 46 42 31 S08 35 34 43 35 42 S09 41 23 51 27 30 S10 39 24 35 26 32 S11 47 25 36 21 27 S12 53 32 48 46 54 S13 38 33 42 48 49 S14 60 41 67 53 50 S15 37 35 45 34 46 S16 59 37 52 36 52 S17 67 33 61 31 50 S18 43 27 36 33 32 S19 64 53 62 40 43 S20 41 34 47 37 46 ; proc print; id subjno group; proc glm; class group; model p1-p5 = group / nouni; repeated position 5 profile /short printh printe; manova h = group / printe printh canonical short; title2 'Two Sample Profile Analysis'; * A similar analysis is produced by CANDISC, but with more detail; proc candisc data=probe2 short /* suppress most printed output */ out=disc /* output disc. function scores */ outstat=canstat; /* output sscp matrices */ class group; var p1-p5; title2 'Output from CANDISC Procedure'; /* CANDISC produces two output data sets: DISC contains the original scores + discriminant function scores CANSTAT contains the various SSCP and other matrices (identified by a _TYPE_ variable */ proc print data=disc; var subjno p1-p5 group can1; title2 'Scores on discriminant function'; proc ttest data=disc; class group; var can1; title3 'Univariate t-test on discriminant scores '; *-- Plotting means; proc transpose data=probe2 out=long; var p1-p5; by group subjno; run; data long; set long; rename col1=Speed _name_=Pos; drop _label_; label _name_='Position'; run; title; %meanplot(data=long, var=Speed, class=Pos Group);