similar to: Optimizing C code

Vince Carey found this (thank you!). Since the fix to the problem is not entirely obvious, I post this to R-devel as RFC: help(dist) says: >> Missing values are allowed, and are excluded from all computations >> involving the rows within which they occur. If some columns are >> excluded in calculating a Euclidean, Manhattan or Canberra >> distance, the sum is

R-2.9.0/include/R_ext/Arith.h has: int R_IsNA(double); /* True for R's NA only */ int R_IsNaN(double); /* True for special NaN, *not* for NA */ int R_finite(double); /* True if none of NA, NaN, +/-Inf */ #define ISNA(x) R_IsNA(x) /* True for *both* NA and NaN. The first and last lines are contradictory - if R_IsNA is true only for NA, not NaN, then ISNA should be the same.

Hello, I'm a new R user and I'm having a little trouble getting started. I'm hoping someone can help me out. I read numbers (integer) from a SQL database. Some calculations as mean()are possible but hist() give an error. Here are the commands: library(RMySQL) Loading required package: DBI > con <- dbConnect(dbDriver("MySQL"), dbname = "test") >

Full_Name: Albrecht Gebhardt Version: 0.65.0 OS: Digital Unix 4.0E Submission from: (NULL) (212.17.104.62) Plottimg on the Alpha stopped working with 0.65.0. The tickmarks have a length of -Inf and go across the whole plot. A first solution is the following patch: --- src/main/plot.c.alpha-patch Mon Sep 13 01:37:11 1999 +++ src/main/plot.c Mon Sep 13 01:58:16 1999 @@ -832,7 +832,7 @@

Hi, Apropos of a recent Inf question, I've previously wondered if dnorm "does the right thing" with dnorm(0, 0, -Inf) which gives zero. Should that be zero or NaN (or NA)? The help says "'sd < 0' is an error and returns 'NaN'" and since -Inf < 0 is TRUE, then... is this a bug? Thank you, Stephen Rochester, MN USA

Yes, that looks like a bug and an easily fixable one too. However, I spy another issue: Why do we check the !R_FINITE(x) && mu == x before checking for sd < 0 ? The difference is whether we return ML_NAN; or ML_ERR_return_NAN; but surely negative sd should always be an error? I'd be inclined to do if (sigma < 0) ML_ERR_return_NAN; if(!R_FINITE(sigma)) return R_D__0;

Hi, poids taille fumeur sexe sport etat 85 184 oui homme 1 malade 65 175 oui homme 1 malade 74 180 oui homme 2 gueri 79 175 oui homme 2 malade 71 165 non homme 3 gueri 80 185 non homme 3 gueri 75 180 non homme 4 malade 69 155 non homme 4 malade 74 168 oui

I have a problem building an extension using ISNAN() on R version 2.0.x. In R 1.9.1 Arith.h and Rmath.h contained code like #ifdef IEEE_754 # define ISNAN(x) (isnan(x)!=0) #else # define ISNAN(x) R_IsNaNorNA(x) #endif #define R_FINITE(x) R_finite(x) int R_IsNaNorNA(double); int R_finite(double); which works. R 2.0.x has # define ISNAN(x) (isnan(x)!=0) unconditionally. This breaks

I have a problem building an extension using ISNAN() on R version 2.0.x. In R 1.9.1 Arith.h and Rmath.h contained code like #ifdef IEEE_754 # define ISNAN(x) (isnan(x)!=0) #else # define ISNAN(x) R_IsNaNorNA(x) #endif #define R_FINITE(x) R_finite(x) int R_IsNaNorNA(double); int R_finite(double); which works. R 2.0.x has # define ISNAN(x) (isnan(x)!=0) unconditionally. This breaks

Hi, Dear R users I have problem with the following code. The matrix result must be a matrix (3x3). But I have obtained a matrix(3x1) and I don't know why. So, I need your help Best regards ##################################################################### taille <- function (delta, t, prob = 0.2) { niv.conf <- c(0.90, 0.95, 0.99) if(niv.conf <- 0.90) { t <- 1.645 }

Hi everyone, Has anyone written a parser in Java for either the ASCII or binary format produced by save()? I need to parse a single large 2D array that is structured like this: list( "32609_1" = c(-9549.39231289146, -9574.07159324482, ... ), "32610_2" = c(-6369.12526971635, -6403.99620977124, ... ), "32618_2" = c(-2138.29095689061, -2057.9229403233, ... ),

Bonjour. Une fois mon application ouverte, je voudrais savoir, si cela est possible, comment dimensionner automatiquement la totalité des colonnes du grid à la taille de mon boxsizer, quand je change la taille de mon application? (puique mon grid est inséré dans mon boxsizer). merci

en fait je vous renvoi ce message car maintenant la taille de la police change au bout de 500 lignes environ, alors que j''avais reussi à garder la taille jusqu''à environ 5000 lignes. Je ne comprend pas. Est-ce que c''est instable? le bout de code que vous m''avez envoyer ne marche pas. merci Sebastien _______________________________________________ wxruby-users

I'm observing that base::sum(x, na.rm=FALSE) for typeof(x) == "double" is much more time consuming when there are missing values versus when there are not. I'm observing this on both Window and Linux, but it's quite surprising to me. Currently, my main suspect is settings in on how R was built. The second suspect is my brain. I hope that someone can clarify the below

Dear R users, When I try to resize a graphics device in R, I often get this warning message (mostly when I already have several other graphics devices open) : "Not enough memory to modify the size. Alpha blending is desactivated" (translation from the french message: "Mémoire insuffisante pour modifier la taille. L'"alpha blending" est désactivé"). Following

Dear r-help list, I'm trying to use multiple imputation for my MSc thesis. Having good exemples using the Hmisc package, I tried the aregImpute function. But with my own dataset, I have the following error : Erreur dans matxv(X, xcof) : columns in a (51) must be <= length of b (50) De plus : Warning message: In f$xcoef[, 1] * f$xcenter :   la taille d'un objet plus long n'est pas

In my package 'glmmML' I'm using old C code and linpack in the optimizing procedure. Specifically, one part of the code looks like this: F77_CALL(dpoco)(*hessian, &bdim, &bdim, &rcond, work, info); if (*info == 0){ F77_CALL(dpodi)(*hessian, &bdim, &bdim, det, &job); ........ This usually works OK, but with an ill-conditioned data

Hi R-devel, I have a question about the differentiation between NA and NaN values as implemented in R. In arithmetic.c, we have int R_IsNA(double x) { if (isnan(x)) { ieee_double y; y.value = x; return (y.word[lw] == 1954); } return 0; } ieee_double is just used for type punning so we can check the final bits and see if they're equal to 1954; if they are, x is NA, if they're

Hi all, I mounted by smbmount a partition on a winNT machine without problem... but when i tried to copy a file whitch have size more than 2 Gb i have the following error : Overflow size allowed for a file (translated from french ;-) :D?bordement de la taille permise pour un fichier ) PS. the file system of the harde drive is NTFS PS. The command of smbmount : /usr/sbin/smbmount

Bonjour, Je veux cherché l’estimateur de vraisemblance maximal (MLE)d’une fonction à 3 paramètre inconue étant donné une échantillon de taille 50 (les observations des valeurs de x) alors comment je peux procédé La fonction de densité est définie par : f(x)= 1/3(g(a1)+g(a2)+g(a3)) avec g(ai)=(exp(ai)*ai^x)/x! pour i=1,2,3. Je vous remercie beaucoup. A. Elhabti