else {
System.out.println( " failed." );
}
- // timeNeighborJoining();
+ // timeNeighborJoining();
}
public static boolean test( final File test_dir ) {
final NeighborJoiningR nj2 = NeighborJoiningR.createInstance( true, 6 );
//nj = NeighborJoining.createInstance( true, 6 );
final Phylogeny p2 = nj2.execute( m );
- Archaeopteryx.createApplication( p2 );
+ Archaeopteryx.createApplication( p2 );
p2.reRoot( p2.getNode( "Bovine" ) );
if ( isUnequal( p2.getNode( "Chimp" ).getDistanceToParent(), 0.151675 ) ) {
return false;
if ( otu1 < i ) {
if ( otu2 > i ) {
_d_values[ m_otu1 ][ m_i ] = ( _d_values[ m_otu1 ][ m_i ] + _d_values[ m_i ][ m_otu2 ] - d ) / 2;
- System.out.print( DF.format( _d_values[ m_otu1 ][ m_i ] ) );
+ System.out.print( DF.format( _d_values[ m_otu1 ][ m_i ] ) );
}
else {
_d_values[ m_otu1 ][ m_i ] = ( _d_values[ m_otu1 ][ m_i ] + _d_values[ m_otu2 ][ m_i ] - d ) / 2;
- System.out.print( DF.format( _d_values[ m_otu1 ][ m_i ] ) );
-
+ System.out.print( DF.format( _d_values[ m_otu1 ][ m_i ] ) );
}
}
else {
if ( otu2 > i ) {
_d_values[ m_i ][ m_otu1 ] = ( _d_values[ m_i ][ m_otu1 ] + _d_values[ m_i ][ m_otu2 ] - d ) / 2;
- System.out.print( DF.format( _d_values[ m_i ][ m_otu1 ] ) );
-
+ System.out.print( DF.format( _d_values[ m_i ][ m_otu1 ] ) );
}
else {
_d_values[ m_i ][ m_otu1 ] = ( _d_values[ m_i ][ m_otu1 ] + _d_values[ m_otu2 ][ m_i ] - d ) / 2;
- System.out.print( DF.format( _d_values[ m_otu1 ][ m_i ] ) );
-
+ System.out.print( DF.format( _d_values[ m_otu1 ][ m_i ] ) );
}
}
System.out.print( " " );
final PhylogenyNode node = new PhylogenyNode();
//final double d = getDvalueUnmapped( otu1, _mappings[ otu2 ] );
final double d = _d_values[ otu1 ][ _mappings[ otu2 ] ];
- final double d1 = ( d / 2 ) + ( ( _r[ _rev_mappings[ otu1 ] ] - _r[ otu2 ] ) / ( 2 * ( _n - 2 ) ) );
+ final double d1 = ( d / 2 ) + ( ( _r[ _rev_mappings[ otu1 ] ] - _r[ otu2 ] ) / ( 2 * ( _n - 2 ) ) );
final double d2 = d - d1;
if ( _df == null ) {
_external_nodes[ otu1 ].setDistanceToParent( d1 );
if ( j == otu2 ) {
continue;
}
- if ( otu1 < _mappings[ j ] ) {
+ if ( otu1 < _mappings[ j ] ) {
updateDvalue( otu1, otu2, j, d );
}
}
private final void updateDvalue( final int otu1, final int otu2, final int j, final double d ) {
final double new_d = ( getDvalueUnmapped( otu1, _mappings[ j ] ) + getDvalue( j, otu2 ) - d ) / 2;
System.out.print( DF.format( new_d ) + " " );
- // System.out.println( "going to remove: " + getDvalueUnmapped( otu1, _mappings[ j ] ) + ", " + otu1 + ", "
- // + _mappings[ j ] );
-
- if ( otu1< _mappings[ j ] ) {
+ // System.out.println( "going to remove: " + getDvalueUnmapped( otu1, _mappings[ j ] ) + ", " + otu1 + ", "
+ // + _mappings[ j ] );
+ if ( otu1 < _mappings[ j ] ) {
_s.removePairing( getDvalueUnmapped( otu1, _mappings[ j ] ), otu1, _mappings[ j ] );
}
else {
- _s.removePairing( getDvalueUnmapped( otu1, _mappings[ j ] ), _mappings[ j ] , otu1 );
+ _s.removePairing( getDvalueUnmapped( otu1, _mappings[ j ] ), _mappings[ j ], otu1 );
}
-
-
- // System.out.println( "going to remove: " + getDvalue( j, otu2 ) + ", " +_mappings[ otu2 ] + ", "
- // + _mappings[ j ] );
-
- if ( _mappings[ otu2 ] < _mappings[ j ] ) {
- _s.removePairing( getDvalue( j, otu2 ), _mappings[ otu2 ] , _mappings[ j ] );
+ // System.out.println( "going to remove: " + getDvalue( j, otu2 ) + ", " +_mappings[ otu2 ] + ", "
+ // + _mappings[ j ] );
+ if ( _mappings[ otu2 ] < _mappings[ j ] ) {
+ _s.removePairing( getDvalue( j, otu2 ), _mappings[ otu2 ], _mappings[ j ] );
}
else {
- _s.removePairing( getDvalue( j, otu2 ) , _mappings[ j ], _mappings[ otu2 ] );
-
+ _s.removePairing( getDvalue( j, otu2 ), _mappings[ j ], _mappings[ otu2 ] );
}
-
_s.addPairing( new_d, otu1, _mappings[ j ] );
setDvalueU( otu1, j, new_d );
}
if ( i < _mappings[ j ] ) {
_d_values[ i ][ _mappings[ j ] ] = d;
}
- _d_values[_mappings[ j] ][ i ] = d;
+ _d_values[ _mappings[ j ] ][ i ] = d;
}
private double getDvalue( final int i, final int j ) {