[R-es] gráfico

Estimados
Quisiéramos que la salida final plot (z) fuera de tres colores, uno para
cada especie, no encuentro el comando para ello. Apreciaría su ayuda
Saludos
Jose


rm(list = ls())
library(MASS)
data(iris)
datos <- data.frame(iris[,1:4],clase=as.vector(iris[,5]))
head(datos)
especies <- iris[,5] # La columna cinco se refiere a las especies(clase)
z <- lda(clase~.,datos)
predict(z,iris[1:4, ])$class #Para clasificar las observaciones
Tabla <- table(iris[,5],predict(z,iris[,1:4])$class)
Tabla
(pred.class<-predict(z)$class) # Para clasificar las observaciones
plot(z)



--
Este mensaje le ha llegado mediante el servicio de correo electronico que
ofrece Infomed para respaldar el cumplimiento de las misiones del Sistema
Nacional de Salud. La persona que envia este correo asume el compromiso de
usar el servicio a tales fines y cumplir con las regulaciones establecidas

Infomed: http://www.sld.cu/

Usa

plot(z, col = as.numeric(pred.class))

Saludos,
Jorge.-


2017-09-03 8:29 GMT-05:00 <jbetancourt en iscmc.cmw.sld.cu>:

Estimados
Quisiéramos que la salida final plot (z) fuera de tres colores, uno
para
cada especie, no encuentro el comando para ello. Apreciaría su ayuda
Saludos
Jose


rm(list = ls())
library(MASS)
data(iris)
datos <- data.frame(iris[,1:4],clase=as.vector(iris[,5]))
head(datos)
especies <- iris[,5] # La columna cinco se refiere a las especies(clase)
z <- lda(clase~.,datos)
predict(z,iris[1:4, ])$class #Para clasificar las observaciones
Tabla <- table(iris[,5],predict(z,iris[,1:4])$class)
Tabla
(pred.class<-predict(z)$class) # Para clasificar las observaciones
plot(z)



--
Este mensaje le ha llegado mediante el servicio de correo electronico
que
ofrece Infomed para respaldar el cumplimiento de las misiones del
Sistema
Nacional de Salud. La persona que envia este correo asume el compromiso
de
usar el servicio a tales fines y cumplir con las regulaciones
establecidas

Infomed: http://www.sld.cu/