From abd1d73a8fc4a5197fc3ffdcb40044547f0dfb0b Mon Sep 17 00:00:00 2001 From: Morten Hjorth-Jensen Date: Wed, 3 Sep 2025 13:40:18 +0200 Subject: [PATCH] update week 36 --- .../_build/.doctrees/environment.pickle | Bin 273657 -> 273661 bytes .../_build/.doctrees/exercisesweek37.doctree | Bin 36403 -> 36790 bytes .../html/_sources/exercisesweek37.ipynb | 59 ++++---- .../_build/html/exercisesweek37.html | 15 +- doc/LectureNotes/_build/html/searchindex.js | 2 +- .../jupyter_execute/exercisesweek37.ipynb | 59 ++++---- doc/LectureNotes/exercisesweek37.ipynb | 59 ++++---- doc/Programs/Regression/ooreg.py | 34 ++--- .../Regression/predictions_kernel_ridge.csv | 101 +++++++++++++ doc/Programs/Regression/predictions_lasso.csv | 101 +++++++++++++ .../Regression/predictions_linear.csv | 101 +++++++++++++ doc/Programs/Regression/predictions_ridge.csv | 101 +++++++++++++ doc/pub/week36/ipynb/week36.ipynb | 7 +- doc/src/week35/programs/lasso.py | 139 ++++++++++++++++++ doc/src/week37/exercisesweek37.do.txt | 9 +- 15 files changed, 668 insertions(+), 119 deletions(-) create mode 100644 doc/Programs/Regression/predictions_kernel_ridge.csv create mode 100644 doc/Programs/Regression/predictions_lasso.csv create mode 100644 doc/Programs/Regression/predictions_linear.csv create mode 100644 doc/Programs/Regression/predictions_ridge.csv create mode 100644 doc/src/week35/programs/lasso.py diff --git a/doc/LectureNotes/_build/.doctrees/environment.pickle b/doc/LectureNotes/_build/.doctrees/environment.pickle index 84a382b7141aeafdca03fe60295e8d525adc19be..98ef94a83541a27eae1bf14ac1d65c5f561f52b7 100644 GIT binary patch delta 10234 zcmZ9Sd0h`0GILVKJ%e$@MSf#a)04)G0(0iL3|C5WPSfI| z3fi-2lj5TCsp&c6GC-G;GBIa(Mn+mTl(DMxIYJ>**|NAOAUh{1r%7sdPHHmv`sbu1 zC67)S2l)S*gs#bBa$D};fHh+ro#`r<8?7p#jP!;}{%{bdOjJvKIrGu&r{?}HnTv3qK_2!LI zGxO4i7rP#zy0TU`v$T(&)%tvS`Dc7vSiTpH4uiuRK6ti7E6Y~y*%EP<_L4DuaJUgO zI542KrMhcu8XV2;8Ak?JwfzHC8BIG^=ND|oo_003>$Va7MooUlW;AVIqrmkL@8nP8 z)@baAV%JraE1BbLRd<t$(Pa2ytG$&(dW>U`R;wi;3uKQqf`9b4c ze<9hCt|JxrBjRAXDj?zbxdM=I%v=?baGYF~@aH1x0+DchT$Pb0!tlBxkuavNC?t%eD;f#o=c*0FT~ECIhXjZMjF89fp5=@@-QJI)W&^-2-o()LCz}YJ@%inObLVRvCr6ON)K@yJoc#Drijt@*w^-rr-B%9kNw55@3B8IObMg!vA5XYmTWiz9{ZJY?MT8=@Yqi_G$kAfk9~l<1KD6Ss(=_% zvonb}kv;YgI=>Wo+PIa}ja@L>4h$wd5 zNPypejnVM?$Z(H{;*V{{x)G6#8wWsj(YOe|4-Eg5C|<6-Q9mV;T{03t5ao?j5bhf8 zlqkp;{{huMHsc_uE_g4dbYtxH0@p-i{irY49fWhEzo?D#a!Rpl40IkBS2|o@>nmIt za80!3X1lH$dsCy>W#cUTUirk;ig8xnh)=7{1p-D)%T9CYj0SrZ&@v4L~eJ80hEzxzp{`CY?OrhF3O@P3&l}OsCNzZ=z9>?KDdB45I?vdFG_~#@bA$EG*$Q zjlTIEjYm_QhRQEC4rV!x^vT1GI=N06v6ELd0;Y5{QZt;!*eUajoE)bSHMQ6{I@u{N zEf?*L(1KpHU|lu+U>#FRi$kK6_7W0hv_43b(UOoTqh%mbMw@~}8Ep;{Wwgael+o5A zQAXR2L>cWXB+6*V{4jA*N;`u}N@-V-aM8IQB2h;3w*yf|3rC`i_8byrw8lu3(b^$V zM(d758EpU(Wwa4Ul+nf^;aY0qDnO)^R)j<;Z7~wAAJ<1fxka{|4A@rC6TrJ*S{)`c zZ?F)VS(*8JOQ$y=(IHFw@E~u*j5Jg9E_%mh8LskJ?~AiWYs1dju>vb_lnBwn8#AYk zcKgUrR56CEtms`>bX~j3cSAY3$`z-DpWGQEg1pD(yR|iKL3`8c-BtXVsVf;Q+Izn) z@o0}>i%LwA*IyPkZ>1%3wAiQ#5A^O>>eYVv%Qq8+dE2in*1F$UnPtPgbKi&hMz-M@ zu!>q3yVfMgtsS`Te#D2+d=SIDsUM{hIlOO?5vRt7(nb^6Ri2+C+YT?wF=&!p=#W#4bO5kBHCbYu)bM?`aygXM(0ndnap} zuuqfR(vJms!@i!V4J-Dq*EIOs4>ScG{7}<3-!0JOI#i(P;o)>m`K4nt-9P%Sro}&e zqG`yF?=xx3hY+4g?`UlPIYAyNv_cs4zt#*zhq3HzfOo_%&^VO5H;+%zH2oxWn{=V4 zzSH#g>7|;MovqGI>$kr_Tgh+gAJAQ(3y=}uki zsCA8RD6PA96WTv%xyXRT*}!P=PlD|Di_lYG;XkbyiVp9oJMRFMw}p9g??G*(^LqdT zN;LD~DNQ+#yKz%j`lN`-VJAhH%3us;2s)|?f_mdy8Vpu5Y3D{T!KL#vo`H{#O!;~#+-zDTdy^qL)U_zQ z1L3;b45l=q9J57KQ*S!ia3+Iq^Dz0q3;q714uer6*Rl23YK_q`Y^yG{r9iY&)nXZ> z81h=!fW57;enY^L#-fVc(vDSBhZ#jabec7r|-Ft<7Z#o#0mRC||yRl>1e_tZ|%|yAf zC#$5se2pE}CXjb%9N7y;uudBE>94HXGN3Q?X=|olaVo7ZYX%E56NU#ezz~g*16T%O zm|8uM!OSNI*I=9$fMF``O*W0O2Mc~qhpj-I9P$=RFnTO9SAqKCEf&u}bSVE}>M#oBQbJMUr}P3?})G3sRBy&H#{u;(7uOXDYd zS+|-OfNqxocXxy`(!NaDw7i2Nmz^2Ha1Ffl5{{U;tP-C0uv`W zWdG{JMyJ^dtv!8)MRKF+nhe}zOQXjgy%HmShs6M*L%sP2%O#Q>e-NSS&Y#c~GIY7X zD!@tvm$xplS*GsJWmp96*C2xA?0+57KGEAZSfE+K>RYUaPWQbSkfuhN^Hohm;1gLk;rh)0y$)rwsj`QCW@-;0ear;~}zO z1%%1!X~&;a3j=uvQ;dt`(aJA~|6uxE4(5d>&2;cJCc!7H_f7gc zoI^6h`sP>VknDd;hbF>njht7V_m>^>pow(f3N<*)B#<4de*}NiMDxzDN9mwVaO*&L zpp#cL|4eJ0&vD2dq&@HCum>XCSC{WGQNpbQ)vD+DEG?gm=M`0#`W#k7L4{)UXMMZ8Di2_2iY*b0p7!`Q2>n|Iwhx-p}%d0WO5Cu1d zMMN0+B8SZl;g=mSVuXR6_)U$$U3eAdNO?(rYK?CGLUe96e9Tm>y7INAPiRrD)A|(K zuz5q%{NyDLD-LUu@VvImt<^<{+VBd8-5cHFW&#i5COQ)NFzvc&8?UC8zG`{NSDmnP zdGH#)p-Vwb7BQI~$ZE;6vAj3Ds>aKH-vZ9So#*v(D?c=2t2#0t|jr|aFxB1 z`4zZKOX1To7cAaU*U<2Z*xVXAEobPS8DW4U~e zN$C^#D3dNs`dIkR{Q&VItw0*n&jztcEw zUC`m1w>dl;2;Y8(f1`2GEPha9b2neA@k}BAK%+5-f3GoPE{}yrAy&429nxX7aGdAE-A*{BE=bI=Y znO9U>H*?rYp-JBVmPe?@TVc+au6;k|e-c(R*Vb?^Z$okwIj09C-u>HnU+vmp2X?a^ zd@I6tmEQQRW~dB4z(PwT#CpX}%sw5DK})La;(^RqaytXQ;PjN|K7$sQ+N>78_Bpl~ zWQQu)%`-K!J$#(TnR{VEnhurrahFE-SD1aN*qyI=md3nqIK0b}7CHu75yI4i9QIIz zXTRg{(nC1-kgwizhjHnVc6cdYr7`Rn&P~)xSizlCz5ntq+62oOwry14{!bj<`3Q~w zaoEui_BzhL)L7>vEF)73>B4I1>!&$vMalQX85}6WjB_|E38R1KB^nR>!H*J#sg39P zUag&ak;9gOe8*kp+cd`i#b;_Py=v8K?3`Uy<^PQdiTv+f$4Q4MmlX4s>d;NhLL@^L zs;geP&AV$EQc(kSC`+p8gsuDrmg>f`xMS15s7gZ)}xP4VzjfmC#%yGnu*xl0wwfzB0X~NRUl- zSxf9hOHq@d?2vm_LVl_dCt%NjvKkgI9%>VyO=YzakO0Y4zMX)XM7XRyhJg?|yjmB) zV3|pRcM`A|$)4F+OxE_H+~%ijU9n|g3RC$nX>=&~k_S%!EPx;N0#r(ZfJf+4`Dix* z3zFo`M{P zb^8fezf|ah{x~dz2?GU8Xv7_Z1*}<8u6sj3+9qr{1m_pw?zeo}Zo|ZEQ;Qoon>wE) zV3Dud3q4Lt7O>8+fH7=QJGpBvEB1oR!BwkRb zvqS-krUO;Z3N{`lIVY1RTXvC!nqK8&a z*)AHZiFsnMmf=%Rtg4(Z-~mI0KH4dsQ^k|T08>_vrwB-}Xj1Sgs90ma>0+Tqc)ip^ zrPjZLYmR(-&l2-A-kL4oVMvw3t7o`8Qv@wqQD}8Tu{plMTVG@ifa>-xu0*Qqz&w$~ zO-x)M=45u;kmXQSCxom=F)W7j3V)?ZnQ?cAu~10HTA z@r+`ooU{{mJquQd2>AZE6Xw7dE5$U8Jy#3ZAfOh$uZZrd#|JnQ$ZORak)%zhQ$@6D z_L1nL7s!;wH{{*nc8l|1x%)H`$VvUnUJGR zX6Rsm01!YQ8!NYzxJv)KxxK_zDf46|gCSm33F_;AIy1y)U^{vu$jL2W;C9CH?7)=6Jy15ROrCGGOlmV0J*rc)6MX^ES@k;`>v*hc(f(d~z_$n?M zx&yrpZM=pnm`oL~V@4$ty%_DeffJuhRc~P)CX}~*;Rw2ed5H>~jD?iT^3Dq& z*SPN+z*!IQQArgveT3PS@Zw`%&y${50o;DxX2tt0ZX2#|@S?xXu6_|VD^k76+tPF` zhy84jDXBv6g8c&FkqWk18dECTW}7H8>%qf~1=^M|IxI{KrIW(MNIE7=452f^#L9F) zm{@^M2NQ$nXfQFD&IJ?w=uj{*fKCLVI35Rr4bK9>6AuF6f`@>K;dBC+$m#eGD#Nos zD2E4sAU%8;uYq*t2U8EHltX=WK|woifs*!iwPpNLXP8EH`&X-OGrSs7_X8EI7+ z>HRX&mZy?jUCY)G-nC0>*_v{iU_**&g24sPq2Ypu&~U*MXt?0Q3g zh6^4!!v!Z8TyTQH1t%C>aDu@FkC@?t=gV-x2?iIOU~r)c23MQ}aKQ-%7o1>l!3hQz zJV=HMPB6IO1cM7sFu35^FxtgO}%haC*0HrH+8{H zJ#bS8+|+(IHQ!CGcT?jFsqJoRx|TB^DavrdAhHqutczLTYj$wb)G! zcEf_$S=%N!yC6Sz+WNykx@1@#TR(O~x$4-$dH9bE&cXkSKS%NB82C5#=Jk)4#;v}sXcz_{eJk>kb`6^0~DOv+41&PvM4O-f2>-m<7L(3rLUf0Zp{ z<-zsFwJj?2x6};GGcwaXNeS7|sQqLcX#am(Yw@tU)kDPBWhIPFO3WFTR8&|YIW0Rg z{gLLI$m4@z3bK<%f=OmpQDMcb?2*~YS=q@6AXu+}b}fKbgR+w*W{*lwPl4M18P_B4 zgg#eH7?+a;o=8pu57@nKUNo;Tba=+d?6JeM(}yRfC!~!`O)APSl5+-!v)l5$!BzQ9 zhdkasf_qNM1B0s$`P=c>V}&8MQ=2zV&X|-os>pKz)m4nTnXP3N9?|-3#iQC5w(V8f zX-H(Fe;#e8m0fr3QEh8mdsa>x5?Sr?qjhcVx{Gq%kUH#^{P&P*4i4{=O`m_9pL58a z9cpmT8CmDmTKusJvT4T}1)lTRKp@RRRMNxO zQ5E`d#%$BOWNJ?;nchKFWteC)eWZLqN{S=DXYM`JwY1y2@=)jZ(tXh7z-*dX#0kZT zh1smO6ODwP^wdE@|9R>HanJp-0YfW_0`!p2 zIk#qF%-)UxF(1AW@ul?yHT`M55|MUzJ(UU&sq36ewgx>Ty&sEJN2|em_9=T*m z=vANdkEbS39eUR1WcikKj_Ts?TM|0h=iECqmn_lAKIhaDOFED90KNIm64Bc}XT4?< z$rL^AbH4J}B$CkUKIhFBCzFJp_c=QcvLy7r&pCSL6sp4+@Hqp%vLu`ZpYxMTmV`6m zb5?IyK$b8YRYAP+`7{zSkbTbjN8gCqE3c01#?Hx(qekO;eQs2={Bl%Rb_EpCvO&VX zAw}IvXv}WPdWq5Ol#GMlzvWo?y&%1b(VSI~A12mjLVgRXtMWAb{v!jEqWLw4te;ey zot1GQTz1H05H3n@QZ%HD4?x8$$UUGsS9Ut78)IkYdnU>?V?M(!Fi);LG4`{%7@GM- zo^jwtTwiHW_J|jrbSUE;Ia!_?a#wORJ1vjH@1>2NR*YYF$Tlfa?5-S;63y<&3Gn+s z7Q^p7xdnc&;7fXr`;6U}!Ku%{z}lsrfPo##dqO^!HdO}Yy5(WSt4!Wry)@IcGwi95vZ3 zn@kLo9dmoi;w)lye4(f!_vdz&cQR2ioSmm0gPHoP9j_t}P3$aVa@@+nSuww}T%F;Tsrhqd(G<6go?0Z6)7|RA646nH7xboe?5Pz1%bPM?Yb45aFCtN* z>x)E*ZX^;Vx^yH;booe>=w>5PqFabWiEcF#CA!T>l<2-dqC|Hn0P`7Ty5CSqneGA- zu0_usBuaFFP9REjkw}#2o~JMfs5u^q#8Hj!DU+(mzjzwOV*o;dDVEOvaZW!n~Fb{n`aeO6t9@Z z)ZG&-9oA6`x$XTpm2ue~Y|WMJ7>cg4F{@Gxom>qsB;A`IHscoUS?e{rgmtqF-CyrB zw5W8Fp(i({8hXFXw(GRH$mssu^17kit&n&oL6pia zADwFG^}(;Lbb|3Da=<|gJB{L=~sd78V*e!c#N zp{$$GW6}lO#!1E6H}CvvDEnSFm^RcMcu;8RbSh&Q5unsTf_-i9goEMyqjqKhGi~Y` zfoy@%ejdbN^ki_M5*n2B=25D1E72Dd{4{lNFb~tPG+Q%7m;NM#fqw~K3S%D{`1h)*XR%L^Xc4&3>zQOPsY^=dKHJRWlx|p|;2d}loO(a8IF1Dssduy|P3^q}H zEE*aIP_tKy=6YXUHkvV}Vkfb1^<0*ys$X-nQDnF)hAk!-`V{+^AhaIaMUVcpU3Ase z8?Y3Tx=j;}^u~rPiKLEii>GzxMy4*-Eg~qa^}fbz7PT_42|Gm4pc&gqaIQJ~mSA&B z1`(RP4%YYZ@paVo)+|mn*=82F9@vIKIs@5CXZoYIY?vkMTJ2fBsltV%VIeg$RaZN* z2L>IT8N@I)po?M|1RY_WE=(Hy>IFY~y4o%3TQ9M_mVunIq>?_~4f=0A`b9iDVWC>y zi&fH}_h9=?708VRV_b#)}Qqx*`;d_ zWa$Q%4Z@HBE0^{R@ng%^*mS1C1EA%}nXT2(VbJp30RFiCY#3_;i8279%{hX-W-xms zgH?rEbzn4vLd68^7+Ez=Vvw;2p$gIo;qb8x(l23vhr!w){8tKkiSVUVh!GX?ICP+P zR%2C=4jm|c(jMCJboMN)*e9XWb;q-8gR3*K7o=^Q&H5Nzl7q{Iv`unZC@gt=tj{Et zXtdu>W|(47J0XwNGX|wo8KxH0wtL-gaC#bpyieL$GuT>#UEgG_&0}+CvKdAjGYg{I zvOVQv)r~>e9DmPp-e#PudMPl*PqN0Vf$zW=N92NMOW$E#8KO%EN(?8$q9O*Fm(W#= ztDA7SV$X03k1d3iVqyQq>?xz2vjnG-Z0jzAWoF^kjrmzg$uD%Lmq7s!296IHc3t0YSLQ$Z&`bWVwld^!{!(a z{|Zfe?EX&ld$(8jC{f?hnb6882JM_kF;xkWJ8U%`cVdZDQWe8 z*cyYG|Mg>|V{Dnx{`wn(kCt2NohX)u=Wiu@mKdkk`~%hnh%WuwNtQ#TUN|hm_4Pkt zG+3bq{4A>JA*a|C%Wl{i_LfE0&%$Oy&n`L7R$DaSB3o`z;3XIr)lR?6f;rK9R~e)= zeA|s{5Sm0e*V#>r3U0C~7FE7&)51F}!_w8f$L3qK^FB*vEMP%>J`h{jJfAPIDq=Z* z!FrK$aEP)ots#5Qnb1sa3gB_-xi;)^9qHtd;;M9kygN9M!!i$?0*lG|{HUVFRs_3{ zJRYNuROF2yo8-Y{x+?Qt7OGhtxLYp>=ADe(sy45y148&=V|pf(!(`J_Z@4(D2Euj` z9L@-Y7b7{$FyY*490DB6&J%MQzwsy>!2FEJp<9D2h`-R5h-5> z{)GOk7W9BBTG!^m+%l=dw_7^*6Z~U~{&MrZ7VUnL?Qan8PxmDmX=NG$wKcmrh5oyG|x+Vko#yE|PPhCaih=;!!5 z7VYlD+gVg4mS46gtP781p$RXVzs1&R{!E>mtR&s$1&%Sg(-oTKw*OAssCgr-5?xKp z*3ozA35fp}6wT=}Iqy?tTxuFYFbrnCkEi zOUF;u*H!4ZEL?To0a3NK53dTtKV+RCLe!CyB1~WH%U`gRx?6t^yA;;z0|PiD0K%j} zJjFb9<}i=d?jgL3kqeGM*uMKJKVVGrhuUqar>^lQ_03_tyD?q;4{xXoM(|9=Zs;o` z`6wt={{(&x#MngsJc>|8s)b4XN1NXq&3`tmv&O)=5R7Mz1&2U+Dw)SJHMuVY(1A(a z)bSJupe6n7!&2-x4u>U>T^hdqATbg~rK4-+3}O*@*bCB!GI*R-wYp&@>~dt(F&hHa z!XI)VNG#0EAL}_!f)4 zo6gr*H1iGql|{qe;@?`-#>*F3baED7VNuC!SWnbS_8eZHBdLrHaK(UM4~$06iz0%%2lM|U#qlWX7&{L6UBUW$z}WQpArtI zCy-sb(;{q!@bF^(zCrI2c(tXCScch$3`&=yp9y=d;2#(adLKqYTHNJ=^w%pn?t0jG z{Ax%fYIHkDhH3Cenlm5543>5>eboKe@H(*Zb%Nu?%(dW23u~?8Q!G@Fx-02T>p32w zz)IcwhS${1OCd@u(?K8eKTQ>VUD6U8T?UKSvL_|2r#+i^Kl9MDTX4Wrmxw5}{9D!s z|4=hrWo(4=0;R_lTQT{7?9w&1@nELrMM3*&*HXhQpy}{ux8-z!0@>K>4Z~KP_sUTd(-Q@#JRp~vypBp^j$(DBWO3?$pQPKGh zA>iAH@RUQqxtb9EiUd0t;_CK+@QF&w9uOX+%Ydwe^*b~gNLOTWch(QsUJKe zDm)^qxwm0QYxJ!BxMkF`X*nbPN^>!ssb3VNkLaoWRAx&Nr*u`+1RPkiQ@x zB1oJ6qJX7|sD{Uhiu!(60mnJ20vSU0Vbzgt0wR;D`p1iz1~0xWVB;Yxe60rm!vZI( zpjSkpsRD<$*2jALI|E*S*I3=|C*aGJ8lcgw8z5k9Py;Uy5)kT$a%OQy9Xmt}F*WO6 z6_D4-sO?aUKElt3`L#Vp*geEEokO1;DIiCYRceBmYoQ!b(m{7m5|b=ZU5sXJbnqAf zc^WO%%J;<+T8R=2U z4?;J=(da_1fWJjheUC{3QU&$DPCeE@P23Ij?kOVHDA#O-L6zi*fktk*O|;Q%r}`T` zwOzcVZx`5|)`O=B_@f!Nw`GQZ($>EriY=uc@RsdPePgD89g|)SAEl8hum`+)#VmU$ zy25M$>6m)CzR*9*J>C|nM!RpW7|$((Uh@Tf3s74NJ{7I>(?xbOs_{=N~KNA;*%+v`6|Y=0h5JV zg~6a97~ot>RuBge;)D6ShYaaJcI^{ArhuM=F`iFp{2DP~6 zr1-#K;$H$z9UFpR84kM08tPuB?JntKXT%a?HRl{=DC#H#_>;Qb1zcPtQ>LaZtMnqS zV0uC5B~0r6dJVgmnND5ns(|wqS-pGB@05t^nB43Ks_FGNaBbMKL&2sMu4~`I1>#p~ z*tuX`Q&X+(ij@X0-1B>5^aFbeH=lFZ{!ieJk;e2#;jsHMpn}6=w1)#6CKh#wipxnQ<)lUBq<70nOUp^i%SkKBNgEzYYGst8 z5qx=nALVGuDZqvnQGh{#*V0hnoir49Aq@rIMni#D(NN$$G!%FV4F%pnLxBMX1qK)t z7+_FffI)#b&QRcWGZYwLP+)*TK>-FO1_2ZpU{GLyL4g4V1>P$|fdK{u1{f3=U{K&y zG8FgKo;r^1G!`!n#Y-dc(m=d44lfPEOQZ18AiOjNFAc#B<9_fqF) zQQy7PbuZnId#T^EsM}u8r3TO$FZFsBb=pgPo<&`rMLqUXhrO^McGPtU&d#Y_ZpQ%l nCk|CT#?hZ$(P=S`NFI5Z!L|4i{P`Y#4&l#1UGkJ8Dd7JADvARn diff --git a/doc/LectureNotes/_build/.doctrees/exercisesweek37.doctree b/doc/LectureNotes/_build/.doctrees/exercisesweek37.doctree index 6eb58b2d1d73eebea855d5adb1c54172e0bbefd2..bd54489041ae98c8ad73221d8c9004a89743b6d9 100644 GIT binary patch delta 1428 zcmb_cU2GIp7~MO|{Q`e;ujJQysV(6LxrEO++?sjJA&aAsL z-5*=~V=&q6A|c{Rrhz3IhJUrcGS+G-E=?^XH=*6pdT9-JBtcl{SD$d7G;H+WyGyW(kW1Koucc%nR4o~HZK zsMG^WZ1?Ye$mx(eq)yuKx*ZN_tK>AIl;69WWiRYF)M|b__7Orms3ZVqgtaJWzAyMV zZomz9`&-Zed>xD0;IJ6C!A()Lfk+;Ba+@4yPXdRw)kU6MO%XK9@k@xCsW{KTxYUVe z;I#A<^eNqNPg+ZnTaGfqEZ+pL$tHT2&7^MF_Lu=(i7*EBhWm<-#&sXRh`2>J{Fr?V zewb)SpHfXfT7;d65^J+rEW+)?E(VxySlx)euX5f|`&Qb&sxfq@+CHg?tj+R|+V2j+ z6Uk}L<$>#MU9hWp4J;*nT+j_8L;dFCFJD@NngC5OomzP%pBrE9xGo6-plr(5Zuh}=k_uGZE!5-qYWwdtj+l%x5oxO`IrqZ z=bvGKkZpqxr*{M)F2AN!KoZqw6ouzves^TE8^Z-?bemA*=kdJ^;^W9uOvJV z8#ua*=P| zPgVcQAAy|O$mco_-KY4=caBXt~QufWa-S++f6SXfaAPsz84+xg>on z{kD`_@Ye$GXX=8bELf_7r6?Fu6XPtw!O~E>W!|!CevOSR^FP^$nJYo~tn>$Jfh#97 Fe*;li=q>;N delta 1065 zcmZuvT}V@57~Xd_J?WXLX#Q1mrO1{VF8|F)11&QNNlLJqbDnjaZpR&4&gm>S%SJ~l z-*q7rUP=&I9_c#45JO143F@LBh(JlRAH0*k?`UpPH_!Jz&-=X3_nwdQVw92|Qeij3 z9_WUihmAxcCyYiqWjrFCWGeMBt_l|`m|HqnZElv<*}=kSo)BaHSFDpd^C#n?1$N7x zotsVErB!^)wZB?0V!!mRY9<0bl5SrpLNx0#+1lG!cwPTRQl9Te|JOZ8vNSwV!i=#3 zhK&Z8HlD&~Hf%Jez;|8)WhpAD#59SJ9zOBeu+JQgGNddsLs5m{piS1GLa1k(g5*7G z36dsDbC7(sgmmmzT1jWsmjevP;K)R70#UZjNEtH1Df4xYuXF82Kja~eB_T|7UO6pBRF1tIZvtJbex_`ux7x7!-IaxE8>Szn?hQW zO;}Y?G!2f1av1HTDHRRck_!0N$xw)tiQq3WzzJ6+S=QBukPd3SsTece1(c44XM8gB zxD8Yh4GU>G@WIXC-m<%V+vo7cgQ+`r2cPXeyMx#C=rLdM8PtL|WbMhGHxx%!ONEQE z2ck}UABQv-JzdBJIf4N^z4A-;W+HTp2nlRgi;QYPo+1~~dsfQuECMYz5h=nHrP z_|Z$@<(&dl0_%4)=qMzJ2hkaZQv935R z9G6-GzrrxA5!g#;5`h_FoW|f3BkU#{d$p^b?Z9{HtP9XPn4e0f{qm$=p7P5RekYeq jJae!(crINI*KCK+lS&R>Clv|DKTXxpZ={CfcRT+8qSS06 diff --git a/doc/LectureNotes/_build/html/_sources/exercisesweek37.ipynb b/doc/LectureNotes/_build/html/_sources/exercisesweek37.ipynb index e6ebd7459..25296c4e0 100644 --- a/doc/LectureNotes/_build/html/_sources/exercisesweek37.ipynb +++ b/doc/LectureNotes/_build/html/_sources/exercisesweek37.ipynb @@ -2,7 +2,7 @@ "cells": [ { "cell_type": "markdown", - "id": "1b638450", + "id": "1b941c35", "metadata": { "editable": true }, @@ -14,7 +14,7 @@ }, { "cell_type": "markdown", - "id": "60060788", + "id": "dc05b096", "metadata": { "editable": true }, @@ -27,7 +27,7 @@ }, { "cell_type": "markdown", - "id": "7cdd88e4", + "id": "2cf07405", "metadata": { "editable": true }, @@ -37,7 +37,7 @@ "After having completed these exercises you will have:\n", "1. Your own code for the implementation of the simplest gradient descent approach applied to ordinary least squares (OLS) and Ridge regression\n", "\n", - "2. Be able to compare the analytical expressions for OLS and Rudge regression with the gradient descent approach\n", + "2. Be able to compare the analytical expressions for OLS and Ridge regression with the gradient descent approach\n", "\n", "3. Explore the role of the learning rate in the gradient descent approach and the hyperparameter $\\lambda$ in Ridge regression\n", "\n", @@ -46,7 +46,7 @@ }, { "cell_type": "markdown", - "id": "0328fa2a", + "id": "3c139edb", "metadata": { "editable": true }, @@ -58,7 +58,7 @@ }, { "cell_type": "markdown", - "id": "ac760265", + "id": "aad4cfac", "metadata": { "editable": true }, @@ -70,7 +70,7 @@ }, { "cell_type": "markdown", - "id": "8d4b0753", + "id": "6682282f", "metadata": { "editable": true }, @@ -83,7 +83,7 @@ }, { "cell_type": "markdown", - "id": "4517d311", + "id": "89e2f4c4", "metadata": { "editable": true }, @@ -99,7 +99,7 @@ }, { "cell_type": "markdown", - "id": "da1834a9", + "id": "b06d4e53", "metadata": { "editable": true }, @@ -120,7 +120,7 @@ { "cell_type": "code", "execution_count": 1, - "id": "590b2fb0", + "id": "63796480", "metadata": { "collapsed": false, "editable": true @@ -140,7 +140,7 @@ }, { "cell_type": "markdown", - "id": "24dd92fc", + "id": "80748600", "metadata": { "editable": true }, @@ -156,7 +156,7 @@ }, { "cell_type": "markdown", - "id": "4ba80e16", + "id": "92751e5f", "metadata": { "editable": true }, @@ -168,18 +168,18 @@ }, { "cell_type": "markdown", - "id": "be65f56f", + "id": "aedfbd7a", "metadata": { "editable": true }, "source": [ - "## Exercise 3, use the analytical formulae for OLS and Ridge regression to find the optimal paramters $\\boldsymbol{\\theta}$" + "## Exercise 3, using the analytical formulae for OLS and Ridge regression to find the optimal paramters $\\boldsymbol{\\theta}$" ] }, { "cell_type": "code", "execution_count": 2, - "id": "c265677e", + "id": "5d1288fa", "metadata": { "collapsed": false, "editable": true @@ -200,7 +200,7 @@ }, { "cell_type": "markdown", - "id": "b989efb9", + "id": "628f5e89", "metadata": { "editable": true }, @@ -214,7 +214,7 @@ }, { "cell_type": "markdown", - "id": "17b08af7", + "id": "f115ba4e", "metadata": { "editable": true }, @@ -226,7 +226,7 @@ }, { "cell_type": "markdown", - "id": "6acd4708", + "id": "a9b5189c", "metadata": { "editable": true }, @@ -238,7 +238,7 @@ }, { "cell_type": "markdown", - "id": "50579422", + "id": "a3969ff6", "metadata": { "editable": true }, @@ -258,7 +258,7 @@ { "cell_type": "code", "execution_count": 3, - "id": "c57cc917", + "id": "34d87303", "metadata": { "collapsed": false, "editable": true @@ -277,13 +277,14 @@ "cost_history = np.zeros(num_iters)\n", "\n", "# Gradient descent loop\n", - "m = n_samples # number of examples\n", + "m = n_samples # number of data points\n", "for t in range(num_iters):\n", " # Compute prediction error\n", " error = X_norm.dot(theta) - y_centered \n", " # Compute cost for OLS and Ridge (MSE + regularization for Ridge) for monitoring\n", " cost_OLS = ?\n", " cost_Ridge = ?\n", + " # You could add a history for both methods (optional)\n", " cost_history[t] = ?\n", " # Compute gradients for OSL and Ridge\n", " grad_OLS = ?\n", @@ -301,7 +302,7 @@ }, { "cell_type": "markdown", - "id": "e2654903", + "id": "989f70bb", "metadata": { "editable": true }, @@ -313,19 +314,19 @@ }, { "cell_type": "markdown", - "id": "7e7e7de6", + "id": "370b2dad", "metadata": { "editable": true }, "source": [ "### 4b)\n", "\n", - "Try to add a stopping parameter as function of the number iterations. How would you define a stopping criterion?" + "Try to add a stopping parameter as function of the number iterations and the difference between the new and old $\\theta$ values. How would you define a stopping criterion?" ] }, { "cell_type": "markdown", - "id": "75542708", + "id": "ef197cd7", "metadata": { "editable": true }, @@ -351,7 +352,7 @@ { "cell_type": "code", "execution_count": 4, - "id": "06077986", + "id": "4ccc2f65", "metadata": { "collapsed": false, "editable": true @@ -380,7 +381,7 @@ }, { "cell_type": "markdown", - "id": "86d46505", + "id": "00e279ef", "metadata": { "editable": true }, @@ -394,7 +395,7 @@ }, { "cell_type": "markdown", - "id": "f1d2f848", + "id": "c910b3f4", "metadata": { "editable": true }, @@ -406,7 +407,7 @@ }, { "cell_type": "markdown", - "id": "f1c91e84", + "id": "89e6e040", "metadata": { "editable": true }, diff --git a/doc/LectureNotes/_build/html/exercisesweek37.html b/doc/LectureNotes/_build/html/exercisesweek37.html index 498d581ce..a366a3c71 100644 --- a/doc/LectureNotes/_build/html/exercisesweek37.html +++ b/doc/LectureNotes/_build/html/exercisesweek37.html @@ -389,7 +389,7 @@ document.write(`
  • Exercise 2, calculate the gradients
  • -
  • Exercise 3, use the analytical formulae for OLS and Ridge regression to find the optimal paramters \(\boldsymbol{\theta}\)
  • Exercise 2, calculate the gradients
  • -
  • Exercise 3, use the analytical formulae for OLS and Ridge regression to find the optimal paramters \(\boldsymbol{\theta}\)