Hola,
Estás cometiendo un error a la hora de utilizar la función "length".
Lo estás usando con *corchetes*, cuando debieras de utilizarlo con
*paréntesis*. Es una función.
Usando tus propios datos:
> length[do$d]
Error in length[do$d] : objeto de tipo 'builtin' no es
subconjunto> length(do$d)
[1] 18
Si corriges este error, tu código funcionará.
El código puedes optimizarlo y hacerlo más compacto: saber de primeras en
qué casos puedes generar el "loess" y también guardar el resultado en
una
lista en vez de en diferentes vectores....
Pero si merece la pena esto o no, depende de lo que quieras hacer, que no
lo he terminado de entender....
Saludos,
Carlos Ortega
www.qualityexcellence.es
El 29 de septiembre de 2016, 19:43, Manuel Máquez <manuelmx11 en
gmail.com>
escribió:
> Me debo explicar, tengo el siguiente script:
> library(ggplot2)
> da <- read.csv('T20.csv')
> un <- da[da[,2] == 1,]
> do <- da[da[,2] == 2,]
> tr <- da[da[,2] == 3,]
> cu <- da[da[,2] == 4,]
> ci <- da[da[,2] == 5,]
> se <- da[da[,2] == 6,]
> si <- da[da[,2] == 7,]
> if (length[un$d] >= 8) {sua1 <- loess(d ~ n, data = un, span = 0.65,
> degree = 2)
> idx1 <- predict(sua1)
> print(tail(idx1))}
> if (length[do$d] >= 8) {sua2 <- loess(d ~ n, data = do, span = 0.65,
> degree = 2)
> idx2 <- predict(sua2)
> print(tail(idx2))}
> if (length[tr$d] >= 8) {sua3 <- loess(d ~ n, data = tr, span = 0.65,
> degree = 2)
> idx3 <- predict(sua3)
> print(tail(idx3))}
> if (length[cu$d] >= 8) {sua4 <- loess(d ~ n, data = cu, span = 0.65,
> degree = 2)
> idx4 <- predict(sua4)
> print(tail(idx4))}
> if (length[ci$d] >= 8) {sua5 <- loess(d ~ n, data = ci, span = 0.65,
> degree = 2)
> idx5 <- predict(sua5)
> print(tail(idx5))}
> if (length[se$d] >= 8) {sua6 <- loess(d ~ n, data = se, span = 0.65,
> degree = 2)
> idx6 <- predict(sua6)
> print(tail(idx6))}
> if (length[si$d] >= 8) {sua7 <- loess(d ~ n, data = si, span = 0.65,
> degree = 2)
> idx7 <- predict(sua7)
> print(tail(idx7))}
>
> Adjunto archivo T20.csv
> El objeto de incluir if (length['x'$d] >= 8) es eliminar
aquellos grupos
> cuyo numero sea de poca, ninguna ayuda o quizá hasta se genere un error.
> No puedo usar la instrucción antes citada porque se genera un warning
> diciendo que sólo se usará el primer elemento.
> RStudio me da la cantidad de observaciones de cada grupo (frecuencia),
> pero si elimino manualmente el/los grupos, es lento, y por ello intento
> almacenar la cantidad de elementos en alguna variable con la cual se evite
> usar esos datos.
> Agradezco a Javier Marcuzzi su rápida respuesta y espero que en esta
> ocasión haya sido explicito.
> Nuevamente doy las gracias anticipadas por la ayuda que me presten.
>
> *MANOLO MÁRQUEZ P.*
>
> _______________________________________________
> 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
[[alternative HTML version deleted]]