From 12971932c54513ceddcf515b2577608714ed1ab9 Mon Sep 17 00:00:00 2001 From: Jan Grewe Date: Tue, 6 Feb 2018 23:20:55 +0100 Subject: [PATCH] [plotting] add code and figure to creta errorbar plots --- plotting/code/errorbarplot.m | 43 ++++++++++++++++++++++++++ plotting/lecture/images/errorbars.pdf | Bin 0 -> 6558 bytes 2 files changed, 43 insertions(+) create mode 100644 plotting/code/errorbarplot.m create mode 100644 plotting/lecture/images/errorbars.pdf diff --git a/plotting/code/errorbarplot.m b/plotting/code/errorbarplot.m new file mode 100644 index 0000000..6910d1e --- /dev/null +++ b/plotting/code/errorbarplot.m @@ -0,0 +1,43 @@ +% fake some data +x_data = 0:2*pi/10:2*pi; +sym_data = sin(repmat(x_data, 10,1)) * 0.75 + (randn(10, length(x_data)) .* cos(repmat(x_data, 10,1))); +asym_data = sin(repmat(x_data, 10,1)) * 0.75 + (randn(10, length(x_data)) .* cos(repmat(x_data, 10,1)) + 0.5); + +% get some data characteristics +avg_sym = mean(sym_data, 1); +err_sym = std(sym_data, [], 1); +avg_asym = median(asym_data, 1); +err_upper_asym = prctile(asym_data, 75, 1); +err_lower_asym = prctile(asym_data, 25, 1); + +fig = figure(); +set(fig, 'paperunits', 'centimeters', 'papersize', [15 6.5], ... + 'paperposition', [0.0 0.0 15, 6.5], 'color', 'white') +subplot(1,3,1) +errorbar(x_data, avg_sym, err_sym, 'marker', 'o') +xlim([-.5, 6.5]) +xlabel('x-data') +ylabel('y-data') +box('off') + +subplot(1,3,2) +errorbar(x_data, avg_asym, err_lower_asym, err_upper_asym, 'marker', 'o') +xlim([-.5, 6.5]) +yticklabels([]) +xlabel('x-data') +box('off') + +subplot(1,3,3) +hold on +p = fill(cat(2, x_data, fliplr(x_data)), ... + cat(2, avg_sym - err_sym, fliplr(avg_sym + err_sym)), ... + 'b'); +p.FaceAlpha = 0.125; +p.EdgeColor = 'w'; +xlim([-.5, 6.5]) +yticklabels([]) +box('off') +xlabel('x-data') +plot(x_data, avg_sym, 'b', 'linewidth', 1.) + +saveas(fig, '../lecture/images/errorbars', 'pdf') diff --git a/plotting/lecture/images/errorbars.pdf b/plotting/lecture/images/errorbars.pdf new file mode 100644 index 0000000000000000000000000000000000000000..49ade7a70d4faa46985f50bc220c17182ecd0909 GIT binary patch literal 6558 zcmbVRc|26>|8F5TywUZ?oQ3ho>V{IlggODX#B_dLiN-B~yYeO(8veV5zcEb0Jdv)F1*X#FtnSbVdw&$}v@AJIod5?&Jnuat=1`8K?mYk9bM*&EH zYUKbY5a0*{8r8xjXvBSv-| zF>>pOkzYp)ejPC=S`Qph>%jwRJ%B*12NCG? zAOgJ}M4Io1fKC-U13G}t!7s!P znJb4$=Kc0+Bz^Sx5ker3$g!8acAMG09jM*x0i(epMAbW`!?S;H@o5k1ZtK4>G-m2; z$mEvi%$%yg>{mb0^M`=*ME$N2*wU-f*84CiGZk7t*(;9zbq|10$A zcRxDC7W-$#rTYhO7kgJO?DCo!_+x2tqA}5D)3oNrrHlBfx#|Z?x$~WWR+GZt#@?fs zB%L1h5z?Jwy8Fxw?%eNvgExu#QvT*GryWCoB&HXaJ|B6|HOzH(akHNL^aHVkwg-6+ ziv_R73+WP5Opl($UA8)VjGSL#*JG*@8XF(%D9H;H7G47aguid1xr5-s2ke zryyBC_?1_9Tpo}irzdbF{+nn?;Hf%p=PR=tX4~D`edb4~EyK0v6*L=VKm8zWecFoN zQA<0|fqlqb*RhQwvEu}uz3hg7)=3roP*Z2qR_v23Wzx3Vk_eiT&1Gw%hSa@nrwSU! z&=>Pg^hWG_sg)fU6@_UXZH`S-%$orXi@U=;j~H{QYj09XJb_ZW9FO~?x$%I?78%se zC`OUa#0~S0Dy5^9+*jxVlhZ{Y{ogbAg zid^&6GJ_pV9*%8l>3tt~;BoK%%0jn+5*HrwLS(d|?*{VgF>r!MSolmZ`O3pWyIiuD2ygRSaBD#pc1yaX%Upe}Zrjs?T|l7sx@ON}{e+~S3|QtnCiI(%?d-@r-mcP8ZN zb)-J@Q&(xYL#~GT1(v(ZBgMO&s%H{(zGqYXPx7S!AIgcnNm{&IaeRCE&CVXXJ7a(jd;6TO`H;FQb4Igby4WaNS{*q&SEq@gTl$7RS{^&dtbaly3}`< zs0sMMshEb`u(aiQhL$Wbg5>qIMg9>&dt1SJ+IszkaV<8ozEu$3aY?u2YqEYXPAh9XwAH|W_~-9G%qce*?#}?yd=@d=grq?6DMl%Iga@uWp~ne?LD%d)?RaD+F0&b zoGT%Y`?&bRyE^8FtkcNF-W~78x)t3QnwftVCQgoQ3#HMf#&~{wIhu~hn``6P*u^AL zNLHR!1dDNL9H_34qW0{yqMfY4W40%~mtz_nKEfO&1b4k_6`ApMU_bM*EXEWkkWy*2 zQv$plzEMKV*yCx^Nn&Mi?)PZ~e*@cb?YEP?7K5sK`zihEL5j!Ts5Z;&@#E7Kx9!=p z9hdq!0CxOp%1t9z+yUuoSQGw8_9sFr`>gVo*UfOHvP+WmzL|_CQ}&qW3IaO9iuGOg zJELQcG?iSvB@x0&PwD&S<{Ozv6>{vmeE&-hm-dbp+J~+>n`?=?Uw0U8SIW4T9Nh=V zn+2QoY(Bcb!M7;#9>R2dX{6n(>ofakpFH-xIcHFpPQVG^=^2Lj^h7GxS*P}x6u|cP z9eg%ve6{oUU%71DI@}UdWg8V|iN7@uR}W5w#Mwkjfv5TlxZh4uY7Op`tB4ej8_Prldnz=tBXg;a7L`58 zS!~3`jJWe^jUH)73v3-_V0hk-LQ#xU6T<=gt;IoQ&(k*sS+RdT7o9g4z~53FB*^J$ z^>h*h8_$2Ad^;(IMv%{xr&+xTY6kgA)cSVU`kj&z1)DNlHtO0~phyO?{BVcIdA-2tSD&iJvNr|6Tc4hd`$z$7QdK1_#% zd4-)L`;EMUUQzD)QAD1)9Usyi8%PjuN(XQL_WhVdm?mEF9cHLlt$FmJDIQ;C7S3z=7y=VQpApx|){6Ok#qR`N`~s zZ@i}@4_+QV#(Bpoo2?I56A_8YcT$o;yextFqNQED2M6D^N@DefIU*ZO zyK0U_^9#JY{K-LAQCG0b^wyq{x^u)yu_bMBeRK5wH`S(&wdIP%V3(`uM9IrA={~0# z;Y0;`>BaJFcJXm=_7x9dOR7L+rZl-e8^zQ~T~xT5RG_RcDNes>w=bsC-bfsAD2HOj z^wisNLzk@Y{Jt|Bm6~-@@xl?b+WVciKY?w$;xAK*9!P}E!7fQk2qF6PRhCl2rL;#1 zl7gZE-FpKY5SKSa0WfwI`al3Bm~E45zav{eFYuKg_VRnc;bSohkEkJ=UzfIX-24%` zvDjNETnK6LdNT6L=5P&Ej>Ld&R=)H$z*1B-AYgn!xNC!Z(#sp8X6}@b%~v8&y*HUg z4eoUQjICT^x^vi*v3nl*BEsL51QTvYoVl{;w%53K%gN!qy^Y>2*9P+vC%oz2L+PjH zy<3=5CN*b^kBG_xbOFSYu^YwN)?|89*0pApVD`hs5MEIzYMabfN;v=y`zkPW2 zNyx?I&}nPwq)y1Y|E6w&yEN z=Ov3BhjF(@q))3HG`umwJ2gOK^b1)6OoDeJ?5t)E?8l2_kqaB%d>)uPm>I@ou9wtvJ~BZ&D~wdU+vNL?)T;} zjQgn+p;+S6u@`oYEYCMFMFr*U*D*3*Tn^OqSXSaBBbe1y(_@4|>3}Zxz4(UoPv7iKdvzVzn#dJ(pJvR3x5-&P`)+1drPXrM2*VN`ENJ(+g8olG5WjCADc2G{l_ILg&nxh-|ZPtQHo>z zf}PLNsm-r~7aft`QPiPZ`NBZy7UzATAKwZTe8j_xx0$L>OdZP}E>JLQ#JC$ydFs}G z|6UU3FiOixp(*ORneVmo-u8XcvwA+q=W0O#%$CneV^g>1pBLdSE{;p%rGJDR^z>AgkXoMUIg;ACEhG-)B{Z1?oswd(jR)lgG^+!Kx34;Hy|jwR&5$^KL6 zqy^kr;?$hyb7uXS6(X0-xCF4gM(tJujO5Hg^g8v}^*$Mc_l^0|zba2N6dCzIDEYc{!Ee~Mur^QEZRqg>nkOS7@^qZg6Vc^*BLd_BztOuFPikK3?Qri@s>cnv^&IbP;qxz+(CJCR~W|IzSh9ON$ zT;t6~k@xJ|i;{xm%&h7gDD_PuAK4LJSz6Lgo)(2f>)TprdmZ?_17*ttLW)UY|z3UMeIOx_Ol7g*i)9T;(bZd)LY*JKVj+ zq-(;fr+*4g(0x?x;&$7+DCXkk|tA^%3 z4XX$P06-XmbgB!@8l(eg7Aar|hjf2YBLxPX*;{cRA)|_Qg zrBToSr$NOD1a1D)Vnnk)_n!s>OBzJl0caM$N)$PhzA7_+UKAY zDt~{gLJfcy7|Z+ilgb1jf)$F>Bti?V%w3a+2A0W;zbqKwukkOtVS(j7|FYnKpMcRM z$^t(d(Im}KXoYVf9UYgxPI|NW1&I)qQfAe`TVNGLEH0-4vSoy z1B+Xm1BYCj11GyShb(K-ozB}3!G8KTF|9V-p)c_I!FDr{xRlz~);<0FT z6&zMx4T;r2qi}c)RaG@SN(uPCU6#)SqP+*m0yYG8I>VC2ScW=I4v&M2h^Qaefd3DN C=Lbpv literal 0 HcmV?d00001