Simon and the other developers of the wonderful R.Gui for the Mac,
The following code (anything using mvtnorm) reliably crashes R (under
development) when running from the GUI, but not when running from
the X11 window or under R-Studio.
I reported this several months ago but was told to wait for a newer
release. I believe I have the latest versions of R (unstable), the
R GUI, and mvtnorm.
library(mvtnorm)
sessionInfo()
R Under development (unstable) (2011-07-08 r56336)
Platform: i386-apple-darwin9.8.0/i386 (32-bit)
locale:
[1] en_US.UTF-8/en_US.UTF-8/en_US.UTF-8/C/en_US.UTF-8/en_US.UTF-8
attached base packages:
[1] stats graphics grDevices utils datasets methods base
other attached packages:
[1] mvtnorm_0.9-9991 psych_1.0.99
qmvnorm(.95,sigma=diag(2),tail="both")
*** caught bus error ***
address 0x0, cause 'non-existent physical address'
Traceback:
1: .Fortran("mvtdst", N = as.integer(n), NU = as.integer(df), LOWER
= as.double(lower), UPPER = as.double(upper), INFIN =
as.integer(infin), CORREL = as.double(corrF), DELTA =
as.double(delta), MAXPTS = as.integer(x$maxpts), ABSEPS =
as.double(x$abseps), RELEPS = as.double(x$releps), error =
as.double(error), value = as.double(value), inform =
as.integer(inform), PACKAGE = "mvtnorm")
2: probval.GenzBretz(algorithm, n, df, lower, upper, infin, corr,
corrF, delta)
3: probval(algorithm, n, df, lower, upper, infin, corr, corrF, delta)
4: mvt(lower = lower, upper = upper, df = 0, corr = corr, delta =
mean, algorithm = algorithm, ...)
5: pmvnorm(lower = low, upper = upp, mean = args$mean, corr =
args$corr, sigma = args$sigma, algorithm = algorithm)
6: f(upper, ...)
7: uniroot(pfct, interval = interval)
8: qmvnorm(0.95, sigma = diag(2), tail = "both")
Good luck and thanks for all the great work.
Bill
--
William Revelle http://personality-project.org/revelle.html
Professor http://personality-project.org
Department of Psychology http://www.wcas.northwestern.edu/psych/
Northwestern University http://www.northwestern.edu/
Use R for psychology http://personality-project.org/r
It is 6 minutes to midnight http://www.thebulletin.org
_______________________________________________
R-SIG-Mac mailing list
[email protected]
https://stat.ethz.ch/mailman/listinfo/r-sig-mac
