Resultados

Estatísticas Básicas:

Para obter as estatísticas básicas, instalou-se o pacote fBasics, além de carregar o datasets, primordialmente.

>data(airquality)

Ozônio:

> basicStats(airquality$Ozone)

            X..airquality.Ozone

nobs                 153.000000

NAs                   37.000000

Minimum                1.000000

Maximum              168.000000

1. Quartile           18.000000

3. Quartile           63.250000

Mean                  42.129310

Median                31.500000

Sum                 4887.000000

SE Mean                3.062848

LCL Mean              36.062398

UCL Mean              48.196223

Variance            1088.200525

Stdev                 32.987885

Skewness               1.209866

Kurtosis               1.112243

Vento:

> basicStats(airquality$Wind)

            X..airquality.Wind

nobs                153.000000

NAs                   0.000000

Minimum               1.700000

Maximum              20.700000

1. Quartile           7.400000

3. Quartile          11.500000

Mean                  9.957516

Median                9.700000

Sum                1523.500000

SE Mean               0.284818

LCL Mean              9.394804

UCL Mean             10.520229

Variance             12.411539

Stdev                 3.523001

Skewness              0.341028

Kurtosis              0.028865

Histogramas:

>caminho="C:/Users/Cheyenne/Documents/UFRJ/geoestat/P3/"

>amazon=read.table(paste(caminho,"ov.txt",sep=""))

Ozônio:

> x.norm<-amazon[,1]

> h<-hist(x.norm,breaks=10)

> xhist<-c(min(h$breaks),h$breaks)

> yhist<-c(0,h$density,0)

> xfit<-seq(0,200,by=1.0)

> yfit<-dnorm(xfit,mean=mean(x.norm),sd=sd(x.norm))

> plot(xhist,yhist,type="s",ylim=c(0,max(yhist,yfit)),xlab="Concentração de Ozônio (ppb)",ylab="Frequência  relativa",main="Histograma de Ozônio",col.axis="blue")

> lines(xfit,yfit,col="red")

Vento:

> x.norm<-amazon[,2]
> h<-hist(x.norm,breaks=10)
> xhist<-c(min(h$breaks),h$breaks)
> yhist<-c(0,h$density,0)
> xfit<-seq(0,200,by=1.0)
> yfit<-dnorm(xfit,mean=mean(x.norm),sd=sd(x.norm))
>plot(xhist,yhist,type="s",ylim=c(0,max(yhist,yfit)),xlab="Velocidade Média (mph)",ylab="Frequência Relativa",main="Histograma de Vento",col.axis="blue")
> lines(xfit,yfit,col="red")

    

Dispersograma:

> plot(Ozone~Wind,data=airquality,main="Dispersograma",xlab="Ozônio",ylab="Vento",col.axis="blue")

Coeficiente de Correlação:

> cor(amazon[,1],amazon[,2])

[1] -0.5319017

Modelo de Regressão Linear:

> g.m1<-lm(Ozone~Wind,data=airquality)

> plot(Ozone~Wind,data=airquality,main="Modelo de Regressão Linear",xlab="Ozônio",ylab="Vento")

> abline(g.m1,col="red")