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)
> basicStats(airquality$Ozone)
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")
> 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")
Nenhum comentário:
Postar um comentário