[R-es] Problema con el subset
Jose Manuel Veiga del Baño
chemavb en um.es
Sab Dic 20 09:23:30 CET 2014
Muchas gracias.
Tienes razón era ese el problema. Voy a mirar la opción alternativa.
Un saludo.
Carlos Ortega <cof en qualityexcellence.es> escribió:
> Hola,
>
> El error lo tienes aqui:
>
> #-------------------
> tree <- dfgrupo[ dfgrupo$param == levels( dfinter$param )[ i ], ]
> lines(as.Date(TREE$fecha), TREE$z, type = "b", lwd = 1.5,
> col = colors[i], pch = char[i])
> #-------------------
> Primero seleccionas un subconjunto "tree", pero luego a la hora de
> pintar las lineas, vuelves a escoger "dfgrupo" y une todos los puntos
> del conjunto dfgroup en vez de solamente los del subconjunto "tree".
>
> Como alternativa, aunque estás empezando a usar R, mira esta otra
> forma que evita todos los cálculos intermedios que has utilizado:
>
> #---------------------
> library(lattice)
> parVal <- unique(as.vector(dfgrupo$param))
> colVal <- rainbow(length(parVal))
>
> xyplot(
> z ~ as.Date(fecha)
> ,data=dfgrupo
> ,group=param
> ,type="b"
> ,main="Evaluación de los parámetros"
> ,xlab=list("Fecha del interlaboratorio", font=2, cex=1.1)
> ,ylab=list("Z-Score", font=2, cex=1.1)
> ,scales=list(col="blue", font=2, cex=1.1)
> ,key=simpleKey(parVal, columns=4)
> ,par.settings=simpleTheme(fill=colVal, col=colVal, pch=19)
> )
> #---------------------
>
> Saludos,
> Carlos Ortega
> www.qualityexcellence.es[1]
>
>
> El 17 de diciembre de 2014, 19:59, Jose Manuel Veiga del Baño
> <chemavb en um.es> escribió:
>> _Hola a todos,
>>
>> Agradeceros de antemano vuestro tiempo y paciencia ya que soy un poco
>> novato y tal vez esto sea un poco trivial.
>>
>> Lo que quiero hacer es que me represente en eje de las x las fechas
>> (columna fecha) y los valores de z (columna z) pero de los datos que he
>> filtrado antes en
>> (dfgrupo<-subset(df,df$parametroslaboratorio=="Aflatoxinas ByG")) y que
>> los parámetros iguales (variables de la columna param) se unan entre
sí
>> y tengan un color cada param. En principio estoy usando un código de un
>> ejemplo que funciona muy bien pero lo único que consigo es que se me
>> unan todos los puntos (todos los parm) y que no lo haga por colores.
>>
>> Si alguien se le ocurre algo se lo agradezco._
>>
>> _df <- read.table(file="xxxx.csv", header=T, sep=";",
dec=",")_
>> _df <-na.omit(df)_
>> _dfgrupo<-subset(df,df$parametroslaboratorio=="Aflatoxinas ByG")_
>> _niveles <- as.numeric(dfgrupo$param)_
>> _ntrees <- max(niveles) _
>> _xrange <- range(as.Date(dfgrupo$fecha))_
>> _yrange <- range(dfgrupo$z)_
>> _plot(xrange, yrange, type = "n", xlab = "Fecha del
>> interlaboratorio", _
>> _ ylab = "Z-score")_
>> _colors <- rainbow(ntrees)_
>> _ltipo <- c(1:ntrees)_
>> _char <- seq(18, 18 + ntrees, 1)_
>> _# Añadir las lineas al grafico_
>> _for (i in 1:ntrees) {_
>> _ tree <- dfgrupo[ dfgrupo$param == levels( dfinter$param )[
>> i ], ] _
>> _ lines(as.Date(dfgrupo$fecha), dfgrupo$z, type = "b", lwd =
>> 1.5, _
>> _ col = colors[i], pch = char[i])_
>> _}_
>> _# Añadir un titulo and subtitulo_
>> _title("Evaluacion de los parametros")_
>> _# Añadir la leyenda_
>> _legend(xrange[1], yrange[2], 1:ntrees, cex = 0.8, col = colors,
>> pch = char, _
>> _ lty = ltipo, title = "Parametro")_
>> _ _
>>
>> _Un saludo._ _Dr. José M. Veiga
>> Dpt. Química Agrícola, Geología y Edafología.
>> Universidad de Murcia._
>>
>> ________________________________________________
>> R-help-es mailing list
>> R-help-es en r-project.org
>> https://stat.ethz.ch/mailman/listinfo/r-help-es_
>>
>
> _-- _
> _Saludos,
> Carlos Ortega
> www.qualityexcellence.es[1]_
Vínculos:
---------
[1] http://www.qualityexcellence.es
Dr. José M. Veiga
Dpt. Química Agrícola, Geología y Edafología.
Universidad de Murcia.
[[alternative HTML version deleted]]
Más información sobre la lista de distribución R-help-es