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.

# 2 std dev<br />

r.lines(x, r.rep(r.mean(mu)-(2*r.sd(mu)), nbr),<br />

col=r.rgb(227/256.,26/256.,28/256.), lty=3, lwd=1)<br />

r.lines(x, r.rep(r.mean(mu)+(2*r.sd(mu)), nbr),<br />

col=r.rgb(227/256.,26/256.,28/256.), lty=3, lwd=1)<br />

# 1 std dev<br />

r.lines(x, r.rep(r.mean(mu)-r.sd(mu), nbr),<br />

col='turquoise', lty=4, lwd=1)<br />

r.lines(x, r.rep(r.mean(mu)+r.sd(mu), nbr),<br />

col='turquoise', lty=4, lwd=1)<br />

# mean<br />

r.lines(x, r.rep(r.mean(mu), nbr), lty=2, lwd=2)<br />

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

(1.5*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 probeset intensities and numbers below 2(1.5)<br />

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

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

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

else:<br />

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

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

description = 'avg array probeset intensities below 2 stdevs"+"' where<br />

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

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

intensities below 2 stdevs')<br />

cur.execute(tmp)<br />

conn.commit()<br />

f<strong>or</strong> j in upprs1:<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 probeset intensities above 2 stdevs and probesets<br />

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

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

intensities above 2 stdevs and probesets below 1.5 stdevs')<br />

else:<br />

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

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

description = 'avg probeset intensities above 2 stdevs"+"' where<br />

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

142

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

Saved successfully!

Ooh no, something went wrong!