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.

#determine outliers<br />

lowrs1=nonzero(less_equal(mu, ceil(r.mean(mu)-<br />

(2*r.sd(mu)))))<br />

upprs1=nonzero(greater_equal(mu,<br />

flo<strong>or</strong>(r.mean(mu)+(2*r.sd(mu)))))<br />

lowrs2=nonzero(less_equal(prbs, ceil(r.mean(prbs)-<br />

(2*r.sd(prbs)))))<br />

upprs2=nonzero(greater_equal(prbs,<br />

flo<strong>or</strong>(r.mean(prbs)+(2*r.sd(prbs)))))<br />

f<strong>or</strong> j in lowrs1:<br />

cmn=nonzero(equal(j,lowrs2))<br />

if len(cmn)==1:<br />

r.points(x[j],mu[j], pch='X')<br />

tmp="update sample_mask set exclude = true,<br />

description = 'array probe intensities and numbers below 2 stdevs"+"'<br />

where mask_id = '"+msk[ptr[j]]+"'"<br />

fp.write('\n'+msk[ptr[j]]+' excluded:\tprobe<br />

intensities and numbers below 2 stdevs')<br />

else:<br />

cmn=nonzero(equal(j,upprs2))<br />

if len(cmn)==1:<br />

r.points(x[j],mu[j], pch='X')<br />

tmp="update sample_mask set exclude =<br />

true, description = 'array probe intensities below 2 stdevs and numbers<br />

above stdevs"+"' where mask_id = '"+msk[ptr[j]]+"'"<br />

fp.write('\n'+msk[ptr[j]]+'<br />

excluded:\tprobe intensities below 2 stdevs and numbers above 2<br />

stdevs')<br />

else:<br />

r.points(x[j],mu[j], pch='X')<br />

tmp="update sample_mask set exclude =<br />

true, description = 'avg array probe intensities below 2 stdevs"+"'<br />

where mask_id = '"+msk[ptr[j]]+"'"<br />

fp.write('\n'+msk[ptr[j]]+'<br />

excluded:\tprobe intensities below 2 stdevs')<br />

cur.execute(tmp)<br />

conn.commit()<br />

f<strong>or</strong> j in upprs1:<br />

cmn=nonzero(equal(j,upprs2))<br />

if len(cmn)==1:<br />

r.points(x[j],mu[j], pch='X')<br />

tmp="update sample_mask set exclude = true,<br />

description = 'array probe intensities and numbers above 2 stdevs"+"'<br />

where mask_id = '"+msk[ptr[j]]+"'"<br />

fp.write('\n'+msk[ptr[j]]+' excluded:\tprobe<br />

intensities and numbers above 2 stdevs')<br />

else:<br />

cmn=nonzero(equal(j,lowrs2))<br />

if len(cmn)==1:<br />

r.points(x[j],mu[j], pch='X')<br />

139

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

Saved successfully!

Ooh no, something went wrong!