@@ -25,4 +25,15 @@ class Gene < ActiveRecord::Base "Load gene: #{gene.name} is not a valid coding ORF") end end + + def average_residue_weight + length = Bio::Sequence::NA.new(self.dna).translate.gsub('*','').length + total_protein_weight / length + end + + def total_protein_weight + protein_less_stops = Bio::Sequence::NA.new(self.dna).translate.gsub('*','') + protein_less_stops.molecular_weight + end + end model/gene.rb b1dcdd2ce0469b4f7b2fdc6e0da3a218f553be3e Michael Barton mail@michaelbarton.me.uk http://github.com/michaelbarton/cost_in_evolution/commit/1cc18d17d56d0d41d3acd0ccfa2c38470e2cbfaf 1cc18d17d56d0d41d3acd0ccfa2c38470e2cbfaf 2008-12-03T06:59:28-08:00 2008-12-03T06:59:28-08:00 Created gene methods for total and average molecular weight 43864dc21db05a0e9affb69c09a6623f4b15e766 Michael Barton mail@michaelbarton.me.uk