Hola,

Otra alternativa a diferentes marcas cada cierta distancia y con otro tipo
de símbolo y color...

###################

T = seq(0, 20*pi, length.out=1000)
#Defino data.frame
dat.df <- data.frame(
                      id=1:length(X)
                     ,X=T*cos(T)
                     ,Y=T*sin(T)
                     ,tipo=rep(19,length(X))
                     ,color=rep("blue", length(X))
                    )

#Pongo marcas cada cierta distancia - Marcas diferentes y colores diferentes
mar.ca <- 15
dat.df$tipo <- ifelse(dat.df$id %% mar.ca == 0, 3, 19 )
dat.df$color<- ifelse(dat.df$id %% mar.ca == 0, "red", "blue" )

plot(dat.df$X, dat.df$Y, pch=dat.df$tipo, col=dat.df$color)

###################

Saludos,
Carlos Ortega
www.qualityexcellence.es



El 14 de marzo de 2013 12:56, Griera <griera@gmail.com> escribió:

> Hola Carlos:
>
> On Thu, 14 Mar 2013 12:46:31 +0100
> "Carlos J. Gil Bellosta " <cgb@datanalytics.com> wrote:
>
> > n <- 20
> >
> > T <- sapply(1:n, function(x) pi * seq(x-1, x, length.out = 5 * x))
> > T <- unlist(T)
> > X = T * cos(T)
> > Y = T * sin(T)
> > plot(X, Y)
> >
> > es una "aproximación (y a trozos)" a la solución de tu problema, que
> > en términos matemáticos precisos es el de encontrar una
> > parametrización por longitud de tu espiral.
> >
> > Lo que hago es muestrear cada ciclo más densamente (directamente
> > proporcional al radio "promedio" de la espiral en dicho ciclo).
>
> Es una solución alternativa a lo que quiero representar, ya que todos los
> círculos están equidistantes. No la entiendo en términos "matemáticos" (me
> lo tengo que revisar con más calma), pero igual me va bien!!
>
> Gracias por la sugerencia. Saludos.
>
> >
> > Un saludo,
> >
> > Carlos J. Gil Bellosta
> > http://www.datanalytics.com
> >
> >
> > El día 14 de marzo de 2013 12:36, Griera <griera@gmail.com> escribió:
> > > On Thu, 14 Mar 2013 12:10:20 +0100
> > > Jorge Ayuso Rejas <jayusor@gmail.com> wrote:
> > >
> > >> Te valdría con:
> > >>
> > >> T = seq(0, 20*pi, length.out=1000)
> > >> X = T * cos(T)
> > >> Y = T * sin(T)
> > >> plot(X, Y, type="l")
> > >>
> > >> Poniendo type="l" te dibuja la curva directamente.
> > >>
> > >> Para los puntos te podría vales algo como esto :
> > >>
> > >> T = seq(0, 20*pi, length.out=1000)
> > >> X = T * cos(T)
> > >> Y = T * sin(T)
> > >> plot(X, Y, type="l")
> > >>
> > >> aux<-seq(0,20*pi)
> > >>
> > >> points(aux * cos(aux), aux * sin(aux),pch=20,col=2)
> > >
> > > Es justo esto, pero los puntos rojos quería que estuviesen
> equidistantes (ahora en la parte interna de la espiral están más juntos y a
> parte de aquí se van separando). Igual esto no es posible?
> > >
> > > Muchas gracias por la ayuda!! Griera.
> > >
> > >>
> > >>
> > >>
> > >> Saludos.
> > >>
> > >>
> > >>
> > >> El 14 de marzo de 2013 11:57, Griera <griera@gmail.com> escribió:
> > >>
> > >> > Hola:
> > >> >
> > >> > Funciona!!! El código es:
> > >> >
> > >> > T = seq(0, 20*pi, length.out=1000)
> > >> > X = T * cos(T)
> > >> > Y = T * sin(T)
> > >> > plot(X, Y, type="n")
> > >> > lines(X, Y)
> > >> >
> > >> > Dibuja la espiral. Ahora voy a ver como le inserto marcas
> equidistantes.
> > >> >
> > >> > Gracias!! Salud!!
> > >> >
> > >> > On Thu, 14 Mar 2013 11:27:35 +0100
> > >> > Jorge Ayuso Rejas <jayusor@gmail.com> wrote:
> > >> >
> > >> > > Mira aquí tienes ejemplos con lo que necesitas:
> > >> > >
> > >> > > http://www.statmethods.net/graphs/line.html
> > >> > >
> > >> > > Saludos.
> > >> > >
> > >> > >
> > >> > > El 14 de marzo de 2013 11:21, Griera <griera@gmail.com> escribió:
> > >> > >
> > >> > > > Hola:
> > >> > > >
> > >> > > > Estoy intentando dibujar una espiral a la que le quería colocar
> unas
> > >> > > > marcas equidistantes.
> > >> > > >
> > >> > > > Estoy utilizando:
> > >> > > >
> > >> > > > T = seq(0, 20*pi, length.out=1000)
> > >> > > > X = T * cos(T)
> > >> > > > Y = T * sin(T)
> > >> > > > plot(X, Y)
> > >> > > >
> > >> > > > que dibuja puntos siguiendo una espiral. Como podría
> convertirlo en una
> > >> > > > espiral continua con marcas equidistantes?
> > >> > > >
> > >> > > > He buscado en r-cran, y no he encontrado nada.
> > >> > > >
> > >> > > > Muchas gracias y saludos!! Griera
> > >> > > >
> > >> > > > _______________________________________________
> > >> > > > R-help-es mailing list
> > >> > > > R-help-es@r-project.org
> > >> > > > https://stat.ethz.ch/mailman/listinfo/r-help-es
> > >> > > >
> > >> > >
> > >> > >
> > >> > >
> > >> > > --
> > >> > > Jorge Ayuso Rejas
> > >> >
> > >>
> > >>
> > >>
> > >> --
> > >> Jorge Ayuso Rejas
> > >
> > > _______________________________________________
> > > R-help-es mailing list
> > > R-help-es@r-project.org
> > > https://stat.ethz.ch/mailman/listinfo/r-help-es
>
> _______________________________________________
> R-help-es mailing list
> R-help-es@r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-help-es
>



-- 
Saludos,
Carlos Ortega
www.qualityexcellence.es

	[[alternative HTML version deleted]]

