[R-es] Script sin resultados (SOLUCIONADO)
Manuel Máquez
manuelmx11 en gmail.com
Mar Mayo 3 02:31:06 CEST 2016
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]]
Más información sobre la lista de distribución R-help-es