Difference between revisions of "Input files"
Line 4: | Line 4: | ||
*".blkcsv" for comma separated value ascii text files in block format | *".blkcsv" for comma separated value ascii text files in block format | ||
*".bin" for binary files in block format | *".bin" for binary files in block format | ||
Note that this mechanism does not apply to files containing genotypes. | |||
== ".csv" files == | == ".csv" files == | ||
Line 62: | Line 64: | ||
T,T | T,T | ||
T,F | T,F | ||
== genotype files == | |||
The genotype files must be in ascii text format where each line contains a single genotype coded 0,1 or 2 for homozygous "aa", heterozygous "Aa" and homozygous "AA", respectively, and 3 for missing. A single genotype must not contain any space. The file must have a many lines as there are genotypes. A valid genotype file maybe | |||
0102210 | |||
1021021 |
Revision as of 09:59, 28 December 2020
File content determination via file name extensions
lmt automatically detects the format of input files by the filename extension. Supported extensions are
- ".csv" for ordinary comma separated values ascii text files
- ".blkcsv" for comma separated value ascii text files in block format
- ".bin" for binary files in block format
Note that this mechanism does not apply to files containing genotypes.
".csv" files
".csv" files may contain several commented lines at the top only where the comment character is "#".
The type of the file content is determined by its prospective use, that is
- the data file is supposed to contain only real/float numbers which are transferred to integer if required,
- a file containing an ordinary pedigree is supposed to contain only integer numbers,
- a file containing a missing value indicator matrix is supposed to contain only character strings.
data file
lmt accepts only a single file containing the actual data. A data file in ".csv" format must follow the following formatting rules:
- file must have at least one commented line where the last commented line must containing the column names separated by comma,
- the column names
- must be alpha-numeric only
- must not be quoted
- must be unique
- below the header the data file must contain only numeric values where the decimal separator is a dot(".").
An example for a data file with three columns is shown below.
#y,mu,id 25.0,1,5 33.1,1,6 36.0,1,7 28.3,1,8
co-variance matrix file
A co-variance matrix files must contain a single full squared symmetric matrix, for instance
1.5,0.8,0.1 0.8,2.1,1.1 0.1,1.1,1.9
co-variance mask files
A co-variance mask file communicates to lmt which co-variances should remain constant when lmt is used to estimate variances. The mask file must contain characters which can be interpreted as a boolean data type, and where "T" codes for co-variances which must remain constant and "F" codes for variances which are allowed to float. The mask file must have the same dimensions as the associated co-variance matrix. For instance the above co-variance matrix file maybe accompanied by a mask file containing
T,F,F F,T,F F,F,F
which communicates to lmt that the first and second diagonal element should remain at their original values.
missing observations indicator file
The pattern of missing observations maybe communicated via an indicator file, where the file must contain characters which can be interpreted as a boolean data type, and where "T" codes for an available observation and "F" codes for a missing observation. Further, similar to the data file, the missing value indicator file must contain a header with the same column names as the observation columns in the data file.
For example a data file
#y1,y2,mu,id 25.0,0.0,1,5 0.0,0.8,1,6 36.0,-1.5,1,7 28.3,0.0,1,8
maybe accompanied by a missing value indicator file
#y1,y2 T,F F,T T,T T,F
genotype files
The genotype files must be in ascii text format where each line contains a single genotype coded 0,1 or 2 for homozygous "aa", heterozygous "Aa" and homozygous "AA", respectively, and 3 for missing. A single genotype must not contain any space. The file must have a many lines as there are genotypes. A valid genotype file maybe
0102210 1021021