---
title: "Matematikai statisztika 1. gyakorlat, 2022/2023 I. félév"
output:
  html_document:
    df_print: paged
    #theme: darkly
  pdf_document: default
  word_document: default
date: '2022-09-12'
self_contained: yes
---

## R bevezető

## 1. feladat
```{r,eval=TRUE}
x=c(-5,1,3.5,7.5)#az osztályközepek
f=c(2,3,5,2) #gyakoriságok
mean=sum(x*f)/sum(f) #átlaghőm.
mean
xx=c(-5,-5,1,1,1,3,3,4,4,4,6,8) #egy minta adatsor
a=hist(xx,breaks=c(-10,0,2,5,10),main="Sűrűséghisztogram",prob="T",ylab="sűrűség",xlab="Cfok")
str(a)
```

## 2. feladat

```{r , echo=TRUE}
x=c(21,35,45,61)
plot(ecdf(x)) #tap.elo.fv, jobbról folytonos verzió (P(X<=z) a def)
for (i in 1:7)
print (quantile(x,0.33,type=i))
```


```{r,eval=TRUE}
#napi kozephomerseklet adatok 1951-1988 között

homb<-read.table("http://zempleni.elte.hu/nyir-51-88m.hom") #beolvasás:Nyíregyháza
homk<-read.table("http://zempleni.elte.hu/karc-51-88.hom")  #Karcag, az eredmény 1 oszlopos mátrix
#nézzük is meg:
head(homb)
#OK-nak tűnik
r=c(0:37) #az egész számok 0 és 37 között

jan1b<-homb[365*r+1,1] #minden év első napja
jan1k<-homk[365*r+1,1]
summary(jan1b)
#mi lenne szept. 12?
jan1b<-homb[365*r+31+28+31+30+31+30+31+31+12,1]  #az év 255.napja
jan1k<-homk[365*r+31+28+31+30+31+30+31+31+12,1]
summary(jan1b)

#mindenki küldje be a chat-be a saját születésnapjára vonatkozó alapstatisztikákat
```
### Adattisztítás




```{r,eval=TRUE}
#keressük meg minden napra a minimum es #maximum hőmérsékleteket, Nyíregyházára
#ábrázoljuk is
par(mfrow=c(1,1)) #grafikus beállítás: 1 ábra
minb=rep(0, times=365)
maxb=minb
for (i in 1:365)
{bpi<-homb[365*r+i,1]
minb[i]=min(bpi)
maxb[i]=max(bpi)
}
plot(minb,type="l",col=4,ylim=c(min(minb),max(maxb))) #ábra készítés. "l": vonalas, 
#az y tengelyt manuálisan állítjuk be, mert másik görbe is lesz
lines(maxb,col=2)
#hoppá! Keressük meg a hibát
homb[order(homb[,1])[dim(homb)[1]],1] #mennyi a max - kicsit bonyolultan, mert így azt is tudjuk...
order(homb[,1])[dim(homb)[1]] #hogy melyik megfigyelés a hibás
k=order(homb[,1])[dim(homb)[1]] 
i=k-365*trunc(k/365) 
r=c(c(1:34),36,37)
bpi<-mean(homb[365*r+i,1])

plot(homb[(k-5):(k+5),1],type="l") #vonaldiagram
lines(homk[(k-5):(k+5),1],col=2) #piros a Karcag

```

### Házi feladat (2 pont): algoritmust adni a hibás adat kijavítására. Canvasban beküldendő a kód és a javasolt javított érték. Határidő: szeptember 12, este 8 óra
##keressük meg a hibás adatot, nézzük meg a környező értékeket, a másik adatsort is és ez alapján javasoljunk javítást
