Diagrama de dispersión 2D

En este script se utilizan los datos quakes del paquete datasets para visualizar la localización de los terremotos ocurridos desde 1964 en el archipiélago Fiji, divididos en dos grupos según su intensidad. 

#INICIO -------------------------

rm(list=ls(all=TRUE))

#Elegimos el marco de datos 'quakes' del paquete 'datasets'

library(datasets)

data(quakes)

attach(quakes)

#Dividimos los seísmos en dos grupos (Débil y Fuerte) según superen o no el valor 5.5 de magnitud en la escala Richter

a<-long[mag<5.5];m<-length(a)

b<-lat[mag<5.5]

s<-long[mag>=5.5];n<-length(s)

t<-lat[mag>=5.5]

d<-data.frame(longitud=c(a,s),latitud=c(b,t),Intensidad=c(rep('Débil',m),rep('Fuerte',n)))

#Cargamos el paquete gráfico 'ggplot2'
#En este caso el símbolo + sirve para añadir opciones al gráfico establecido por defecto con 'ggplot()'

library(ggplot2)

ggplot() + 

geom_point(data=d,mapping=aes(x=longitud,y=latitud,fill=Intensidad),size=5,shape=21,color='white')

#FIN -------------------------