similar to: "evaluating expression 0" bug in R-1.4.1 (PR#1385)

I'm using the 'pnt' C function from Rmath library in some C-code. How can I catch the warning message: "full precision was not achieved in 'pnt'" in R. I call the function using the .C(). (options(warn=-1) didn't work) Thanks in advance -- Maarten van Iterson Center for Human and Clinical Genetics Leiden University Medical Center (LUMC) Research Building,

Dear list, If have some repeated measurement data which looks something like: time <- rep(1:5 , each=2*4) groups <- rep(c("Case", "Control"), each=4) subjects <- factor(rep(1:(2*4), 5)) responses <- time + rnorm(5*2*4) + as.integer(factor(groups)) data <- data.frame(responses, time, groups, subjects) Now I want to plot each subject in a separate panel:

Dear all, I would like to use predict.lm() with an existing lm object but with new arbitrary coefficients. I modify 'fit$coef' (see example below) "by hand" but the actual model in 'fit' used for prediction does not seem to be altered (although fit$coef is!). Can anyone please help me do this properly? Thanks in advance, J?r?mie > dat <-

Hi all, I should have given a better explanation of my problem. Here it is. I extracted from my code the bit that gives the error. Place this in a file called test.c #include <math.h> #include <R.h> #include <Rmath.h> #include <float.h> #include <R_ext/Print.h> int main(){ double prob[3] = {0.0, 0.0, 0.0}; double prob_tot = 0.; prob[0] = 0.3*dnorm(2, 0,

https://bugzilla.netfilter.org/show_bug.cgi?id=1385 Bug ID: 1385 Summary: Incorrectly evaluated expression with negated ip saddr and negated ip daddr Product: nftables Version: unspecified Hardware: x86_64 OS: Debian GNU/Linux Status: NEW Severity: normal Priority: P5

Full_Name: Jelle Goeman Version: 2.0.0 OS: windows Submission from: (NULL) (145.88.209.33) I get a strange error: > is.vector(1:10) Error: recursive default argument reference What's recursive about is.vector? Kind regards, Jelle Goeman My R: platform i386-pc-mingw32 arch i386 os mingw32 system i386, mingw32 status major 2

> On 10 Mar 2016, at 21:25 , M.van_Iterson at lumc.nl wrote: > > Hi all, > > I should have given a better explanation of my problem. Here it is. > > I extracted from my code the bit that gives the error. Place this in a file called test.c Aha. Missing info #1, C not R... > > #include <math.h> > #include <R.h> > #include <Rmath.h> >

Full_Name: Jelle Goeman Version: 2.4.0 OS: windows XP Submission from: (NULL) (145.88.209.33) Hi, I like to use optimize() to optimize functions whose evaluation is costly in terms of computation time. The Brent algorithm which is implemented in optimize was designed to optimize a function with as few function evaluations as possible. Therefore it bothers me that optimize() always evaluates

Full_Name: Jelle Goeman Version: 1.9.0 OS: mingw32, windows 2000 Submission from: (NULL) (145.88.209.33) Hi Everyone, I get very strange results using the function hasArg with the ... function argument. In my own function: > gt <- globaltest(X,Y) > sampling(gt) works fine, but > sampling(globaltest(X,Y)) results in: Error in eval(expr, envir, enclos) : "missing"

Dear all, I have a questions regarding using the c function rmultinom.c. I got the following error message "rbinom: probability sum should be 1, but is 0.999264" Which is thrown by: if(fabs((double)(p_tot - 1.)) > 1e-7) MATHLIB_ERROR(_("rbinom: probability sum should be 1, but is %g"), (double) p_tot); I understand my probabilities do not sum to one close enough. I

Hi, I want to create a data frame similar to the following, but greatly scaled up: df <- data.frame(aaa= c("a","b","c"), integer(3), integer(3)) names(df)[2:3] <- paste("var",1:2,sep="") which yields aaa var1 var2 1 a 0 0 2 b 0 0 3 c 0 0 I would not relish having to paste 'integer(3)' 5000 times :(

I am not sure if this is already known. I checked BUGS and found some references to similar behavior when R_HOME is set? Here is what I get: (1) Fresh build of 1.4.1 under Linux/RH-7.2. All tests pass. > R CMD check ctest * checking for working latex ... OK * using log directory `/home/andyj/stat/R-1.4.1/library/ctest.Rcheck' Installing *source* package `ctest'

Hello! I've looked at R help and previous postings, but I am not sure I completely understand the mechanism of evaluating R expressions from C++. They have function SEXP eval(SEXP expr, SEXP rho), but I can't get it to work. Could anyone who used it give any comments on how efficient it is and, if possible, give example of how to use it. My main problem is that I don't really

Hi, developers, I notice following asm symbol definitions are eveluated improperly: A = 9 B = 3*A - 2 * A + 1 the value of B should be 10 instead of 8. It seems asm parser treat "B = 3*A - 2 * A + 1" as "B = 3*A - (2 * A + 1)" because the second "*" has higher precedence than "-", and evaluation just performs simple recursion within

On 3 May 2013 01:37, 杨勇勇 <triple.yang at gmail.com> wrote: > Hi, developers, > > I notice following asm symbol definitions are eveluated improperly: > > A = 9 > B = 3*A - 2 * A + 1 > > the value of B should be 10 instead of 8. > > It seems asm parser treat "B = 3*A - 2 * A + 1" as "B = 3*A - (2 * A + 1)" > because the second "*"

All, Thank you for the prompt and useful answers to my questions. I had missed references in 5.7.6 which would have answered some of the points. As Bill pointed out a newer version of acrobat would help, but the Sun system here is still running 5.0. (An oversubscribed sysadmin). Then I could have searched and at at least avoided the most trivial. All three comments were different,

Dear R experts! I want to evaluate expressions from within C code. For instance, I want this code: /* ============ -- C code -- ============== */ #include <R.h> #include <Rinternals.h> SEXP foo(SEXP expr, SEXP rho) { SEXP x; PROTECT(x = allocVector(REALSXP, 1)); REAL(x)[0] = 1; UNPROTECT(1); return

Hallo I'm having trouble figuring out how to evaluate an expression when one of the variables in the expression is defined separately as a sub expression. Here's a simplified example mat <- expression(0, f1*s1*g1) # vector of formulae g1 <- expression(1/Tm) # expansion of the definition of g1 vals <- data.frame(f1=1, s1=.5, Tm=2) # one set of possible values for

Full_Name: Martin Lenze Version: Version 1.9.0 (2004-04-12), ISBN 3-900051-00-3 OS: Microsoft Windows 2000 [Version 5.00.2195] SP4 Submission from: (NULL) (82.83.167.79) Hello, now I switched to R 1.9.0 and did get a problem with subset, see sample: x <- data.frame(a=as.integer(round(runif(5),0)), b=as.integer(round(runif(5),0)),c=as.integer(round(runif(5),0))) # correct results: for (c

I am currently puzzled by a passage in the R Extensions manual, section 5.10: SEXP lapply(SEXP list, SEXP expr, SEXP rho) { R_len_t i, n = length(list); SEXP ans; if(!isNewList(list)) error("`list' must be a list"); if(!isEnvironment(rho)) error("`rho' should be an environment"); PROTECT(ans = allocVector(VECSXP, n));