dictionary - Python script to use coordinates from one file and add values from matching coordinates in another file -

-

i have original set of genomic coordinates (chrom, start, end) in tab delimited bed file. have additional tab delimited bed files contain of original genomic coordinates plus numerical value associated each of these coordinates. these coordinates can show multiple times in bed file different numerical value each time. need final bed file contains each of original genomic coordinates summed number of values found associated specific coordinate. examples of files i'm working below.

original file:

chr1    2100    2300  chr2    3300    3600  chr1    2560    2800

other bed file:

chr1    2100    2300    6  chr2    3300    3600    56  chr1    2100    2300    10

needed output file:

chr1    2100    2300    16  chr2    3300    3600    56  chr1    2560    2800    0

i need write python script this, i'm not sure best way is.

def fetch_data(filename1, filename2):     lines = []     data = {}     open (filename1) f:         lines = f.readlines()      line in lines:         if not line.strip():             continue         data[' '.join(line.split())] = 0      open (filename2) f:         lines = f.readlines()      line in lines:         if not line.strip():             continue         arr = line.split()         data[' '.join(arr[:-1])] += int(arr[3])      return data   open ('output.txt', 'w') f:     key,value in fetch_data('original.txt','data.txt').iteritems():         f.write('{0}    {1} \n\n'.format('    '.join(key.split()), str(value)))

Comments

Post a Comment

Popular posts from this blog

c# - Unity IoC Lifetime per HttpRequest for UserStore -

-
i'm trying clean default implementation of accountcontroller.cs comes out of box in new mvc5/owin security implementation. have modified constructor this: private usermanager<applicationuser> usermanager; public accountcontroller(usermanager<applicationuser> usermanager) { this.usermanager = usermanager; } also, have created lifetime manager unity looks this: public class httpcontextlifetimemanager<t> : lifetimemanager, idisposable { private httpcontextbase _context = null; public httpcontextlifetimemanager() { _context = new httpcontextwrapper(httpcontext.current); } public httpcontextlifetimemanager(httpcontextbase context) { if (context == null) throw new argumentnullexception("context"); _context = context; } public void dispose() { this.removevalue(); } public override object get...
Read more

Change the color of an oval at click in Java AWT -

-
i have draw ovals in java, , @ click change color. beginning tried change color after 20 ms, doesn't work. my code is: public class mycomponentnew extends frame { public graphics2d g2d; public mycomponentnew(string title) { super(title); setsize(400, 550); } @override public void paint(graphics g) { this.g2d = (graphics2d) g; this.g2d.setcolor(color.red); this.g2d.filloval(10, 55, 50, 100); } public void changecolor () { this.g2d.setcolor(color.blue); this.g2d.filloval(10, 55, 50, 100); } } and in class main method have: mycomponentnew m; m = new mycomponentnew("fereastra cu baloane"); m.setvisible(true); m.addwindowlistener(new windowadapter() { @override public void windowclosing(windowevent we) { system.exit(0); } }); try { thread.sleep(20); } catch(interruptedexception e) {} m.changecolor(); the color of oval remains red. ...
Read more

I am trying to solve the error message 'incompatible ranks 0 and 1 in assignment' in a fortran 95 program. -

-
i writing program in fortran find velocity of parachuting person in relation time. keep getting error can't fix. new programming , appreciated. the error is v(i+1)=v(i)+[32-((c*v(i)*v(i))/m)]*(h) 1 error: incompatible ranks 0 , 1 in assignment @ (1) and program is program para integer :: real :: v(11) !velocity real :: q !initial velocity real :: h !time step real :: c !drag coefficient real :: m !mass ! gravity equal 32 ft/s^2 write (*,*)'enter time step' read(*,*)h write(*,*)'enter initial velocity' read(*,*)q write(*,*)'enter drag coefficient' read(*,*)c write(*,*)'enter mass' read(*,*)m i=1,10 ! 1 10, 1 being interval. end v(i+1)=v(i)+[32-((c*v(i)*v(i))/m)]*(h) q=v(1) end program you cannot use [] normal parenthesis in expressions. array constructor, [ items ] means array items elements. end do should after line.
Read more