family ind father mother sex generation affection risk genetic_risk locus_risk_1 locus_risk_2 locus_risk_3 genotype_1 genotype_2 genotype_3
1 111 97 106 -1 9 1 1003 1003 1001 1 1 0/1 1/1 1/1
1 97 93 96 1 8 0 3 3 1 1 1 1/1 1/1 1/1
1 106 0 0 -1 8 1 1003 1003 1001 1 1 0/1 1/1 1/1
1 93 89 59 1 7 0 4 4 1 1 2 1/1 1/1 0/1
1 96 0 0 -1 7 0 3 3 1 1 1 1/1 1/1 1/1
1 107 97 106 -1 9 0 3 3 1 1 1 1/1 1/1 1/1
1 108 97 106 -1 9 0 3 3 1 1 1 1/1 1/1 1/1
1 109 97 106 1 9 1 1003 1003 1001 1 1 0/1 1/1 1/1
1 110 97 106 1 9 0 3 3 1 1 1 1/1 1/1 1/1
1 59 50 57 -1 6 0 3 3 1 1 1 1/1 1/1 1/1
1 89 0 0 1 6 0 4 4 1 1 2 1/1 1/1 0/1
1 98 93 96 1 8 0 3 3 1 1 1 1/1 1/1 1/1
1 99 93 96 1 8 0 4 4 1 1 2 1/1 1/1 0/1
1 100 93 96 -1 8 0 4 4 1 1 2 1/1 1/1 0/1
1 101 93 96 1 8 0 4 4 1 1 2 1/1 1/1 0/1
1 102 93 96 1 8 0 4 4 1 1 2 1/1 1/1 0/1
1 103 93 96 1 8 0 3 3 1 1 1 1/1 1/1 1/1
1 104 93 96 1 8 0 4 4 1 1 2 1/1 1/1 0/1
1 105 93 96 1 8 0 3 3 1 1 1 1/1 1/1 1/1
1 215 109 214 -1 10 0 3 3 1 1 1 1/1 1/1 1/1
1 216 109 214 -1 10 0 3 3 1 1 1 1/1 1/1 1/1
1 50 0 0 1 5 0 4 4 1 1 2 1/1 1/1 0/1
1 57 0 0 -1 5 0 3 3 1 1 1 1/1 1/1 1/1
1 90 89 59 -1 7 0 4 4 1 1 2 1/1 1/1 0/1
1 91 89 59 1 7 0 4 4 1 1 2 1/1 1/1 0/1
1 92 89 59 -1 7 0 4 4 1 1 2 1/1 1/1 0/1
1 113 112 100 -1 9 0 4 4 1 1 2 1/1 1/1 0/1
1 114 112 100 -1 9 1 1003 1003 1001 1 1 0/1 1/1 1/1
1 115 112 100 -1 9 1 1004 1004 1001 1 2 0/1 1/1 0/1
1 214 0 0 -1 9 0 3 3 1 1 1 1/1 1/1 1/1
1 112 0 0 1 8 1 1003 1003 1001 1 1 0/1 1/1 1/1
2 37 0 0 1 -1 1 1003 1003 1001 1 1 0/1 1/1 1/1
2 6 37 38 -1 0 1 1003 1003 1001 1 1 0/1 1/1 1/1
2 7 2 6 -1 1 0 3 3 1 1 1 1/1 1/1 1/1
2 8 2 6 1 1 0 3 3 1 1 1 1/1 1/1 1/1
2 9 2 6 1 1 0 3 3 1 1 1 1/1 1/1 1/1
2 10 2 6 -1 1 1 1003 1003 1001 1 1 0/1 1/1 1/1
2 11 2 6 1 1 0 3 3 1 1 1 1/1 1/1 1/1
2 38 75 76 -1 -1 0 3 3 1 1 1 1/1 1/1 1/1
2 2 67 68 1 0 0 3 3 1 1 1 1/1 1/1 1/1
2 20 19 7 -1 2 0 3 3 1 1 1 1/1 1/1 1/1
2 21 19 7 1 2 0 3 3 1 1 1 1/1 1/1 1/1
2 22 19 7 -1 2 0 3 3 1 1 1 1/1 1/1 1/1
2 23 19 7 1 2 0 3 3 1 1 1 1/1 1/1 1/1
2 24 19 7 -1 2 0 3 3 1 1 1 1/1 1/1 1/1
2 25 19 7 1 2 0 3 3 1 1 1 1/1 1/1 1/1
2 26 19 7 -1 2 0 3 3 1 1 1 1/1 1/1 1/1
2 40 9 39 1 2 0 3 3 1 1 1 1/1 1/1 1/1
2 41 9 39 -1 2 0 3 3 1 1 1 1/1 1/1 1/1
2 42 9 39 -1 2 0 3 3 1 1 1 1/1 1/1 1/1
2 43 9 39 1 2 0 3 3 1 1 1 1/1 1/1 1/1
2 56 55 10 1 2 1 1003 1003 1001 1 1 0/1 1/1 1/1
2 57 55 10 1 2 1 1003 1003 1001 1 1 0/1 1/1 1/1
2 58 55 10 1 2 1 1003 1003 1001 1 1 0/1 1/1 1/1
2 59 55 10 -1 2 1 1003 1003 1001 1 1 0/1 1/1 1/1
2 75 0 0 1 -2 1 1003 1003 1001 1 1 0/1 1/1 1/1
2 76 0 0 -1 -2 0 3 3 1 1 1 1/1 1/1 1/1
2 67 0 0 1 -1 0 3 3 1 1 1 1/1 1/1 1/1
2 68 0 0 -1 -1 0 3 3 1 1 1 1/1 1/1 1/1
2 19 0 0 1 1 0 3 3 1 1 1 1/1 1/1 1/1
2 108 23 107 -1 3 0 3 3 1 1 1 1/1 1/1 1/1
2 109 23 107 -1 3 0 3 3 1 1 1 1/1 1/1 1/1
2 110 23 107 -1 3 0 3 3 1 1 1 1/1 1/1 1/1
2 70 69 24 -1 3 0 3 3 1 1 1 1/1 1/1 1/1
2 71 69 24 1 3 0 3 3 1 1 1 1/1 1/1 1/1
2 72 69 24 -1 3 0 3 3 1 1 1 1/1 1/1 1/1
2 73 69 24 -1 3 0 3 3 1 1 1 1/1 1/1 1/1
2 74 69 24 1 3 0 3 3 1 1 1 1/1 1/1 1/1
2 141 25 140 1 3 0 3 3 1 1 1 1/1 1/1 1/1
2 142 25 140 1 3 0 3 3 1 1 1 1/1 1/1 1/1
2 143 25 140 1 3 0 3 3 1 1 1 1/1 1/1 1/1
2 130 129 26 -1 3 0 3 3 1 1 1 1/1 1/1 1/1
2 131 129 26 -1 3 0 3 3 1 1 1 1/1 1/1 1/1
2 132 129 26 -1 3 0 3 3 1 1 1 1/1 1/1 1/1
2 39 0 0 -1 1 0 3 3 1 1 1 1/1 1/1 1/1
2 47 46 41 -1 3 0 3 3 1 1 1 1/1 1/1 1/1
2 48 46 41 -1 3 0 4 4 1 1 2 1/1 1/1 0/1
2 49 46 41 -1 3 0 4 4 1 1 2 1/1 1/1 0/1
2 50 46 41 1 3 0 3 3 1 1 1 1/1 1/1 1/1
2 51 46 41 1 3 0 3 3 1 1 1 1/1 1/1 1/1
2 52 46 41 1 3 0 4 4 1 1 2 1/1 1/1 0/1
2 53 46 41 1 3 0 3 3 1 1 1 1/1 1/1 1/1
2 54 46 41 -1 3 0 4 4 1 1 2 1/1 1/1 0/1
2 55 0 0 1 1 0 3 3 1 1 1 1/1 1/1 1/1
2 69 0 0 1 2 0 3 3 1 1 1 1/1 1/1 1/1
2 129 0 0 1 2 0 3 3 1 1 1 1/1 1/1 1/1
2 46 0 0 1 2 0 4 4 1 1 2 1/1 1/1 0/1
2 107 0 0 -1 2 0 3 3 1 1 1 1/1 1/1 1/1
2 140 0 0 -1 2 0 3 3 1 1 1 1/1 1/1 1/1
3 3 0 0 -1 0 1 1003 1003 1001 1 1 0/1 1/1 1/1
3 46 45 3 1 1 0 3 3 1 1 1 1/1 1/1 1/1
3 47 45 3 -1 1 1 1003 1003 1001 1 1 0/1 1/1 1/1
3 45 125 126 1 0 0 3 3 1 1 1 1/1 1/1 1/1
3 267 46 266 -1 2 0 3 3 1 1 1 1/1 1/1 1/1
3 125 0 0 1 -1 0 3 3 1 1 1 1/1 1/1 1/1
3 126 0 0 -1 -1 0 3 3 1 1 1 1/1 1/1 1/1
3 266 0 0 -1 1 0 4 4 1 1 2 1/1 1/1 0/1
4 190 0 0 1 4 1 1003 1003 1001 1 1 0/1 1/1 1/1
4 191 190 137 1 5 0 3 3 1 1 1 1/1 1/1 1/1
4 192 190 137 1 5 0 3 3 1 1 1 1/1 1/1 1/1
4 193 190 137 -1 5 0 3 3 1 1 1 1/1 1/1 1/1
4 194 190 137 1 5 0 3 3 1 1 1 1/1 1/1 1/1
4 137 133 107 -1 4 0 3 3 1 1 1 1/1 1/1 1/1
4 107 102 36 -1 3 0 3 3 1 1 1 1/1 1/1 1/1
4 133 0 0 1 3 0 3 3 1 1 1 1/1 1/1 1/1
4 36 0 0 -1 2 0 3 3 1 1 1 1/1 1/1 1/1
4 102 0 0 1 2 0 3 3 1 1 1 1/1 1/1 1/1
4 134 133 107 -1 4 0 3 3 1 1 1 1/1 1/1 1/1
4 135 133 107 1 4 0 3 3 1 1 1 1/1 1/1 1/1
4 136 133 107 1 4 0 3 3 1 1 1 1/1 1/1 1/1
4 138 133 107 1 4 0 3 3 1 1 1 1/1 1/1 1/1
5 42 38 34 1 4 1 1003 1003 1001 1 1 0/1 1/1 1/1
5 34 33 20 -1 3 0 3 3 1 1 1 1/1 1/1 1/1
5 38 137 138 1 3 1 1003 1003 1001 1 1 0/1 1/1 1/1
5 20 10 19 -1 2 0 3 3 1 1 1 1/1 1/1 1/1
5 33 0 0 1 2 0 3 3 1 1 1 1/1 1/1 1/1
5 39 38 34 1 4 1 1003 1003 1001 1 1 0/1 1/1 1/1
5 40 38 34 1 4 1 1003 1003 1001 1 1 0/1 1/1 1/1
5 41 38 34 1 4 0 3 3 1 1 1 1/1 1/1 1/1
5 137 0 0 1 2 0 3 3 1 1 1 1/1 1/1 1/1
5 138 0 0 -1 2 1 1003 1003 1001 1 1 0/1 1/1 1/1
5 10 6 4 1 1 0 3 3 1 1 1 1/1 1/1 1/1
5 19 0 0 -1 1 0 3 3 1 1 1 1/1 1/1 1/1
5 35 33 20 1 3 0 3 3 1 1 1 1/1 1/1 1/1
5 36 33 20 -1 3 0 3 3 1 1 1 1/1 1/1 1/1
5 37 33 20 1 3 0 3 3 1 1 1 1/1 1/1 1/1
5 4 0 0 -1 0 0 3 3 1 1 1 1/1 1/1 1/1
5 6 0 0 1 0 0 3 3 1 1 1 1/1 1/1 1/1
5 142 37 141 -1 4 0 3 3 1 1 1 1/1 1/1 1/1
5 143 37 141 -1 4 0 3 3 1 1 1 1/1 1/1 1/1
5 144 37 141 1 4 0 3 3 1 1 1 1/1 1/1 1/1
5 145 37 141 1 4 0 3 3 1 1 1 1/1 1/1 1/1
5 146 37 141 -1 4 0 3 3 1 1 1 1/1 1/1 1/1
5 147 37 141 1 4 0 3 3 1 1 1 1/1 1/1 1/1
5 141 0 0 -1 3 0 3 3 1 1 1 1/1 1/1 1/1
6 225 219 204 -1 7 1 1003 1003 1001 1 1 0/1 1/1 1/1
6 204 203 152 -1 6 0 3 3 1 1 1 1/1 1/1 1/1
6 219 0 0 1 6 1 1003 1003 1001 1 1 0/1 1/1 1/1
6 152 88 147 -1 5 0 3 3 1 1 1 1/1 1/1 1/1
6 203 0 0 1 5 0 3 3 1 1 1 1/1 1/1 1/1
6 220 219 204 1 7 0 3 3 1 1 1 1/1 1/1 1/1
6 221 219 204 1 7 1 1003 1003 1001 1 1 0/1 1/1 1/1
6 222 219 204 -1 7 0 3 3 1 1 1 1/1 1/1 1/1
6 223 219 204 1 7 0 3 3 1 1 1 1/1 1/1 1/1
6 224 219 204 -1 7 0 3 3 1 1 1 1/1 1/1 1/1
6 226 219 204 1 7 0 3 3 1 1 1 1/1 1/1 1/1
6 88 86 60 1 4 0 3 3 1 1 1 1/1 1/1 1/1
6 147 0 0 -1 4 0 4 4 1 1 2 1/1 1/1 0/1
6 205 203 152 -1 6 0 3 3 1 1 1 1/1 1/1 1/1
6 206 203 152 1 6 0 3 3 1 1 1 1/1 1/1 1/1
6 207 203 152 1 6 0 3 3 1 1 1 1/1 1/1 1/1
6 208 203 152 1 6 0 3 3 1 1 1 1/1 1/1 1/1
6 236 235 224 -1 8 0 3 3 1 1 1 1/1 1/1 1/1
6 237 235 224 1 8 0 3 3 1 1 1 1/1 1/1 1/1
6 238 235 224 -1 8 0 3 3 1 1 1 1/1 1/1 1/1
6 239 235 224 -1 8 0 3 3 1 1 1 1/1 1/1 1/1
6 240 235 224 1 8 0 3 3 1 1 1 1/1 1/1 1/1
6 241 235 224 1 8 0 3 3 1 1 1 1/1 1/1 1/1
6 60 0 0 -1 3 0 3 3 1 1 1 1/1 1/1 1/1
6 86 0 0 1 3 0 3 3 1 1 1 1/1 1/1 1/1
6 148 88 147 1 5 0 3 3 1 1 1 1/1 1/1 1/1
6 149 88 147 -1 5 0 4 4 1 1 2 1/1 1/1 0/1
6 150 88 147 -1 5 0 4 4 1 1 2 1/1 1/1 0/1
6 151 88 147 1 5 0 3 3 1 1 1 1/1 1/1 1/1
6 235 0 0 1 7 0 3 3 1 1 1 1/1 1/1 1/1
6 256 241 255 1 9 0 3 3 1 1 1 1/1 1/1 1/1
6 257 241 255 1 9 0 3 3 1 1 1 1/1 1/1 1/1
6 255 0 0 -1 8 0 3 3 1 1 1 1/1 1/1 1/1
7 162 160 127 -1 4 1 1003 1003 1001 1 1 0/1 1/1 1/1
7 127 124 28 -1 3 1 2003 2003 2001 1 1 0/0 1/1 1/1
7 160 0 0 1 3 0 3 3 1 1 1 1/1 1/1 1/1
7 286 285 162 1 5 1 1004 1004 1001 1 2 0/1 1/1 0/1
7 287 285 162 -1 5 1 1003 1003 1001 1 1 0/1 1/1 1/1
7 288 285 162 1 5 1 1004 1004 1001 1 2 0/1 1/1 0/1
7 289 285 162 -1 5 1 1004 1004 1001 1 2 0/1 1/1 0/1
7 290 285 162 1 5 0 4 4 1 1 2 1/1 1/1 0/1
7 28 15 26 -1 2 1 1003 1003 1001 1 1 0/1 1/1 1/1
7 124 197 198 1 2 1 1003 1003 1001 1 1 0/1 1/1 1/1
7 161 160 127 -1 4 1 1003 1003 1001 1 1 0/1 1/1 1/1
7 163 160 127 -1 4 1 1003 1003 1001 1 1 0/1 1/1 1/1
7 164 160 127 -1 4 1 1003 1003 1001 1 1 0/1 1/1 1/1
7 285 0 0 1 4 1 1004 1004 1001 1 2 0/1 1/1 0/1
7 15 3 14 1 1 0 3 3 1 1 1 1/1 1/1 1/1
7 26 167 168 -1 1 1 1003 1003 1001 1 1 0/1 1/1 1/1
7 125 124 28 1 3 1 1003 1003 1001 1 1 0/1 1/1 1/1
7 126 124 28 -1 3 1 1003 1003 1001 1 1 0/1 1/1 1/1
7 128 124 28 1 3 1 1003 1003 1001 1 1 0/1 1/1 1/1
7 129 124 28 -1 3 1 1003 1003 1001 1 1 0/1 1/1 1/1
7 130 124 28 1 3 1 2003 2003 2001 1 1 0/0 1/1 1/1
7 131 124 28 -1 3 1 1003 1003 1001 1 1 0/1 1/1 1/1
7 197 219 220 1 1 0 3 3 1 1 1 1/1 1/1 1/1
7 198 0 0 -1 1 1 1003 1003 1001 1 1 0/1 1/1 1/1
7 216 215 161 1 5 1 1003 1003 1001 1 1 0/1 1/1 1/1
7 217 215 161 1 5 0 3 3 1 1 1 1/1 1/1 1/1
7 218 215 161 1 5 1 1003 1003 1001 1 1 0/1 1/1 1/1
7 195 194 163 -1 5 1 1003 1003 1001 1 1 0/1 1/1 1/1
7 196 194 163 -1 5 0 3 3 1 1 1 1/1 1/1 1/1
7 3 0 0 1 0 0 3 3 1 1 1 1/1 1/1 1/1
7 14 0 0 -1 0 1 1003 1003 1001 1 1 0/1 1/1 1/1
7 27 15 26 1 2 0 3 3 1 1 1 1/1 1/1 1/1
7 29 15 26 -1 2 0 3 3 1 1 1 1/1 1/1 1/1
7 30 15 26 -1 2 1 1003 1003 1001 1 1 0/1 1/1 1/1
7 167 0 0 1 0 1 2003 2003 2001 1 1 0/0 1/1 1/1
7 168 0 0 -1 0 0 4 4 1 2 1 1/1 0/1 1/1
7 219 0 0 1 0 0 3 3 1 1 1 1/1 1/1 1/1
7 220 0 0 -1 0 0 3 3 1 1 1 1/1 1/1 1/1
7 215 0 0 1 4 0 3 3 1 1 1 1/1 1/1 1/1
7 236 216 235 -1 6 1 1003 1003 1001 1 1 0/1 1/1 1/1
7 237 216 235 -1 6 1 1003 1003 1001 1 1 0/1 1/1 1/1
7 238 216 235 1 6 1 1003 1003 1001 1 1 0/1 1/1 1/1
7 194 0 0 1 4 0 3 3 1 1 1 1/1 1/1 1/1
7 235 0 0 -1 5 0 3 3 1 1 1 1/1 1/1 1/1
8 79 0 0 1 0 1 1003 1003 1001 1 1 0/1 1/1 1/1
8 70 79 80 1 1 0 3 3 1 1 1 1/1 1/1 1/1
8 71 70 30 -1 2 0 3 3 1 1 1 1/1 1/1 1/1
8 72 70 30 -1 2 0 3 3 1 1 1 1/1 1/1 1/1
8 80 0 0 -1 0 0 3 3 1 1 1 1/1 1/1 1/1
8 30 25 2 -1 1 0 3 3 1 1 1 1/1 1/1 1/1
8 2 0 0 -1 0 0 3 3 1 1 1 1/1 1/1 1/1
8 25 0 0 1 0 0 3 3 1 1 1 1/1 1/1 1/1
9 226 5 221 -1 1 1 1003 1003 1001 1 1 0/1 1/1 1/1
9 5 0 0 1 0 1 1003 1003 1001 1 1 0/1 1/1 1/1
9 221 277 278 -1 0 0 3 3 1 1 1 1/1 1/1 1/1
9 222 5 221 1 1 1 1003 1003 1001 1 1 0/1 1/1 1/1
9 223 5 221 1 1 1 1003 1003 1001 1 1 0/1 1/1 1/1
9 224 5 221 -1 1 1 1003 1003 1001 1 1 0/1 1/1 1/1
9 225 5 221 -1 1 0 3 3 1 1 1 1/1 1/1 1/1
9 277 0 0 1 -1 0 3 3 1 1 1 1/1 1/1 1/1
9 278 0 0 -1 -1 0 3 3 1 1 1 1/1 1/1 1/1
10 69 66 48 1 4 1 1003 1003 1001 1 1 0/1 1/1 1/1
10 48 26 47 -1 3 1 1003 1003 1001 1 1 0/1 1/1 1/1
10 66 208 209 1 3 0 3 3 1 1 1 1/1 1/1 1/1
10 26 24 12 1 2 0 3 3 1 1 1 1/1 1/1 1/1
10 47 104 105 -1 2 1 1003 1003 1001 1 1 0/1 1/1 1/1
10 67 66 48 1 4 0 3 3 1 1 1 1/1 1/1 1/1
10 68 66 48 1 4 0 3 3 1 1 1 1/1 1/1 1/1
10 208 0 0 1 2 0 3 3 1 1 1 1/1 1/1 1/1
10 209 0 0 -1 2 0 4 4 1 1 2 1/1 1/1 0/1
10 12 10 3 -1 1 0 3 3 1 1 1 1/1 1/1 1/1
10 24 0 0 1 1 0 4 4 1 1 2 1/1 1/1 0/1
10 49 26 47 -1 3 0 3 3 1 1 1 1/1 1/1 1/1
10 104 121 122 1 1 1 1003 1003 1001 1 1 0/1 1/1 1/1
10 105 0 0 -1 1 0 3 3 1 1 1 1/1 1/1 1/1
10 3 0 0 -1 0 0 4 4 1 2 1 1/1 0/1 1/1
10 10 0 0 1 0 0 3 3 1 1 1 1/1 1/1 1/1
10 25 24 12 -1 2 0 4 4 1 1 2 1/1 1/1 0/1
10 27 24 12 1 2 0 3 3 1 1 1 1/1 1/1 1/1
10 28 24 12 -1 2 0 4 4 1 1 2 1/1 1/1 0/1
10 29 24 12 -1 2 0 3 3 1 1 1 1/1 1/1 1/1
10 30 24 12 -1 2 0 4 4 1 1 2 1/1 1/1 0/1
10 121 0 0 1 0 1 2003 2003 2001 1 1 0/0 1/1 1/1
10 122 0 0 -1 0 0 3 3 1 1 1 1/1 1/1 1/1
11 136 103 135 -1 5 1 1003 1003 1001 1 1 0/1 1/1 1/1
11 103 99 77 1 4 0 3 3 1 1 1 1/1 1/1 1/1
11 135 0 0 -1 4 1 1003 1003 1001 1 1 0/1 1/1 1/1
11 249 248 136 1 6 0 3 3 1 1 1 1/1 1/1 1/1
11 250 248 136 -1 6 1 1003 1003 1001 1 1 0/1 1/1 1/1
11 251 248 136 -1 6 1 1003 1003 1001 1 1 0/1 1/1 1/1
11 252 248 136 1 6 1 1003 1003 1001 1 1 0/1 1/1 1/1
11 77 61 76 -1 3 0 3 3 1 1 1 1/1 1/1 1/1
11 99 159 160 1 3 0 3 3 1 1 1 1/1 1/1 1/1
11 248 0 0 1 5 0 3 3 1 1 1 1/1 1/1 1/1
11 61 58 7 1 2 0 3 3 1 1 1 1/1 1/1 1/1
11 76 79 80 -1 2 0 3 3 1 1 1 1/1 1/1 1/1
11 100 99 77 -1 4 0 3 3 1 1 1 1/1 1/1 1/1
11 101 99 77 -1 4 0 3 3 1 1 1 1/1 1/1 1/1
11 102 99 77 -1 4 0 3 3 1 1 1 1/1 1/1 1/1
11 104 99 77 -1 4 0 3 3 1 1 1 1/1 1/1 1/1
11 159 0 0 1 2 0 3 3 1 1 1 1/1 1/1 1/1
11 160 0 0 -1 2 0 3 3 1 1 1 1/1 1/1 1/1
11 7 0 0 -1 1 0 4 4 1 1 2 1/1 1/1 0/1
11 58 0 0 1 1 0 3 3 1 1 1 1/1 1/1 1/1
11 78 61 76 -1 3 0 3 3 1 1 1 1/1 1/1 1/1
11 79 0 0 1 1 0 3 3 1 1 1 1/1 1/1 1/1
11 80 0 0 -1 1 0 3 3 1 1 1 1/1 1/1 1/1
11 199 198 100 1 5 0 4 4 1 1 2 1/1 1/1 0/1
11 200 198 100 -1 5 0 3 3 1 1 1 1/1 1/1 1/1
11 201 198 100 1 5 0 3 3 1 1 1 1/1 1/1 1/1
11 202 198 100 1 5 0 4 4 1 1 2 1/1 1/1 0/1
11 203 198 100 1 5 0 4 4 1 1 2 1/1 1/1 0/1
11 204 198 100 -1 5 0 4 4 1 1 2 1/1 1/1 0/1
11 198 0 0 1 4 0 4 4 1 1 2 1/1 1/1 0/1
12 170 250 251 1 3 1 1003 1003 1001 1 1 0/1 1/1 1/1
12 148 170 171 1 4 0 3 3 1 1 1 1/1 1/1 1/1
12 250 319 320 1 2 0 3 3 1 1 1 1/1 1/1 1/1
12 251 0 0 -1 2 1 1003 1003 1001 1 1 0/1 1/1 1/1
12 149 148 89 1 5 0 3 3 1 1 1 1/1 1/1 1/1
12 150 148 89 1 5 0 3 3 1 1 1 1/1 1/1 1/1
12 151 148 89 -1 5 0 3 3 1 1 1 1/1 1/1 1/1
12 152 148 89 1 5 0 3 3 1 1 1 1/1 1/1 1/1
12 153 148 89 1 5 0 3 3 1 1 1 1/1 1/1 1/1
12 154 148 89 1 5 0 3 3 1 1 1 1/1 1/1 1/1
12 171 384 385 -1 3 1 1003 1003 1001 1 1 0/1 1/1 1/1
12 319 0 0 1 1 0 3 3 1 1 1 1/1 1/1 1/1
12 320 0 0 -1 1 0 4 4 1 1 2 1/1 1/1 0/1
12 89 42 87 -1 4 0 3 3 1 1 1 1/1 1/1 1/1
12 365 152 364 1 6 0 3 3 1 1 1 1/1 1/1 1/1
12 366 152 364 1 6 0 3 3 1 1 1 1/1 1/1 1/1
12 367 152 364 1 6 0 3 3 1 1 1 1/1 1/1 1/1
12 368 152 364 1 6 0 3 3 1 1 1 1/1 1/1 1/1
12 369 152 364 1 6 0 3 3 1 1 1 1/1 1/1 1/1
12 370 152 364 -1 6 0 3 3 1 1 1 1/1 1/1 1/1
12 384 0 0 1 2 0 3 3 1 1 1 1/1 1/1 1/1
12 385 0 0 -1 2 1 1004 1004 1001 1 2 0/1 1/1 0/1
12 42 0 0 1 3 0 3 3 1 1 1 1/1 1/1 1/1
12 87 0 0 -1 3 0 3 3 1 1 1 1/1 1/1 1/1
12 364 0 0 -1 5 0 3 3 1 1 1 1/1 1/1 1/1
13 60 19 54 1 3 1 1003 1003 1001 1 1 0/1 1/1 1/1
13 19 10 18 1 2 1 1003 1003 1001 1 1 0/1 1/1 1/1
13 54 0 0 -1 2 0 4 4 1 1 2 1/1 1/1 0/1
13 10 1 9 1 1 0 3 3 1 1 1 1/1 1/1 1/1
13 18 252 253 -1 1 1 1004 1004 1001 1 2 0/1 1/1 0/1
13 55 19 54 1 3 1 1003 1003 1001 1 1 0/1 1/1 1/1
13 56 19 54 -1 3 0 3 3 1 1 1 1/1 1/1 1/1
13 57 19 54 -1 3 0 3 3 1 1 1 1/1 1/1 1/1
13 58 19 54 -1 3 0 4 4 1 1 2 1/1 1/1 0/1
13 59 19 54 -1 3 0 3 3 1 1 1 1/1 1/1 1/1
13 61 19 54 -1 3 1 1004 1004 1001 1 2 0/1 1/1 0/1
13 62 19 54 -1 3 1 1003 1003 1001 1 1 0/1 1/1 1/1
13 63 19 54 1 3 0 4 4 1 1 2 1/1 1/1 0/1
13 1 153 154 1 0 0 3 3 1 1 1 1/1 1/1 1/1
13 9 24 25 -1 0 0 3 3 1 1 1 1/1 1/1 1/1
13 20 10 18 1 2 1 1003 1003 1001 1 1 0/1 1/1 1/1
13 252 283 284 1 0 0 4 4 1 1 2 1/1 1/1 0/1
13 253 0 0 -1 0 1 1003 1003 1001 1 1 0/1 1/1 1/1
13 69 68 56 1 4 0 3 3 1 1 1 1/1 1/1 1/1
13 70 68 56 1 4 0 3 3 1 1 1 1/1 1/1 1/1
13 71 68 56 -1 4 0 3 3 1 1 1 1/1 1/1 1/1
13 72 68 56 1 4 0 3 3 1 1 1 1/1 1/1 1/1
13 73 68 56 -1 4 0 3 3 1 1 1 1/1 1/1 1/1
13 74 68 56 -1 4 0 3 3 1 1 1 1/1 1/1 1/1
13 75 68 56 1 4 0 3 3 1 1 1 1/1 1/1 1/1
13 115 114 59 1 4 1 1003 1003 1001 1 1 0/1 1/1 1/1
13 217 63 216 -1 4 0 3 3 1 1 1 1/1 1/1 1/1
13 218 63 216 1 4 0 4 4 1 1 2 1/1 1/1 0/1
13 219 63 216 -1 4 0 3 3 1 1 1 1/1 1/1 1/1
13 220 63 216 1 4 0 4 4 1 1 2 1/1 1/1 0/1
13 221 63 216 -1 4 0 4 4 1 1 2 1/1 1/1 0/1
13 222 63 216 -1 4 0 3 3 1 1 1 1/1 1/1 1/1
13 11 1 9 1 1 0 3 3 1 1 1 1/1 1/1 1/1
13 12 1 9 -1 1 0 3 3 1 1 1 1/1 1/1 1/1
13 13 1 9 1 1 0 3 3 1 1 1 1/1 1/1 1/1
13 14 1 9 -1 1 0 3 3 1 1 1 1/1 1/1 1/1
13 153 0 0 1 -1 0 3 3 1 1 1 1/1 1/1 1/1
13 154 0 0 -1 -1 0 3 3 1 1 1 1/1 1/1 1/1
13 24 0 0 1 -1 0 3 3 1 1 1 1/1 1/1 1/1
13 25 0 0 -1 -1 0 3 3 1 1 1 1/1 1/1 1/1
13 27 20 26 1 3 1 1003 1003 1001 1 1 0/1 1/1 1/1
13 28 20 26 1 3 0 3 3 1 1 1 1/1 1/1 1/1
13 29 20 26 -1 3 1 1003 1003 1001 1 1 0/1 1/1 1/1
13 30 20 26 1 3 0 3 3 1 1 1 1/1 1/1 1/1
13 31 20 26 1 3 0 3 3 1 1 1 1/1 1/1 1/1
13 283 0 0 1 -1 0 4 4 1 1 2 1/1 1/1 0/1
13 284 0 0 -1 -1 0 3 3 1 1 1 1/1 1/1 1/1
13 68 0 0 1 3 0 3 3 1 1 1 1/1 1/1 1/1
13 96 70 95 -1 5 0 3 3 1 1 1 1/1 1/1 1/1
13 97 70 95 -1 5 0 3 3 1 1 1 1/1 1/1 1/1
13 98 70 95 -1 5 0 3 3 1 1 1 1/1 1/1 1/1
13 99 70 95 1 5 0 3 3 1 1 1 1/1 1/1 1/1
13 100 70 95 -1 5 0 3 3 1 1 1 1/1 1/1 1/1
13 231 230 73 -1 5 0 3 3 1 1 1 1/1 1/1 1/1
13 232 230 73 1 5 0 3 3 1 1 1 1/1 1/1 1/1
13 233 230 73 1 5 0 3 3 1 1 1 1/1 1/1 1/1
13 234 230 73 1 5 0 3 3 1 1 1 1/1 1/1 1/1
13 235 230 73 -1 5 0 3 3 1 1 1 1/1 1/1 1/1
13 236 230 73 1 5 0 3 3 1 1 1 1/1 1/1 1/1
13 92 75 91 1 5 0 3 3 1 1 1 1/1 1/1 1/1
13 93 75 91 -1 5 0 3 3 1 1 1 1/1 1/1 1/1
13 94 75 91 -1 5 0 3 3 1 1 1 1/1 1/1 1/1
13 114 0 0 1 3 1 1003 1003 1001 1 1 0/1 1/1 1/1
13 216 0 0 -1 3 0 3 3 1 1 1 1/1 1/1 1/1
13 230 0 0 1 4 0 3 3 1 1 1 1/1 1/1 1/1
13 26 0 0 -1 2 0 3 3 1 1 1 1/1 1/1 1/1
13 95 0 0 -1 4 0 3 3 1 1 1 1/1 1/1 1/1
13 91 0 0 -1 4 0 3 3 1 1 1 1/1 1/1 1/1
14 279 278 147 -1 5 1 1003 1003 1001 1 1 0/1 1/1 1/1
14 147 63 143 -1 4 0 3 3 1 1 1 1/1 1/1 1/1
14 278 0 0 1 4 1 1003 1003 1001 1 1 0/1 1/1 1/1
14 285 284 279 -1 6 0 3 3 1 1 1 1/1 1/1 1/1
14 286 284 279 1 6 1 1003 1003 1001 1 1 0/1 1/1 1/1
14 287 284 279 1 6 1 1003 1003 1001 1 1 0/1 1/1 1/1
14 63 14 59 1 3 0 3 3 1 1 1 1/1 1/1 1/1
14 143 0 0 -1 3 0 3 3 1 1 1 1/1 1/1 1/1
14 280 278 147 -1 5 0 3 3 1 1 1 1/1 1/1 1/1
14 281 278 147 1 5 1 1003 1003 1001 1 1 0/1 1/1 1/1
14 284 0 0 1 5 0 3 3 1 1 1 1/1 1/1 1/1
14 14 12 8 1 2 0 3 3 1 1 1 1/1 1/1 1/1
14 59 0 0 -1 2 0 3 3 1 1 1 1/1 1/1 1/1
14 144 63 143 1 4 0 3 3 1 1 1 1/1 1/1 1/1
14 145 63 143 1 4 0 3 3 1 1 1 1/1 1/1 1/1
14 146 63 143 -1 4 0 3 3 1 1 1 1/1 1/1 1/1
14 8 0 0 -1 1 0 3 3 1 1 1 1/1 1/1 1/1
14 12 0 0 1 1 0 3 3 1 1 1 1/1 1/1 1/1
14 60 14 59 1 3 0 3 3 1 1 1 1/1 1/1 1/1
14 61 14 59 -1 3 0 3 3 1 1 1 1/1 1/1 1/1
14 62 14 59 1 3 0 3 3 1 1 1 1/1 1/1 1/1
14 180 145 179 -1 5 0 3 3 1 1 1 1/1 1/1 1/1
14 181 145 179 1 5 0 3 3 1 1 1 1/1 1/1 1/1
14 182 145 179 -1 5 0 3 3 1 1 1 1/1 1/1 1/1
14 183 145 179 1 5 0 3 3 1 1 1 1/1 1/1 1/1
14 179 0 0 -1 4 0 3 3 1 1 1 1/1 1/1 1/1
15 219 218 37 -1 3 1 1003 1003 1001 1 1 0/1 1/1 1/1
15 37 31 36 -1 2 0 4 4 1 1 2 1/1 1/1 0/1
15 218 0 0 1 2 1 1003 1003 1001 1 1 0/1 1/1 1/1
15 31 28 2 1 1 0 3 3 1 1 1 1/1 1/1 1/1
15 36 0 0 -1 1 0 4 4 1 1 2 1/1 1/1 0/1
15 220 218 37 1 3 0 4 4 1 1 2 1/1 1/1 0/1
15 221 218 37 -1 3 1 1004 1004 1001 1 2 0/1 1/1 0/1
15 222 218 37 1 3 1 1003 1003 1001 1 1 0/1 1/1 1/1
15 223 218 37 -1 3 1 1003 1003 1001 1 1 0/1 1/1 1/1
15 224 218 37 1 3 0 4 4 1 1 2 1/1 1/1 0/1
15 225 218 37 -1 3 0 3 3 1 1 1 1/1 1/1 1/1
15 226 218 37 -1 3 0 3 3 1 1 1 1/1 1/1 1/1
15 2 0 0 -1 0 0 3 3 1 1 1 1/1 1/1 1/1
15 28 34 35 1 0 0 3 3 1 1 1 1/1 1/1 1/1
15 38 31 36 1 2 0 4 4 1 1 2 1/1 1/1 0/1
15 39 31 36 1 2 0 3 3 1 1 1 1/1 1/1 1/1
15 40 31 36 1 2 0 3 3 1 1 1 1/1 1/1 1/1
15 41 31 36 1 2 0 4 4 1 1 2 1/1 1/1 0/1
15 42 31 36 1 2 0 3 3 1 1 1 1/1 1/1 1/1
15 43 31 36 -1 2 0 4 4 1 1 2 1/1 1/1 0/1
15 228 227 223 1 4 0 3 3 1 1 1 1/1 1/1 1/1
15 229 227 223 -1 4 0 3 3 1 1 1 1/1 1/1 1/1
15 230 227 223 1 4 0 3 3 1 1 1 1/1 1/1 1/1
15 231 227 223 -1 4 1 1003 1003 1001 1 1 0/1 1/1 1/1
15 232 227 223 1 4 0 3 3 1 1 1 1/1 1/1 1/1
15 29 28 2 -1 1 0 3 3 1 1 1 1/1 1/1 1/1
15 30 28 2 1 1 0 3 3 1 1 1 1/1 1/1 1/1
15 32 28 2 -1 1 0 3 3 1 1 1 1/1 1/1 1/1
15 33 28 2 -1 1 0 3 3 1 1 1 1/1 1/1 1/1
15 34 0 0 1 -1 0 4 4 1 1 2 1/1 1/1 0/1
15 35 0 0 -1 -1 0 3 3 1 1 1 1/1 1/1 1/1
15 102 38 101 1 3 0 3 3 1 1 1 1/1 1/1 1/1
15 103 38 101 1 3 0 4 4 1 1 2 1/1 1/1 0/1
15 104 38 101 1 3 0 3 3 1 1 1 1/1 1/1 1/1
15 163 39 162 1 3 0 3 3 1 1 1 1/1 1/1 1/1
15 164 39 162 1 3 0 4 4 1 2 1 1/1 0/1 1/1
15 165 39 162 1 3 0 3 3 1 1 1 1/1 1/1 1/1
15 166 39 162 -1 3 0 4 4 1 2 1 1/1 0/1 1/1
15 98 42 97 -1 3 0 3 3 1 1 1 1/1 1/1 1/1
15 99 42 97 -1 3 0 3 3 1 1 1 1/1 1/1 1/1
15 100 42 97 -1 3 0 3 3 1 1 1 1/1 1/1 1/1
15 227 0 0 1 3 0 3 3 1 1 1 1/1 1/1 1/1
15 101 0 0 -1 2 0 3 3 1 1 1 1/1 1/1 1/1
15 162 0 0 -1 2 0 4 4 1 2 1 1/1 0/1 1/1
15 97 0 0 -1 2 0 3 3 1 1 1 1/1 1/1 1/1
16 124 81 123 1 5 1 1003 1003 1001 1 1 0/1 1/1 1/1
16 81 76 25 1 4 0 3 3 1 1 1 1/1 1/1 1/1
16 123 0 0 -1 4 1 1003 1003 1001 1 1 0/1 1/1 1/1
16 25 20 15 -1 3 0 3 3 1 1 1 1/1 1/1 1/1
16 76 0 0 1 3 1 1003 1003 1001 1 1 0/1 1/1 1/1
16 125 81 123 1 5 1 1003 1003 1001 1 1 0/1 1/1 1/1
16 126 81 123 -1 5 0 3 3 1 1 1 1/1 1/1 1/1
16 127 81 123 1 5 0 3 3 1 1 1 1/1 1/1 1/1
16 15 8 11 -1 2 0 3 3 1 1 1 1/1 1/1 1/1
16 20 37 38 1 2 1 1003 1003 1001 1 1 0/1 1/1 1/1
16 77 76 25 -1 4 0 3 3 1 1 1 1/1 1/1 1/1
16 78 76 25 -1 4 0 3 3 1 1 1 1/1 1/1 1/1
16 79 76 25 1 4 1 1003 1003 1001 1 1 0/1 1/1 1/1
16 80 76 25 1 4 0 3 3 1 1 1 1/1 1/1 1/1
16 82 76 25 1 4 1 1003 1003 1001 1 1 0/1 1/1 1/1
16 83 76 25 1 4 1 1003 1003 1001 1 1 0/1 1/1 1/1
16 84 76 25 1 4 0 3 3 1 1 1 1/1 1/1 1/1
16 8 0 0 1 1 0 3 3 1 1 1 1/1 1/1 1/1
16 11 0 0 -1 1 1 1003 1003 1001 1 1 0/1 1/1 1/1
16 21 20 15 -1 3 1 1003 1003 1001 1 1 0/1 1/1 1/1
16 22 20 15 1 3 1 1003 1003 1001 1 1 0/1 1/1 1/1
16 23 20 15 -1 3 1 1003 1003 1001 1 1 0/1 1/1 1/1
16 24 20 15 -1 3 1 1003 1003 1001 1 1 0/1 1/1 1/1
16 26 20 15 1 3 0 3 3 1 1 1 1/1 1/1 1/1
16 37 0 0 1 1 1 1003 1003 1001 1 1 0/1 1/1 1/1
16 38 0 0 -1 1 1 1003 1003 1001 1 1 0/1 1/1 1/1
16 350 349 78 -1 5 0 3 3 1 1 1 1/1 1/1 1/1
16 351 349 78 -1 5 0 3 3 1 1 1 1/1 1/1 1/1
16 352 349 78 -1 5 0 4 4 1 1 2 1/1 1/1 0/1
16 353 349 78 1 5 0 3 3 1 1 1 1/1 1/1 1/1
16 354 349 78 -1 5 0 4 4 1 1 2 1/1 1/1 0/1
16 120 79 119 1 5 1 1003 1003 1001 1 1 0/1 1/1 1/1
16 121 79 119 -1 5 0 3 3 1 1 1 1/1 1/1 1/1
16 122 79 119 -1 5 1 1003 1003 1001 1 1 0/1 1/1 1/1
16 226 83 225 1 5 0 3 3 1 1 1 1/1 1/1 1/1
16 227 83 225 1 5 0 3 3 1 1 1 1/1 1/1 1/1
16 112 84 111 1 5 0 3 3 1 1 1 1/1 1/1 1/1
16 113 84 111 1 5 0 4 4 1 1 2 1/1 1/1 0/1
16 114 84 111 -1 5 0 4 4 1 1 2 1/1 1/1 0/1
16 115 84 111 -1 5 0 3 3 1 1 1 1/1 1/1 1/1
16 116 84 111 -1 5 0 3 3 1 1 1 1/1 1/1 1/1
16 349 0 0 1 4 0 4 4 1 1 2 1/1 1/1 0/1
16 119 0 0 -1 4 0 3 3 1 1 1 1/1 1/1 1/1
16 225 0 0 -1 4 0 3 3 1 1 1 1/1 1/1 1/1
16 111 0 0 -1 4 0 4 4 1 1 2 1/1 1/1 0/1
17 108 103 58 -1 6 1 1003 1003 1001 1 1 0/1 1/1 1/1
17 58 48 57 -1 5 0 3 3 1 1 1 1/1 1/1 1/1
17 103 344 345 1 5 1 1003 1003 1001 1 1 0/1 1/1 1/1
17 355 354 108 1 7 1 1003 1003 1001 1 1 0/1 1/1 1/1
17 356 354 108 -1 7 1 1003 1003 1001 1 1 0/1 1/1 1/1
17 357 354 108 1 7 0 3 3 1 1 1 1/1 1/1 1/1
17 358 354 108 1 7 0 3 3 1 1 1 1/1 1/1 1/1
17 359 354 108 1 7 1 1003 1003 1001 1 1 0/1 1/1 1/1
17 48 27 45 1 4 0 3 3 1 1 1 1/1 1/1 1/1
17 57 99 100 -1 4 0 3 3 1 1 1 1/1 1/1 1/1
17 104 103 58 -1 6 0 3 3 1 1 1 1/1 1/1 1/1
17 105 103 58 1 6 1 1003 1003 1001 1 1 0/1 1/1 1/1
17 106 103 58 1 6 0 3 3 1 1 1 1/1 1/1 1/1
17 107 103 58 1 6 0 3 3 1 1 1 1/1 1/1 1/1
17 109 103 58 1 6 1 1003 1003 1001 1 1 0/1 1/1 1/1
17 110 103 58 1 6 0 3 3 1 1 1 1/1 1/1 1/1
17 344 0 0 1 4 1 1003 1003 1001 1 1 0/1 1/1 1/1
17 345 0 0 -1 4 1 1003 1003 1001 1 1 0/1 1/1 1/1
17 354 0 0 1 6 0 3 3 1 1 1 1/1 1/1 1/1
17 27 26 22 1 3 0 3 3 1 1 1 1/1 1/1 1/1
17 45 0 0 -1 3 0 3 3 1 1 1 1/1 1/1 1/1
17 59 48 57 1 5 0 3 3 1 1 1 1/1 1/1 1/1
17 60 48 57 1 5 0 3 3 1 1 1 1/1 1/1 1/1
17 61 48 57 1 5 0 3 3 1 1 1 1/1 1/1 1/1
17 62 48 57 1 5 0 3 3 1 1 1 1/1 1/1 1/1
17 99 0 0 1 3 0 4 4 1 1 2 1/1 1/1 0/1
17 100 0 0 -1 3 0 3 3 1 1 1 1/1 1/1 1/1
17 251 105 250 1 7 0 3 3 1 1 1 1/1 1/1 1/1
17 252 105 250 -1 7 0 3 3 1 1 1 1/1 1/1 1/1
17 22 0 0 -1 2 0 3 3 1 1 1 1/1 1/1 1/1
17 26 0 0 1 2 0 3 3 1 1 1 1/1 1/1 1/1
17 46 27 45 -1 4 0 3 3 1 1 1 1/1 1/1 1/1
17 47 27 45 1 4 0 3 3 1 1 1 1/1 1/1 1/1
17 49 27 45 -1 4 0 3 3 1 1 1 1/1 1/1 1/1
17 50 27 45 -1 4 0 3 3 1 1 1 1/1 1/1 1/1
17 258 60 257 -1 6 0 3 3 1 1 1 1/1 1/1 1/1
17 259 60 257 -1 6 0 3 3 1 1 1 1/1 1/1 1/1
17 260 60 257 -1 6 0 3 3 1 1 1 1/1 1/1 1/1
17 261 60 257 1 6 0 3 3 1 1 1 1/1 1/1 1/1
17 262 60 257 -1 6 0 3 3 1 1 1 1/1 1/1 1/1
17 263 60 257 1 6 0 3 3 1 1 1 1/1 1/1 1/1
17 264 60 257 1 6 0 3 3 1 1 1 1/1 1/1 1/1
17 265 60 257 1 6 0 3 3 1 1 1 1/1 1/1 1/1
17 250 0 0 -1 6 0 3 3 1 1 1 1/1 1/1 1/1
17 257 0 0 -1 5 0 3 3 1 1 1 1/1 1/1 1/1
18 64 0 0 1 1 1 1003 1003 1001 1 1 0/1 1/1 1/1
18 65 64 19 1 2 0 3 3 1 1 1 1/1 1/1 1/1
18 66 64 19 -1 2 0 3 3 1 1 1 1/1 1/1 1/1
18 19 3 16 -1 1 0 3 3 1 1 1 1/1 1/1 1/1
18 68 67 66 1 3 0 4 4 1 1 2 1/1 1/1 0/1
18 69 67 66 -1 3 0 3 3 1 1 1 1/1 1/1 1/1
18 70 67 66 1 3 1 1003 1003 1001 1 1 0/1 1/1 1/1
18 71 67 66 1 3 1 1003 1003 1001 1 1 0/1 1/1 1/1
18 72 67 66 1 3 1 1003 1003 1001 1 1 0/1 1/1 1/1
18 73 67 66 -1 3 0 3 3 1 1 1 1/1 1/1 1/1
18 74 67 66 1 3 1 1004 1004 1001 1 2 0/1 1/1 0/1
18 75 67 66 -1 3 1 1003 1003 1001 1 1 0/1 1/1 1/1
18 3 0 0 1 0 0 3 3 1 1 1 1/1 1/1 1/1
18 16 25 26 -1 0 0 3 3 1 1 1 1/1 1/1 1/1
18 67 0 0 1 2 1 1004 1004 1001 1 2 0/1 1/1 0/1
18 107 68 106 1 4 0 5 5 1 1 3 1/1 1/1 0/0
18 108 68 106 -1 4 0 5 5 1 1 3 1/1 1/1 0/0
18 109 68 106 -1 4 0 4 4 1 1 2 1/1 1/1 0/1
18 110 68 106 1 4 0 5 5 1 1 3 1/1 1/1 0/0
18 111 68 106 1 4 0 4 4 1 1 2 1/1 1/1 0/1
18 112 68 106 1 4 0 4 4 1 1 2 1/1 1/1 0/1
18 96 95 69 -1 4 0 3 3 1 1 1 1/1 1/1 1/1
18 97 95 69 1 4 0 4 4 1 2 1 1/1 0/1 1/1
18 151 71 150 1 4 0 3 3 1 1 1 1/1 1/1 1/1
18 152 71 150 -1 4 0 3 3 1 1 1 1/1 1/1 1/1
18 17 3 16 1 1 0 3 3 1 1 1 1/1 1/1 1/1
18 18 3 16 1 1 0 3 3 1 1 1 1/1 1/1 1/1
18 20 3 16 -1 1 0 3 3 1 1 1 1/1 1/1 1/1
18 25 0 0 1 -1 0 3 3 1 1 1 1/1 1/1 1/1
18 26 0 0 -1 -1 0 3 3 1 1 1 1/1 1/1 1/1
18 106 0 0 -1 3 0 5 5 1 1 3 1/1 1/1 0/0
18 254 107 253 -1 5 0 4 4 1 1 2 1/1 1/1 0/1
18 255 107 253 -1 5 0 4 4 1 1 2 1/1 1/1 0/1
18 256 107 253 1 5 0 4 4 1 1 2 1/1 1/1 0/1
18 95 0 0 1 3 0 4 4 1 2 1 1/1 0/1 1/1
18 150 0 0 -1 3 0 3 3 1 1 1 1/1 1/1 1/1
18 157 151 156 -1 5 0 3 3 1 1 1 1/1 1/1 1/1
18 158 151 156 -1 5 0 3 3 1 1 1 1/1 1/1 1/1
18 159 151 156 1 5 0 3 3 1 1 1 1/1 1/1 1/1
18 160 151 156 -1 5 0 3 3 1 1 1 1/1 1/1 1/1
18 253 0 0 -1 4 0 3 3 1 1 1 1/1 1/1 1/1
18 156 0 0 -1 4 0 3 3 1 1 1 1/1 1/1 1/1
19 288 162 286 -1 5 1 1003 1003 1001 1 1 0/1 1/1 1/1
19 162 75 156 1 4 0 3 3 1 1 1 1/1 1/1 1/1
19 286 0 0 -1 4 1 1003 1003 1001 1 1 0/1 1/1 1/1
19 75 22 71 1 3 0 3 3 1 1 1 1/1 1/1 1/1
19 156 282 283 -1 3 0 3 3 1 1 1 1/1 1/1 1/1
19 287 162 286 -1 5 1 1003 1003 1001 1 1 0/1 1/1 1/1
19 289 162 286 1 5 0 3 3 1 1 1 1/1 1/1 1/1
19 22 20 12 1 2 1 1003 1003 1001 1 1 0/1 1/1 1/1
19 71 0 0 -1 2 0 3 3 1 1 1 1/1 1/1 1/1
19 157 75 156 -1 4 0 3 3 1 1 1 1/1 1/1 1/1
19 158 75 156 -1 4 0 3 3 1 1 1 1/1 1/1 1/1
19 159 75 156 -1 4 0 3 3 1 1 1 1/1 1/1 1/1
19 160 75 156 1 4 0 3 3 1 1 1 1/1 1/1 1/1
19 161 75 156 1 4 0 3 3 1 1 1 1/1 1/1 1/1
19 282 0 0 1 2 0 3 3 1 1 1 1/1 1/1 1/1
19 283 0 0 -1 2 0 4 4 1 1 2 1/1 1/1 0/1
19 12 0 0 -1 1 1 1003 1003 1001 1 1 0/1 1/1 1/1
19 20 0 0 1 1 0 3 3 1 1 1 1/1 1/1 1/1
19 72 22 71 1 3 1 1003 1003 1001 1 1 0/1 1/1 1/1
19 73 22 71 -1 3 1 1003 1003 1001 1 1 0/1 1/1 1/1
19 74 22 71 -1 3 0 3 3 1 1 1 1/1 1/1 1/1
19 76 22 71 -1 3 1 1003 1003 1001 1 1 0/1 1/1 1/1
19 77 22 71 -1 3 0 3 3 1 1 1 1/1 1/1 1/1
20 214 185 213 -1 5 1 1004 1004 1001 2 1 0/1 0/1 1/1
20 185 181 89 1 4 1 1005 1005 1001 2 2 0/1 0/1 0/1
20 213 281 282 -1 4 0 3 3 1 1 1 1/1 1/1 1/1
20 89 43 88 -1 3 0 4 4 1 2 1 1/1 0/1 1/1
20 181 0 0 1 3 1 1005 1005 1001 1 3 0/1 1/1 0/0
20 215 185 213 -1 5 0 4 4 1 1 2 1/1 1/1 0/1
20 216 185 213 1 5 0 5 5 1 2 2 1/1 0/1 0/1
20 217 185 213 1 5 1 1004 1004 1001 2 1 0/1 0/1 1/1
20 218 185 213 -1 5 0 3 3 1 1 1 1/1 1/1 1/1
20 219 185 213 -1 5 1 1004 1004 1001 1 2 0/1 1/1 0/1
20 281 0 0 1 3 0 3 3 1 1 1 1/1 1/1 1/1
20 282 0 0 -1 3 0 3 3 1 1 1 1/1 1/1 1/1
20 43 14 40 1 2 0 3 3 1 1 1 1/1 1/1 1/1
20 88 0 0 -1 2 0 5 5 1 3 1 1/1 0/0 1/1
20 182 181 89 -1 4 1 1005 1005 1001 2 2 0/1 0/1 0/1
20 183 181 89 -1 4 0 5 5 1 2 2 1/1 0/1 0/1
20 184 181 89 -1 4 1 1005 1005 1001 2 2 0/1 0/1 0/1
20 14 0 0 1 1 0 3 3 1 1 1 1/1 1/1 1/1
20 40 0 0 -1 1 0 3 3 1 1 1 1/1 1/1 1/1
20 90 43 88 -1 3 0 4 4 1 2 1 1/1 0/1 1/1
20 91 43 88 -1 3 0 4 4 1 2 1 1/1 0/1 1/1
20 92 43 88 1 3 0 4 4 1 2 1 1/1 0/1 1/1
20 93 43 88 1 3 0 4 4 1 2 1 1/1 0/1 1/1
20 94 43 88 1 3 0 4 4 1 2 1 1/1 0/1 1/1
pedigree
{
delimiters = "\t"
delimiter_mode = "single"
individual_missing_value = "0"
sex_code, male = "1", female = "-1", missing = ".", trait
pedigree_id = "family"
individual_id = "ind"
parent_id = "father"
parent_id = "mother"
sex_field = "sex"
trait = "affection", binary, affected = "1", unaffected = "0", missing = "."
}
segreg, out = "SEGREG2"
{
title = "y"
trait = "affection"
}
my.pedigree <- function(n.init=3,n.iter=20,mean.kid=2){
trio <- matrix(0,n.init,5)
m <- length(trio[1,])
trio[1:n.init,1] <- 1:n.init
trio[1:n.init,4] <- sample(c(-1,1),n.init,replace=TRUE)
trio[,5] <- 0
cnt <- n.init
tmp.f <- function(m,k){
(m/(1-exp(-m))-k)^2
}
xmin <- optimize(tmp.f, c(0, mean.kid), tol = 0.0001, k = mean.kid)
mean.kid. <- xmin[[1]]
for(i in 1:n.iter){
s <- sample(trio[,1],1)
if(length(which(trio[,2:3]==s))==0){
trio <- rbind(trio,c(cnt+1,rep(0,m-1)))
dp <- dpois(1:100,mean.kid.)
dp <- dp/sum(dp)
num.kids <- sample(1:100,1,prob=dp)
tmp <- matrix(0,num.kids,m)
tmp[,1] <- cnt+1+(1:num.kids)
if(trio[s,4]==1){
tmp[,2] <- s
tmp[,3] <- cnt+1
}else{
tmp[,2] <- cnt+1
tmp[,3] <- s
}
tmp[,4] <- sample(c(-1,1),num.kids,replace=TRUE)
tmp[,5] <- trio[s,5] + 1
trio <- rbind(trio,tmp)
trio[cnt+1,4] <- -trio[s,4]
trio[cnt+1,5] <- trio[s,5]
cnt <- cnt+1+num.kids
}else if(trio[s,2] == 0 & trio[s,3] == 0){
trio <- rbind(trio,c(cnt+1,rep(0,m-1)))
trio <- rbind(trio,c(cnt+2,rep(0,m-1)))
trio[cnt+1,4] <- 1
trio[cnt+2,4] <- -1
trio[c(cnt+1,cnt+2),5] <- trio[s,5]-1
trio[s,2] <- cnt+1
trio[s,3] <- cnt+2
cnt <- cnt+2
}
}
trio
}
make.genotype <- function(trio,prob){
ret <- rep(0,length(trio[,1]))
ord <- order(trio[,5])
for(i in 1:length(ord)){
if(trio[ord[i],2]==0 & trio[ord[i],3]==0){
ret[ord[i]] <- sample(c(0,1,2),1,prob=prob)
}else{
parents <- trio[ord[i],2:3]
pat <- ret[trio[ord[i],2]]/2
mat <- ret[trio[ord[i],3]]/2
if(pat==0.5){
pat <- sample(c(0,1),1)
}
if(mat==0.5){
mat <- sample(c(0,1),1)
}
ret[ord[i]] <- pat+mat
}
}
ret
}
make.affected.pedigree <- function(trio,p.g,r.g,h,prev){
genotype <- matrix(0,length(trio[,1]),n.loci)
for(i in 1:n.loci){
genotype[,i] <- make.genotype(trio,p.g[i,])
}
R.mat <- genotype
for(i in 1:n.loci){
tmp <- r.g[i,]
R.mat[,i] <- tmp[genotype[,i]+1]
}
R.sum <- apply(R.mat,1,sum)
var.R <- var(R.sum)
var.E <- var.R*(1-herit)/herit
E <- rnorm(length(trio[,1]),0,sqrt(var.E))
R.E <- R.sum + E
affected <- rep(0,length(trio[,1]))
affected[which(R.E>= prev)] <- 1
return(list(pedigree=trio,genotype=genotype,gen.risk.mat=R.mat,gen.risk=R.sum,risk=R.E,affected=affected))
}
n.init <- 10
n.iter <- 50
mean.kid <-2
n.loci <- 10
p.a <- runif(n.loci)*0.1
p.g <- cbind(p.a^2,2*p.a*(1-p.a),(1-p.a)^2)
r.a <- c(100,rep(1,n.loci-1))
r.g <- cbind(r.a*2+1,r.a+1,rep(1,n.loci))
herit <- 0.9
prev <- r.a[1]*0.9
my.ped <- my.pedigree(n.init=n.init,n.iter=n.iter,mean.kid=mean.kid)
aff.pedigree <- make.affected.pedigree(my.ped,p.g,r.g,herit,prev)
library(kinship2)
plot.affected.pedigree <- function(ap){
sex <- ap$pedigree[,4]
sex[which(sex==-1)] <- 2
my.ped <- pedigree(ap$pedigree[,1],ap$pedigree[,2],ap$pedigree[,3],sex,ap$affected)
plot(my.ped,symbolosize=0.001)
}
plot.affected.pedigree(aff.pedigree)
library(igraph)
affected.pedigree <- function(pedigree,genotype,gen.risk.mat,gen.risk,risk,affected){
return(list(pedigree=pedigree,genotype=genotype,gen.risk.mat=gen.risk.mat,gen.risk=gen.risk,risk=risk,affected=affected))
}
el.from.trio <- function(trio){
non.zero <- which(apply((trio[,2:3])^2,1,sum)!=0)
el <- rbind(cbind(trio[non.zero,2],trio[non.zero,1]),cbind(trio[non.zero,3],trio[non.zero,1]))
el
}
sample.aff.sub <- function(g,aff.pedigree,k,n=NULL){
if(is.null(n)){
n <- sample(which(aff.pedigree$affected==1),1)
}
nb <- neighborhood(g,k,n,mode="all")[[1]]
ped <- aff.pedigree$pedigree[nb,]
tmp <- ped[,2:3]
tmp[!is.element(tmp,ped[,1])] <- 0
no.parents <- which(apply(tmp,1,sum)==0)
ped[no.parents,2:3] <- 0
tmp2 <- ped[,2:3]
tobe.added <- unique(tmp2[which(!is.element(tmp2,ped[,1]))])
zero <- which(tobe.added==0)
tobe.added <- tobe.added[-zero]
nb <- c(nb,tobe.added)
ped <- aff.pedigree$pedigree[nb,]
tmp <- ped[,2:3]
tmp[!is.element(tmp,ped[,1])] <- 0
ped[,2:3] <- tmp
affected.pedigree(ped,aff.pedigree$genotype[nb,],aff.pedigree$gen.risk.mat[nb,],aff.pedigree$gen.risk[nb],aff.pedigree$risk[nb],aff.pedigree$affected[nb])
}
n.loci <- 3
p.a <- runif(n.loci)*0.1
p.g <- cbind(p.a^2,2*p.a*(1-p.a),(1-p.a)^2)
r.a <- c(1000,rep(1,n.loci-1))
r.g <- cbind(r.a*2+1,r.a+1,rep(1,n.loci))
herit <- 1
prev <- r.a[1]*0.9
n.init <- 5
n.iter <- 100
mean.kid <-4
n.family <- 20
k<- 4
F <- list()
for(i in 1:n.family){
my.ped <- my.pedigree(n.init=n.init,n.iter=n.iter,mean.kid=mean.kid)
aff.pedigree <- make.affected.pedigree(my.ped,p.g,r.g,herit,prev)
el <- el.from.trio(aff.pedigree$pedigree)
g <- graph.edgelist(el)
smpl.aff.ped <- sample.aff.sub(g,aff.pedigree,k)
F[[i]] <- smpl.aff.ped
}
plot.affected.pedigree(F[[1]])
genotype.change <- function(g.mat){
ret <- g.mat
ret[which(ret==0)] <- "0/0"
ret[which(ret==1)] <- "0/1"
ret[which(ret==2)] <- "1/1"
ret
}
make.multi.family.table <- function(F){
family.id <- rep(1,length(F[[1]]$pedigree[,1]))
pedigree <- F[[1]]$pedigree
genotype <- genotype.change(F[[1]]$genotype)
gen.risk.mat <- F[[1]]$gen.risk.mat
gen.risk <- F[[1]]$gen.risk
risk <- F[[1]]$risk
affected <- F[[1]]$affected
if(length(F)>1){
for(i in 2:length(F)){
family.id <- c(family.id,rep(i,length(F[[i]]$pedigree[,1])))
pedigree <- rbind(pedigree,F[[i]]$pedigree)
genotype <- rbind(genotype,genotype.change(F[[i]]$genotype))
gen.risk.mat <- rbind(gen.risk.mat,F[[i]]$gen.risk.mat)
gen.risk <- c(gen.risk,F[[i]]$gen.risk)
risk <- c(risk,F[[i]]$risk)
affected <- c(affected,F[[i]]$affected)
}
}
tmp <- cbind(family.id,pedigree,affected,risk,gen.risk,gen.risk.mat,genotype)
dimnames(tmp)[[2]] <- c("family","ind","father","mother","sex","generation","affection","risk","genetic_risk",paste("locus_risk",1:length(genotype[1,]),sep="_"),paste("genotype",1:length(genotype[1,]),sep="_"))
tmp
}
tbl <- make.multi.family.table(F)
library(MASS)
write.matrix(tbl, file = "out.txt", sep = "\t")