[R] tapply and error bars

[Jim Lemon]

Hi Ogbos,
The problem is almost certainly with the data. I get the plot I expect
with the sample data that you first posted, so I know that the code
works. If you try thIs what do you get?

oodf<-read.table(text="S/N      A        B
1    -5  64833
2    -4  95864
3    -3  82322
4    -2  95591
5    -1  69378
6     0  74281
7     1 103261
8     2  92473
9     3  84344
10    4 127415
11    5 123826
12    6 100029
13    7  76205
14    8 105162
15    9 119533
16   10 106490
17   -5  82322
18   -4  95591
19   -3  69378
20   -2  74281
21   -1 103261
22    0  92473
23    1  84344
24    2 127415
25    3 123826
26    4 100029
27    5  76205
28    6 105162
29    7 119533
30    8 106490
31    9 114771
32   10  55593
33   -5  85694
34   -4  65205
35   -3  80995
36   -2  51723
37   -1  62310
38    0  53401
39    1  65677
40    2  76094
41    3  64035
42    4  68290
43    5  73306
44    6  82176
45    7  75566
46    8  89762
47    9  88063
48   10  94395
49   -5  80651
50   -4  81291
51   -3  63702
52   -2  70297
53   -1  64117
54    0  71219
55    1  57354
56    2  62111
57    3  42252
58    4  35454
59    5  33469
60    6  38899
61    7  64981
62    8  85694
63    9  79452
64   10  85216
65   -5  71219
66   -4  57354
67   -3  62111
68   -2  42252
69   -1  35454
70    0  33469
71    1  38899
72    2  64981
73    3  85694
74    4  79452
75    5  85216
76    6  81721
77    7  91231
78    8 107074
79    9 108103
80   10  7576",
header=TRUE)
library(plotrix)
std.error<-function(x) return(sd(x)/(sum(!is.na(x))))
oomean<-as.vector(by(oodf$B,oodf$A,mean))
oose<-as.vector(by(oodf$B,oodf$A,std.error))
plot(-5:10,oomean,type="b",ylim=c(50000,110000),
 xlab="days (epoch is the day of Fd)",ylab="strikes/km2/day")
dispersion(-5:10,oomean,oose)

I get the attached plot;

Jim

On Mon, Jun 25, 2018 at 1:58 AM, Ogbos Okike <giftedlife2014 using gmail.com> wrote:
> Hi Jim
>
> Thanks again for returning to this.
> please not that the line "oomean<-as.vector(by(oodf$B,oodf$A,mean))" was
> omitted (not sure whether deliberate)  after you introduced the standard
> error function.
> When I used it, empty plot window with the correct axes were generated but
> no data was displayed. No error too.
>
> library(plotrix)
> std.error<-function(x) return(sd(x)/(sum(!is.na(x))))
> oose<-as.vector(by(oodf$B,oodf$A,std.error))
> plot(-5:10,oomean,type="b",ylim=c(50000,110000),
>  xlab="days (epoch is the day of Fd)",ylab="strikes/km2/day")
> dispersion(-5:10,oomean,oose)
>
> When I included the line, the same empty graph window was generated but with
> the former error "Error in FUN(X[[1L]], ...) : could not find function
> "FUN""
> library(plotrix)
> std.error<-function(x) return(sd(x)/(sum(!is.na(x))))
> oomean<-as.vector(by(oodf$B,oodf$A,mean))
> oose<-as.vector(by(oodf$B,oodf$A,std.error))
> plot(-5:10,oomean,type="b",ylim=c(50000,110000),
>  xlab="days (epoch is the day of Fd)",ylab="strikes/km2/day")
> dispersion(-5:10,oomean,oose)
>
> I am sure am missing something but can't place it. Please have a look again
> to track my mistake.
>
> Warmest regards
> Ogbos
>

-------------- next part --------------
A non-text attachment was scrubbed...
Name: ooplot.pdf
Type: application/pdf
Size: 5776 bytes
Desc: not available
URL: <https://stat.ethz.ch/pipermail/r-help/attachments/20180625/fe014680/attachment-0002.pdf>