[R] plot average from aggregated data

[Duncan Mackay]

Hi Luigi

It always helps to know what your data is

str(AVG)
'data.frame':   8 obs. of  2 variables:
 $ stimulation: Factor w/ 8 levels "Unstimulated",..: 1 2 3 4 5 6 7 8
 $ copy       : num  1006 1750 7439 1058 2952 ...
> AVG
   stimulation       copy
1 Unstimulated 1005.80527
2        ESAT6 1750.04691
3        CFP10 7438.63889
4      Rv3615c 1057.95323
5       Rv2654 2952.39095
6       Rv3879   79.69747
7       Rv3873  319.35418
8          PHA 4677.58578

Now if you wanted the x variable as numeric you could do 

plot(1:8, AVG$copy)

if it is not as simple as that and you have groups with multiple values then
this is a starting point

library(lattice)

xyplot(copy ~ 1:dim(AVG)[1], AVG)

Regards

Duncan 

Duncan Mackay
Department of Agronomy and Soil Science
University of New England
Armidale NSW 2351
Email: home: mackay at northnet.com.au


-----Original Message-----
From: r-help-bounces at r-project.org [mailto:r-help-bounces at r-project.org] On
Behalf Of Luigi Marongiu
Sent: Monday, 17 March 2014 04:27
To: r-help at r-project.org
Subject: [R] plot average from aggregated data

Dear all,
I was trying to plot the average of different variables at the same time; I
used the aggregate function and I obtained correctly the averages I wanted.
However when I have tried to plot the data the outlook of the data is a bar
rather than by a symbol; even using the "pch" or the "col" arguments, the
bar have remained unchanged.

Is there a way to modify those bar into symbols or to eliminate them
(perhaps drawing them in white) so that I can superimpose the symbols using
the "points" function (as I have done in the following example)?

best regards
Luigi



my.data<-structure(list(
column_1 = 1:32,
column_2 = structure(c(    1, 2, 3, 4, 5, 6, 7, 8,
        1, 2, 3, 4, 5, 6, 7, 8,
        1, 2, 3, 4, 5, 6, 7, 8,
        1, 2, 3, 4, 5, 6, 7, 8),.Label = c("Unstimulated", "ESAT6",
"CFP10", "Rv3615c", "Rv2654", "Rv3879", "Rv3873", "PHA"),
class = "factor"),
column_3 = structure(c(1, 1, 1, 1, 1, 1, 1, 1, 2, 2, 2, 2, 2, 2, 2, 2, 3,
3, 3,
3, 3, 3, 3, 3, 4, 4, 4, 4, 4, 4, 4, 4), .Label = c("bd", "2m", "1m",
"0.5m"),
class = "factor"),
column_4 = c(71.62097178, 2.892223468, 144.8618621, 40.90079262,
37.0516856,
32.78206822, 72.44424152, 1512.516293, 3164.748333, 1092.634557,
28733.20269,
1207.87783, 729.6090973, 151.8706088, 241.2466141, 9600.963594, 682.865522,
5869.375226, 554.8924896, 2982.759858, 7690.028092, 40.37298719,
942.3549279,
3403.679967, 103.9862389, 35.28562613, 321.5985377, 0.274424607,
3352.874906,
93.76421187, 21.3709382, 4193.183281)),
.Names = c("row", "stimulation", "type", "copy"), row.names = c(NA, -32L),
class = "data.frame")
attach(my.data)

AVG<-aggregate(copy ~ stimulation  , my.data, mean, na.rm = T)

plot(AVG, pch=16, col="white")

points(AVG)

	[[alternative HTML version deleted]]

______________________________________________
R-help at r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.