Tak jsme včera konečně odevzdali ten strašný Ass no. 1 (samozřejmě, jde o Assignment z ekonometrie). Byla to hrůza a strávili jsme s tím několik desítek drahých hodin, plných stresu, zoufalství a lámání si hlavy nad nepochopitelnými vzorci. Na konci příspěvku přikládám výňatek ze skriptu R programu, ve kterém to všechno programujeme...
Ovšem hned zítra dostaneme další Ass, takže není moc důvod k radosti. Ale dnes aspoň na "oslavu" jdeme využít naše volňásky do Tivoli (volňásky jsme sice dostali zadarmo na IGOPu, ale celý IGOP byl placený, takže jsme si je vlastně poctivě zaplatili...), protože jim končí platnost v neděli. Tivoli totiž na pár měsíců zavírá, otevře až na konci listopadu a po Vánocích zase zavře. Není se čemu divit, asi ani Dánům se nechce v mrazivém vzduchu fičet na horské dráze a omrzat na všech možných místech... stačí, že budou jezdit na kole! Však včera večer už docela studeně pofukovalo (krásné letní dny jsou asi tatam) a člověku je tak... nepříjemno... už přesně vím, které části těla budou na kole v zimě nejvíc trpět. Prsty u nohou, ruce, uši a krk. I když, tím šlapáním se člověk aspoň zahřeje... Jsem zvědavá, kolik lidí budu v zimě na kole potkávat (jestli to taky nezabalím a nekoupím si lítačku na metro). Mimochodem, tohle větrno- studeno- občas slunečné počasí nejspíš mate i Dány, protože každé ráno, když se stovkami ostatních frčím z Amageru přes mosty do centra, potkávám lidi oblečené velice různorodě- někdo je v krátké sukni a silonkách, někdo v dlouhých kalhotech a přes to má ještě plastové "rain pants", někdo je zahalený v šále a někdo naopak odhalený a má jen vlající rozepnutý svetr, někdo má vysoké gumové boty do deště a někdo naopak sandály, někdo má na hlavě sluneční brýle a jiný ušanku... jako by si ani oni nebyli jistí, jak dneska bude...
... I když, ruku na srdce- kdo si tady může být jistý? Zvlášť když jde o počasí.
Jinak se Kodaň teď na pár dní stala dějištěm velice důležité sportovní události- koná se tu UCI Road World Championships 2011 (http://www.cyclingnews.com/races/uci-road-world-championships-2011), což je taková cyklistická obdoba Formule 1 v Monte Carlu, protože závody na kolech se jezdí přímo na ulicích v centru města, každý den několik. Spousta lidí z toho byla nadšená už dopředu, evidentně je to velice populární závod. Ovšem, museli kvůli tomu zahradit , zablokovat nebo zúžit spoustu jinak dost rušných ulic, včetně hlavních tahů centrem města, kudy jezdí každý den tisíce lidí na kole i v autě.
No a že řidiči nemůžou projet autem, přesedlala jich většina na kolo. Ale na kole to taky projet nejde a je to dost záludné, protože člověk si frčí ráno hezky načas svou obvyklou cestou do práce nebo na univerzitu, a najednou- ha! Cesta je zahrazená, policie mává a směruje proud cyklistů do vedlejších ulic. Pro člověka, který v Kodani vyrostl, určitě žádný problém. Kdybych nemohla projet Národní třídou a Na Příkopech od Vltavy k Náměstí Republiky, taky si poradím... jenže kdo to ve městě nezná, ten se lehko ztratí... Což se mi sice zatím nestalo, ale ještě se klidně stát může... závody končí až zítra... a vůbec, pokud vás to nedonutí se ztratit, tak vás to minimálně dost zpomalí a zdrží... Čili je nutno ráno vyrazit ještě o pět minut dřív.
Ještě jednu perličku na závěr si dovolím- náš houmles profesor Frank z předmětu Natural Resources Economics (to je ten, kterého sebrali z Hlaváku) je sice vágus, ale je docela přátelský a ukázal i smysl pro humor. Když se mu včera nechtělo pokračovat ve výkladu, pokoušel se dva spolužáky přemluvit, aby zatančili nahatí, a pak ještě pověděl sprostý vtip. A pak dal jednu z mnoha přestávek "aby si tady kluci z Řecka mohli jít zakouřit". Haha, stejně chce kouřit sám...
Tak ještě na závěr slibovaný výtažek R skriptu z Ass no.1:
par(mfrow=c(1,2))
plot(density(volpers))
plot(density(log(volpers)))
eq1<-lm(log(volpers)~log(price)+age+I(age*age)+high_edu+kid06+kid714+kid1520+capital+urban+single+log(totfoodpers))
eq1
summary(eq1)
lower<-summary(eq1)$coef["single",1]-1.96*summary(eq1)$coef["single",2]
lower
upper<-summary(eq1)$coef["single",1]+1.96*summary(eq1)$coef["single",2]
upper
white<-eq1$res^2
white
fit<-eq1$fit
fit
fit2<-eq1$fit^2
fit2
whitetest<-lm(white~fit+fit2)
whitetest
summary(whitetest)
library(sandwich)
robstderr<-sqrt(diag(vcovHC(eq1)))
robstderr
stderr<-sqrt(diag(vcov(eq1)))
stderr
cbind(coef(eq1),stderr,robstderr)
Kdo tomu porozumí, má u mě lízátko! (A může mi rovnou pomoci se skriptem na Ass no.2)
No comments:
Post a Comment