! From FQuehenberger@lumc.nl Fri Feb 14 11:04:08 2003
! Date: Thu, 13 Feb 2003 17:58:05 +0100
! From: "Quehenberger,  F. (MStat)" <FQuehenberger@lumc.nl>
! To: "'davidD@qimr.edu.au'" <davidD@qimr.edu.au>
! Subject: sib-pair write mendel bug
! 
! Dear David,
! 
! I have encounters problems with sib-pair writing mendel pedigree files
! (Windows version).
! 
! It seems that write mendel does not export the affection status. It is left
! out from the file if a variable is delclared as affection, but included if
! it is declared to be quantitative. However, affection variables are exported
! properly by write sage, for example. Furthermore space is reserved for the
! affection variable in the fortran format statement in the second line of the
! output file. So it seems that the affection variable was just omitted in the
! write statement.
! 
! I will give a small test file that could demonstrate the problem.
! 
! Thank you very much for providing this great program!
! 
! Best,
! Franz
! 
! 
! -- 
! mailto:fquehenberger@lumc.nl
! Department of Medical Statistics, University of Leiden, The Netherlands
! tel: ++31 71 527 6736 (office), ++31 71 5324019 (private)
! http://www.medstat.medfac.leidenuniv.nl/ms/index.html
!
! Demonstrate that the affection variable is not written to mendel output
!
set loc ca affection
set loc age quantitative
set loc m1 marker
read ped inline
A 0020    x    x m 2. 72 x x                    
A 0019    x    x f 2. 48 x x                    
A 0017    x    x m 2. 47 x x                    
A 0012 0020 0019 f 1. 78 2 1                    
A 0030    x    x m 2.  x x x                    
A 0015    x    x m 2.  x x x                    
A 0022 0020 0019 m 2. 41 x x                    
A 0051    x    x f 2. 44 x x                    
A 0023 0020 0019 f 2. 79 2 1                    
A 0024 0020 0019 m 2. 68 1 1                    
A 0028 0020 0019 f 2. 55 2 1                    
A 0029 0020 0019 f 2. 66 2 1                    
A 0018    x    x m 2. 58 x x                    
;;;;
run
write mendel write_mendel_out.ped
write locus mendel write_mendel_out.lom trait
# write sage write_sage_out.ped

