Errors in PREPLINK on SUNs

[Daniel E. Weeks]

Heikki Lehvaslaiho  <Heikki.Lehvaslaiho at Helsinki.FI> writes:
***********************************************************************
Running the Sun version of preplink (part of the LINKAGE package v. 5.1)
I constantly get error messages when trying to read from or write to a
file. The error message is: 

Label of 3 not found in case
*** Program execution terminated by trace trap 
???() 
In readdata.readloci.getlocus (system = 2), line 175 in "./preplink.p"
called by readdata.readloci (), line 219 in "./preplink.p"

Examining the code I found out that procedure getlocus first determines
value "numbers" (3) to the variable "which" and then a few lines later
doesn't check that possibility in a case statement. The same logical flaw
is the corresponding procedure writing to a file, writedata.

Could someone please tell me what is the right way of correcting the
error?
***********************************************************************

I too ran into problems with preplink on the SUN. My 
solution to the "numbers" problem described by Heikki
Lehvaslaiho is to add two lines (one around 176, and one
around 369):
            CASE which OF
              binary : getbin(thislocus[system]);
              affection : getpen(thislocus[system]);
              quantitative : getquan(thislocus[system]);
              numbers:; {Added for SUN Pascal}
            END;


                FOR i:=1 TO ntrait DO
                  FOR j:=i TO ntrait DO
                    write(outdata,vmat[i,j]:7:3);
                writeln(outdata,'  << VARIANCE - COVARIANCE MATRIX');
                write(outdata,contrait:7:3);
                writeln(outdata,'  << MULTIPLIER FOR VARIANCE IN HETEROZYGOTES');
              END;
            numbers:; {Added for SUN Pascal}
          END;
          IF risk AND (risklocus=system)
          THEN
            writeln(outdata,riskall:3,' << ALLELE AT RISK');
        END;


Heikki Lehvaslaiho  <Heikki.Lehvaslaiho at Helsinki.FI> writes:
***********************************************************************
There seems to be other bugs in preplink, too. In the affection status
phenotype, some of the several penetrances that I type in are not
recorded. In the next screen there is a zero in their place. I haven't
look into that part of the code yet.
***********************************************************************

I don't seem to have that problem on my SUN.  Could you let me
know exactly when it occurs, so I can try to see if I can make
it happen here?

  -- Dan Weeks --
********************************
| Daniel E. Weeks              | 
| Department of Human Genetics |  Internet: 
| University of Pittsburgh     |    dweeks at watson.hgen.pitt.edu
| Crabtree Hall, Room A310     |
| 130 DeSoto Street            |  Bitnet:
| Pittsburgh, PA 15261         |    weeks at pittvms.bitnet
|                              |
| (412) 624-3066               |
| FAX: (412) 624-3020          |
********************************