;procedure evalcoran does a standart evaluation of the results of
;correspondence analysis for 8 factors:
;
;a) create differential reconstituted images for both ends of all 8
;   factors. these are arranged in a montage with 4 images/row. 
;   order of images: image 1 = factor 1, value -1
;                    image 2 = factor 1, value +1
;                    image 3 = factor 2, value -1
;                    image 4 = factor 2, value +1 ..etc.
;
;b) create 7 visual maps for all 8 factors. maps are:
;                             image file name:   postscript file names:
;                    1 versus 2    vismap#12     vismappost#12 ..etc.
;                    2 versus 3    vismap#23
;                    3 versus 4    vismap#34   ..etc.
; (# stands for any number between 1 and 9 and depends on the value
; entered to the query: ?vismap number offset (e.g. 100,200 etc)?
; 
; author:
; m.radermacher, 2/95
rr x20
?coord. file code?
;[enter file code like used in ca s]
;------------------------------------------------------------------1
rr x30
?vismap number offset (e.g. 100,200 etc)?
;[see above, this number is used to construct file numbers
; for the visual maps. the numbers will be 12+offset for map 1v2,
; 23+offset for map 2v3 etc.] 
;------------------------------------------------------------------2
fr
?filename for montage of reconstituted images?
;[enter name of file that will contain the 4x4 montage of the differ-
; ential reconstituted images.
;------------------------------------------------------------------3
fr
?filename of postscript file of reconst. images?
;[enter name of file that will receive the postscript version of the
; montage of differential reconstituted images]
;------------------------------------------------------------------4
x21=8
sd c
x20,x21
?imc document file name?
;[enter name of document file that is created as a document file 
; format copy of the file imc***, the image coordinate file]
;------------------------------------------------------------------5
fr
?one of the image files (must exist)?
;[enter the name of one file in the image series that will be used
; to create the visual maps]
;------------------------------------------------------------------6
rr x22
?upper threshold in vismap as multiple of sigma?
;[enter threshold used for each individual average in the visual 
; map. a good starting value is 2]
;------------------------------------------------------------------7
rr x23
?lower threshold in vismap?
;[see 7. again a good start value is 2]
;------------------------------------------------------------------8
rr x24
?first image number?
;[enter first image number in image file series]
;------------------------------------------------------------------9
rr x25
?last image number?
;[enter last image number of file series (gaps in series are 
; acceptable)]
;-----------------------------------------------------------------10
ca sr
screconst0001
5011,5012
x20
1
-1
1
ca sr
screconst0001
5021,5022
x20
2
-1
1
ca sr
screconst0001
5031,5032
x20
3
-1
1
ca sr
screconst0001
5041,5042
x20
4
-1
1
ca sr
screconst0001
5051,5052
x20
5
-1
1
ca sr
screconst0001
5061,5062
x20
6
-1
1
ca sr
screconst0001
5071,5072
x20
7
-1
1
ca sr
screconst0001
5081,5082
x20
8
-1
1
mn s
screconst0001
5011,5012,5021,5022,5031,5032,5041,5042,5051,5052,
5061,5062,5071,5072,5081,5082
4,1
2
<3>
cp po
<3>
<4>
0,2
18
-1
n
x40=1+x30
do lb1 i=1,7
x30=x0
x31=x30+1
x40=x40+11
ca vis
(700,700)
<5>
visdocx40
<6>
(10,10)
(2.3)
x22,x23
vismapx40
x24,x25
x30,x31
fs x50,x51
vismapx40
cp po
vismapx40
vismappostx40
x51,x50
(18)
(-1)
n
lb1
en