8.3: Trazado...

Al inicio, verificar visualmente la distribución de los datos. Hacer histograma:

Código\(\PageIndex{1}\) (R):

data <- read.table("data/bugs.txt", h=TRUE)
hist(data$WEIGHT, breaks=3)

(Para ver un histograma más detallado, aumente el número de descansos).

Si para el histograma desea dividir los datos de la manera específica (por ejemplo, por 20 unidades, comenzando desde 0 y terminando en 100), escriba:

Código\(\PageIndex{2}\) (R):

data <- read.table("data/bugs.txt", h=TRUE)
hist(data$WEIGHT, breaks=c(seq(0, 100, 20))

Las parcelas de caja muestran valores atípicos, máximo, mínimo, rango de cuartil y mediana para cualquier variable de medición:

Código\(\PageIndex{3}\) (R):

data <- read.table("data/bugs.txt", h=TRUE)
boxplot(data$LENGTH)

... ahora para machos y hembras por separado, usando interfaz de fórmula:

Código\(\PageIndex{4}\) (R):

data <- read.table("data/bugs.txt", h=TRUE)
boxplot(data$LENGTH ~ data$SEX)

Hay dos comandos que juntos ayudan a verificar la normalidad del personaje:

Código\(\PageIndex{5}\) (R):

data <- read.table("data/bugs.txt", h=TRUE)
qqnorm(data$WEIGHT); qqline(data$WEIGHT)

(Estos dos comandos separados trabajan juntos para hacer una sola trama, es por eso que usamos punto y coma. Cuantos más puntos en la gráfica resultante se desvíen de la línea, más no normales son los datos.)

Hacer diagrama de dispersión donde todos los errores representados con pequeños círculos. El eje X representará la longitud mientras que el eje Y: el peso:

Código\(\PageIndex{6}\) (R):

data <- read.table("data/bugs.txt", h=TRUE)
plot(data$LENGTH, data$WEIGHT, type="p")

(type="p” es el valor por defecto para plot (), por lo tanto suele omitirse.)

Es posible cambiar el tamaño de los puntos variando el parámetro cex. Comparar

Código\(\PageIndex{7}\) (R):

data <- read.table("data/bugs.txt", h=TRUE)
plot(data$LENGTH, data$WEIGHT, type="p", cex=0.5)

con

Código\(\PageIndex{8}\) (R):

data <- read.table("data/bugs.txt", h=TRUE)
plot(data$LENGTH, data$WEIGHT, type="p", cex=2)

¿Cómo comparar? La mejor manera es tener más de una ventana gráfica en el escritorio. Para iniciar nueva ventana, escriba dev.new ().

También es posible cambiar el tipo de símbolo de trazado. La figura\(\PageIndex{1}\) muestra sus números. Si quieres esta tabla en la computadora, puedes ejecutar:

Código\(\PageIndex{9}\) (R):

Ex.points() # gmoon.r

Figura Tipos de\(\PageIndex{1}\) puntos en gráficas estándar R. Para los tipos 21—25, es posible especificar diferentes colores de fondo y marco.

Para obtener ejemplos gráficos similares sobre tipos de líneas, colores predeterminados, caras de fuente y tipos de trazado, cargue gmoon.r\(^{[1]}\) y ejecute:

Código\(\PageIndex{10}\) (R):

Ex.lines() # gmoon.r
Ex.cols() # gmoon.r
Ex.fonts() # gmoon.r
Ex.types() # gmoon.r

Utilice el símbolo 2 (triángulo vacío):

Código\(\PageIndex{11}\) (R):

data <- read.table("data/bugs.txt", h=TRUE)
plot(data$LENGTH, data$WEIGHT, type="p", pch=2)

Use códigos de texto (0/1) para el SEXO en lugar del símbolo gráfico:

Código\(\PageIndex{12}\) (R):

data <- read.table("data/bugs.txt", h=TRUE)
plot(data$LENGTH, data$WEIGHT, type="n")
text(data$LENGTH, data$WEIGHT, labels=data$SEX)

(Aquí ambos comandos hacen una trama juntos. El primero traza el campo vacío con ejes, el segundo agrega allí símbolos de texto).

La misma trama es posible hacer con el comando único, pero esto solo funciona para etiquetas de una letra:

Código\(\PageIndex{13}\) (R):

data <- read.table("data/bugs.txt", h=TRUE)
plot(data$LENGTH, data$WEIGHT, pch=as.character(data$SEX))

Si queremos que estos números tengan diferentes colores, escriba:

Código\(\PageIndex{14}\) (R):

data <- read.table("data/bugs.txt", h=TRUE)
plot(data$LENGTH, data$WEIGHT, type="n")
text(data$LENGTH, data$WEIGHT, labels=data$SEX, col=data$SEX+1)

(Nuevamente, ambos comandos hacen una trama. Agregamos +1 porque de lo contrario los signos femeninos serían de 0 color, que es “invisible”.)

Diferentes símbolos para machos y hembras:

Código\(\PageIndex{15}\) (R):

data <- read.table("data/bugs.txt", h=TRUE)
plot(data$LENGTH, data$WEIGHT, type="n")
points(data$LENGTH, data$WEIGHT, pch=data$SEX)

La variante más complicada: usa símbolos de fuentes Hershey\(^{[2]}\) que son internas en R (Figura\(\PageIndex{2}\)):

Código\(\PageIndex{16}\) (R):

data <- read.table("data/bugs.txt", h=TRUE)
plot(data$LENGTH^3, data$WEIGHT, type="n", xlab=expression("Volume (cm"^3*")"), ylab="Weight")
text(data$LENGTH^3, data$WEIGHT, labels=ifelse(data$SEX, "\MA", "\VE"), vfont=c("serif","plain"), cex=1.5)

(Obsérvese también cómo se empleó expression () para hacer etiquetas de ejes avanzados. Dentro de expression (), diferentes partes se unen con estrella *. Para saber más, ¿correr? plotmath.)

Figura\(\PageIndex{2}\) Distribución de bichos machos y hembras por tamaño y peso (fuentes Hershey utilizadas)

Podemos pintar símbolos con diferentes colores:

Código\(\PageIndex{17}\) (R):

data <- read.table("data/bugs.txt", h=TRUE)
plot(data$LENGTH, data$WEIGHT, type="n")
points(data$LENGTH, data$WEIGHT, pch=data$SEX*3, col=data$SEX+1)

Por último, es bueno tener una leyenda:

Código\(\PageIndex{18}\) (R):

legend("bottomright", legend=c("male", "female"), pch=c(0, 3), col=1:2)

Y luego guardar la trama como archivo PDF:

Código\(\PageIndex{19}\) (R):

dev.copy(pdf, "graph.pdf")
dev.off()

¡Atención! Guardando en el archivo externo, ¡nunca olvides escribir dev.off ()!

Si no desea ninguna de las etiquetas de eje y principal, inserte las opciones main= "”, xlab= "”, ylab= "” en su comando plot ().

También hay una mejor manera de guardar parcelas porque no se duplica a la pantalla y por lo tanto funciona mejor en los scripts R:

Código\(\PageIndex{20}\) (R):

data <- read.table("data/bugs.txt", h=TRUE)
pdf("graph.pdf")
plot(data$LENGTH, data$WEIGHT, type="n")
points(data$LENGTH, data$WEIGHT, pch=data$SEX*3, col=data$SEX+1)
legend("bottomright", legend=c("male", "female"), pch=c(0, 3), col=1:2)
dev.off()

(Tenga en cuenta aquí que R no emite ninguna advertencia si el archivo con el mismo nombre ya existe en el disco, simplemente lo borra y guarda el nuevo. ¡Ten cuidado!)

Referencias

1. Con fuente de comando (“ashipunov.info/r/gmoon.r”).

2. Para saber qué símbolos están disponibles, ejecute demo (Hershey).