fneighbor Error distance file...

Bonjour :coucou:

Ne trouves-tu pas que tes distances sont vachement grandes pour des distances phylogénétiques ?

Quand je divise tes valeurs par 1000, il n'y a pas d'erreur.

Code:

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
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
 
 
Output file   .stdout
 
 
    Cycle   9: species 1 (   0.51140) joins species 8 (   0.33560)
    Cycle   8: node 1 (   0.01628) joins species 9 (   0.43722)
    Cycle   7: species 2 (   0.43466) joins species 7 (   0.30034)
    Cycle   6: node 2 (   0.01986) joins species 4 (   0.50664)
    Cycle   5: species 3 (   0.52823) joins species 11 (   0.33077)
    Cycle   4: node 3 (   0.02159) joins species 5 (   0.45291)
    Cycle   3: node 1 (   0.04444) joins species 12 (   0.45056)
    Cycle   2: species 6 (  -0.21344) joins species 10 (   0.53044)
    Cycle   1: node 1 (   0.04434) joins node 2 (   0.07528)
    last cycle:
     node 1  (   0.00425) joins node 3  (   0.04691) joins node 6  (   0.02209)
 
    Output written on file "outfile"
 
    Tree written on file "outtreefile"
 
    Done.
 
Output file   outfile
 
    Neighbor-Joining/UPGMA method version 3.69.650
 
 
      12 Populations
 
     Neighbor-joining method
 
     Negative branch lengths allowed
 
 
      +-------------------sul       
      ! 
      ! +-------------------------hho       
      ! ! 
      ! !  +--------------------------aae       
      1-2  ! 
      ! !  !        +-------------------------hys       
      ! !  !      +-3 
      ! +--7 +----4 +-----------------pmx       
      !    ! !    ! 
      !    ! !    +-----------------------------hte       
      !    ! ! 
      !    +-9      +------------------------------tal       
      !      !    +-5 
      !      !  +-6 +-------------------saf       
      !      !  ! ! 
      !      +-10 +---------------------------hya       
      !         !  
      !         ! +ttk       
      !         +-8 
      !           +-------------------------------trd       
      ! 
      +------------------------------hth       
 
 
    remember: this is an unrooted tree!
 
    Between        And            Length
    -------        ---            ------
       1          sul             0.33560
       1             2            0.01628
       2          hho             0.43722
       2             7            0.04444
       7          aae             0.45056
       7             9            0.04434
       9             4            0.07528
       4             3            0.01986
       3          hys             0.43466
       3          pmx             0.30034
       4          hte             0.50664
       9            10            0.00425
      10             6            0.04691
       6             5            0.02159
       5          tal             0.52823
       5          saf             0.33077
       6          hya             0.45291
      10             8            0.02209
       8          ttk            -0.21344
       8          trd             0.53044
       1          hth             0.51140
 
 
Output file   outtreefile
 
    (sul:0.33560,(hho:0.43722,(aae:0.45056,(((hys:0.43466,pmx:0.30034):0.01986,
    hte:0.50664):0.07528,(((tal:0.52823,saf:0.33077):0.02159,
    hya:0.45291):0.04691,(ttk:-0.21344,trd:0.53044):0.02209):0.00425):0.04434):0.04444):0.01628,hth:0.51140);

Pour transformer, on peut utiliser awk :

Code:

awk '(NR==1){print;next;} {printf("%s\t",$1); for (i=2;i<14;i++) printf("%01.3f  ",$i/1000);print "";}' fichier_casse.dat >fichier_bon.dat

[edit]
Après vérification, c'est le 0 qui ne lui plaît pas.
Il veut toujours des nombres à virgule.
Pas besoin de diviser par 1000.
[/edit]