# COMPUTATION OF JACKSON INDEX OF INDIVIDUAL RELIABILITY
# CREATE FACET ITEM VECTORS;
$K=0;
for [$J = 1; $J < 11; $J++] {
$ANX[$J]=ITEMS[1+$K];
$FRI[$J]=ITEMS[2+$K];
$IMA[$J]=ITEMS[3+$K];
$TRU[$J]=ITEMS[4+$K];
$SFEFF[$J]=ITEMS[5+$K];
$ANG[$J]=ITEMS[6+$K];
$GRE[$J]=ITEMS[7+$K];
$ART[$J]=ITEMS[8+$K];
$MOR[$J]=ITEMS[9+$K];
$ORD[$J]=ITEMS[10+$K];
$DEP[$J]=ITEMS[11+$K];
$ASS[$J]=ITEMS[12+$K];
$EMO[$J]=ITEMS[13+$K];
$ALT[$J]=ITEMS[14+$K];
$DUT[$J]=ITEMS[15+$K];
$SFCON[$J]=ITEMS[16+$K];
$ACT[$J]=ITEMS[17+$K];
$ADV[$J]=ITEMS[18+$K];
$COO[$J]=ITEMS[19+$K];
$ACH[$J]=ITEMS[20+$K];
$IMM[$J]=ITEMS[21+$K];
$EXC[$J]=ITEMS[22+$K];
$INT[$J]=ITEMS[23+$K];
$MOD[$J]=ITEMS[24+$K];
$SFDIS[$J]=ITEMS[25+$K];
$VUL[$J]=ITEMS[26+$K];
$CHE[$J]=ITEMS[27+$K];
$LIB[$J]=ITEMS[28+$K];
$SYM[$J]=ITEMS[29+$K];
$CAU[$J]=ITEMS[30+$K];
$K=$K+30;
}
# CREATING TWO HALF-SCALE SCORES FOR JACKSON RELIABILTY;
for [$I = 1; $I < 61; $I++] {
$SH[$I]=0;
} 
for [$I = 2; $I < 11; $I++] {
$SH[1]=$SH[1]+$ANX[$I-1];
$SH[2]=$SH[2]+$FRI[$I-1];
$SH[3]=$SH[3]+$IMA[$I-1];
$SH[4]=$SH[4]+$TRU[$I-1];
$SH[5]=$SH[5]+$SFEFF[$I-1];
$SH[6]=$SH[6]+$ANG[$I-1];
$SH[7]=$SH[7]+$GRE[$I-1];
$SH[8]=$SH[8]+$ART[$I-1];
$SH[9]=$SH[9]+$MOR[$I-1];
$SH[10]=$SH[10]+$ORD[$I-1];
$SH[11]=$SH[11]+$DEP[$I-1];
$SH[12]=$SH[12]+$ASS[$I-1];
$SH[13]=$SH[13]+$EMO[$I-1];
$SH[14]=$SH[14]+$ALT[$I-1];
$SH[15]=$SH[15]+$DUT[$I-1];
$SH[16]=$SH[16]+$SFCON[$I-1];
$SH[17]=$SH[17]+$ACT[$I-1];
$SH[18]=$SH[18]+$ADV[$I-1];
$SH[19]=$SH[19]+$COO[$I-1];
$SH[20]=$SH[20]+$ACH[$I-1];
$SH[21]=$SH[21]+$IMM[$I-1];
$SH[22]=$SH[22]+$EXC[$I-1];
$SH[23]=$SH[23]+$INT[$I-1];
$SH[24]=$SH[24]+$MOD[$I-1];
$SH[25]=$SH[25]+$SFDIS[$I-1];
$SH[26]=$SH[26]+$VUL[$I-1];
$SH[27]=$SH[27]+$CHE[$I-1];
$SH[28]=$SH[28]+$LIB[$I-1];
$SH[29]=$SH[29]+$SYM[$I-1];
$SH[30]=$SH[30]+$CAU[$I-1];
$SH[31]=$SH[31]+$ANX[$I];
$SH[32]=$SH[32]+$FRI[$I];
$SH[33]=$SH[33]+$IMA[$I];
$SH[34]=$SH[34]+$TRU[$I];
$SH[35]=$SH[35]+$SFEFF[$I];
$SH[36]=$SH[36]+$ANG[$I];
$SH[37]=$SH[37]+$GRE[$I];
$SH[38]=$SH[38]+$ART[$I];
$SH[39]=$SH[39]+$MOR[$I];
$SH[40]=$SH[40]+$ORD[$I];
$SH[41]=$SH[41]+$DEP[$I];
$SH[42]=$SH[42]+$ASS[$I];
$SH[43]=$SH[43]+$EMO[$I];
$SH[44]=$SH[44]+$ALT[$I];
$SH[45]=$SH[45]+$DUT[$I];
$SH[46]=$SH[46]+$SFCON[$I];
$SH[47]=$SH[47]+$ACT[$I];
$SH[48]=$SH[48]+$ADV[$I];
$SH[49]=$SH[49]+$COO[$I];
$SH[50]=$SH[50]+$ACH[$I];
$SH[51]=$SH[51]+$IMM[$I];
$SH[52]=$SH[52]+$EXC[$I];
$SH[53]=$SH[53]+$INT[$I];
$SH[54]=$SH[54]+$MOD[$I];
$SH[55]=$SH[55]+$SFDIS[$I];
$SH[56]=$SH[56]+$VUL[$I];
$SH[57]=$SH[57]+$CHE[$I];
$SH[58]=$SH[58]+$LIB[$I];
$SH[59]=$SH[59]+$SYM[$I];
$SH[60]=$SH[60]+$CAU[$I];
$I++;
}
# COMPUTE SPLIT-HALF RELIABILITY COEFFICIENT;
$SUMODD=0;
$SUMEVEN=0;
$SPRDOE=0;
$SUMOSQ=0;
$SUMESQ=0;
for [$I = 1; $I < 31; $I++] {
$SUMODD=$SUMODD+$SH[$I];
$SUMEVEN=$SUMEVEN+$SH[$I+30];
$SPRDOE=$SPRDOE+[[$SH[$I]]*[$SH[$I+30]]];
$SUMOSQ=$SUMOSQ+[$SH[$I]]**2;
$SUMESQ=$SUMESQ+[$SH[$I+30]]**2;
}
$SPLTHALF=[$SPRDOE-[$SUMODD*$SUMEVEN/30]]/[[$SUMOSQ-[$SUMODD**2]/30]*[$SUMESQ-[$SUMEVEN**2]/30]]**.5;
$SPLTHALF=[2*SPLTHALF]/[1+SPLTHALF];