Re: [R] use different symbols for frequency in a plot

From: Marc Schwartz (via MN) <mschwartz_at_mn.rr.com>
Date: Tue 09 Aug 2005 - 06:44:10 EST

On Mon, 2005-08-08 at 11:57 -0700, Kerry Bush wrote:
> suppose I have the following data
>
> x<-c(rep(.1,5),rep(.2,6),rep(.4,10),rep(.5,20))
> y<-c(rep(.5,3),rep(.6,8),rep(1.2,8),rep(2.5,18),rep(3,4))
>
> If I plot(x,y) in R, I will only get seven distinct
> points. What I want to do is to use different symbols
> to show the frequency at each point.
>
> e.g. if the frequncey is between 1 and 5, then I plot
> the point as a circle; if the frequency is between 6
> and 10, then I plot the point as a square; if the
> frequency is above 10, then I plot the point as a
> triangle.
>
> I am not sure how to do this in R. Can anybody help me?

You might want to review this recent post by Deepayan Sarkar:

https://stat.ethz.ch/pipermail/r-help/2005-July/074042.html

with modest modification you can replace his example, which plots the frequencies with:

x <- c(rep(.1,5),rep(.2,6),rep(.4,10),rep(.5,20)) y <- c(rep(.5,3),rep(.6,8),rep(1.2,8),rep(2.5,18),rep(3,4))

temp <- data.frame(x, y)

foo <- subset(as.data.frame(table(temp)), Freq > 0)

> foo

     x y Freq

1  0.1 0.5    3
5  0.1 0.6    2
6  0.2 0.6    6
11 0.4 1.2    8
15 0.4 2.5    2
16 0.5 2.5   16
20 0.5   3    4

# Use cut() to create the bins and specify the plotting symbols # for each bin, which are the 'label' values foo$sym <- with(foo, cut(Freq, c(0, 5, 10, Inf),

                         labels = c(21, 22, 24))) 


# convert 'foo' to all numeric from factors above for plotting foo <- apply(foo, 2, function(x) as.numeric(as.character(x)))

> foo

     x y Freq sym

1  0.1 0.5    3  21
5  0.1 0.6    2  21
6  0.2 0.6    6  22
11 0.4 1.2    8  22
15 0.4 2.5    2  21
16 0.5 2.5   16  24
20 0.5   3    4  21


# Now do the plot. Keep in mind that 'foo' is now # a matrix, rather than a data frame
plot(foo[, "x"], foo[, "y"], pch = foo[, "sym"])

HTH, Marc Schwartz



R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html Received on Tue Aug 09 06:47:14 2005

This archive was generated by hypermail 2.1.8 : Sun 23 Oct 2005 - 15:10:08 EST