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

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

Appendix D<br />

# Aggregate BaFL cleansed samples from the Aggregation.py file<br />

# Individuals within the same disease class, if the sample size is large the population is<br />

randomized and divided into sub-populations bef<strong>or</strong>e intersecting the entire population, f<strong>or</strong><br />

computational efficiency.<br />

def StateIntersect(usr, pswd, db, logfile, rgr=4):<br />

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

states=get_unique_states(cur)<br />

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

notes='intersection <strong>of</strong> '<br />

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

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

k=i<br />

pt=find(i,' ')<br />

if pt>0:<br />

k=i[:pt]+i[pt+1:]<br />

msk, Lmsk, state, cel=get_inc_state(cur, i)<br />

if len(msk)>50:<br />

div=2<br />

while len(msk)/div>30:<br />

div+=1<br />

cur, conn =splitIntersect(cur, conn, k, msk, div)<br />

else:<br />

tmp= 'create table Intersect_'+k+' as select '<br />

+msk[0]+'_sr'+str(rgr)+'.probeset_id, '<br />

+msk[0]+'_sr'+str(rgr)+'.pm_mm_other, '+msk[0]+'_sr'+str(rgr)+'.x,<br />

'+msk[0]+'_sr'+str(rgr)+'.y, '+msk[0]+'_sr'+str(rgr)+'.probe_index, '<br />

tmp= tmp+ msk[0]+'_sr'+str(rgr)+'.signalrawintensity<br />

as '+msk[0]<br />

f<strong>or</strong> j in range(1, len(msk)):<br />

tmp= tmp+ ',<br />

'+msk[j]+'_sr'+str(rgr)+'.signalrawintensity as '+msk[j]<br />

tmp= tmp +' from ' +msk[0]+'_sr'+str(rgr)+' inner<br />

join '+msk[1]+'_sr'+str(rgr)+' on ('+msk[0]+'_sr'+str(rgr)+'.x =<br />

'+msk[1]+'_sr'+str(rgr)+'.x and '+msk[0]+'_sr'+str(rgr)+'.y =<br />

'+msk[1]+'_sr'+str(rgr)+'.y ) '<br />

f<strong>or</strong> j in range(2, len(msk)):<br />

tmp=tmp+ 'inner join '+msk[j]+'_sr'+str(rgr)+'<br />

on (' +msk[j]+'_sr'+str(rgr)+'.x = '+msk[j-1]+'_sr'+str(rgr)+'.x and<br />

'+msk[j]+'_sr'+str(rgr)+'.y = '+msk[j-1]+'_sr'+str(rgr)+'.y ) '<br />

tmp=tmp +' <strong>or</strong>der by<br />

'+msk[0]+'_sr'+str(rgr)+'.probeset_id,<br />

'+msk[0]+'_sr'+str(rgr)+'.probe_index'<br />

cur.execute(tmp)<br />

conn.commit()<br />

144

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

Saved successfully!

Ooh no, something went wrong!