Moja databáza obsahuje údaje o počasí zozbierané v meste Dillí za
obdobie 4 rokov (od roku 2013 do roku 2017).
Štatistický popis denného klimatického testu v Delhí
Teplota: V Delhí sa priemerná denná teplota pohybuje
medzi 11 °C (najchladnejší deň) a 34.5 °C (najteplejší deň), pričom
väčšina dní má teplotu okolo 20 °C (medián 19.88 °C). To naznačuje
mierne teplé počasie s občasnými teplými dňami.
Vlhkosť: Hodnoty vlhkosti ukazujú, že vzduch je
pomerne suchý (min 17.75 %) až veľmi vlhký (max 95.83 %), s mediánom
okolo 57.75 %, čo znamená strednú úroveň vlhkosti.
Rýchlosť vetra: Vietor je pomerne mierny, väčšina
dní má rýchlosť medzi 1.4 až 19.3 jednotkami, s mediánom okolo 8.07.
Tlak vzduchu: Hodnota minimálneho tlaku 59 vyzerá
nezvyčajne nízka, pravdepodobne chyba v dátach alebo iná jednotka.
Ostatné hodnoty (medzi 1007 a 1023) sú štandardné pre atmosférický tlak,
kde medián je okolo 1013 hPa, čo je blízko normálneho tlaku na úrovni
mora.
Boxplot
df <- read.csv("DailyDelhiClimateTest.csv")
boxplot(df$meantemp,
main = "Boxplot priemernej dennej teploty",
ylab = "Teplota (°C)",
col = "yellow",
border = "purple")

Väčšina denných teplôt sa pohybuje medzi 17 °C a 27 °C. Medián
(typická hodnota) je okolo 20 °C, čo naznačuje, že polovica dní má
teplotu pod 20 °C a polovica nad. Rozptyl je pomerne veľký – od cca 13
°C do 35 °C, čo znamená, že sú dni s výrazne nižšími aj vyššími
teplotami. Žiadne extrémne odľahlé hodnoty (outliers) nie sú
zobrazené.
Lineárna regresia
model <- lm(meantemp ~ humidity, data = df)
summary(model)
Call:
lm(formula = meantemp ~ humidity, data = df)
Residuals:
Min 1Q Median 3Q Max
-7.1322 -2.2974 0.2871 2.6521 7.7067
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) 37.8081 0.9621 39.30 <2e-16 ***
humidity -0.2861 0.0162 -17.66 <2e-16 ***
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
Residual standard error: 3.284 on 112 degrees of freedom
Multiple R-squared: 0.7357, Adjusted R-squared: 0.7333
F-statistic: 311.8 on 1 and 112 DF, p-value: < 2.2e-16
Korelácia
cor(data$wind_speed, data$meanpressure, use = "complete.obs")
[1] 0.1303523
plot(data$wind_speed, data$meanpressure,
main = "Bodový graf: Rýchlosť vetra vs. Priemerný tlak",
xlab = "Rýchlosť vetra (wind_speed)",
ylab = "Priemerný tlak (meanpressure)",
pch = 19, # Nastaví plné bodky
col = "pink") # Nastaví farbu bodiek

Lineárny vzťah medzi rýchlosťou vetra a priemerným tlakom je veľmi
slabý (r=0.13$).
Graf na regresiu
library(ggplot2)
ggplot(df, aes(x = humidity, y = meantemp)) +
geom_point(color = "blue") +
geom_smooth(method = "lm", color = "red", se = TRUE)

labs(title = "Lineárna regresia: Teplota vs. Vlhkosť",
x = "Vlhkosť (%)",
y = "Priemerná teplota (°C)") +
theme_minimal()
NULL
Keď vlhkosť stúpa, teplota klesá. Vysoká vlhkosť je spojená s nižšími
teplotami. Model je veľmi silný (R² ≈ 0.74), takže vlhkosť je
dobrý prediktor teploty.
T-test
t.test(data$meantemp, data$humidity, paired = TRUE)
Paired t-test
data: data$meantemp and data$humidity
t = -14.909, df = 113, p-value < 2.2e-16
alternative hypothesis: true mean difference is not equal to 0
95 percent confidence interval:
-39.13596 -29.95461
sample estimates:
mean difference
-34.54528
Nulová hypotéza (H₀): medzi meantemp a humidity nie je
rozdiel (priemerný rozdiel = 0) Alternatívna hypotéza (H₁):
medzi nimi je rozdiel
Výsledok párového t-testu (t = –14.909, df = 113, p < 0.001)
ukazuje, že rozdiel medzi premennými meantemp a humidity je
štatisticky vysoko významný. Nulová hypotéza o rovnosti priemerov
sa zamieta. Priemerná teplota je teda významne nižšia ako priemerná
vlhkosť, pričom 95 % interval spoľahlivosti pre rozdiel sa pohybuje
medzi –39.14 a –29.95.
Heatmap korelačnej matice numerických veličín
library(corrplot)
# Krok 2: Vytvorte korelačnú maticu (iba z numerických stĺpcov)
# Vylúčime stĺpec 'date' (ktorý je na pozícii 1)
cor_matrix <- cor(data[, -1], use = "pairwise.complete.obs")
# Krok 3: Vytvorte Heatmap
corrplot(cor_matrix,
method = "color",
type = "upper",
addCoef.col = "black",
tl.col = "black",
tl.srt = 45,
diag = FALSE)

Dominantný vzťah: V dátach dominuje veľmi silná inverzná
závislosť medzi priemernou teplotou a vlhkosťou (r=-0.86). Slabé
vzťahy: Väčšina ostatných vzťahov je slabá (v rozsahu od 0.03 do
-0.34). Takmer žiadny vzťah: Korelácie zahŕňajúce priemerný
tlak (meanpressure) sú extrémne blízke nule (0.03, -0.10, 0.13), čo
naznačuje, že priemerný tlak v rámci tohto súboru dát je nezávislý od
ostatných sledovaných premenných. ```
LS0tCnRpdGxlOiAiRGF0YWLDoXphIgphdXRob3I6ICJNw6FyaWEgTWF0w7rFoW92w6EiCmRhdGU6ICIyMDI1LTEwLTE2IgpvdXRwdXQ6IAogIGh0bWxfbm90ZWJvb2s6CiAgICB0b2M6IHRydWUKICAgIHRvY19mbG9hdDogdHJ1ZQogICAgdGhlbWU6IGNvc21vCiAgICBoaWdobGlnaHQ6IHplbmJ1cm4KZWRpdG9yX29wdGlvbnM6IAogIG1hcmtkb3duOiAKICAgIHdyYXA6IDcyCgoKLS0tCk1vamEgZGF0YWLDoXphIG9ic2FodWplIMO6ZGFqZSBvIHBvxI1hc8OtIHpvemJpZXJhbsOpIHYgbWVzdGUgRGlsbMOtIHphIG9iZG9iaWUgNCByb2tvdiAob2Qgcm9rdSAyMDEzIGRvIHJva3UgMjAxNykuCmBgYHtyfQpkYXRhIDwtIHJlYWQuY3N2KCJEYWlseURlbGhpQ2xpbWF0ZVRlc3QuY3N2IikKaGVhZChkYXRhKQpgYGAKCioqU3VtbWFyeSoqCmBgYHtyfQpzdW1tYXJ5KGRhdGEpCmBgYAoKIyDFoHRhdGlzdGlja8O9IHBvcGlzIGRlbm7DqWhvIGtsaW1hdGlja8OpaG8gdGVzdHUgdiBEZWxow60KCioqVGVwbG90YSoqOiBWIERlbGjDrSBzYSBwcmllbWVybsOhIGRlbm7DoSB0ZXBsb3RhIHBvaHlidWplIG1lZHppIDExIMKwQyAobmFqY2hsYWRuZWrFocOtIGRlxYgpIGEgMzQuNSDCsEMgKG5hanRlcGxlasWhw60gZGXFiCksIHByacSNb20gdsOkxI3FoWluYSBkbsOtIG3DoSB0ZXBsb3R1IG9rb2xvIDIwIMKwQyAobWVkacOhbiAxOS44OCDCsEMpLiBUbyBuYXpuYcSNdWplIG1pZXJuZSB0ZXBsw6kgcG/EjWFzaWUgcyBvYsSNYXNuw71taSB0ZXBsw71taSBkxYhhbWkuCgoqKlZsaGtvc8WlKio6IEhvZG5vdHkgdmxoa29zdGkgdWthenVqw7osIMW+ZSB2emR1Y2ggamUgcG9tZXJuZSBzdWNow70gKG1pbiAxNy43NSAlKSBhxb4gdmXEvm1pIHZsaGvDvSAobWF4IDk1LjgzICUpLCBzIG1lZGnDoW5vbSBva29sbyA1Ny43NSAlLCDEjW8gem5hbWVuw6Egc3RyZWRuw7ogw7pyb3ZlxYggdmxoa29zdGkuCgoqKlLDvWNobG9zxaUgdmV0cmEqKjogVmlldG9yIGplIHBvbWVybmUgbWllcm55LCB2w6TEjcWhaW5hIGRuw60gbcOhIHLDvWNobG9zxaUgbWVkemkgMS40IGHFviAxOS4zIGplZG5vdGthbWksIHMgbWVkacOhbm9tIG9rb2xvIDguMDcuCgoqKlRsYWsgdnpkdWNodSoqOiBIb2Rub3RhIG1pbmltw6FsbmVobyB0bGFrdSA1OSB2eXplcsOhIG5lenZ5xI1ham5lIG7DrXprYSwgcHJhdmRlcG9kb2JuZSBjaHliYSB2IGTDoXRhY2ggYWxlYm8gaW7DoSBqZWRub3RrYS4gT3N0YXRuw6kgaG9kbm90eSAobWVkemkgMTAwNyBhIDEwMjMpIHPDuiDFoXRhbmRhcmRuw6kgcHJlIGF0bW9zZsOpcmlja8O9IHRsYWssIGtkZSBtZWRpw6FuIGplIG9rb2xvIDEwMTMgaFBhLCDEjW8gamUgYmzDrXprbyBub3Jtw6FsbmVobyB0bGFrdSBuYSDDunJvdm5pIG1vcmEuCgoKKipCb3hwbG90KioKYGBge3J9CmRmIDwtIHJlYWQuY3N2KCJEYWlseURlbGhpQ2xpbWF0ZVRlc3QuY3N2IikKYm94cGxvdChkZiRtZWFudGVtcCwKICAgICAgICBtYWluID0gIkJveHBsb3QgcHJpZW1lcm5laiBkZW5uZWogdGVwbG90eSIsIAogICAgICAgIHlsYWIgPSAiVGVwbG90YSAowrBDKSIsIAogICAgICAgIGNvbCA9ICJ5ZWxsb3ciLAogICAgICAgIGJvcmRlciA9ICJwdXJwbGUiKQpgYGAKClbDpMSNxaFpbmEgZGVubsO9Y2ggdGVwbMO0dCBzYSBwb2h5YnVqZSBtZWR6aSAxNyDCsEMgYSAyNyDCsEMuCk1lZGnDoW4gKHR5cGlja8OhIGhvZG5vdGEpIGplIG9rb2xvIDIwIMKwQywgxI1vIG5hem5hxI11amUsIMW+ZSBwb2xvdmljYSBkbsOtIG3DoSB0ZXBsb3R1IHBvZCAyMCDCsEMgYSBwb2xvdmljYSBuYWQuClJvenB0eWwgamUgcG9tZXJuZSB2ZcS+a8O9IOKAkyBvZCBjY2EgMTMgwrBDIGRvIDM1IMKwQywgxI1vIHpuYW1lbsOhLCDFvmUgc8O6IGRuaSBzIHbDvXJhem5lIG5pxb7FocOtbWkgYWogdnnFocWhw61taSB0ZXBsb3RhbWkuCsW9aWFkbmUgZXh0csOpbW5lIG9kxL5haGzDqSBob2Rub3R5IChvdXRsaWVycykgbmllIHPDuiB6b2JyYXplbsOpLgoKCioqTGluZcOhcm5hIHJlZ3Jlc2lhKioKYGBge3J9Cm1vZGVsIDwtIGxtKG1lYW50ZW1wIH4gaHVtaWRpdHksIGRhdGEgPSBkZikKc3VtbWFyeShtb2RlbCkKYGBgCgoKKipLb3JlbMOhY2lhKioKYGBge3J9CmNvcihkYXRhJHdpbmRfc3BlZWQsIGRhdGEkbWVhbnByZXNzdXJlLCB1c2UgPSAiY29tcGxldGUub2JzIikKcGxvdChkYXRhJHdpbmRfc3BlZWQsIGRhdGEkbWVhbnByZXNzdXJlLAogICAgIG1haW4gPSAiQm9kb3bDvSBncmFmOiBSw71jaGxvc8WlIHZldHJhIHZzLiBQcmllbWVybsO9IHRsYWsiLAogICAgIHhsYWIgPSAiUsO9Y2hsb3PFpSB2ZXRyYSAod2luZF9zcGVlZCkiLAogICAgIHlsYWIgPSAiUHJpZW1lcm7DvSB0bGFrIChtZWFucHJlc3N1cmUpIiwKICAgICBwY2ggPSAxOSwgICMgTmFzdGF2w60gcGxuw6kgYm9ka3kKICAgICBjb2wgPSAicGluayIpICMgTmFzdGF2w60gZmFyYnUgYm9kaWVrCmBgYApMaW5lw6FybnkgdnrFpWFoIG1lZHppIHLDvWNobG9zxaVvdSB2ZXRyYSBhIHByaWVtZXJuw71tIHRsYWtvbSBqZSB2ZcS+bWkgc2xhYsO9IChyPTAuMTMkKS4KCgoKKipHcmFmIG5hIHJlZ3Jlc2l1KioKYGBge3J9CmxpYnJhcnkoZ2dwbG90MikKCmdncGxvdChkZiwgYWVzKHggPSBodW1pZGl0eSwgeSA9IG1lYW50ZW1wKSkgKwogIGdlb21fcG9pbnQoY29sb3IgPSAiYmx1ZSIpICsgICAgICAgICAgICAgICAgICAgICAgICAgICAKICBnZW9tX3Ntb290aChtZXRob2QgPSAibG0iLCBjb2xvciA9ICJyZWQiLCBzZSA9IFRSVUUpIAogIGxhYnModGl0bGUgPSAiTGluZcOhcm5hIHJlZ3Jlc2lhOiBUZXBsb3RhIHZzLiBWbGhrb3PFpSIsCiAgICAgICB4ID0gIlZsaGtvc8WlICglKSIsCiAgICAgICB5ID0gIlByaWVtZXJuw6EgdGVwbG90YSAowrBDKSIpICsKICB0aGVtZV9taW5pbWFsKCkKYGBgCktlxI8gdmxoa29zxaUgc3TDunBhLCB0ZXBsb3RhIGtsZXPDoS4KVnlzb2vDoSB2bGhrb3PFpSBqZSBzcG9qZW7DoSBzIG5pxb7FocOtbWkgdGVwbG90YW1pLgpNb2RlbCBqZSAqdmXEvm1pIHNpbG7DvSogKFLCsiDiiYggMC43NCksIHRha8W+ZSB2bGhrb3PFpSBqZSBkb2Jyw70gcHJlZGlrdG9yIHRlcGxvdHkuCgpgYGB7cn0KCmBgYAoKKipULXRlc3QqKgpgYGB7cn0KdC50ZXN0KGRhdGEkbWVhbnRlbXAsIGRhdGEkaHVtaWRpdHksIHBhaXJlZCA9IFRSVUUpCmBgYAoqTnVsb3bDoSBoeXBvdMOpemEgKEjigoApKjogbWVkemkgbWVhbnRlbXAgYSBodW1pZGl0eSBuaWUgamUgcm96ZGllbCAocHJpZW1lcm7DvSByb3pkaWVsID0gMCkKKkFsdGVybmF0w612bmEgaHlwb3TDqXphIChI4oKBKSo6IG1lZHppIG5pbWkgamUgcm96ZGllbAoKVsO9c2xlZG9rIHDDoXJvdsOpaG8gdC10ZXN0dSAodCA9IOKAkzE0LjkwOSwgZGYgPSAxMTMsIHAgPCAwLjAwMSkgdWthenVqZSwgxb5lICpyb3pkaWVsIG1lZHppIHByZW1lbm7DvW1pIG1lYW50ZW1wIGEgaHVtaWRpdHkgamUgxaF0YXRpc3RpY2t5IHZ5c29rbyB2w716bmFtbsO9LioKTnVsb3bDoSBoeXBvdMOpemEgbyByb3Zub3N0aSBwcmllbWVyb3Ygc2EgemFtaWV0YS4gUHJpZW1lcm7DoSB0ZXBsb3RhIGplIHRlZGEgdsO9em5hbW5lIG5pxb7FoWlhIGFrbyBwcmllbWVybsOhIHZsaGtvc8WlLCBwcmnEjW9tIDk1ICUgaW50ZXJ2YWwgc3BvxL5haGxpdm9zdGkgcHJlIHJvemRpZWwgc2EgcG9oeWJ1amUgbWVkemkg4oCTMzkuMTQgYSDigJMyOS45NS4KYGBgCmBgYAoKKipIZWF0bWFwIGtvcmVsYcSNbmVqIG1hdGljZSBudW1lcmlja8O9Y2ggdmVsacSNw61uKioKYGBge3J9CmxpYnJhcnkoY29ycnBsb3QpCgojIEtyb2sgMjogVnl0dm9ydGUga29yZWxhxI1uw7ogbWF0aWN1IChpYmEgeiBudW1lcmlja8O9Y2ggc3TEunBjb3YpCiMgVnlsw7rEjWltZSBzdMS6cGVjICdkYXRlJyAoa3RvcsO9IGplIG5hIHBvesOtY2lpIDEpCmNvcl9tYXRyaXggPC0gY29yKGRhdGFbLCAtMV0sIHVzZSA9ICJwYWlyd2lzZS5jb21wbGV0ZS5vYnMiKQoKIyBLcm9rIDM6IFZ5dHZvcnRlIEhlYXRtYXAKY29ycnBsb3QoY29yX21hdHJpeCwKICAgICAgICAgbWV0aG9kID0gImNvbG9yIiwgICAgICAgCiAgICAgICAgIHR5cGUgPSAidXBwZXIiLCAgICAgICAgIAogICAgICAgICBhZGRDb2VmLmNvbCA9ICJibGFjayIsICAKICAgICAgICAgdGwuY29sID0gImJsYWNrIiwgICAgICAgCiAgICAgICAgIHRsLnNydCA9IDQ1LCAgICAgICAgICAgIAogICAgICAgICBkaWFnID0gRkFMU0UpCmBgYAoqRG9taW5hbnRuw70gdnrFpWFoKjogViBkw6F0YWNoIGRvbWludWplIHZlxL5taSBzaWxuw6EgaW52ZXJ6bsOhIHrDoXZpc2xvc8WlIG1lZHppIHByaWVtZXJub3UgdGVwbG90b3UgYSB2bGhrb3PFpW91IChyPS0wLjg2KS4KKlNsYWLDqSB2esWlYWh5KjogVsOkxI3FoWluYSBvc3RhdG7DvWNoIHZ6xaVhaG92IGplIHNsYWLDoSAodiByb3pzYWh1IG9kIDAuMDMgZG8gLTAuMzQpLgoqVGFrbWVyIMW+aWFkbnkgdnrFpWFoKjogS29yZWzDoWNpZSB6YWjFlcWIYWrDumNlIHByaWVtZXJuw70gdGxhayAobWVhbnByZXNzdXJlKSBzw7ogZXh0csOpbW5lIGJsw616a2UgbnVsZSAoMC4wMywgLTAuMTAsIDAuMTMpLCDEjW8gbmF6bmHEjXVqZSwgxb5lIHByaWVtZXJuw70gdGxhayB2IHLDoW1jaSB0b2h0byBzw7pib3J1IGTDoXQgamUgbmV6w6F2aXNsw70gb2Qgb3N0YXRuw71jaCBzbGVkb3ZhbsO9Y2ggcHJlbWVubsO9Y2guCmBgYA==