Carlos: Muchísimas gracias nuevamente, en mi anterior había dicho que te comentaría, y, ahora lo hago. Quería 'comerme el pastel con una mordida', era necesario tomar un vector y 'exprimirlo' quedando finalmente así: library(ggplot2) ogl <- read.csv('T06.csv') re <- ogl[ogl[,3] <= 9,] # Gracias por tu post *Carlos* qu <- ogl[ogl[,3] >= 10,] sua <- loess(d ~ n, data = re, span = 0.65, degree = 2) idx <- predict(sua) sv <- loess(d ~ n, data = qu, span = 0.65, degree = 2) ix <- predict(sv) ogl$det <-ifelse(ogl$d >= 10, 'Myr', 'Nor') gr <- ggplot(data = ogl, aes(x = n, y = d, col = det)) + geom_point(aes(shape = det))+ geom_smooth(span = 0.65) + scale_colour_manual(values = c('red', 'blue')) + theme(legend.position = c(0.17, 0.85)) + ggtitle('Norm, c/y>=10; Curvas Sep e Iterv d Conf [06]') print (gr) # tail(idx,20) ix también # sink( ? ?.txt) para sacar valores de idx e ix
Bueno pues, gracias y buenos días. Atentamente; *MANOLO MÁRQUEZ P.* [[alternative HTML version deleted]] _______________________________________________ R-help-es mailing list R-help-es@r-project.org https://stat.ethz.ch/mailman/listinfo/r-help-es