1#NEXUS 2 3BEGIN TAXA; 4 DIMENSIONS NTAX = 13; 5 TAXLABELS 6 'U68496' 'U68497' 'U68498' 'U68499' 'U68500' 'U68501' 'U68502' 'U68503' 'U68504' 'U68505' 'U68506' 'U68507' 'U68508' ; 7END; 8 9BEGIN CHARACTERS; 10 DIMENSIONS NCHAR = 273; 11 FORMAT 12 DATATYPE = DNA 13 GAP=- 14 MISSING=? 15 ; 16 17MATRIX 18 'U68496' GTAGTAATTAGATCTGAAAACTTCTCGAACAATGCTAAAACCATAATAGTACAGCTAAATAAATCTGTAGAAATTAATTGTACAAGACCCAACAACAATACAAGAAGAAGTATACATTTTGGACCAGGGAAAGCATTTTATGCAGGAGAAATAATAGGAGATATAAGACAAGCATATTGTACCCTTAATGGAACAGAATGGAATAACACTTTAAAACAGGTAGCTGAAAAATTAAGAGAACAATTTATTAAAACAATAGTTTTTAATCAATCC 19 'U68497' ATAGTAATTAGATCTGAAAACTTCTCGAACAATGCTAAAACCATAATAGTACAGCTAAATAAATCTGTAGAAATTAATTGTACAAGACCCAACAACAATACAAGAAGAAGTATACATTTCGGACCAGGGAAAGCATTTTATGCAGGAGAAATAATAGGAGATATAAGACAAGCATATTGTACTCTTAATGGAGCAGAATGGAATAACACTGTAAAACAGGTAGCTGCAAAATTAAGAGAACAATTTAATAAAACAATAATCTTTAATCAATCC 20 'U68498' GTAGTAATTAGATCTGAAAACTTCACGAACAATGCTAAAACCATAATAGTACAGCTAAATAAATCTGTAGAAATTAATTGTGTAAGACCCGGCAACAATACAAGAAGAAGTATACATATAGGACCAGGGAGAGCATATTATACAGGAGAAGTAATAGGAGATATAAGACAAGCACATTGTAACCTTAGTAGAACAGACTGGAATAAAACTTTAAAACAGGTAGCTGAAAAATTAAGAGAACAATTTAATACAACAATAGTCTTTAATCAATCC 21 'U68499' ATAGTAATTAGATCTGAAAACTTCACGAACAATGCTAAAACCATAATAGTACAGCTAAATAAATCTGTAGAAATTAAGTGTGAAAGACCCAACAACAATACAAGAAAAAGTGTACATATAGGACCAGGGAGAGCATATTATACAGGAGAAATAATAGGAGATATAAGACAAGCACATTGTAACCTTAGTGGAACAGAATGGAGGGAAACTTTAAAACAGGTAGCTGAAAAATTAAGAGAACAATTTAATAAAACAATAGTCTTTAATCAATCC 22 'U68500' ATAGTAATTAGATCTGAAAACTTCACGAACAATGCTAAAACCATAATAGTACATCTAAATGAATCTGTAGAAATTATTTGTGAAAGACCCAACAACAATACAAGAAAAAGTGTACATATGGGACCAGGGAGAGCATATTACACAGGAGAAATAATAGGAGATATAAGACAAGCACATTGTAACATTAGTAGAACAAATTGGACGGAAACTTTAAAACAGGTAGCTGAAAAATTAAGAGAACAATTTAATAAAACAATAGTCTTTAATCAATCC 23 'U68501' GTAGTAATTAGATCTGAAAACTTCACGAACAATGCTAAGACCATAATAGTACAGCTAAATAAACCTGTAAAAATTAATTGTACAAGACCCAACAACAATGCAAAAATAAGAATACATATAGGACCAGGGAGACCATTTTATACAGCAGGAGAAATAGGAAATATAAGACAAGCACATTGTAACCTTAGTAGAACAGACTGGAATAACACTTTAAAACTGGTAGCTGAAAAATTAAGAGAACAATTTAATAAAACAATAGTCTTTAATCAATCC 24 'U68502' GTAGTAATTAGATCTGAAAACTTCACGAACAATGCTAAGACCATAATAGTACAGCTAAATAACTCTGTAGCAATTAAGTGTGAAAGACCCAACAACAATACAAGAAAAAGTATACCTATAGGACCAGGGAGAGCCTTTTATACAACAGGAGACATAGGAGATATAAGACAAGCACATTGTAACCTTAGTAGAAAAGACTGGAATGACACTTTAAGACAGGTAGTTGGAAAGTTAAGAGAACAATTTGGAAGAACAATAATCTTTAATCAATCC 25 'U68503' ATAGTAATTAGATCTGAAAACTTCACGAACAATGCTAAAACCATAATAGTACAGCTAAAGGAACCTGTAGACATTACTTGTGAAAGACCCAGCAACAATACAAGAAAAAGTATACATATAGGACCAGGAAAAGCATTTTATGCAACAGGAGAAATAGGAGATATAAGACGAGCACATTGTAACCTTAATAGAACAGCATGGAATAAAACTTTAAAACAGGTAGTTGAAAAATTAAGAGAACAATTTAAGAAAACAATAACCTTTAACCAATCC 26 'U68504' ATAGTAATTAGATCTGAAAACTTCACAGACAATGCTAAAACCATAATAGTACAGCTAAATAAATCTGTAGAAATTAATTGTACAAGACCCAACAACAATACAAGAAAAAGTATACATATAGGACCAGGGAGAGCTTTTTATACAACAGGTGAAATAGGAGATTTAAGACAAGCACATTGTAACCTTAGTAGAACAGCATGGAATGAAACTTTAAGACAAGTAGCTCAAAAATTAAAAGAACAATTTAATAGAACAATAGTTTTTAATCAATCC 27 'U68505' ATAGTAATTAGATCTGAAAACTTCACAGACAATGCTAAAACCATAATAGTACAGCTAAACAAATCTGTAGAAATTAATTGTATAAGACCCAACAACAATACAAGAAAAAGTATACATATAGGACCAGGGAGAGCATTTTATACAACAGGAGACATAGGAGATATAAGACAAGCATATTGTAACCTTAGTAGAACAGCATGGAATGAAACTTTAAGACAAGTAGCTCAAAAATTAAAAGAACAATTTAATAGAACAATAGTCTTTAATCAATCC 28 'U68506' ATAGTAATTAGATCTGAAAACTTCACGGACAATGCTAAAACCATAATAGTACAGCTAAATAAATCTGTAGAAATTAATTGTACAAGACCCAACAACAATACAAGAAAAAGTATAAATTTTGGACCAGGGAGAGCATTTTATACAACAGGTGAAATAGGAAATTTAAGACAAGCACATTGTAACCTTAGTAGAACAGCATGGAATGAAACTTTAAGACAAGTAGCTAAAAAACTAAAAGAACAATTTAATAGAACAATAGTTTTTAATCAATCC 29 'U68507' GTAGTAATTCGATCTGAAAACTTCACGGACAATGCTAAAACCATAATAGTACAGCTAAACAAATCTGTAGAAATTACTTGTGTAAGACCCAACAACAATACAAGAAAAAGTATAAATATAAGACCAGGGAGAGCATTTTATACAACAGGAGAAATAGGAGATATAAGACAAGCACATTGTAACCTTAGTAGAACAGCATGGAATGAAGCTTTAAGACAAGTAGCTAAAAAATTAAAAGAACAATTTAATAGAACAATAGTCTTTAATCAATCC 30 'U68508' ATAGTAATCAGATCTGAAAACTTCTCGGACAATGCTAAAACCATAATAGTACAGCTAAACAACACTGTAAACATTACTTGTGAAAGACCCAACAACAATACAAGAAAAAGGATACATATAGGACCAGGGAGAGCAGTTTATACAACAGGACAAATAGGAGATATAAGAAAAGCACATTGTAACCTTAGTAGAACAAATTGGACTGAAACTTTAAGACAAGTAGCTGAAAAATTAAAAGAACAATTTAATAAAACAATAATCTTTAATAATTCC; 31END; 32 33BEGIN TREES; 34 TREE tree = (U68496: 0.023746, U68497: 0.079178, ((U68498: 0.050407, (U68499: 0.028429, U68500: 0.099890): 0.083460): 0.023647, (U68501: 0.131681, (U68502: 0.191150, (U68503: 0.193336, ((((U68504: 0.018984, U68506: 0.070370): 0.055603, U68505: 0.030301): 0.036688, U68507: 0.061120): 0.034773, U68508: 0.198346): 0.050591): 0.064013): 0.033695): 0.048101): 0.152813); 35END; 36 37BEGIN HYPHY; 38 39 40global P2=0.7104121522639678; 41P2:<1; 42global W_1=0.05768869986246565; 43W_1:<1; 44global W_2=3.641844577696067; 45W_2:>1; 46global P1=0.3739141342952348; 47P1:<1; 48global kappa=0.3597183898000834; 49 50c.weights={1,3}; 51c.weights[0][0]:=P1; 52c.weights[0][1]:=(1-P1)*P2; 53c.weights[0][2]:=(1-P1)*(1-P2); 54 55 56c.points={1,3}; 57c.points[0][0]:=W_1; 58c.points[0][1]:=1; 59c.points[0][2]:=W_2; 60 61 62category c=(3,c.weights,MEAN,,c.points,0,1,); 63 64modelMatrix={61,61}; 65modelMatrix[0][1]:=kappa*c*t; 66modelMatrix[0][2]:=t; 67modelMatrix[0][3]:=kappa*c*t; 68modelMatrix[0][4]:=kappa*c*t; 69modelMatrix[0][8]:=c*t; 70modelMatrix[0][12]:=kappa*c*t; 71modelMatrix[0][16]:=kappa*c*t; 72modelMatrix[0][32]:=c*t; 73modelMatrix[1][0]:=kappa*c*t; 74modelMatrix[1][2]:=kappa*c*t; 75modelMatrix[1][3]:=t; 76modelMatrix[1][5]:=kappa*c*t; 77modelMatrix[1][9]:=c*t; 78modelMatrix[1][13]:=kappa*c*t; 79modelMatrix[1][17]:=kappa*c*t; 80modelMatrix[1][33]:=c*t; 81modelMatrix[1][48]:=kappa*c*t; 82modelMatrix[2][0]:=t; 83modelMatrix[2][1]:=kappa*c*t; 84modelMatrix[2][3]:=kappa*c*t; 85modelMatrix[2][6]:=kappa*c*t; 86modelMatrix[2][10]:=c*t; 87modelMatrix[2][14]:=kappa*c*t; 88modelMatrix[2][18]:=kappa*c*t; 89modelMatrix[2][34]:=c*t; 90modelMatrix[3][0]:=kappa*c*t; 91modelMatrix[3][1]:=t; 92modelMatrix[3][2]:=kappa*c*t; 93modelMatrix[3][7]:=kappa*c*t; 94modelMatrix[3][11]:=c*t; 95modelMatrix[3][15]:=kappa*c*t; 96modelMatrix[3][19]:=kappa*c*t; 97modelMatrix[3][35]:=c*t; 98modelMatrix[3][49]:=kappa*c*t; 99modelMatrix[4][0]:=kappa*c*t; 100modelMatrix[4][5]:=kappa*t; 101modelMatrix[4][6]:=t; 102modelMatrix[4][7]:=kappa*t; 103modelMatrix[4][8]:=kappa*c*t; 104modelMatrix[4][12]:=c*t; 105modelMatrix[4][20]:=kappa*c*t; 106modelMatrix[4][36]:=c*t; 107modelMatrix[4][50]:=kappa*c*t; 108modelMatrix[5][1]:=kappa*c*t; 109modelMatrix[5][4]:=kappa*t; 110modelMatrix[5][6]:=kappa*t; 111modelMatrix[5][7]:=t; 112modelMatrix[5][9]:=kappa*c*t; 113modelMatrix[5][13]:=c*t; 114modelMatrix[5][21]:=kappa*c*t; 115modelMatrix[5][37]:=c*t; 116modelMatrix[5][51]:=kappa*c*t; 117modelMatrix[6][2]:=kappa*c*t; 118modelMatrix[6][4]:=t; 119modelMatrix[6][5]:=kappa*t; 120modelMatrix[6][7]:=kappa*t; 121modelMatrix[6][10]:=kappa*c*t; 122modelMatrix[6][14]:=c*t; 123modelMatrix[6][22]:=kappa*c*t; 124modelMatrix[6][38]:=c*t; 125modelMatrix[6][52]:=kappa*c*t; 126modelMatrix[7][3]:=kappa*c*t; 127modelMatrix[7][4]:=kappa*t; 128modelMatrix[7][5]:=t; 129modelMatrix[7][6]:=kappa*t; 130modelMatrix[7][11]:=kappa*c*t; 131modelMatrix[7][15]:=c*t; 132modelMatrix[7][23]:=kappa*c*t; 133modelMatrix[7][39]:=c*t; 134modelMatrix[7][53]:=kappa*c*t; 135modelMatrix[8][0]:=c*t; 136modelMatrix[8][4]:=kappa*c*t; 137modelMatrix[8][9]:=kappa*c*t; 138modelMatrix[8][10]:=t; 139modelMatrix[8][11]:=kappa*c*t; 140modelMatrix[8][12]:=kappa*c*t; 141modelMatrix[8][24]:=kappa*t; 142modelMatrix[8][40]:=c*t; 143modelMatrix[9][1]:=c*t; 144modelMatrix[9][5]:=kappa*c*t; 145modelMatrix[9][8]:=kappa*c*t; 146modelMatrix[9][10]:=kappa*c*t; 147modelMatrix[9][11]:=t; 148modelMatrix[9][13]:=kappa*c*t; 149modelMatrix[9][25]:=kappa*c*t; 150modelMatrix[9][41]:=c*t; 151modelMatrix[9][54]:=kappa*c*t; 152modelMatrix[10][2]:=c*t; 153modelMatrix[10][6]:=kappa*c*t; 154modelMatrix[10][8]:=t; 155modelMatrix[10][9]:=kappa*c*t; 156modelMatrix[10][11]:=kappa*c*t; 157modelMatrix[10][14]:=kappa*c*t; 158modelMatrix[10][26]:=kappa*t; 159modelMatrix[10][42]:=c*t; 160modelMatrix[10][55]:=kappa*c*t; 161modelMatrix[11][3]:=c*t; 162modelMatrix[11][7]:=kappa*c*t; 163modelMatrix[11][8]:=kappa*c*t; 164modelMatrix[11][9]:=t; 165modelMatrix[11][10]:=kappa*c*t; 166modelMatrix[11][15]:=kappa*c*t; 167modelMatrix[11][27]:=kappa*c*t; 168modelMatrix[11][43]:=c*t; 169modelMatrix[11][56]:=kappa*c*t; 170modelMatrix[12][0]:=kappa*c*t; 171modelMatrix[12][4]:=c*t; 172modelMatrix[12][8]:=kappa*c*t; 173modelMatrix[12][13]:=kappa*t; 174modelMatrix[12][14]:=c*t; 175modelMatrix[12][15]:=kappa*t; 176modelMatrix[12][28]:=kappa*c*t; 177modelMatrix[12][44]:=c*t; 178modelMatrix[12][57]:=kappa*c*t; 179modelMatrix[13][1]:=kappa*c*t; 180modelMatrix[13][5]:=c*t; 181modelMatrix[13][9]:=kappa*c*t; 182modelMatrix[13][12]:=kappa*t; 183modelMatrix[13][14]:=kappa*c*t; 184modelMatrix[13][15]:=t; 185modelMatrix[13][29]:=kappa*c*t; 186modelMatrix[13][45]:=c*t; 187modelMatrix[13][58]:=kappa*c*t; 188modelMatrix[14][2]:=kappa*c*t; 189modelMatrix[14][6]:=c*t; 190modelMatrix[14][10]:=kappa*c*t; 191modelMatrix[14][12]:=c*t; 192modelMatrix[14][13]:=kappa*c*t; 193modelMatrix[14][15]:=kappa*c*t; 194modelMatrix[14][30]:=kappa*c*t; 195modelMatrix[14][46]:=c*t; 196modelMatrix[14][59]:=kappa*c*t; 197modelMatrix[15][3]:=kappa*c*t; 198modelMatrix[15][7]:=c*t; 199modelMatrix[15][11]:=kappa*c*t; 200modelMatrix[15][12]:=kappa*t; 201modelMatrix[15][13]:=t; 202modelMatrix[15][14]:=kappa*c*t; 203modelMatrix[15][31]:=kappa*c*t; 204modelMatrix[15][47]:=c*t; 205modelMatrix[15][60]:=kappa*c*t; 206modelMatrix[16][0]:=kappa*c*t; 207modelMatrix[16][17]:=kappa*c*t; 208modelMatrix[16][18]:=t; 209modelMatrix[16][19]:=kappa*c*t; 210modelMatrix[16][20]:=kappa*c*t; 211modelMatrix[16][24]:=c*t; 212modelMatrix[16][28]:=kappa*c*t; 213modelMatrix[16][32]:=kappa*c*t; 214modelMatrix[17][1]:=kappa*c*t; 215modelMatrix[17][16]:=kappa*c*t; 216modelMatrix[17][18]:=kappa*c*t; 217modelMatrix[17][19]:=t; 218modelMatrix[17][21]:=kappa*c*t; 219modelMatrix[17][25]:=c*t; 220modelMatrix[17][29]:=kappa*c*t; 221modelMatrix[17][33]:=kappa*c*t; 222modelMatrix[17][48]:=c*t; 223modelMatrix[18][2]:=kappa*c*t; 224modelMatrix[18][16]:=t; 225modelMatrix[18][17]:=kappa*c*t; 226modelMatrix[18][19]:=kappa*c*t; 227modelMatrix[18][22]:=kappa*c*t; 228modelMatrix[18][26]:=c*t; 229modelMatrix[18][30]:=kappa*c*t; 230modelMatrix[18][34]:=kappa*c*t; 231modelMatrix[19][3]:=kappa*c*t; 232modelMatrix[19][16]:=kappa*c*t; 233modelMatrix[19][17]:=t; 234modelMatrix[19][18]:=kappa*c*t; 235modelMatrix[19][23]:=kappa*c*t; 236modelMatrix[19][27]:=c*t; 237modelMatrix[19][31]:=kappa*c*t; 238modelMatrix[19][35]:=kappa*c*t; 239modelMatrix[19][49]:=c*t; 240modelMatrix[20][4]:=kappa*c*t; 241modelMatrix[20][16]:=kappa*c*t; 242modelMatrix[20][21]:=kappa*t; 243modelMatrix[20][22]:=t; 244modelMatrix[20][23]:=kappa*t; 245modelMatrix[20][24]:=kappa*c*t; 246modelMatrix[20][28]:=c*t; 247modelMatrix[20][36]:=kappa*c*t; 248modelMatrix[20][50]:=c*t; 249modelMatrix[21][5]:=kappa*c*t; 250modelMatrix[21][17]:=kappa*c*t; 251modelMatrix[21][20]:=kappa*t; 252modelMatrix[21][22]:=kappa*t; 253modelMatrix[21][23]:=t; 254modelMatrix[21][25]:=kappa*c*t; 255modelMatrix[21][29]:=c*t; 256modelMatrix[21][37]:=kappa*c*t; 257modelMatrix[21][51]:=c*t; 258modelMatrix[22][6]:=kappa*c*t; 259modelMatrix[22][18]:=kappa*c*t; 260modelMatrix[22][20]:=t; 261modelMatrix[22][21]:=kappa*t; 262modelMatrix[22][23]:=kappa*t; 263modelMatrix[22][26]:=kappa*c*t; 264modelMatrix[22][30]:=c*t; 265modelMatrix[22][38]:=kappa*c*t; 266modelMatrix[22][52]:=c*t; 267modelMatrix[23][7]:=kappa*c*t; 268modelMatrix[23][19]:=kappa*c*t; 269modelMatrix[23][20]:=kappa*t; 270modelMatrix[23][21]:=t; 271modelMatrix[23][22]:=kappa*t; 272modelMatrix[23][27]:=kappa*c*t; 273modelMatrix[23][31]:=c*t; 274modelMatrix[23][39]:=kappa*c*t; 275modelMatrix[23][53]:=c*t; 276modelMatrix[24][8]:=kappa*t; 277modelMatrix[24][16]:=c*t; 278modelMatrix[24][20]:=kappa*c*t; 279modelMatrix[24][25]:=kappa*t; 280modelMatrix[24][26]:=t; 281modelMatrix[24][27]:=kappa*t; 282modelMatrix[24][28]:=kappa*c*t; 283modelMatrix[24][40]:=kappa*c*t; 284modelMatrix[25][9]:=kappa*c*t; 285modelMatrix[25][17]:=c*t; 286modelMatrix[25][21]:=kappa*c*t; 287modelMatrix[25][24]:=kappa*t; 288modelMatrix[25][26]:=kappa*t; 289modelMatrix[25][27]:=t; 290modelMatrix[25][29]:=kappa*c*t; 291modelMatrix[25][41]:=kappa*c*t; 292modelMatrix[25][54]:=c*t; 293modelMatrix[26][10]:=kappa*t; 294modelMatrix[26][18]:=c*t; 295modelMatrix[26][22]:=kappa*c*t; 296modelMatrix[26][24]:=t; 297modelMatrix[26][25]:=kappa*t; 298modelMatrix[26][27]:=kappa*t; 299modelMatrix[26][30]:=kappa*c*t; 300modelMatrix[26][42]:=kappa*c*t; 301modelMatrix[26][55]:=c*t; 302modelMatrix[27][11]:=kappa*c*t; 303modelMatrix[27][19]:=c*t; 304modelMatrix[27][23]:=kappa*c*t; 305modelMatrix[27][24]:=kappa*t; 306modelMatrix[27][25]:=t; 307modelMatrix[27][26]:=kappa*t; 308modelMatrix[27][31]:=kappa*c*t; 309modelMatrix[27][43]:=kappa*c*t; 310modelMatrix[27][56]:=c*t; 311modelMatrix[28][12]:=kappa*c*t; 312modelMatrix[28][16]:=kappa*c*t; 313modelMatrix[28][20]:=c*t; 314modelMatrix[28][24]:=kappa*c*t; 315modelMatrix[28][29]:=kappa*t; 316modelMatrix[28][30]:=t; 317modelMatrix[28][31]:=kappa*t; 318modelMatrix[28][44]:=kappa*c*t; 319modelMatrix[28][57]:=t; 320modelMatrix[29][13]:=kappa*c*t; 321modelMatrix[29][17]:=kappa*c*t; 322modelMatrix[29][21]:=c*t; 323modelMatrix[29][25]:=kappa*c*t; 324modelMatrix[29][28]:=kappa*t; 325modelMatrix[29][30]:=kappa*t; 326modelMatrix[29][31]:=t; 327modelMatrix[29][45]:=kappa*c*t; 328modelMatrix[29][58]:=c*t; 329modelMatrix[30][14]:=kappa*c*t; 330modelMatrix[30][18]:=kappa*c*t; 331modelMatrix[30][22]:=c*t; 332modelMatrix[30][26]:=kappa*c*t; 333modelMatrix[30][28]:=t; 334modelMatrix[30][29]:=kappa*t; 335modelMatrix[30][31]:=kappa*t; 336modelMatrix[30][46]:=kappa*c*t; 337modelMatrix[30][59]:=t; 338modelMatrix[31][15]:=kappa*c*t; 339modelMatrix[31][19]:=kappa*c*t; 340modelMatrix[31][23]:=c*t; 341modelMatrix[31][27]:=kappa*c*t; 342modelMatrix[31][28]:=kappa*t; 343modelMatrix[31][29]:=t; 344modelMatrix[31][30]:=kappa*t; 345modelMatrix[31][47]:=kappa*c*t; 346modelMatrix[31][60]:=c*t; 347modelMatrix[32][0]:=c*t; 348modelMatrix[32][16]:=kappa*c*t; 349modelMatrix[32][33]:=kappa*c*t; 350modelMatrix[32][34]:=t; 351modelMatrix[32][35]:=kappa*c*t; 352modelMatrix[32][36]:=kappa*c*t; 353modelMatrix[32][40]:=c*t; 354modelMatrix[32][44]:=kappa*c*t; 355modelMatrix[33][1]:=c*t; 356modelMatrix[33][17]:=kappa*c*t; 357modelMatrix[33][32]:=kappa*c*t; 358modelMatrix[33][34]:=kappa*c*t; 359modelMatrix[33][35]:=t; 360modelMatrix[33][37]:=kappa*c*t; 361modelMatrix[33][41]:=c*t; 362modelMatrix[33][45]:=kappa*c*t; 363modelMatrix[33][48]:=kappa*c*t; 364modelMatrix[34][2]:=c*t; 365modelMatrix[34][18]:=kappa*c*t; 366modelMatrix[34][32]:=t; 367modelMatrix[34][33]:=kappa*c*t; 368modelMatrix[34][35]:=kappa*c*t; 369modelMatrix[34][38]:=kappa*c*t; 370modelMatrix[34][42]:=c*t; 371modelMatrix[34][46]:=kappa*c*t; 372modelMatrix[35][3]:=c*t; 373modelMatrix[35][19]:=kappa*c*t; 374modelMatrix[35][32]:=kappa*c*t; 375modelMatrix[35][33]:=t; 376modelMatrix[35][34]:=kappa*c*t; 377modelMatrix[35][39]:=kappa*c*t; 378modelMatrix[35][43]:=c*t; 379modelMatrix[35][47]:=kappa*c*t; 380modelMatrix[35][49]:=kappa*c*t; 381modelMatrix[36][4]:=c*t; 382modelMatrix[36][20]:=kappa*c*t; 383modelMatrix[36][32]:=kappa*c*t; 384modelMatrix[36][37]:=kappa*t; 385modelMatrix[36][38]:=t; 386modelMatrix[36][39]:=kappa*t; 387modelMatrix[36][40]:=kappa*c*t; 388modelMatrix[36][44]:=c*t; 389modelMatrix[36][50]:=kappa*c*t; 390modelMatrix[37][5]:=c*t; 391modelMatrix[37][21]:=kappa*c*t; 392modelMatrix[37][33]:=kappa*c*t; 393modelMatrix[37][36]:=kappa*t; 394modelMatrix[37][38]:=kappa*t; 395modelMatrix[37][39]:=t; 396modelMatrix[37][41]:=kappa*c*t; 397modelMatrix[37][45]:=c*t; 398modelMatrix[37][51]:=kappa*c*t; 399modelMatrix[38][6]:=c*t; 400modelMatrix[38][22]:=kappa*c*t; 401modelMatrix[38][34]:=kappa*c*t; 402modelMatrix[38][36]:=t; 403modelMatrix[38][37]:=kappa*t; 404modelMatrix[38][39]:=kappa*t; 405modelMatrix[38][42]:=kappa*c*t; 406modelMatrix[38][46]:=c*t; 407modelMatrix[38][52]:=kappa*c*t; 408modelMatrix[39][7]:=c*t; 409modelMatrix[39][23]:=kappa*c*t; 410modelMatrix[39][35]:=kappa*c*t; 411modelMatrix[39][36]:=kappa*t; 412modelMatrix[39][37]:=t; 413modelMatrix[39][38]:=kappa*t; 414modelMatrix[39][43]:=kappa*c*t; 415modelMatrix[39][47]:=c*t; 416modelMatrix[39][53]:=kappa*c*t; 417modelMatrix[40][8]:=c*t; 418modelMatrix[40][24]:=kappa*c*t; 419modelMatrix[40][32]:=c*t; 420modelMatrix[40][36]:=kappa*c*t; 421modelMatrix[40][41]:=kappa*t; 422modelMatrix[40][42]:=t; 423modelMatrix[40][43]:=kappa*t; 424modelMatrix[40][44]:=kappa*c*t; 425modelMatrix[41][9]:=c*t; 426modelMatrix[41][25]:=kappa*c*t; 427modelMatrix[41][33]:=c*t; 428modelMatrix[41][37]:=kappa*c*t; 429modelMatrix[41][40]:=kappa*t; 430modelMatrix[41][42]:=kappa*t; 431modelMatrix[41][43]:=t; 432modelMatrix[41][45]:=kappa*c*t; 433modelMatrix[41][54]:=kappa*c*t; 434modelMatrix[42][10]:=c*t; 435modelMatrix[42][26]:=kappa*c*t; 436modelMatrix[42][34]:=c*t; 437modelMatrix[42][38]:=kappa*c*t; 438modelMatrix[42][40]:=t; 439modelMatrix[42][41]:=kappa*t; 440modelMatrix[42][43]:=kappa*t; 441modelMatrix[42][46]:=kappa*c*t; 442modelMatrix[42][55]:=kappa*c*t; 443modelMatrix[43][11]:=c*t; 444modelMatrix[43][27]:=kappa*c*t; 445modelMatrix[43][35]:=c*t; 446modelMatrix[43][39]:=kappa*c*t; 447modelMatrix[43][40]:=kappa*t; 448modelMatrix[43][41]:=t; 449modelMatrix[43][42]:=kappa*t; 450modelMatrix[43][47]:=kappa*c*t; 451modelMatrix[43][56]:=kappa*c*t; 452modelMatrix[44][12]:=c*t; 453modelMatrix[44][28]:=kappa*c*t; 454modelMatrix[44][32]:=kappa*c*t; 455modelMatrix[44][36]:=c*t; 456modelMatrix[44][40]:=kappa*c*t; 457modelMatrix[44][45]:=kappa*t; 458modelMatrix[44][46]:=t; 459modelMatrix[44][47]:=kappa*t; 460modelMatrix[44][57]:=kappa*c*t; 461modelMatrix[45][13]:=c*t; 462modelMatrix[45][29]:=kappa*c*t; 463modelMatrix[45][33]:=kappa*c*t; 464modelMatrix[45][37]:=c*t; 465modelMatrix[45][41]:=kappa*c*t; 466modelMatrix[45][44]:=kappa*t; 467modelMatrix[45][46]:=kappa*t; 468modelMatrix[45][47]:=t; 469modelMatrix[45][58]:=kappa*c*t; 470modelMatrix[46][14]:=c*t; 471modelMatrix[46][30]:=kappa*c*t; 472modelMatrix[46][34]:=kappa*c*t; 473modelMatrix[46][38]:=c*t; 474modelMatrix[46][42]:=kappa*c*t; 475modelMatrix[46][44]:=t; 476modelMatrix[46][45]:=kappa*t; 477modelMatrix[46][47]:=kappa*t; 478modelMatrix[46][59]:=kappa*c*t; 479modelMatrix[47][15]:=c*t; 480modelMatrix[47][31]:=kappa*c*t; 481modelMatrix[47][35]:=kappa*c*t; 482modelMatrix[47][39]:=c*t; 483modelMatrix[47][43]:=kappa*c*t; 484modelMatrix[47][44]:=kappa*t; 485modelMatrix[47][45]:=t; 486modelMatrix[47][46]:=kappa*t; 487modelMatrix[47][60]:=kappa*c*t; 488modelMatrix[48][1]:=kappa*c*t; 489modelMatrix[48][17]:=c*t; 490modelMatrix[48][33]:=kappa*c*t; 491modelMatrix[48][49]:=t; 492modelMatrix[48][51]:=kappa*c*t; 493modelMatrix[48][54]:=c*t; 494modelMatrix[48][58]:=kappa*c*t; 495modelMatrix[49][3]:=kappa*c*t; 496modelMatrix[49][19]:=c*t; 497modelMatrix[49][35]:=kappa*c*t; 498modelMatrix[49][48]:=t; 499modelMatrix[49][53]:=kappa*c*t; 500modelMatrix[49][56]:=c*t; 501modelMatrix[49][60]:=kappa*c*t; 502modelMatrix[50][4]:=kappa*c*t; 503modelMatrix[50][20]:=c*t; 504modelMatrix[50][36]:=kappa*c*t; 505modelMatrix[50][51]:=kappa*t; 506modelMatrix[50][52]:=t; 507modelMatrix[50][53]:=kappa*t; 508modelMatrix[50][57]:=c*t; 509modelMatrix[51][5]:=kappa*c*t; 510modelMatrix[51][21]:=c*t; 511modelMatrix[51][37]:=kappa*c*t; 512modelMatrix[51][48]:=kappa*c*t; 513modelMatrix[51][50]:=kappa*t; 514modelMatrix[51][52]:=kappa*t; 515modelMatrix[51][53]:=t; 516modelMatrix[51][54]:=kappa*c*t; 517modelMatrix[51][58]:=c*t; 518modelMatrix[52][6]:=kappa*c*t; 519modelMatrix[52][22]:=c*t; 520modelMatrix[52][38]:=kappa*c*t; 521modelMatrix[52][50]:=t; 522modelMatrix[52][51]:=kappa*t; 523modelMatrix[52][53]:=kappa*t; 524modelMatrix[52][55]:=kappa*c*t; 525modelMatrix[52][59]:=c*t; 526modelMatrix[53][7]:=kappa*c*t; 527modelMatrix[53][23]:=c*t; 528modelMatrix[53][39]:=kappa*c*t; 529modelMatrix[53][49]:=kappa*c*t; 530modelMatrix[53][50]:=kappa*t; 531modelMatrix[53][51]:=t; 532modelMatrix[53][52]:=kappa*t; 533modelMatrix[53][56]:=kappa*c*t; 534modelMatrix[53][60]:=c*t; 535modelMatrix[54][9]:=kappa*c*t; 536modelMatrix[54][25]:=c*t; 537modelMatrix[54][41]:=kappa*c*t; 538modelMatrix[54][48]:=c*t; 539modelMatrix[54][51]:=kappa*c*t; 540modelMatrix[54][55]:=kappa*c*t; 541modelMatrix[54][56]:=t; 542modelMatrix[54][58]:=kappa*c*t; 543modelMatrix[55][10]:=kappa*c*t; 544modelMatrix[55][26]:=c*t; 545modelMatrix[55][42]:=kappa*c*t; 546modelMatrix[55][52]:=kappa*c*t; 547modelMatrix[55][54]:=kappa*c*t; 548modelMatrix[55][56]:=kappa*c*t; 549modelMatrix[55][59]:=kappa*c*t; 550modelMatrix[56][11]:=kappa*c*t; 551modelMatrix[56][27]:=c*t; 552modelMatrix[56][43]:=kappa*c*t; 553modelMatrix[56][49]:=c*t; 554modelMatrix[56][53]:=kappa*c*t; 555modelMatrix[56][54]:=t; 556modelMatrix[56][55]:=kappa*c*t; 557modelMatrix[56][60]:=kappa*c*t; 558modelMatrix[57][12]:=kappa*c*t; 559modelMatrix[57][28]:=t; 560modelMatrix[57][44]:=kappa*c*t; 561modelMatrix[57][50]:=c*t; 562modelMatrix[57][58]:=kappa*c*t; 563modelMatrix[57][59]:=t; 564modelMatrix[57][60]:=kappa*c*t; 565modelMatrix[58][13]:=kappa*c*t; 566modelMatrix[58][29]:=c*t; 567modelMatrix[58][45]:=kappa*c*t; 568modelMatrix[58][48]:=kappa*c*t; 569modelMatrix[58][51]:=c*t; 570modelMatrix[58][54]:=kappa*c*t; 571modelMatrix[58][57]:=kappa*c*t; 572modelMatrix[58][59]:=kappa*c*t; 573modelMatrix[58][60]:=t; 574modelMatrix[59][14]:=kappa*c*t; 575modelMatrix[59][30]:=t; 576modelMatrix[59][46]:=kappa*c*t; 577modelMatrix[59][52]:=c*t; 578modelMatrix[59][55]:=kappa*c*t; 579modelMatrix[59][57]:=t; 580modelMatrix[59][58]:=kappa*c*t; 581modelMatrix[59][60]:=kappa*c*t; 582modelMatrix[60][15]:=kappa*c*t; 583modelMatrix[60][31]:=c*t; 584modelMatrix[60][47]:=kappa*c*t; 585modelMatrix[60][49]:=kappa*c*t; 586modelMatrix[60][53]:=c*t; 587modelMatrix[60][56]:=kappa*c*t; 588modelMatrix[60][57]:=kappa*c*t; 589modelMatrix[60][58]:=t; 590modelMatrix[60][59]:=kappa*c*t; 591 592 593vectorOfFrequencies={ 594{0.09810421316102973} 595{0.02278249499948084} 596{0.0103838582650695} 597{0.05207427428452762} 598{0.05450234064501652} 599{0.01265694166637824} 600{0.005768810147260831} 601{0.02893015238029313} 602{0.05009576842265348} 603{0.01163361446782} 604{0.005302395709822721} 605{0.02659111878358857} 606{0.0716647798268515} 607{0.01664253180813139} 608{0.007585371640440838} 609{0.03804007270430033} 610{0.02422953796550475} 611{0.005626764740804422} 612{0.002564579847849634} 613{0.01286117655041011} 614{0.01346085442528042} 615{0.003125980411558012} 616{0.001424766582138685} 617{0.00714509808356117} 618{0.0123725300249386} 619{0.00287324156977247} 620{0.001309572688263643} 621{0.006567409302337075} 622{0.01769959156345383} 623{0.004110331690091173} 624{0.001873416484599378} 625{0.009395043863065538} 626{0.04862851325943961} 627{0.01129287748678929} 628{0.005147093820509405} 629{0.02581229139837552} 630{0.02701584069968867} 631{0.006273820825994051} 632{0.002859496566949669} 633{0.0143401618879864} 634{0.02483158123886278} 635{0.005766575737934958} 636{0.002628303227494164} 637{0.01318074454385133} 638{0.03552295649448425} 639{0.00824940695843473} 640{0.003759933783776374} 641{0.0188557873335651} 642{0.006846552901398386} 643{0.01564926377462488} 644{0.01637894174824331} 645{0.003803640500776881} 646{0.001733632064979939} 647{0.008694035430347156} 648{0.003496112119863006} 649{0.00159346606823688} 650{0.007991113416829728} 651{0.02153656595832843} 652{0.005001382615915133} 653{0.002279541736505537} 654{0.0114317316935203} 655}; 656Model theModel=(modelMatrix,vectorOfFrequencies); 657TRY_NUMERIC_SEQUENCE_MATCH=0; 658ACCEPT_ROOTED_TREES=0; 659 660UseModel (theModel); 661Tree givenTree=(U68496,U68497,((U68498,(U68499,U68500)Node7)Node5,(U68501,(U68502,(U68503,((((U68504,U68506)Node19,U68505)Node18,U68507)Node17,U68508)Node16)Node14)Node12)Node10)Node4); 662 663givenTree.U68496.t=0.1820765235473329; 664givenTree.U68497.t=0.6321073455214453; 665givenTree.U68498.t=0.3343564501791256; 666givenTree.U68499.t=0.1876640197504523; 667givenTree.U68500.t=0.8005801879569905; 668givenTree.Node7.t=0.750108304945318; 669givenTree.Node5.t=0.189875002788432; 670givenTree.U68501.t=1.055571283712865; 671givenTree.U68502.t=1.747134065548263; 672givenTree.U68503.t=1.779098089271837; 673givenTree.U68504.t=0.1466122835681893; 674givenTree.U68506.t=0.5212546540008767; 675givenTree.Node19.t=0.4171959202637034; 676givenTree.U68505.t=0.2414131900046768; 677givenTree.Node18.t=0.2762126262689943; 678givenTree.U68507.t=0.4784362608407944; 679givenTree.Node17.t=0.006818214720109361; 680givenTree.U68508.t=1.945808887957573; 681givenTree.Node16.t=0.5086054961633685; 682givenTree.Node14.t=0.508023542650578; 683givenTree.Node12.t=0.1816851776949857; 684givenTree.Node10.t=0.394526379685196; 685givenTree.Node4.t=1.247106653589127;SetParameter (DEFER_CONSTRAINT_APPLICATION, 1, 0); 686SetParameter (DEFER_CONSTRAINT_APPLICATION, 0, 0); 687 688DataSet ds = ReadDataFile(USE_NEXUS_FILE_DATA); 689DataSetFilter filteredData = CreateFilter(ds,3,"0-272","0-8,10,9,11,12","TAA,TAG,TGA"); 690ASSUME_REVERSIBLE_MODELS=0; 691USE_LAST_RESULTS=1; 692LikelihoodFunction lf = (filteredData,givenTree); 693 694END;