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
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