02.08.2013 Views

Sample A: Cover Page of Thesis, Project, or Dissertation Proposal

Sample A: Cover Page of Thesis, Project, or Dissertation Proposal

Sample A: Cover Page of Thesis, Project, or Dissertation Proposal

SHOW MORE
SHOW LESS

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

fp.write('\tTable '), fp.write(msk[k]+ '_BioP '),<br />

fp.write(' created.\n')<br />

pid_tbles.append(msk[k]+ '_BioP')<br />

trl='drop table tmp'<br />

cur.execute(trl)<br />

conn.commit()<br />

cur, conn=UpdateW<strong>or</strong>kReg(cur, conn, msk[k]+'_BioP', exp,<br />

'biophysical_filter', usr, msk[k]+'_XH', 'known_seq_biophysical',<br />

notes)<br />

k+=1<br />

onto='select * from '+msk[k-1]+'_BioP'<br />

cur.execute(onto)<br />

desc=cur.description<br />

conn.close()<br />

fp=UpdateOntology2(usr, pswd, db, desc, exp,<br />

'biophysical_filter', '_BioP', notes, '', fp)<br />

fp.close()<br />

def DriverBioPhy(usr, pswd, db, logfile):<br />

# first step when combines chip and cel info<br />

cur, conn= make_connect(usr, pswd, db)<br />

msk, msk_tab, Lmsk, state=get_mask2(cur)<br />

#print len(msk), len(tables), len(msk_tab)<br />

t = datetime.datetime.now()<br />

EpochSeconds=time.mktime(t.timetuple())<br />

now = datetime.datetime.fromtimestamp(EpochSeconds)<br />

fp=open(logfile, 'a')<br />

fp.write('\n')<br />

fp.write('Biophysical filters:\n')<br />

conn.close()<br />

BioPhy_filter(usr, pswd, db, msk, msk_tab, fp)<br />

def Fnl_filter(usr, pswd, db, msk, msk_tab, fp, rgr):<br />

cur, conn= make_connect(usr, pswd, db)<br />

exp=get_exp(usr, pswd, db)<br />

notes='probesets w/ statistical rig<strong>or</strong>'<br />

k=0<br />

#print 'Map length ', len(map), map[k]<br />

pid_tbles=[]<br />

#fp=open(logfile,'a')<br />

f<strong>or</strong> i in msk:<br />

trl= 'create table tmp as select count(*), probeset_id from<br />

'<br />

trl=trl+ msk[k]+ '_BioP group by probeset_id having<br />

count(*) >='+str(rgr)<br />

cur.execute(trl)<br />

conn.commit()<br />

trl= 'create table ' +msk[k]+'_SR'+str(rgr)+' as select '<br />

trl=trl+msk[k]+'_BioP.* from '+msk[k]+'_BioP inner join tmp<br />

on (tmp.probeset_id = ' +msk[k]+'_BioP.probeset_id) <strong>or</strong>der by '+<br />

msk[k]+'_BioP.probeset_id, '+msk[k]+'_BioP.probe_index'<br />

cur.execute(trl)<br />

conn.commit()<br />

136

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!