David Rio drio

🐢

I don't know

Builder of things.

drio / gist:111723

Created May 14, 2009 15:43

	I am trying to get my head around the sorting example. I use this as a input:

	$ hadoop fs -cat /input/small*
	9971681
	9686036
	2592322
	4518219
	1467363
	607354
	....

drio / gist:120707

Created May 31, 2009 01:49

	+ cartera DNI
	+ iphone
	+ pasaporte + documentacion
	+ laptop, power adaptor, international adaptor
	+ book -- the economist
	+ papers

	+ gomina
	+ deodorant
	+ cuchilla + crema ?

drio / gist:133299

Created June 20, 2009 21:47

drio / gist:133785

Created June 22, 2009 03:01

	DATA.each do \|line\|
	i = 0
	# For each match of "any character followed by the same character" do:
	marked = line.gsub(/(.)\1/) do \|match\|
	i += 1
	# return the match surrounded by asterisks from the block, which will
	# tell gsub to replace the match with that
	"#{ match }"
	end
	puts "{#{ i }} #{marked}"

drio / gist:138951

Created July 1, 2009 18:20

	#!/usr/bin/env ruby
	#
	# vim: tw=80 ts=2 sw=2
	#
	# merge_eg_output.rb: Merges multiple egenotype outputs
	#
	# Usage:
	#
	# 1. Concat all your eg outputs in 1 file:
	#

drio / gist:138962

Created July 1, 2009 18:34

	#!/usr/bin/env ruby

	in_file = ARGV[0]
	sqtocs = { "AA" => 0, "AC" => 1, "AG" => 2, "AT" => 3,
	"CA" => 1, "CC" => 0, "CG" => 3, "CT" => 2,
	"GA" => 2, "GC" => 3, "GG" => 0, "GT" => 1,
	"TA" => 3, "TC" => 2, "TG" => 1, "TT" => 0 }

	# seq space Old/NEW format:
	# 1 775852 rs2980300 A/G A GACTTCACTAACTCANAGAGACACAGTCATT

drio / gist:139529

Created July 2, 2009 15:35

	#!/usr/bin/env ruby
	#
	# Author:: David Rio Deiros (mailto:[email protected])
	# Copyright:: Copyright (c) 2009 David Rio Deiros
	# License:: BSD
	#
	# vim: tw=80 ts=2 sw=2
	#
	# genome_hasher.rb: Hashes a whole genome and allows you to query it for k-mers
	# sequence

drio / gist:139530

Created July 2, 2009 15:35

	#!/usr/bin/env ruby
	#
	# Author:: David Rio Deiros (mailto:[email protected])
	# Copyright:: Copyright (c) 2009 David Rio Deiros
	# License:: BSD
	#
	# vim: tw=80 ts=2 sw=2
	#
	# genome_hasher.rb: Hashes a whole genome and allows you to query it for k-mers
	# sequence

drio / gist:141567

Created July 6, 2009 17:31

	class EgOutput
	attr_accessor :chrm, :pos, :id, :al_cs, :al_seq, :counters

	def initialize(others, counters)
	@chrm, @pos, @id, @ref, @var = others
	@counters = counters.map{\|c\| c.to_i }
	end

	def +(ec)
	new_obj = self.clone

drio / gist:141576

Created July 6, 2009 17:47