I don't recompile R but use the command

ln -sfn libRblas.0.dylib libRblas.dylib

to use the R BLAS

The problem is still there.

Huang

On Thu, Oct 21, 2010 at 9:14 PM, Marc Schwartz <marc_schwartz@me.com> wrote:

> Hi all,
>
> To the extent that it may be helpful here and I can do more if need be, I
> built 32 bit R 2.12.0 patched on Snow Leopard (10.6.4), using the R BLAS
> rather than Apple's veclib. This is on an early 2009 17" MBP with a 2.93 Ghz
> Core 2 Duo (MacBookPro5,2) and 4Gb of RAM.
>
> Based upon Doug's comment in this thread that the issue may be related to
> the use of Apple's veclib BLAS, as opposed to R's reference BLAS, I ran some
> tests.
>
> My config includes:
>
>  --without-blas --without-lapack
>
> just to be sure that the above is the correct invocation, based upon what I
> found online.
>
> Using this build, with all CRAN packages freshly installed using this
> build, I ran the example used here with lme4 0.999375-35. I get:
>
> library(lme4)
> y <- (1:20)*pi; x <- (1:20)^2;group <- gl(2,10)
> M2. <- lmer (y ~ 1 + x + (1 + x | group))
> M2 <- lmer (y ~ x + ( x | group))
>
> > identical(fixef(M2), fixef(M2.))
> [1] TRUE
>
>
>
> I then created a function so that I could use replicate() to run this test
> a "larger" number of times:
>
> testlme4 <- function()
> {
>   y <- (1:20)*pi; x <- (1:20)^2;group <- gl(2,10)
>  M2. <- lmer (y ~ 1 + x + (1 + x | group))
>  M2 <- lmer (y ~ x + ( x | group))
>  identical(fixef(M2), fixef(M2.))
> }
>
>
> RES <- replicate(1000, testlme4())
>
> > all(RES)
> [1] TRUE
>
> > table(RES)
> RES
> TRUE
> 1000
>
> Does the example need to be run a "very large" number of times to be sure
> that it does not fail, or is the above a reasonable indication that the use
> of R's BLAS is a more appropriate default option for R on OSX?  If I am not
> mistaken (and somebody correct me if wrong), R's BLAS is the default on
> Windows and Linux (from my recollections on Fedora). Why should OSX be
> different in that regard?
>
> Also, as an aside to Federico, I use 32 bit R on OSX largely because I have
> to interact with an Oracle server via RODBC. The only ODBC drivers available
> for Oracle on OSX are 32 bit and they are not compatible with 64 bit R. It
> would be rather cumbersome when running reports (via Sweave) to first
> extract the data in 32 bit R and then switch to 64 bit R to run the reports.
> I can run it all in a single step using 32 bit R. I also do not have a need
> for the larger memory address space afforded by 64 bit R.
>
> HTH,
>
> Marc Schwartz
>
>

	[[alternative HTML version deleted]]


