### Korrelációs mátrix készítése ###

# Az ehhez szükséges csomag telepítése az R-ben rendszergazdaként:

install.packages("corrplot")

# Telepítés után bezárjuk és újranyitjuk az R-t, majd betöltjük a csomagot

library(corrplot)

# Korrelációkat tartalmazó adatmátrix betöltése statisztika mappából

load("ratings.cor.RData")

# vagy közvetlenül a honlapról

load(url("https://phon.nytud.hu/mady/courses/statistics/2025/ratings.cor.RData"))

# correlációs mátrix létrehozása

m = cor(ratings.cor)

# mátrix ábrázolása

corrplot(m)

corrplot(m,method="ellipse")
corrplot(m,method="number")


### Házi feladat ###

plot(ratings$FreqSingular,ratings$FreqPlural,col="blue",main="Növény- és állatnevek gyakorisága",xlab="Egyes szám", ylab = "Többes szám")

# Ábrázolás a két kilógó érték nélkül

plot(ratings$FreqSingular,ratings$FreqPlural,xlim=c(0,1000),ylim=c(0,600))

# Dobozdiagram készítése becsült testsúlyra osztályonként

boxplot(meanWeightRating~Class, ratings)

### Logikai vektorok ###

# a Class oszlopban "plant" betűsort tartalmazó sorok azonosítása
noveny = ratings$Class == "plant"

# ugyanígy az állatok kiválasztása (szűrése)
allat = ratings$Class == "animal"

# ábrázolás azonos skálázású tengelyeken

plot(ratings$Length[noveny],ratings$Frequency[noveny], col="green", xlim = c(3,10), ylim=c(1.5,8)
par(new=T)
plot(ratings$Length[allat],ratings$Frequency[allat], col="brown", xlim = c(3,10), ylim=c(1.5,8)


###### Normális eloszlás sűrűségfüggvénye ######

# Hozzunk létre egy normális eloszlású halmazt 1000 elemmel.
# Az átlag legyen 100, a szórás 20.

norm1 = rnorm(1000,mean=100,sd=20)

# Az átlag legyen 200, a szórás 20.

norm2 = rnorm(1000,mean=200,sd=20)

# Most az átlag legyen 200, a szórás 40.

norm2 = rnorm(1000,mean=200,sd=40)

# Standardizáljuk az eloszlást a z-transzformáció segítségével.

znorm1 = scale(norm1)
znorm2 = scale(norm2)
znorm3 = scale(norm3)

# Ábrázoljuk egymás mellett az eloszláshoz tartozó normális és standard normális eloszlást.
# X tengely határértékei: normál: 0 és 300, standard normál: -3,+3

par(mfrow=c(3,2))
plot(density(norm1),xlim=c(0,500))
plot(density(znorm1),xlim=c(-3,3))
plot(density(norm2),xlim=c(0,500))
plot(density(znorm2),xlim=c(-3,3))
plot(density(norm3),xlim=c(0,500))
plot(density(znorm3),xlim=c(-3,3))