data HRSCPV;
infile './HRSCPV.csv' dlm=',';
input subject $ cond $ period $ HR SC PV;

/* 課題、期間ごとに平均値を求める */
proc sort;
by cond period;
proc means data=HRSCPV;
by cond period;

/* 2要因とも対応がある場合の分散分析:HR */
proc glm data=HRSCPV;
class cond period subject;
model HR= subject cond cond*subject period period*subject
cond*period cond*period*subject /ss3;

test h= cond e=cond*subject;
test h= period e=period*subject;
test h= cond*period e=cond*period*subject;

/* 多重比較 */
means cond / tukey e = period*subject alpha=0.05;
means period / tukey e = period*subject alpha=0.05;

/* 2要因とも対応がある場合の分散分析:SC */
proc glm data=HRSCPV;
class cond period subject;
model SC= subject cond cond*subject period period*subject
cond*period cond*period*subject /ss3;

test h= cond e=cond*subject;
test h= period e=period*subject;
test h= cond*period e=cond*period*subject;

/* 多重比較 */
means cond / tukey e = period*subject alpha=0.05;
means period / tukey e = period*subject alpha=0.05;

/* 2要因とも対応がある場合の分散分析:PV */
proc glm data=HRSCPV;
class cond period subject;
model PV= subject cond cond*subject period period*subject
cond*period cond*period*subject /ss3;

test h= cond e=cond*subject;
test h= period e=period*subject;
test h= cond*period e=cond*period*subject;

/* 多重比較 */
means cond / tukey e = period*subject alpha=0.05;
means period / tukey e = period*subject alpha=0.05;