Estadisticas Del Gasoducto En Manzanillo

ESTADISTICAS DEL GASODUCTO MANZANILLO-GUADALAJARA
> install.packages(c("maptools","PBSmapping"))
> library(PBSmapping)
>nombrearchivo<-importShapefile("nombredarchivoshape",readDBF=TRUE)
Mueven los archivos a la ruta que les aparece
>plotPolys(nombrearchivo,axes=TRUE,bg="beige",main="tituloarchivo",xlab="Longitud",ylab="Altitud")

Para saber cuantas variables contieneel archivo shape usamos:
> ls(nombrearchivo)
[1] "PID" "POS" "SID" "X" "Y"
En los cálculos estadísticos tomamos [[1]] como PID, [[2]] como POS, [[3]] como SID, [[4]] como X, y [[5]] como Y.PROMEDIOS
> mean(MiArchivoShape[[1]])
[1] 3.195759
> mean(MiArchivoShape[[2]])
[1] 1
> mean(MiArchivoShape[[3]])
[1] 69.89396
> mean(MiArchivoShape[[4]])
[1] 642882
>mean(MiArchivoShape[[5]])
[1] 2165225
MEDIANA
> median(MiArchivoShape[[1]])
[1] 3
> median(MiArchivoShape[[2]])
[1] 1
> median(MiArchivoShape[[3]])
[1] 56
> median(MiArchivoShape[[4]])[1] 639194.1
> median(MiArchivoShape[[5]])
[1] 2148315

MAXIMOS
> max(MiArchivoShape[[1]])
[1] 6
> max(MiArchivoShape[[2]])
[1] 1
> max(MiArchivoShape[[3]])
[1] 210
>max(MiArchivoShape[[4]])
[1] 713657.7
> max(MiArchivoShape[[5]])
[1] 2293408
MINIMOS
> min(MiArchivoShape[[1]])
[1] 1
> min(MiArchivoShape[[2]])
[1] 1
> min(MiArchivoShape[[3]])
[1] 1> min(MiArchivoShape[[4]])
[1] 546869.5
> min(MiArchivoShape[[5]])
[1] 2066026

VARIANZA
> var(MiArchivoShape[[1]])
[1] 2.373381
> var(MiArchivoShape[[2]])
[1] 0
>var(MiArchivoShape[[3]])
[1] 2926.732
> var(MiArchivoShape[[4]])
[1] 2159261606
> var(MiArchivoShape[[5]])
[1] 4770993709

DESVIACIÓN ESTANDAR
> sd(MiArchivoShape[[1]])
[1] 1.540578
>sd(MiArchivoShape[[2]])
[1] 0
> sd(MiArchivoShape[[3]])
[1] 54.09928
> sd(MiArchivoShape[[4]])
[1] 46467.86
> sd(MiArchivoShape[[5]])
[1] 69072.38

INTEGRANDO GOOGLEMAPS EN R
USA(TX1)...