python - curve fitting by a sum of gaussian with scipy -


i'm doing bioinformatics , map small rna on mrna. have mapping coordinate of protein on each mrna , calculate relative distance between place protein bound mrna , site bound small rna.

i obtain following dataset :

dist    eff -69 3 -68 2 -67 1 -66 1 -60 1 -59 1 -58 1 -57 2 -56 1 -55 1 -54 1 -52 1 -50 2 -48 3 -47 1 -46 3 -45 1 -43 1 0   1 1   2 2   12 3   18 4   18 5   13 6   9 7   7 8   5 9   3 10  1 13  2 14  3 15  2 16  2 17  2 18  2 19  2 20  2 21  3 22  1 24  1 25  1 26  1 28  2 31  1 38  1 40  2

when plot data, have 3 pics : 1 @ around 3 -4 1 around 20 , last 1 around -50.

i try cubic spline interpolation, does'nt work data.

my idea curve fitting sum of gaussians. example in case, estimate 3 gaussian curve @ point 5,20 , -50.

how can ?

i looked @ scipy.optimize.curve_fit(), how can fit curve @ precise intervalle ? how can add curve have 1 single curve ?

import numpy np import matplotlib.pyplot plt import scipy.stats import scipy.optimize  data = np.array([-69,3, -68, 2, -67, 1, -66, 1, -60, 1, -59, 1,                  -58, 1, -57, 2, -56, 1, -55, 1, -54, 1, -52, 1,                  -50, 2, -48, 3, -47, 1, -46, 3, -45, 1, -43, 1,                  0, 1, 1, 2, 2, 12, 3, 18, 4, 18, 5, 13, 6, 9,                  7, 7, 8, 5, 9, 3, 10, 1, 13, 2, 14, 3, 15, 2,                  16, 2, 17, 2, 18, 2, 19, 2, 20, 2, 21, 3, 22, 1,                  24, 1, 25, 1, 26, 1, 28, 2, 31, 1, 38, 1, 40, 2]) x, y = data.reshape(-1, 2).t  def tri_norm(x, *args):     m1, m2, m3, s1, s2, s3, k1, k2, k3 = args     ret = k1*scipy.stats.norm.pdf(x, loc=m1 ,scale=s1)     ret += k2*scipy.stats.norm.pdf(x, loc=m2 ,scale=s2)     ret += k3*scipy.stats.norm.pdf(x, loc=m3 ,scale=s3)     return ret   params = [-50, 3, 20, 1, 1, 1, 1, 1, 1]  fitted_params,_ = scipy.optimize.curve_fit(tri_norm,x, y, p0=params)  plt.plot(x, y, 'o') xx = np.linspace(np.min(x), np.max(x), 1000) plt.plot(xx, tri_norm(xx, *fitted_params)) plt.show()

enter image description here

>>> fitted_params array([ -60.46845528,    3.801281  ,   13.66342073,   28.26485602,           1.63256981,   10.31905367,  110.51392765,   69.11867159,          63.2545624 ])

so can see idea of 3 peaked function doesn't agree real data.


Comments

Post a Comment

Popular posts from this blog

Why does Ruby on Rails generate add a blank line to the end of a file? -

running rails generate controller foo home results in: class foocontroller < application controller def home end end (there's nothing on line in actual file, blank line) is there purpose blank line? it's common use newline indicate end of line, last line. many editors (like vi, use) add in newline silently. consensus text files (especially in unix world) should end in newline, , there have been problems historically if wasn't present. why should text files end newline? the tool used count lines in file "wc" counts newlines in file, without trailing newline show 3 instead of 4. it improves readability of templates used in generator. consider: https://github.com/rails/rails/blob/master/railties/lib/rails/generators/rails/controller/templates/controller.rb to remove trailing newline, template have last line of: end<% end -%> instead of: end <% end -%> that seems less readable me. the discussion below refers ...
Read more

keyboard - Smiles and long press feature in Android -

i've created android keyboard. how add features inserting smilies, , using long keypresses display symbols? to create long-press behavior allows select several options pop-up menu, follow guide (specifically pop-up menu section): http://developer.android.com/guide/topics/ui/menus.html smileys (also known emoticons) combination of keys, understood programs , parsed same way. have do, when clicks smiley face on keyboard, simulate pressing of key combination matching emoticon. example, when clicks smiley, should simulate key presses ':' , ')' :). here's list of emoticons: http://en.wikipedia.org/wiki/list_of_emoticons
Read more

node.js - Bad Request - node js ajax post -

i'm using express 3.0 , node v0.11.0. have button submits form , i'm using ajax post data json object node server. client code is: $('#contact').submit(function(e) { console.log('submit called'); var formdata = json.stringify($('form').serializeobject()); console.log(formdata); $.ajax({ url: "http://localhost:3000/savedata/", type: "post", datatype: 'json', data: {objectdata: formdata}, contenttype: "application/json", complete: function() { console.log('process complete'); }, success: function(data) { console.log('process sucess'); }, error: function() { console.log('process error'); }, }); return false; }); then on server: var express = require('express'); var app = express(); app.set('views', __dirname + '/views'); app.set('view engine', 'jade'); app.use(express.bodyparser());...
Read more