From cb2bd448f4c28e1e752238bc01d8a8d73a78c5a9 Mon Sep 17 00:00:00 2001 From: Jakub Adamski Date: Mon, 5 Apr 2021 14:03:14 +0200 Subject: [PATCH] zajecia4-cd --- zajecia4/.RData | Bin 0 -> 4172 bytes zajecia4/.Rhistory | 109 +++++++++++++++++++++++++++++++++++++++++++++ 2 files changed, 109 insertions(+) create mode 100644 zajecia4/.RData create mode 100644 zajecia4/.Rhistory diff --git a/zajecia4/.RData b/zajecia4/.RData new file mode 100644 index 0000000000000000000000000000000000000000..739f2085bb8f94046df80cc82d6572f10fbcb8cd GIT binary patch literal 4172 zcmZA1`9Bi?;pi85|lZD)UCF8ti*;+vAVGrY4;+2mfiZ*gG% zek}ncn)Hgre?{b}nZx@`hAbyYB%T^3tPx7#y41U9^*wCv#sJ2+3UI782VpdXF+4+I zkv7NvZF!)sdp)fQ5{-Lw%Gv{qt!x;}j!Y<$-b7#8+T`b|qC$_aa1DeH!CMj;Hd}cY zjva#<7Ixd92##Lf@>dV|vLzYSq5g+tg!uaA`r2iA{@$&Cxapv#*rp~n6qbesu%dpc z+FZZ_$twd1qDj*2tmkBxMVh%t7e#lNE`hF6nAkx1J4`@OK7|PYa)^4E>2^j;^xb0_ z(atwaVz=M9$$aj7b6U(TDPKJbzzP)qA6y1Hxg`~-#{f>?Nxpg{fCb3ob`rcmGq)!N z>d}A`7|vIJ1vn3scRL9o;0w3m0`*wH3CQHPM*x_C9BwDU4K#3*DQJHU_;2J!UI$(- zNP44jiLLKuVSZAa29WK~O{e^e?vDz&lC>WJx&8+lKvsA0!sI0F6L2a_PS(~2UUokT zbs&p7PoYVY_6e93nj~v$0eRg|LKAr2UA}N6N&5s|6pkco>i~t^PoU=^klCH1P&Vl# z3<_nFdmaI~-Tw;(A*}9XyQub1-k|)8DysW={|Gv76~^UT6>53a<4>ap=AskOdUjKg zJXv-@DCa-PoV@wb>CuJz!X$o{qfhzvr$a4S&`P}@&d%doALJ_Z?@rC+UynF?#c{jCZMZZC4sNd}Q*&@dq7j)w=_05;mO^zb1Z6*e;;0 zgz$tg=i&7`1_&wWmAaaIHP$Mlg@G#Gh%Z$#LSRMV1C4w#Jm32Y#MOCN8o;UkV)0R1 zfW?_3Ea4C`mXG6g5aRP467&Qpspt_rwc)J$stf8lpFR@4`+6fuMY8N?Yl_+r*LQII zn<_UUnF3Hu6tq&h1e&}Sp?aR&_*{zo(lO^0#bWgc@bK7l!6PCdc@0}0-bNR14;#Gi zaP$y~C8#Xyem=K8X}nUAB@97mY|wpXryxE;gF8LUXw}SCl;r+l3QNUJ0n6ABWOdKO zR=E=6==Ym533mMO_n^1RQp7pvUI3InVR@-7F(T<`05BDhx#9rSK>=RPGDx& zO7mWi9DKi@#rgf)&iGWtzhT2q+OW)IRjuIJm5(|+Rv#S#aJ$+_YH7Aaz9aZMU54<$ zX3K84>kMO9HlrgXKIN>k9J^Gbd(a_Eu&)NOldh?rxU4>O>iwF!H~&4WcL??sVmBji z(}98H<_n;a_s)dXG3lDd%-S4*ndls^RKeiJT~>b?0gXKyI*iUr6LzIV4AQWIR(s2N zt{BhPhI@wcUzjL2g56C``|V++}m||0$(z| zTo0b|e@L}{g`hfqEPwZx??!)zl96JOlq5A0-`|Oe3zqT>xJNqBR~ZS>8O_}@rt&Sn zCDByY@BaK1d~2%M#FDHSuGqH7ZA`ep5RDJv%^MD^zf<3aB0!30(8Z%^A*0>I&{w>v zWxVveX7e;aj;EdEo;jk(_N`+)Ir{Y1Y*5aq`$Re%dcjwT!Q-y{anB;{kOzd-qO$b!>N(wFH%VLA{&WVP|LeN?B}Xnb+8W!9uctA35d<5uWpm8t8&i=CSD z+zhKxlcL7P(%Dk@Jc*Z}bHB#d*KI=Lz25q9s00h_Z&CdHo60k#NXJG@>q{1qYWZNk zEY$Aj-h=(#-%|<-ajnZuMlK?r1>`Jwy7LSz?Pp*|>;VUA={WI+`S(jC%gieooWh3h zGx{O-wk({#IBAtsohIHUnOKsnox62@26X6;&)rl;j`k9-s(IT((PKcKC%&rW;b3LL zgt0Q3ZF4WD^XYB2U#d6Z72PNiLGR6s<1Lav#lpwTkdh$3)+ouBd*d||Yoo|ye!V;4 zhnKvw5oojVKd;u9;e#<>OtCr=4XF|!e&VEqskDytV-rY=in>={R*c>LK$NZH3*Q*J zd7pqjbzcxl(<|Q?A_cT8UrqJhTfWSBYy5jn$TY;u?qA`STmKJAB|PWxm~vQWlgf)f z&76J(ng=MF2@wsJhQG1SAg_M{Be%NT5IdcNxq_9c?-DQn!+B%#Gk@+5uh7Jo=BdCET8w1G4ChnRh7qv42B^kIXuGp}(L%-<2TjWaZd$`;sX}cXSW7 zb$(>_O)ryZV_ZXU{>4FyIoStJ^;amulTX6_P}OE#v>R3<*zkuR`E}P5YC6dqal&5R zj`&P(tx{WpOM<$re&1KFmF0BnK|DZ5h}eVrp)Nx%k?yXW2GOP5xA11h`T9QO;bl$c1_v!wp`LOlX9%KlwKqiifT^K2V6*uVr_#ApoG*OqT&k~4F?vp*y)mT4 zq)EV$FzRgKy-sNW1#F^setXYU-|kaDqxy_dDTo9jo=3MMejk2C)mEhhJ&N$P)q4r? za!O6-2;)~zi!khyPe-~AYbfq)Ij>=e0sW>KV{e0pKVTT@V>Ml?x$G{J@|AEd(Gs+? z9*YNA;d~QbP7aCs>(3Y@?TTfU31yQP-buw`)i#?VKo9zJ+kUnLmP~9DHg`(gRXY&d zkoby4EkhD&-DWLYtC*vyseXDo08!jnT+lZ{_hUO~Q=_uS% zn0`?>=&1AgConcFOK-Y%J?kj$3sSxr0KFs_Lwd*gk#h#Vnpp%nvzPvSe9mPMb};@8 zl{fUk&0RlF3hwJ+a6YK`(otb^>+LGZS+nPO$WHeicEiC(?%PLrOd;D5*_L2Fe|@d= zq6qa{7K`1rom1+M;k^e>s>0%Xy*OTaV2(`d)DKJkCSOhVInOd0UMS&f8bSj-$N(y~(i@YqDeIQyqPO^HPm9XP%O6KBr2 z(ALOTJ$!6m=xu#s<&Ui~Htq>-F|v|F4}L(do*uelFw0T@}SO)#iJ6*P3^*12y_k8n3A9R{SI}O}xnxLE9v$W`{tQ zx`$IDH9=v&%eeno1N)WxcH&vXUgFj|Ui*~_HLQG(T47EDVj5-V`KsUaf!Z?+|B1+0 zOE(yZAD5aTD=sZbfqab_DC>TR&)CsMIY$=JSMoub;!Hna{x~Zs%NI6UdYfJDl?T!U4~MrVx%ehp;V& zc4Qllhp?L*%HW~Ryo;_Ws|b~`U!3cN9&JBg$ywb4J@vRTC@^I4-dkb2$B|j5)-I3}+(t*r zm3Er^T~3z>@f2g!x)E;R=bN6A%=#7v6Wys3<)r4fuBr%hT71}b;?YZN*G4Fpx8`>% zO2awxqFh1goK<3YC1w}s)Fosari0t}i+T}V&>Tv-|JFZ}k3T$RLJ1E0jKBC53ts}Y zx>!RvjI^&~jz^rI`b!$O@L5`QnEXwdx~TRjO(A$ZWp!6J$Urk(^Bs!bReM`9zDJOz zvlqi5klnPpq^4$vpI2Nt)<8DAxfWSnCPTgW>;*Q_s^Gbfr?IKNytUN9!viAi_9j|= z>ssLA*OBR=3t%P1Y(8O!Og2uwURED=Cf#RS*RuDw&Wm8%I~6yt(dI=ex@}pWSeSgl W^8Fp93(ie4n9njkM0j2~b?QH^DOM8z literal 0 HcmV?d00001 diff --git a/zajecia4/.Rhistory b/zajecia4/.Rhistory new file mode 100644 index 0000000..40c297d --- /dev/null +++ b/zajecia4/.Rhistory @@ -0,0 +1,109 @@ +load(url("http://ls.home.amu.edu.pl/data_sets/czas_oczek_tramwaj.RData")) +a_est_mm <- mean(czas_oczek_tramwaj) - +sqrt(3) * sd(czas_oczek_tramwaj) * sqrt((length(czas_oczek_tramwaj) - 1) / (length(czas_oczek_tramwaj))) +b_est_mm<- mean(czas_oczek_tramwaj) + +sqrt(3) * sd(czas_oczek_tramwaj) * sqrt((length(czas_oczek_tramwaj) - 1) / (length(czas_oczek_tramwaj))) +load(url("http://ls.home.amu.edu.pl/data_sets/czas_oczek_tramwaj.RData")) +a_est_mm <- mean(czas_oczek_tramwaj) - +sqrt(3) * sd(czas_oczek_tramwaj) * sqrt((length(czas_oczek_tramwaj) - 1) / (length(czas_oczek_tramwaj))) +b_est_mm<- mean(czas_oczek_tramwaj) + +sqrt(3) * sd(czas_oczek_tramwaj) * sqrt((length(czas_oczek_tramwaj) - 1) / (length(czas_oczek_tramwaj))) +library(EnvStats) +EnvStats::eunif(czas_oczek_tramwaj, method = "mme") +a_est_mm <- mean(czas_oczek_tramwaj) - +sqrt(3) * sd(czas_oczek_tramwaj) +b_est_mm<- mean(czas_oczek_tramwaj) + +sqrt(3) * sd(czas_oczek_tramwaj) +# ZAD 1 +load(url("http://ls.home.amu.edu.pl/data_sets/czas_oczek_tramwaj.RData")) +a_est_mm <- mean(czas_oczek_tramwaj) - sqrt(3) * sd(czas_oczek_tramwaj) +b_est_mm<- mean(czas_oczek_tramwaj) + sqrt(3) * sd(czas_oczek_tramwaj) +library(EnvStats) +EnvStats::eunif(czas_oczek_tramwaj, method = "emm") +# ZAD 1 +load(url("http://ls.home.amu.edu.pl/data_sets/czas_oczek_tramwaj.RData")) +a_est_mm <- mean(czas_oczek_tramwaj) - sqrt(3) * sd(czas_oczek_tramwaj) +b_est_mm<- mean(czas_oczek_tramwaj) + sqrt(3) * sd(czas_oczek_tramwaj) +library(EnvStats) +EnvStats::eunif(czas_oczek_tramwaj, method = "mle") +# ZAD 1 +load(url("http://ls.home.amu.edu.pl/data_sets/czas_oczek_tramwaj.RData")) +a_est_mm <- mean(czas_oczek_tramwaj) - sqrt(3) * sd(czas_oczek_tramwaj) +b_est_mm<- mean(czas_oczek_tramwaj) + sqrt(3) * sd(czas_oczek_tramwaj) +library(EnvStats) +EnvStats::eunif(czas_oczek_tramwaj, method = "mme") +load(url("http://ls.home.amu.edu.pl/data_sets/czas_oczek_tramwaj.RData")) +a_est_mm <- mean(czas_oczek_tramwaj) - sqrt(3) * sd(czas_oczek_tramwaj) +b_est_mm<- mean(czas_oczek_tramwaj) + sqrt(3) * sd(czas_oczek_tramwaj) +library(EnvStats) +EnvStats::eunif(czas_oczek_tramwaj, method = "mme") +hist(czas_oczek_tramwaj, +xlab = "Czas oczekiwania na tramwaj", +main = "Rozklad empiryczny czasu oczekiwania na tramwaj", +probability = TRUE) +lines(density(czas_oczek_tramwaj), col = "red", lwd = 2) +curve(dunif(x, a_est, b_est), +add = TRUE, col = "blue", lwd = 2) +curve(dunif(x, a_est_mm, b_est_mm), +add = TRUE, col = "green", lwd = 2) +legend(x = 5, y = 0.04, legend = c("empiryczny", "teoretyczny ENW", "teoretyczny EMM"), +col = c("red", "blue", "green"), lwd = 2) +# ZAD 1 +load(url("http://ls.home.amu.edu.pl/data_sets/czas_oczek_tramwaj.RData")) +a_est_mm <- mean(czas_oczek_tramwaj) - sqrt(3) * sd(czas_oczek_tramwaj) +b_est_mm<- mean(czas_oczek_tramwaj) + sqrt(3) * sd(czas_oczek_tramwaj) +a_est <- min(czas_oczek_tramwaj) +b_est <- max(czas_oczek_tramwaj) +library(EnvStats) +EnvStats::eunif(czas_oczek_tramwaj, method = "mme") +hist(czas_oczek_tramwaj, +xlab = "Czas oczekiwania na tramwaj", +main = "Rozklad empiryczny czasu oczekiwania na tramwaj", +probability = TRUE) +lines(density(czas_oczek_tramwaj), col = "red", lwd = 2) +curve(dunif(x, a_est, b_est), +add = TRUE, col = "blue", lwd = 2) +curve(dunif(x, a_est_mm, b_est_mm), +add = TRUE, col = "green", lwd = 2) +legend(x = 5, y = 0.04, legend = c("empiryczny", "teoretyczny ENW", "teoretyczny EMM"), +col = c("red", "blue", "green"), lwd = 2) +load(url("http://ls.home.amu.edu.pl/data_sets/Centrala.RData")) +lambda_est <- mean(Centrala$Liczba) +probs <- dpois(sort(unique(Centrala$Liczba)), lambda = lambda_est) +sum(probs) +counts <- matrix(c(prop.table(table(Centrala$Liczba)), probs), nrow = 2, byrow = TRUE) +rownames(counts) <- c("empiryczny", "teoretyczny") +colnames(counts) <- sort(unique(Centrala$Liczba)) +barplot(counts, +xlab = "Liczba zgloszen", ylab = "Prawdopodobienstwo", +main = "Rozklady empiryczny i teoretyczny liczby zgloszen", +col = c("red", "blue"), legend = rownames(counts), beside = TRUE) +qqplot(rpois(length(Centrala$Liczba), lambda = lambda_est), Centrala$Liczba, +xlab = "Kwantyle teoretyczne", ylab = "Kwantyle empiryczne", +main = "Wykres kwantyl-kwantyl dla liczby zgloszen") +qqline(Centrala$Liczba, distribution = function(probs) { qpois(probs, lambda = lambda_est) }) +mean(Centrala$Liczba < 4) +ppois(3, lambda = lambda_est) +x <- c(0.9, 6.2, 2.1, 4.1, 7.3, +1.0, 4.6, 6.4, 3.8, 5.0, +2.7, 9.2, 5.9, 7.4, 3.0, +4.9, 8.2, 5.0, 1.2, 10.1, +12.2, 2.8, 5.9, 8.2, 0.5) +lambda_est <- mean(x^2) +hist(x, +xlab = "Srednia szybkosci wiatru", +main = "Rozklady empiryczny i teoretyczny sredniej szybkosc wiatru", +probability = TRUE, +col = "lightgreen") +lines(density(x), col = "red", lwd = 2) +curve(VGAM::drayleigh(x, sqrt(lambda_est / 2)), +add = TRUE, col = "blue", lwd = 2) +legend("topright", legend = c("empiryczny", "teoretyczny"), col = c("red", "blue"), lwd = 2) +qqplot(VGAM::rrayleigh(length(x), sqrt(lambda_est / 2)), x, +xlab = "Kwantyle teoretyczne", ylab = "Kwantyle empiryczne", +main = "Wykres kwantyl-kwantyl dla sredniej szybkosci wiatru") +qqline(x, distribution = function(probs) { VGAM::qrayleigh(probs, sqrt(lambda_est / 2)) }) +mean((x > 4) & (x < 8)) +VGAM::prayleigh(8, sqrt(lambda_est / 2)) - VGAM::prayleigh(4, sqrt(lambda_est / 2)) +0.5 * sqrt(pi * lambda_est) +lambda_est * (4 - pi) / 4