[R] Possible bug in apply()

[Clive Jenkins]

Thanks for confirming this on i586 linux-mandrake R-1.1.1 .

Sorry about the typo: the file is named 'distr.dat', but in the code all
references to it (except the first) are shown as 'dist.dat'. I re-named
the file while reducing the example. And yes, I forgot to mention
load(ctest).

Clive Jenkins.

mikalzet at libero.it wrote:
> 
> At my first attempt I got an early syntax error code:
> very silly of me - a misplaced bracket.
> 
> On Sat, 7 Oct 2000, Clive Jenkins wrote:
> 
> > Error in names(x) == ans.names : comparison (1) is possible only for
> > vector types
> 
> If I load c.test library before running your code and check to see it's
> 'distr.dat' everytime instead of 'dist.dat' I get the same reply.
> 
> > If we look at the structure of the value returned by shap() for the
> > first sample, we see that W has a "names" attribute whereas p has
> > not (or its "names" attribute is the empty string):
> >
> > > str(shap(dist.dat[1,]))
> >  Named num [1:2] 0.887519 0.000622
> >  - attr(*, "names")= chr [1:2] "W" ""
> 
> Same reply
> 
> > but in the case of the second sample (all NA except one), shapiro.test()
> > was not called, and the preset default value c(NA, NA) returned:
> >
> > > str(shap(dist.dat[2,]))
> >  logi [1:2] NA NA
> 
> Same reply
> 
> > Looking directly at the result of shapiro.test(), we get:
> >
> > > str(shapiro.test(dist.dat[1,]))
> > List of 4
> >  $ statistic: Named num 0.888
> >   ..- attr(*, "names")= chr "W"
> >  $ p.value  : num 0.000622
> >  $ method   : chr "Shapiro-Wilk normality test"
> >  $ data.name: chr "dist.dat[1, ]"
> >  - attr(*, "class")= chr "htest"
> 
> Idem
> 
> > If I remove the "names" attribute from W by changing the line
> > of code marked "# Problem line" to
> >
> >         result.W <- as.numeric(shap.res$statistic)
> >
> > the error disappears.
> 
> Right. I get the following output:
> 
>         [,1]            [,2]
> [1,]    0.8875191054    NA
> [2,]    0.0006219975    NA
> 
> If it is a bug it's replicable on i586 linux-mandrake R-1.1.1 .
> 
> --
> Dr. Michele Alzetta
-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-
r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html
Send "info", "help", or "[un]subscribe"
(in the "body", not the subject !)  To: r-help-request at stat.math.ethz.ch
_._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._