[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