From 30f9b0548da10ca6375506f2abcb82ba1eec86fb Mon Sep 17 00:00:00 2001 From: Shuyue Li Date: Fri, 26 Oct 2018 19:37:13 -0400 Subject: [PATCH 1/2] Add script for question 1 and output figure --- Q1fig.png | Bin 0 -> 13709 bytes script.py | 44 ++++++++++++++++++++++++++++++++++++++++++++ 2 files changed, 44 insertions(+) create mode 100644 Q1fig.png create mode 100644 script.py diff --git a/Q1fig.png b/Q1fig.png new file mode 100644 index 0000000000000000000000000000000000000000..c0b2eb2d79319d8e7676fc2bb706a33f6a8541de GIT binary patch literal 13709 zcma*ObyStn*Dea8bSd3#=`NK9*@&bdA)O-K-Jxu{q(Qo*yIUG*kdW@~Za8cIzH{$4 zzH`qXmoXUn!eYJeyXIPRKJ$5=9rQ^~0t=lK9RUFWOG@&+A_4+p0Qkm1MFu~~zenu{ z|Df4Os@o$VV7!2TBaY2`-GLtpxGAeSC|MghIP2RPA|Rk9_J3qy?ES<3=Z~x+gA6Ov z(7><~#~Ak*_oPBE10!SE*wnBe^GH9t^b`t0E)VW$Dy63<0>VoKsrRDFF3J1zu4>!N`8JY?C97w{r9@O-M;d2#+{X0Tmc`33k~$<w}Yy^?1Z)gY`2M0&2heN&AdhwUj{)~UP{#Z5^v;0h% z9n+XsP*BjJp9k0Dg2IOCIao-#kPS*bX zo1TtLNzrZ!Yi<@&S69!uz{0_y7ZfD-5qZ3e4kP1taCtNZcuC;B@xGp3=T zAp&8q&w>V5wFM8SXpm4}osf?G+1h%?U6tcY_lv!Zk4bMaNw{Evfv5L3^uFisBMIgPKcO>2C0zyMYiYdg~Rrk-!~+bu~J=;jiHqOg=Y8Tv$L;< z3(W&50;*}{RZN_mk#8+$V$^DEmZ5~K8nDt*>*ET?-ARQ)6;^a~bOc_{o7ej@)ql-Z zAh7Kq>A;ngw8hcu&rh2dipMj>F~GuoHmlFdqQ1R6oHE95t}HKK?u%ufYj6w!d&cf; zOF^UBOx(hPMZ3Xad22N5Y-_Y*gwqfA0~<3t1_~KUN^685q?mqVW5eNKE>YL(`k4ejJOz#_E2RbCwFOy-ex4A7gNwqG@(@(_u z;ZMqU0RsNX#9d<(t3|DTG?~JNuC{jv&J-g91EG{cZbt`=hhv4R?EKCLp^c3KL~L3w zg@jU)lgTxzOfojBps1Fg2=XQCV-@d~z1M#=Alx+)&#~1NS6`A%nk7tzJBoa+0BKrSZBjy8_0a4^QJG3@RA>q z4u^#eq<~Rep|`J}>6+-V`S9`KNWL{%=L?lRxBc8s{b*+mno!bNcD|+3&zXFYGw|zD zOUs#4Wl7tSZ>23a!TCzLT=pBX2!t$ZK475`gxtI zmAsK{kEyz;Q~h`Immkr(cD$ryGD#~qx7b+xF*irn%B9caHjrutc-kjG$1Zn8|DZ@UZkWKpmE674;yUIYrdmTc>mpY8?JH{B^atv`kOj z(`rdj83srfh$#Mc8$FCTE{(B2< z9o*>hJk{mN%Ib28^HQwp*;wl*!uq0YLn9+2!SivYu2KvP3|l+9FC*!qK5qX`M17G* zTRh!8JUoaw^bin6M@PG6?=0J4Xch}i>fXzK9M3tWINv9SzFhhGv~=G)j9t#F&71Nb z=U(uJgkPg;e;O%8K+KM7Z{``dF=-KpPcom0Hb@^qAtBT}JVaz!;|ii+9VZ|qcO?pHou zZEbBVI!y^+HIBif;gOQ|Eq4SUCh$8l8FYv0ulL0YK3uKJCi9lN=vYEv8y4pOb_}=7 zQlw;I`W;K-j4r#Z3EVyD@sSU$Oq%4F`Pz+6G<0+oR`U%j!hWsi`?K*f@oeShQ*2HL zb5a4=6oDZj?n|E*p{T6ayX@zuPyGD+2W?O329k1z$_-G?4z2$Q8s=4WiAJ{`IhywV zfDHV66SB2s!o|%kQ)F*rqc>M)8%o9#4u9O$2$8W0!+x%Sn-bmD7GVPn{j=~pSJ>=z zZF#Y&W($RU%J)d><@RUuPnp@Rtph`!%FBzm(^FEEe9Rd$>tsccbZM)&B0UZ>T08t^_G7qoc2Bu`A3NO z%cV$Dm8o;vIE-6vF%|R?pol9Job9J^nJ4EJd{Hvke*EZ-9)iwybcWc_;OKcJm*$O) zMb@4zrKi3wJcL;wEd|gEoMgkgP zqC%{wIPOg82b0tv`NlMnX=P;v0vnl{YPH3{uT*I?_@cuaRez_Z(^`MddG*NVNbfSb z(Y7#m`S5b4ddD@EUe4pHCUdSp*|w~rBGc6#&nsCc(9B4)@2xQJ8R9&d^JYM~$dcY_ z`gC}3W=r;XVf?7@d!wC#aYI@N%=8n!pZI8qi^#GkvJ1A}p_2O>c1#3`h=NTOY^(|G zUb7PE6h2f)Xj*|}myCe!m-FL^Gso@zOM{l4zV@AE0i7n){QR9mkt>9qaRr^*^OYU8MlxJBD*NNv z|5lr;3JM94a2sI)eC8iW3mYjv9P;WpNsHsIU;P?~!a8+VbJ`Ad)OXn+`PZl_RI%fJ zxc%1Rjq#6_vM*ieP0h_Ay2FXw0bqt7KYsWC6ly+M#QB~D0>i$KJVwE*kxQ4!V)OX- z*bz+8`4^AUoBRpVdXZ?)nGShWR1|=JD8quMnCS-z{5DGn)Q8#vm}qGFyOYH?;owm+ zlJCgBa`LFQIvN?#IXOA0Ew~4E4GxB_ zt$p?Ec>&!D%w&+Ac^hO;viXn04k^r%EdKE&8Bam$WJ&S7%}UGJ45Zof)+sq9<@ofJ zIdwI?C^zxYnSx7DF!|`n2BD&&LV6(jWgk5~y~IpuR5sR|N=8#WGIp!4cBg+5bUEp7 zq#cg=KhBgL-<;QZpP=#F-SPg-e8yyY`i1|VLix<(V6I*or1=z{tSx313_R;JX6w8k-qCCnqOkj&n@=puqcw>z!m{Q5}3*fpmR|^y9L)=g7&+ z%NxHzsmYMeVW>1>L}%zr;0%t6&uvVnD|prF>*eJ|YMDF59ZjJiv#DCDB)*8B{a zp03QP$&XDyJw7(10?WHxSgJSkTCM!azrBljUL)V_ZdJOzJ;qw`GZ-{{m9gm*M;tY;dO$uKa(9U||28*|;f!JWHFE=G^9@v2KS z&}{9)0oDOw+`U%27$Z8%BTkjY(+?4wp{`kxfZVR>}e zTEg#^RKLXvNFKWDikPW1F^u!XT;4X9ULVc2(@5IlX8z=Lxqi7}v{AIz`og{v-GQg& zl1;0%q($kEqL64rbC?iylZ5h~OM9Nwe+EaoJ)xqfN^%`jerdoqRs=>X-0`e|fFKHD#@vtJER4KFcXYk?O4@ z^k#E!`cq+?)f4nA>iAe#9J{l`%Fn;pOMZ;Py;tSglB^ z8##DrGj*Y*ijnQVM}5p?)>Kty&lD6!T3ZRn(3NPqL{agwm=u>EG6YnubQq_+@I2^a znqLWN3+-Ldr}-4KdR#$d6Q8K4#Nsf^#(@eYE32;n0Yft)+^O!8Zp^D95)(-bG$I)5 zj~*|XWm5{+UWA9w;IPVu!w1KBD|v66|5SaTAb3_|t$N8wL-VfH>uAA!X}|dXNBdW) zfGXSxDJ}A|eGD=Gynxy8;O{8{IUIqCa7c@1&!R4rQ{q?O)xjetrEp2>al@GXs~~=Q zw(wBc?ZI_h759=B{sJe1hl{uktoS%6INdcQnin5}pGbt~-Qu;}VkMN$+WU9J{?MaAYKY=)enLQ+fHQ?m@s5IX|d-DP1`oZ1C-dAotIZG533!ImB zgty7lKIRmvZall%W+}N@^hyHqo=0LYx$jy-hRIoCzSuj^rB$TesSfW4q}F>nkY%x&oh1vE znLIo_ZJnIx7#L8Uot^jg_R5-^_s9RcdSBYv1(K0TMk>M%+avVH8BFGg!>Lr|C)Bh2 zGw{$AWDlFw2?+@S$riP@XW@@U)(zmTjHQ87+HFloM+X`X4$4brW_X97{h@fw&dx6T zo+-b{a_e|hqEqLZWahi3wmh3G#mAE}eFk;)Pahw3dk^*3w<_*2oeF(ZQYZi=*xA|n zBPHb{h(%mn8~}!+W%?almx+R|T+#IM$mxdtae!gj<+Po6K|#T~uPU`}){*A*{jha& ztWqLp0Z)uzXKz1UYfZJVuyA%Tp9IJ_B;4oE z+eR|P&PilY=llAVME(Glf?jjH9du+HMkj%h99C{Z_L3hz(gMOTR$5#^;dz!sgb`?3 zudc3e$;k(HOIzhYBbD&lGzu_c8IB_N*Ld8Q18Yugmgn4P&3`l0xrcJ@yD?I9I2;5a zlm%K{%Lns~9w6#~3QaClEv+z{U;;>Rx=bGl^s61Sv$x-FYt5&^|NZ-?*sob-(pGNJ z(-*IWN`2lzCBor$YC<}0`c5K=(B=U1wds>V{uRr~#1ZTa|8Ry*F;gw_p)?WsPoHpk ztQR_mha*5FK=1eLH)tZK0p!#5JdYjvj;@d15_tdJ*xMy8yo4f%v@-0b83kIw-PuzR z(1_dq-o#`1G$&psQ!1Jc-aNWrTE%z^Ej<1+h_GGjeV+0Zza1YinJ3i!X!C^9dLi$y zVdUu#J-Yd6+G^F(=)N3$q-lKQ!X62LGt@UW4v@}aDj(|V@*5vlE)KwJXQ1Fkux=Lb z$)?HEjDLl1ya7#;rO- z2}1z_o(Q!3fq^>rOExOcv*$?2m?P?v5X_wmOSKmA1uspZrTMWO z4tu&DFaF-(ae}>;g!UvvYp&AGrd#3PjIjLCI^Mg>*2f%)2uc?h7hHUN7${gFMC=;_ ziQH?Oo2T>6i`!hexHq{BAf^MUSW;3PR>=*Lb@h?1^wEjAMkf+OR37|SN~zmdp&5VC zGd`0N1pvCrerKFm`~pGinaXm8*{>3bE z*(^VqyTeG$L!xZaOL}_Sjby{P`@)pZ)Pkhhk6uG4V{;t3mHMG8BR$;*&=jM9zK&y2 zZ;zstaY&ItO_6Di?F0Ra>E~7gtew}N*jX5SWa1}Inbg$3#S_&v#Qn^zuU0ZQ;z?um z6IE6xB|rgO&66ij0)vAa4(lYK2|snCyb+g9A!&b`Md_~=ANAUv)H!14B$qrqAo6JR z`q*YN{mtV|NJ?7W-j>qP_*rs_n!(dbbN@ZIo`7P=NK`~OFd`uR<<%0>xgZ8q*?@Gh z?`9*(y29o7i(Hy;alasROULQ$K$^VlIQ`W^*)OBTZ#!+rHCBG}ji4Zj^6_WKvuR6% z{-dR(<+a&3DiKkfbE`__>IDS#B)w@kFiTdr0x zv^Q{}JQs%ns=A-SA=t%mzj*fvmy=%jS$dKt4-yJ~JMU)`;egOkOwhmB*x3BoB!k&r zA4{l}e}D2B@$U{5^}Z8nRqeNOuaUgxUrWB97xG4XAxn@+%Ly|wt~ zQ(MaenwBcF33Ov)W6G9gZwq}h{t&_Xet?6F~xJRkw{^(3n zl^l?x;S`9RQK@zQX4DtM2xxWXU|){V%vY#r-$WjKgxo*wxi7e|V#|mLRsPW(AIqjy_dD=5 zR7T`WCo1@Vdi5u0Ouwde+vQ@1_ag=(VrOoTJ>M;tuOJEuer^&-bOjCQcQ@kwi+zY-lZ#E${$}mp(Q4|joi|)I?;x`A-CfIJ zJPYaIoRpPu01o+6(_DAAzj~z+UF8yqnJkIBM15g%ige@n=9@b$=kY`AG* z*AX80Dv)yr1_xbF2YC|Tm=oIC+GhXywX(A#o%GhS45(9vA0WM|>)g(ZS{0`ih|jRe zpAIu~6tgF5VCIg`(k@E>l4Sg1)8B?X$PCM8z7(^$xVSjmnUEn;ht-vGPYoY=4czk1 z@G3Ib;<8vlB?Njy0%VS9Yp~3ozn7Q_N0w%KKap;0n;_yXOhQ89aq@>gXYykr*E7&Q zCq9+jnj@_xBEDHuBo@G4xNjJz5rm@dk6c;h0GVa}GBiWARL28!n;ICjou19NlvU;^ zX*d;jR@XcgZeGZtaWjR8bW&9WNgC7weqRD`Lic(KzB zWNe1_6S>ERIY1<-1jpWK@1MdKS7l{DD~-cGdPxaFa+6b( z!Ji(rQER`@OeR<=IM8O3MSoI--2XMZ^Pu2R(r|O*fdu;m2?-4yea!HKLJAFiM`>$8 zk4yvo-Fx0QtmY@#1OMn?3o*BR0WF@lW3|@0+FzePsj*b@rr^i-0R-pu&1EUcIyX9y zm5`9}GTj)bzKVD^+HYdR%ZpO3oZ+7+TCI9J@vQq2on~pEfgN95`~op30v)~&8IKjb zgyiPZi;9YVH8lK86PJfxWA%)k-DCpY2g~@a4f7*m^gSIJZvyT_KT=2ss4~YVCrssg zpgDKg`xm+WH@8y11DTLrNA~t&pU?Z@PHb^tY%I!o(K8;w93OeC+uMA1D*fF(F_!7uk8MvfBvvVrvFYERpZs32Ia(4; z^{7~LTyS1+4LIE#J_YALDl$@3S{i-u=bOKIVl!ar!zm^@NXLsciL9-yqoSjOLPLS_ z`qbyg58Rpu&)VQT9j61>_b2B*Zol)*bFXNxiRhe@Stv%KThmClOi;N+e;g|u zQ<0LQ`$R|MaT^b@a4HtkfnNuOh2@w}m1O7SEV~^J;j5~u0)|;kmViF9D?D7gblb>UxqQS735(c$+%itLyw26>b8-T>BOu3z)C& z%3*K;OTAYQH%| z3n#2Dtqkzp@Clxkmm|BHi!1!|vO{V`Q4=j*Jz8aK;k(CLIB(9O)<8va)w50D(0>B3 zN|8o&FhJP+PJ6vzw*hLFX@L%a?%rO>up4r(n{E2Z5^b5;8q4GJ^JjX~6^1|tj)$Lp zqe0T2uT9(GbltTM4ce0;Y1GjGBcij0R-8G;oz8bGZxcfia*vbEi5KblZY&*iK*K|gi5-E>xSa^M3x zTak;!^JoG2wsx(1W>X_0{~Xu)v{8JhNzmOc-qGpsPN|T4d+QLoqv2s!vuxN;>GrUv z-Bpww?w*TT`t$z$-nu|@0i2ejt&QzBqd|w8w~V!lE{2!P7QB^+fGfBVtPgnUOcF1H zY5;hUyb!U=1O6oS;mwtDsn-vvA+oUl3`1sN%}}B*CYMvdpL}SEQ+GFtOicRIudQpOa3~Ix zTY2UlBIs5C7<;gpmf`81_qlK^AXq+P(~pAPW~&d<@OGy!HTPm^_O~hq#1HpICF7d> zJ%NH6YEO8+Ko;||TpPv4ev-s(#O$g2%%KO@7GAgfru>qG4=C@Md*(V{8%!3TWaJTW}y zS?9Iy?Hl|ediDe`7lS|I#2jK|QS~vid6ku_*>#@U+aFy_x}?#Hut7kOE^k`sJM>T$ zx}mXVRv`NQL7?&FJV37XP#k1%9M(kg|lDp=Knq8dQDGWhvDHSW(L#@HaqyP8hP2#)TcO=4??_|^MSf&~iw zCb;=2_WqhaSxD3Ex zLFje&{;Xk1%ld9ev+Fc{IN7kpSHRr5@f;4ujcVrUtBb`;0n3GCu-NBp#Won)4#~pB zf!|EEjHMl&P)e%d9;qu^)bV22W@L9n8LJD~Ta};lJ4~=V4X~8rhpjjg;qTr*6?||> zOv3*l5wB5m+ivy2nldb{kV>F?KfT6rJM$GRj}i=iV?6Zx?)s<`56I_+R#AX1Qku@? zXgwJ=%qiq4JsH3mURqu27&;a(95RL!@G)Y4mJ&r2ynp>~rky@*&Y9srpUmHOZJ64o zVC`g5u6X)(TxZX6?;KYJReu0XPb(AOZr$1z$GW6-p_~$5$z~76Vw&@YnnuVo>sa}d zq>wrO3U{RN%prghbAs8sN?Y5LtlOjI+mDtn=Pf+v+z>wSP=0)Oir|@p2TdSm#P&~T z*g+XRIjM4P3c||}MY`I`Lb1BBoU0WayLwJK9;43zOilD%ZMXnuI1`)C9?a_*J1Kp~<3@hGYQw#1GkOFc*T%qC7aahj0P6bJO z#ZX^gA2>O{bOEsm$iK#5&^yo0G6lDTx!}k@*<1iG0?@Rb z_uGz3%@S0e_$Vj=FZzM$3Jnbn5ty!E#hSG{+))qvbNOtzSrWVFmIhH@Km+qtYu47y z6QoyDuB*HvV0x-+Xn3-;v(tCJJB5`j+4I&)Sy$2ROrg0s2S{_7Vj+2^>e2KX+v%lZ zs0!-|$^i~8F1!E`AOHmQT)P_=w|hX4J%Vz!#IF8{x1q@iy($c{z>+3nAwy)=h z@+=mL7}6PTFL(D@MFfr$+o3V00>5I_yZ{VnHq0ydZ{OR0b|vnLQh|;rP!yseme6o<#M^{93r7; zxqGsNvIH!9fE`tE z&;ab6LuOkkDFysja^f7UW~6Ue%#ZF94%ZW(L106`-IJ5~DuJKCY~uUlJ&xOD$%wwy z#{Ad5aK@qt&0#>YukJ4EuBr@w?AQ0hU7=DWci1j0UD{H(y0bsmmL>rFc*j3hgL=19 zCK-2jL#yF>wFFHQzC2uL45-Yi6_vn?k>T!K$K#~znis4Lqd+qN>6~=-?Z#1OhlBN; zHj1bpZsu!7?1q!H$9G1*+`^xop5W-8&s}V);Whd%y>0UlJTkB_pV@ht4dxtNLOT2` zMK9mJGr?ilCz3L{UIwv3GXFBgk=sMwLh(?Hj|_n=NZF3mvHEi!)>0vV=|EepRc9Wp*$6Af{zD@WU%I_{QeWi+TKjwZ z#`?vEx5LMR840beB_%(3FBX%@CVV;&W8tfUWvND-3?5mV=r2MWuR&zUvw2PWb7Xcl ziIK5kwyS;fXr5Pn`f-#xeSi1U;-+U{xXAQfqDhplPmP9HF0It-)k7bx1PAwjR#rSp zQ#{23@?7MKo?dV!hwJTEMSDd)Qn#3fH#?;KQMgbLsNg*d_>#M%E}#Azo+fWiEc?yD z&$CnHYXb~VUUtldfXTrzGCF#HjodB}ZyO%()chV~LZi!s;;jD+Q6y*Dozk7~-T9b= zOd(*_G&|qk7&^;^O(sJIFH!Z0ccr*cVx0OSWfDKX9X?Xh*FLl|ebu|>x4P698Z41* zk*|wpKS>~hh_m)`~ERROYuT;u_%^Av@fi1Bh5X$$rRrrp^Z;PQnie{>h1S za;iVU73SiULmVy zyTzcD2&-bS1GzLRq)LX3IW|5q(QPl=2O&>EEK7UHI29$`wz<1O7EGoAiffbEB+I2` zAzc33USf>C00t|6EtN4}JVyE>^DW~oVkl|`C!hVuD}abJ$T}0xxM|~b`|JJF3gMQI zZQw!?_x}3O_RxQe2ABf-#*ME2V%8#(MW-gjXf(+CYY?@}txgjJ6&$Uw0N()p4phi{ zTFEqDLUH;JfAg%U=l0j9!9=FutkMQ&;-(OXd2@RYy!QB_0n8Dxqxb>9STtH;C?hB+ z2#Nx?L$cj1^(oGnRqpcQz~YDk*mGl{nY?+*Ouyt(`i4>jf+8Y}LA%x2(?iCuQhztU z2vkf|1p;DmTVAI3HTB*DLTI9 zm0Hv@C-tNr+d1j!E{y2Wa_e`>1t7k(~O`F|YPxpFB$oss&lkrMIy z^guHP+Z<@)U@tHdVgSdf$rdib`sSwb7Alyr6h0%|)a@|gMK=q4gt$=WwF-~{0u^OT zc=u+iLqL2;~?S4ux)fZXKHr8?4F%X0z^h7 z@U*;}DH8;ijIFU84B$P^1q2C;S{Z`R|HmzZ!>4%uOEDjmQY2qwtS1h$UUZ+ym6-_O zOv~m_3fN?oR`dKTtE*TP{aYUsIKF|+8mk6$s}Fj5F98k#fJQpQ6G28sX2J7Pdni>% zMphORtl=wRVV#*Wasb8l`w1FnYVLU``G>{+)g+2ic9jLzDBMMIeTJX~5-^jjuF&ZB2Y6kv@Co<4mF zV8@{ltxop#)O>u?K0&goz=B-@crSpT>;M;vM9+l< zEcl(>-DALUWZGc|Tu4_pH!FL4vOpy|I5+@`ioR*S7oqBdkjOtk68vv>kGF-fkSg3YKwPnPDir+hyRb0m4S$t$VNm&1UL;_M@QN>ZwLV6EhfQb1{BUvV$Sb? z1?Kj+vH^P{6UdZDmzVT(bpC)~K*zu+ud0gi9>zKPD=$)RIr}y}XZ zeSQ7jY%LC8wE(R-2>d*VrX47*yT!GV($cShe;r^`JY3x0Jc}M#vlo9g4*gR1{*{XX zpSv>Lqzx$Yaj-y88=JB*z>ao}jVYRr=S7}O8>cY>vR3-tJ0EadGT?y)wckJ7v=UHhkVUitFFaIPnD&E1D!9|i^<8!XK^y1Sr>7KG2}jat21(P zKIi7;|onXsg$KU}z{X8MM;#@)G-qhY*f}&;m1UP<8bi1fbt$feiazTs(H`5Eq&Tx16w| zVPdj8wCFUu@$&H<0}T~SqtGG6@17GYOgk1_VR4{{DVuuawBhHqfV(1B&f<5iGg%fc7g`3qXw!%-*)PZ-A~p z9^f91(NJym@u^df)#Sf77g$ literal 0 HcmV?d00001 diff --git a/script.py b/script.py new file mode 100644 index 0000000..ff2b629 --- /dev/null +++ b/script.py @@ -0,0 +1,44 @@ +#!/usr/bin/env python2 +# -*- coding: utf-8 -*- +""" +Created on Fri Oct 26 10:47:03 2018 + +Biocomputing Exercise 8 + +@author: syli +""" + +#%% Challenge 1 + +import pandas as pd +import numpy as np +import matplotlib.pyplot as plt + +dat = pd.read_table("UWvMSU_1-22-13.txt") +dat['UWscore'] = 0 +dat['MSUscore'] = 0 + +for i in range(0,dat.shape[0]): + if i == 0: + if dat.loc[i,'team'] == "UW": + dat.loc[i,'UWscore'] = dat.loc[i,'score'] + else: + dat.loc[i,'MSUscore'] = dat.loc[i,'score'] + else: + if dat.loc[i,'team'] == "UW": + dat.loc[i,'UWscore'] = dat.loc[i-1,'UWscore']+ dat.loc[i,'score'] + dat.loc[i,'MSUscore'] = dat.loc[i-1,'MSUscore'] + else: + dat.loc[i,'MSUscore'] = dat.loc[i-1,'MSUscore']+ dat.loc[i,'score'] + dat.loc[i,'UWscore'] = dat.loc[i-1,'UWscore'] + +plt.plot(dat.loc[:,"time"],dat.loc[:,"UWscore"],'r-',dat.loc[:,"time"],dat.loc[:,"MSUscore"],'g-') +plt.legend(['UW','MSU'],loc='upper left') +plt.xticks(np.arange(0, 41, 10)) +plt.xlabel('Time', fontsize=15) +plt.ylabel('Score', fontsize=15) +plt.grid() +plt.show() + + +#%% Challenge 2 From 2d55f296b9919b5b34faa87a2bc015cd4f1f8dcc Mon Sep 17 00:00:00 2001 From: Shuyue Li Date: Tue, 30 Oct 2018 10:50:53 -0400 Subject: [PATCH 2/2] Add answer for Question 2 --- script.py | 33 ++++++++++++++++++++++++++++++--- 1 file changed, 30 insertions(+), 3 deletions(-) diff --git a/script.py b/script.py index ff2b629..8a2423a 100644 --- a/script.py +++ b/script.py @@ -8,12 +8,12 @@ @author: syli """ -#%% Challenge 1 - import pandas as pd import numpy as np import matplotlib.pyplot as plt +#%% Q1 + dat = pd.read_table("UWvMSU_1-22-13.txt") dat['UWscore'] = 0 dat['MSUscore'] = 0 @@ -41,4 +41,31 @@ plt.show() -#%% Challenge 2 +#%% Q2 guess my number + +print "I'm thinking of a number 1-100..." + +from random import randint +number = randint(0, 100) + +try: + guess = int(raw_input('Guess:')) +except ValueError: + print "Not a number" + +while guess != number: + if guess > number: + print "Lower!" + try: + guess = int(raw_input('Guess:')) + except ValueError: + print "Not a number" + else: + print "Higher!" + try: + guess = int(raw_input('Guess:')) + except ValueError: + print "Not a number" + +print "Correct!" +