similar to: Re: [R-SIG-Mac] NaN and linear algebra

On 21/03/2005, at 10:09 PM, David Firth wrote: > I am sorry that I wasn't clear. All that I meant was that *this* > problem can result in different behaviour in "ordinary" statistical > applications. For example, if the objective function in a call to > optim() involves calling one of these linear algebra routines, the > result may be NaN (on systems other than Mac

Dear R-Help List, I'm trying to implement Firth's (1993) bias correction for log-linear models. Firth (1993) states that such a correction can be implemented by supplementing the data with a function of h_i, the diagonals from the hat matrix, but doesn't provide further details. I can see that for a saturated log-linear model, h_i=1 for all i, hence one just adds 1/2 to each count,

Is there a reliable way to link to a vignette or other overview document, which has been included in the "inst/doc" directory of a source package, within a package help file in Rd format? or maybe link to an index of such documents? David Professor David Firth Dept of Statistics University of Warwick Coventry CV4 7AL United Kingdom Voice: +44 (0)247 657 2581 Fax: +44 (0)247 652

This bug is not quite fixed - the example from my original report now = works using R-2.2.1, but plot(Uniform, 6) does not. The bug is due to if (show[6]) { ymx <- max(cook, na.rm =3D TRUE) * 1.025 g <- hatval/(1 - hatval) # Potential division by zero here # plot(g, cook, xlim =3D c(0, max(g)), ylim =3D c(0, ymx),=20 main =3D main, xlab =3D

This concerns behaviour which has been noted previously, for example see http://tolstoy.newcastle.edu.au/R/help/05/03/0844.html http://tolstoy.newcastle.edu.au/R/help/03b/6873.html and especially the thread http://tolstoy.newcastle.edu.au/R/help/04/08/3025.html I find that it affects both my Linux and Mac OS X setups. It is a feature which has caused some (at least John Fox in the

Full_Name: David Firth Version: 2.0.0 OS: Mac OS 10.3.5 Submission from: (NULL) (81.178.234.156) When running R in R Cocoa GUI 1.0 (v2004-10-14), I have > Sys.getenv()[["R_LIBS"]] [1] "/Users/david/Library/R" > .libPaths() [1] "/Users/david/Library/R/library" [2] "/Library/Frameworks/R.framework/Resources/library" The pasting of

Ashutosh, Thanks for the repy. Related earlier topic on this appears in the review of the SVML patch (@mmasten). Adding few names from there. https://reviews.llvm.org/D19544 There, I see Hal's review comment "let's start only with the directly-legal calls". Apparently, what we have right now in the trunk is "not legal enough". I'll work on the patch to stop

It may not be a full solution for the problems you're trying to solve, but I don't know why adding to include/llvm/CodeGen/RuntimeLibcalls.def is a problem in itself. Certainly, it's a mess that could be organized, especially so we're not repeating everything for each data type as we do right now. So yes, I think that would allow us to remove the VecLib mappings because we are

Hi All, When I tried to do logistic regression (with high maximum number of iterations) I got the following warning message Warning message: fitted probabilities numerically 0 or 1 occurred in: (if (is.empty.model(mt)) glm.fit.null else glm.fit)(x = X, y = Y, As I checked from the Archive R-Help mails, it seems that this happens when the dataset exhibits complete separation. However, p-values

Adding to Ashutosh's comments, We are also interested in making LLVM generate vector math library calls that are available with glibc (version > 2.22). reference: https://sourceware.org/glibc/wiki/libmvec Using the example case given in the reference, we found there are 2 vector versions for "sin" (4 X double) with same VF namely _ZGVcN4v_sin (avx) version and _ZGVdN4v_sin

Full_Name: David Firth Version: 2.1.1 OS: Mac OS X Submission from: (NULL) (137.205.240.25) A small inconsistency in the "Writing R Extensions" manual: at http://cran.r-project.org/doc/manuals/R-exts.html#The-DESCRIPTION-file the displayed example has Title: My first collection of functions but the text that follows it says "The Title field...should be capitalized",

Full_Name: David Firth Version: 1.8.0 OS: Mac OS 10.2.8 Submission from: (NULL) (80.177.12.48) In .Rd files, we use \url{stuff} to refer to a URL on the web. Just as the documentation says. Perhaps this isn't a bug, then, but certainly it's a trap for the unwary. If the construction \link{url{something}} appears in a .Rd file, then Rdconv --type=html ... gets stuck and appears

Hi, Having searched google '[R] aov effect size' without any results I wonder if I not completely miss something. Is there any R function that calculates the effect size of an AOV's main effect or interaction effect? It should be related to the F's and the degree of freedom of the error, but the multitude in numbers in aov() baffle me a bit. Thanks, ---david

On 07/02/2018 04:33 PM, Saito, Hideki wrote: > >   > > >It may not be a full solution for the problems you're trying to solve > >   > > If we are inventing a new solution, I’d like it also to solve OpenMP > declare simd legalization issue. If a small extension of existing scheme > > works for mathlib only, I’m happy to take that and discuss OpenMP >

R-devel has the --with-lapack flag for configure. If you build R with --with-blas="-framework vecLib" --with-lapack="--framework vecLib" then the build goes through, using the native optimized BLAS and Lapack in /System/Library/Frameworks/vecLib.framework. Also, it works, in the sense that it does eigenvalue problems correctly. I don't have any timings yet, because I

Illustrative Example: clang -fveclib=SVML -O3 svml.c -mavx #include <math.h> void foo(double *a, int N){ int i; #pragma clang loop vectorize_width(8) for (i=0;i<N;i++){ a[i] = sin(i); } } Currently, this results in a call to <8 x double> __svml_sin8(<8 x double>) after the vectorizer. This is 8-element SVML sin() called with 8-element argument. On the surface,

Is there a straightforward way to convert character information from decimal representation to the octal one used by R? I'd like something like a function ascii(number,base=10), such that > ascii(91) [1] "\133" I can easily do the mapping from 91 to 133, but what is a good way to operate on 133 to deliver "\133"? Would a lookup table be a better solution? David

I am trying to use the dmt function in the package {mnormt}. Throughout my algorithm, the covariance matrix is sometime calculated to be singular. When attempting to calculate the dmt function with a covariance that is not positive definite, I would like it to return Inf or NaN instead of an error message. I have been using the try function, however it is not yeilding the desired result. (I did

OS X 10.2 and higher comes standard with highly optimized versions of BLAS and Lapack in /Systems/Library/Frameworks/vecLib.framework. It seems that even for double precision they do much better. See http://sthmac.magnet.fsu.edu/benchmarks/ I am not sure how these numbers would look on G3 Macs, but obviously for double precision there is not much reliance on Altivec. So I tried to configure

I am making my first attempts at using some Fortran code with R, and so far it's going OK. To print from my Fortran programs, it seems I need subroutines dblepr, realpr and intpr. From the excellent "Writing R Extensions" document: "Three subroutines are provided to ease the output of information from FORTRAN code. subroutine dblepr(label, nchar, data, ndata)