zeros in table() output in version 1.3.1 - R-help

Tue, Oct 23, 2001 12:18 PM #

Are the following differences in table() output
between versions 1.3.0 and 1.3.1 a feature, a bug
or a problem with my implementation?

(I am running under Solaris 2.7 SPARC, built using the
Sun C and FORTRAN compilers.)

Version 1.3.1:

1
  1 1

Version 1.3.0:

0 1
  0 1 0
  1 0 1


-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-
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
_._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._

Peter Dalgaard

Tue, Oct 23, 2001 12:35 PM #

Steve Shiboski <steve at biostat.ucsf.edu> writes:

That'll be a bug, but it is not there when R is built with GNU tools:

0 1
  0 1 0
  1 0 1

_                   
platform sparc-sun-solaris2.7
arch     sparc               
os       solaris2.7          
system   sparc, solaris2.7   
status                       
major    1                   
minor    3.1                 
year     2001                
month    08                  
day      31                  
language R

O__  ---- Peter Dalgaard             Blegdamsvej 3  
  c/ /'_ --- Dept. of Biostatistics     2200 Cph. N   
 (*) \(*) -- University of Copenhagen   Denmark      Ph: (+45) 35327918
~~~~~~~~~~ - (p.dalgaard at biostat.ku.dk)             FAX: (+45) 35327907
-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-
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
_._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._

Mark Myatt

Wed, Oct 24, 2001 1:54 AM #

Steve Shiboski <steve at biostat.ucsf.edu> writes:

Odd. I don't get this behaviour on my Windows box. I think that,
perhaps, it is down to how the binary was compiled. Did you use the Sun
tools or the GNU tools?

Mark


--
Mark Myatt


-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-
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
_._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._

Brian Ripley

Thu, Oct 25, 2001 1:56 PM #

On Wed, 24 Oct 2001, Mark Myatt wrote:

I don't see it on Solaris 2.7 built with Sun's C and Fortran using either
Sun or GNU tools.   So it seems definitely a priblem with one specific
installation.

Brian D. Ripley,                  ripley at stats.ox.ac.uk
Professor of Applied Statistics,  http://www.stats.ox.ac.uk/~ripley/
University of Oxford,             Tel:  +44 1865 272861 (self)
1 South Parks Road,                     +44 1865 272860 (secr)
Oxford OX1 3TG, UK                Fax:  +44 1865 272595

-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-
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
_._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._

Stephen Shiboski

Fri, Oct 26, 2001 9:34 AM #

It turns out to be associated with the compilation flags
I put in the config.site file for the Sun C compiler.
Recompiling with the default flags chosen by R eliminates
the problem.

I put the following flags in config.site before building:

 CFLAGS=-fast 
 SHLIB_CFLAGS=-lsunmath 
 CPICFLAGS=-xcode=pic32

Not sure which caused the problem with table(), but
R's defaults eliminate -fast, -xcode=pic32 and -lsunmath.

Thanks for the responses.

-S

Stephen Shiboski  -------------------  steve at biostat.ucsf.edu
Division of Biostatistics                 voice: 415-476-0533
University of California, San Francisco     fax: 415-476-6014
500 Parnassus Avenue, MU 420-W;  San Francisco, CA 94143-0560

On Thu, Oct 25, 2001 at 09:56:40PM +0100, Prof Brian D Ripley wrote:
> On Wed, 24 Oct 2001, Mark Myatt wrote:
> 
> > Steve Shiboski <steve at biostat.ucsf.edu> writes:
> > >Are the following differences in table() output
> > >between versions 1.3.0 and 1.3.1 a feature, a bug
> > >or a problem with my implementation?
> > >
> > >(I am running under Solaris 2.7 SPARC, built using the
> > >Sun C and FORTRAN compilers.)
> > >
> > >Version 1.3.1:
> > >
> > >> table(c(0,1),c(0,1))
> > >
> > >    1
> > >  1 1
> > >
> > >Version 1.3.0:
> > >
> > >> table(c(0,1),c(0,1))
> > >
> > >    0 1
> > >  0 1 0
> > >  1 0 1
> >
> > Odd. I don't get this behaviour on my Windows box. I think that,
> > perhaps, it is down to how the binary was compiled. Did you use the Sun
> > tools or the GNU tools?
> 
> I don't see it on Solaris 2.7 built with Sun's C and Fortran using either
> Sun or GNU tools.   So it seems definitely a priblem with one specific
> installation.
> 
> -- 
> Brian D. Ripley,                  ripley at stats.ox.ac.uk
> Professor of Applied Statistics,  http://www.stats.ox.ac.uk/~ripley/
> University of Oxford,             Tel:  +44 1865 272861 (self)
> 1 South Parks Road,                     +44 1865 272860 (secr)
> Oxford OX1 3TG, UK                Fax:  +44 1865 272595
> 

-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-
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
_._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._

Brian Ripley

Fri, Oct 26, 2001 9:52 AM #

On Fri, 26 Oct 2001, Steve Shiboski wrote:

-fast disables IEEE 754 arithmetic. For me R fails make check with that.
(You did run make check, didn't you ...?)

I use -xO5 -xlibmil -dalign (some libraries need -dalign), and f90 which
picks up -lsunmath too.

Brian D. Ripley,                  ripley at stats.ox.ac.uk
Professor of Applied Statistics,  http://www.stats.ox.ac.uk/~ripley/
University of Oxford,             Tel:  +44 1865 272861 (self)
1 South Parks Road,                     +44 1865 272860 (secr)
Oxford OX1 3TG, UK                Fax:  +44 1865 272595

-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-
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
_._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._

Wolfgang Viechtbauer

Fri, Oct 26, 2001 1:04 PM #

I hope somebody can help me out with this. Let's say I have a matrix with
i rows and j columns and I want to find the index of the row AND column of
the maximum value in the matrix. I tried to play around with the which()
command, but couldn't find an easy solution that didn't involve loops and
various other inelegant approaches. I hope that there is an easy way of
doing this. Thanks in advance for any suggestions!

------------------------------------------------------------------------
Wolfgang Viechtbauer                       wviechtb at s.psych.uiuc.edu
Department of Quantitative Psychology      viechtba at students.uiuc.edu

University of Illinois, Champaign-Urbana   [Virus  based  on  an honors-
Room 409, Psychology Building              system:  Forward  this  to 10
603 E. Daniel Street                       random  people  and  delete a
Champaign, IL 61820                        bunch of your files. Thanks.]

-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-
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
_._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._

Ben Bolker

Fri, Oct 26, 2001 1:15 PM #

which(x,arr.ind=TRUE)

On Fri, 26 Oct 2001, Wolfgang Viechtbauer wrote:

318 Carr Hall                                bolker at zoo.ufl.edu
Zoology Department, University of Florida    http://www.zoo.ufl.edu/bolker
Box 118525                                   (ph)  352-392-5697
Gainesville, FL 32611-8525                   (fax) 352-392-3704

-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-
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
_._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._

Kjetil Halvorsen

Mon, Oct 29, 2001 8:31 AM #

[,1] [,2] [,3] [,4] [,5]
[1,]    1    4    3    2    4
[2,]    1    4    4    2    1
[3,]    4    3    3    2    3
[4,]    3    2    2    2    1
[5,]    2    2    3    2    2

[1] 4

row col
[1,]   3   1
[2,]   1   2
[3,]   2   2
[4,]   2   3
[5,]   1   5


Kjetil Halvorsen

Wolfgang Viechtbauer wrote:

-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-
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
_._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._