thanks david,
yes, you are right PART of the confusion is because of what you mentioned
(sorry for that) but going back to my own data this is JUST PART of the
problem…
...see my reproducible example
teq<-structure(list(site = structure(c(4L, 2L, 2L, 4L, 2L, 4L, 4L,
3L, 1L, 3L, 1L, 1L, 3L, 4L, 5L, 4L, 3L, 1L, 3L, 1L, 4L, 4L, 2L,
4L, 4L, 2L, 1L, 3L, 4L, 4L, 2L, 4L, 4L, 2L, 4L, 4L, 2L, 4L, 4L,
2L, 4L, 2L, 4L, 2L, 4L, 2L, 4L, 2L, 2L, 2L, 2L, 4L), .Label = c("A",
"B", "C", "D", "E"), class = "factor"), lower_b_i = c(0.302,
0.956, 0.72, 1.21, 0.887, 0.728, 1.294, 20.493, 0.902, 0.031,
0.468, 2.318, 4.795, 89.581, 4.59, 3.366, 32.786, 5.506, 61.038,
1.485, 42.25, 31.279, 21.375, 96.62, 36.099, 30.173, 0.532, 1.005,
0.697, 6.248, 2.121, 0.576, 0.335, 0.442, 22.704, 18.244, 16.53,
626.664, 42.714, 59.87, 9.983, 110.158, 4.763, 3.363, 2.608,
2.41, 23.904, 73.281, 3618.353, 109.663, 21.797, 3.39), date_sampl_time_recs
= structure(c(11L,
10L, 9L, 8L, 12L, 12L, 13L, 7L, 6L, 14L, 14L, 4L, 3L, 5L, 2L,
1L, 18L, 19L, 24L, 23L, 17L, 15L, 16L, 20L, 21L, 22L, 26L, 25L,
29L, 28L, 27L, 32L, 31L, 30L, 34L, 33L, 33L, 35L, 36L, 37L, 38L,
39L, 40L, 41L, 42L, 43L, 46L, 47L, 48L, 49L, 45L, 44L), .Label =
c("2008-07-15, 2 h, 16/17",
"2008-07-16, 2 h, 14/17", "2008-12-03, 4 h, 13/17", "2008-12-03, 4 h,
15/17",
"2009-01-29, 24 h, 0/17", "2009-03-17, 24 h, 14/17", "2009-03-17, 24 h,
8/17",
"2009-04-17, 135 h, 13/17", "2009-04-20, 96 h, 14/17", "2009-04-21, 24 h,
13/17",
"2009-04-21, 24 h, 14/17", "2009-07-17, 168 h, 13/17", "2009-07-21, 24 h,
12/17",
"2009-08-18, 24 h, 16/17", "2009-10-27, 168 h, 3/17", "2009-10-27, 168 h,
5/17",
"2009-11-01, 24 h, 4/17", "2009-12-15, 24 h, 6/17", "2009-12-15, 24 h,
9/17",
"2010-01-05, 12 h, 3/17", "2010-01-15, 168 h, 3/17", "2010-01-15, 168 h,
5/17",
"2010-01-19, 24 h, 12/17", "2010-01-19, 24 h, 2/17", "2010-03-30, 24 h,
13/17",
"2010-03-30, 24 h, 14/17", "2010-04-22, 168 h, 11/17", "2010-04-22, 168 h,
9/17",
"2010-04-22, 24 h, 14/17", "2010-07-22, 168 h, 14/17", "2010-07-22, 168 h,
15/17",
"2010-07-22, 24 h, 15/17", "2010-10-26, 168 h, 7/17", "2010-10-26, 24 h,
8/17",
"2011-01-05, 12 h, 0/17", "2011-01-11, 168 h, 0/17", "2011-01-19, 168 h,
0/17",
"2011-03-15, 168 h, 2/17", "2011-03-23, 168 h, 0/17", "2011-05-10, 168 h,
5/17",
"2011-05-18, 168 h, 7/17", "2011-07-05, 168 h, 7/17", "2011-07-13, 168 h,
7/17",
"2011-09-13, 168 h, 6/17", "2011-09-21, 168 h, 1/17", "2011-11-08, 168 h,
0/17",
"2011-11-16, 168 h, 0/17", "2012-01-05, 12 h, 1/17", "2012-01-09, 12 h,
2/17"
), class = "factor")), .Names = c("site", "lower_b_i",
"date_sampl_time_recs"
), class = "data.frame", row.names = c(NA, -52L))
...and my code
dotplot(date_sampl_time_recs ~ lower_b_i | site, data=teq,
scales=list(x=list(log=TRUE)),
xscale.components = xscale.components.logpower,
layout=c(5,1),
panel = function(x,y,...) {
panel.grid(h=53, v=-1, lty="dotted", col="gray")
panel.dotplot(x,y,...)
medians <- median(x)
panel.abline(v=medians, col.line="red", lty="dotted")
means<-mean(x)
panel.abline(v=means, col.line="blue", lty="dotted")
}
)
attach(teq)
#check median values
tapply(lower_b_i,site,median)
#check mean values
tapply(lower_b_i,site,mean)
detach(teq)
...and now the question is why the plotting of means (but not medians) is
wrong (check against results of tapply)?
thanks
--
View this message in context:
http://r.789695.n4.nabble.com/how-to-add-a-vertical-line-for-each-panel-in-a-lattice-dotplot-with-log-scale-tp4632513p4632760.html
Sent from the R help mailing list archive at Nabble.com.
______________________________________________
R-help@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.
