@@ -11,55 +11,55 @@ subfolders inside the project folder: ms, phylonet, astral, NJst
11
11
12
12
####SIMULATIONS:
13
13
14
- 1.run
15
- ```
14
+ 1 . run
15
+ ```
16
16
perl hgt_all.pl gamma=0.3 nloci=10 from=1 to=10
17
- ```
18
- This perl script
19
- calls hgt.pl which simulates gene trees with ms, and analyzes them
20
- with astral and phylonet. Careful: no spaces between the variable, =
21
- and the number. That is, gamma=0.3 is correct, but gamma = 0.3 is
22
- wrong.
17
+ ```
18
+ This perl script
19
+ calls hgt.pl which simulates gene trees with ms, and analyzes them
20
+ with astral and phylonet. Careful: no spaces between the variable, =
21
+ and the number. That is, gamma=0.3 is correct, but gamma = 0.3 is
22
+ wrong.
23
23
24
24
2 . run
25
- ```
25
+ ```
26
26
R CMD BATCH getSpeTree.R
27
- ```
28
- in the project folder which will
29
- create subfolders for each scenario in NJst directory and one file
30
- (per folder) with the list of all estimated trees(in order)in NJst:
31
- ` estTrees.out `
27
+ ```
28
+ in the project folder which will
29
+ create subfolders for each scenario in NJst directory and one file
30
+ (per folder) with the list of all estimated trees(in order)in NJst:
31
+ ` estTrees.out `
32
32
33
33
34
34
####SUMMARIZE OUTPUT:
35
35
36
36
1 . run
37
- ```
38
- perl summarizeOutputAstralAll.pl
39
- perl summarizeOutputPhylonetAll.pl
40
- ```
41
- which will create one file (per
37
+ ```
38
+ perl summarizeOutputAstralAll.pl
39
+ perl summarizeOutputPhylonetAll.pl
40
+ ```
41
+ which will create one file (per
42
42
folder) with the list of all estimated trees in astral:
43
43
` estTrees.out ` and all the estimated networks in phylonet:
44
44
` estNetworks.out ` for that folder (given value of gamma, nloci)
45
45
46
46
2 . run
47
- ```
48
- perl PhylonetCmpnetsAll.pl
49
- perl PhylonetSymmDistAll.pl
50
- ```
51
- which
47
+ ```
48
+ perl PhylonetCmpnetsAll.pl
49
+ perl PhylonetSymmDistAll.pl
50
+ ```
51
+ which
52
52
will calculate the distance of estimated networks/trees and
53
53
true network/tree (need to match the one used for simulations
54
54
in ms: ` hgt.pl ` ) in every subfolder in phylonet/astral: this will
55
55
create a file ` distances.out ` inside each subfolder
56
56
57
57
3 . run
58
- ```
59
- R CMD BATCH summarizeAstralDistances.r
60
- R CMD BATCH summarizeNJstDistances.r
61
- ```
62
- (` distances.R ` is script
58
+ ```
59
+ R CMD BATCH summarizeAstralDistances.r
60
+ R CMD BATCH summarizeNJstDistances.r
61
+ ```
62
+ (` distances.R ` is script
63
63
for only one specific case) that will get into each
64
64
subfolder in astral/NJst and will summarize results
65
65
based on distances.out into two files:
@@ -68,11 +68,11 @@ R CMD BATCH summarizeNJstDistances.r
68
68
mean/sd distance)
69
69
70
70
4 . run
71
- ```
72
- R CMD BATCH summarizeAstralAll.R
73
- R CMD BATCH summarizeNJstAll.R
74
- ```
75
- that will
71
+ ```
72
+ R CMD BATCH summarizeAstralAll.R
73
+ R CMD BATCH summarizeNJstAll.R
74
+ ```
75
+ that will
76
76
summarize all summary tables and will summarize it
77
77
into two files:
78
78
` 0.1_alltableastral.txt ` ,` 0.1_alltableNJst.txt ` (that
0 commit comments