1

Estoy realizando un Analisis de componentes principales (PCA) en R. Para ver si una especie se segrega en diferentes localidades a partir de las variables obtenida en cada localidad. Pero para poder ilustrar eso mejor tengo que agregar los poligonos de ordenacion pero el codigo que ofrece el paquete vegan no es muy explicativo. Este s el codigo que uso:

data(dune)
mod <- rda(dune, scale = TRUE)
biplot(mod, scaling = "symmetric")

## different type for species and site scores
biplot(mod, scaling = "symmetric", type = c("text", "points"))

El grafico que obtengo es este: Imagen que obtengo

Y esta es la imagen q la que se quiere llegar: Esto es lo que se desea. Tomado de https://www.benjaminbell.co.uk/2018/02/principal-components-analysis-pca-in-r.html

COMO PUEDO PONER COLORES A LOS PUNTOS INTERIORES EN CORRESPONDENCIA A LAS ELIPSES Muchas gracias

1 respuesta 1

2

Si quieres usar Vegan, para lo que buscas, tienes ordiellipse(), dónde claramente necesitaras una variable que determine el área. Tomando iris como ejemplo, emulamos tu área con la variable Species:

mod <-  rda(iris[,-5])

biplot(mod,
       display = c("sites", 
                   "species"),
       type = c("text",
                "points"),
       scaling = "symmetric")

ordiellipse(mod,  
            group = iris$Species, 
            kind = "ehull", 
            conf = 0.95,
            draw = "polygon",
            alpha = 20,
            col = c("red", "blue", "green"), 
            border = c("red", "blue", "green"), 
            lwd = 2)

legend("bottomright",
       col = c("red", "blue", "green"), 
       lty = 1,
       legend = levels(iris$Species))

introducir la descripción de la imagen aquí

3
  • Muchas gracias ya pude poner las elipses. una pregunta mas: como puedo poner el nombre a cada punto para saber q especie es la que esta en cada punto eso es importante. Commented el 17 dic. 2020 a las 21:33
  • 1
    El data.frame necesita rownames con la etiquetas a motrar, biplot usa esta información automáticamente, siempre que no sea un numero demasiado grande Commented el 17 dic. 2020 a las 23:26
  • Como puedo agregar numeros a la distribucion de puntos? Eso me permite saber si hay alguna especie q no se est'a segregando bien. Commented el 23 dic. 2020 a las 21:15

Tu Respuesta

By clicking “Publica tu respuesta”, you agree to our terms of service and acknowledge you have read our privacy policy.

¿No es la respuesta que buscas? Examina otras preguntas con la etiqueta o formula tu propia pregunta.