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