my $self = shift;
  my $node = shift;
  my $p = shift;

  # If we've already burned enough nodes, return
  if ($self->{num_nodes} >= $self->{n}) {
    return;
  }

  $self->{vertices}{$node} = 1;

  my $dist = new Clair::Statistics::Distributions::Geometric;
  my $net = $self->{net};
  my $oldgraph = $self->{oldnet}->{graph};

  # First, add the node to the network
  $self->{num_nodes}++;

  my $num_out_links = $dist->get_random_value($p);
  my @successors = $oldgraph->successors($node);
  my @nodes = $self->get_random_nodes($num_out_links, subset =>\@ successors);
  foreach my $new_node (@nodes) {
    if (not defined $self->{vertices}{$new_node}) {
      $self->burn_node($new_node, $p);
    }
  }

  my $self = shift;
  my $n = shift;
  my $pf = shift;

  my $net = $self->{oldnet};
  $self->{net} = new Clair::Network(directed =>
                                    $self->get_old_net()->{directed});
  $self->{n} = $n;
  $self->{num_nodes} = 0;
  $self->{vertices} = ();

#  my $mean = $pf / (1 - $pf);
#  my $p = 1 / $mean;
  my $p = 1 - $pf;
#  print "p: $p\n";

  while ($self->{num_nodes} < $self->{n}) {
    my @nodes = $self->get_random_nodes();
    my $node = $nodes[0];
    if (not defined $self->{vertices}{$node}) {
      $self->burn_node($node, $p);
    }
  }

  my @sample_vertices = keys %{$self->{vertices}};
  $net = $self->{oldnet}->create_subset_network(\@sample_vertices);

  return $net;