1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
| #NEXUS
begin trees;
tree tree_1 = [&R] ((((((((((((((((((((Leratiomyces_erythrocephalus:0.063048,Hypholoma_sublateritium:0.061545):0.018118,(Pholiota_conissans:0.042413,Pholiota_molesta:0.044289):0.022555):0.026812,Flammula_alnicola:0.047857):0.018264,Hebeloma_cylindrosporum_:0.084451):0.008944,((Agrocybe_pediades:0.091516,Psilocybe_cubensis:0.082387):0.01529,(Galerina_marginata_:0.053369,Gymnopilus_junonius:0.089199):0.013718):0.013594):0.018388,(Panaeolus_papilionaceus:0.138813,Crepidotus_variabilis:0.159754):0.022418):0.011953,(Cortinarius_glaucopus:0.088188,Crassisporium_funariophilum:0.064286):0.022517):0.04491,(((Coprinellus_micaceus:0.095202,Coprinellus_angulatus:0.072316):0.06939,(Coprinopsis_cinerea:0.085357,Coprinopsis_marcescibilis:0.109342):0.049711):0.084357,(Laccaria_amethystina_:0.017276,Laccaria_bicolor:0.019346):0.089572):0.012834):0.015456,(Crucibulum_laeve:0.084316,Cyathus_striatus:0.136396):0.024199):0.010255,((Macrolepiota_fuliginosa:0.061169,Leucoagaricus_leucothites:0.093137):0.01342,Agaricus_bisporus:0.107783):0.128488):0.016501,(((Amanita_rubescens:0.103903,Amanita_muscaria:0.110267):0.040762,Amanita_thiersii:0.103503):0.086768,(Volvariella_volvacea:0.148124,Pluteus_cervinus:0.14107):0.050209):0.01672):0.011013,((((Hypsizygus_marmoreus:0.068769,Lyophyllum_atratum:0.080005):0.024694,Lepista_nuda:0.087625):0.00857,(wgs23947:0.058519,Tricholoma_matsutake:0.049872):0.061731):0.026186,Infundibulicybe_gibba:0.128545):0.013097):0.01774,((((((((((Lentinula_edodes:0.070125,Rhodocollybia_butyracea:0.073752):0.007909,Gymnopus_luxurians_:0.074138):0.01007,Gymnopus_androsaceus:0.072701):0.073823,Omphalotus_olearius:0.089145):0.02597,Dendrothele_bispora:0.125909):0.015077,((Moniliophthora_perniciosa:0.066363,Moniliophthora_roreri:0.023024):0.059978,Marasmius_fiardii:0.104452)ArchaemarasmiuslegettiiFossil:0.063687):0.036952,(Megacollybia_platyphylla:0.007001,Megacollybia_marginata:0.006065):0.211567):0.010313,((((Armillaria_cepistipes:0.005463,Armillaria_gallica:0.005259):0.005418,Armillaria_solidipes:0.005362):0.015606,Guyanagaster_necrorhizus:0.029474):0.078045,((Hymenopellis_radicata:0.018069,Mucidula_mucida:0.018324):0.072702,Cylindrobasidium_torrendii:0.224841):0.059234):0.065523):0.0197,(((Schizophyllum_radiatum:0.003772,Schizophyllum_commune:0.003633):0.051017,Auriculariopsis_ampla:0.062294):0.164559,Fistulina_hepatica:0.245262):0.070772):0.010382,(Panellus_stipticus:0.079049,Mycena_galopus:0.090902):0.107204):0.010147):0.013642,((Pleurotus_eryngii:0.047323,Pleurotus_ostreatus:0.003405):0.166683,Pterula_gracilis:0.315452):0.029131):0.009785,Sarcomyxa_edulis:0.121765):0.007389,((((((ding384:0.013763,Qin538:0.023404):0.008152,zhao2118:0.050305):0.019386,zhu734:0.038099):0.015096,GT1067:0.069567):0.017235,((WGS1656:0.002411,ZK306:0.003593):0.053003,zhu398:0.032346):0.028474):0.163016,((TLGE34509:0.004295,TLGE20417:0.002647):9.53E-4,TLGE59471:0.001861):0.111779):0.026855):0.014635,(Hygrophorus_russula:0.043658,Hygrophorus_pudorinus:0.034936):0.09572):0.012172,Typhula_sp:0.22921):0.012522,Clavaria_fumosa:0.272602)Palaeoagaricitesantiquusfossil:0.00796,(((Clavulina_sp:0.518244,Ramaria_rubella:0.247376):0.100545,Cristinia_sonorae:0.208585):0.031971,(Xerocomus_badius:0.038878,Boletus_edulis:0.051155):0.188232):0.00796);
end;
begin r8s;
blformat lengths=persite nsites=498262 ulrametric=no;
constrain taxon=Palaeoagaricitesantiquusfossil minage=105 maxage=180;
constrain taxon=ArchaemarasmiuslegettiiFossil minage=92 maxage=210;
divtime method=PL crossv=yes cvstart=0 cvinc=1 cvnum=18；
#set smoothing=100;
#divtime method=PL algorithm=TN;
#showage;
#describe plot=phylogram;
#describe plot=chronogram;
#describe plot=tree_description;
#describe plot=chrono_description;
#describe plot=phylo_description;
#describe plot=rato_description;
end;
nsites: 只碱基数量
taxon是标注的化石位置，minage和maxage分别是估算的最大或最小时间
若化石的时间不是范围而是一个点，则写为：
fixage taxon=Asterids age=128;
下面这行是用来估算smoothing的，估算完了以后就可以删掉了
divtime method=PL crossv=yes cvstart=0 cvinc=1 cvnum=18;
set smoothing=100; 这一行是刚刚估算的最优值，其实应该都差不多
然后把下边的都选上就行
|
