similar to: Problems with Rcmd check on Win 2000 & rw1062

I'm having trouble figuring out how to create a function using "body<-" (). The help file for body() says that the argument should be a list of R expressions. However if I try that I get an error: > tmpfun <- function(a, b=2){} > body(tmpfun) <- list(expression(z <- a + b),expression(z^2)) Error in as.function.default(c(formals(f), value), envir) :

Dear R Gang, I'm interested in using R and the nls package for fitting kinetic models. I'm having some difficulty getting a model specified for nls though. The math for the model that I want to fit is dg(t)/dt = K1 f(t) - k2 g(t) where g(t) and f(t) are measured data at a sequence of times t. K1 and k2 are the parameters of the model. If I solve this, the solution is g(t) = K1

Greetings, I have a meta-analysis problem in which I have fixed effects regression coefficients (and estimated standard errors) from identical models fit to different data sets. I would like to use these results to create pooled estimated regression coefficients and estimated standard errors for these pooled coefficients. In particular, I would like to estimate the model \beta_{i} = \mu +

When I try to open a server socket on the system setup listed below, I get the following error: > a <- make.socket(port=8091,server=TRUE) Error in make.socket(port = 8091, server = TRUE) : sockets are not available on this system I just wanted to make sure that sockets were generally available for this platform (or, perhaps, I need to upgrade to 1.2.2 on this system), before I got

Hi, I am trying to solve a model that consists of rather stiff ODEs in R. I use the package ODEsolve (lsoda) to solve these ODEs. To speed up the integration, the jacobian is also specified. Basically, the model is a one-dimensional advection-diffusion problem, and thus the jacobian is a tridiagonal matrix. The size of this jacobian is 100*100. In the original package

I have a package that solves systems of ordinary differential equations coded as an R function (odesolve, on CRAN). The function is passed to R code, and c and Fortran code called by it uses lang4() and eval() to evaluate the R function inside a compiled c function to use in a compiled ODE solver. This works quite well, but can be slow. I'd like to be able to supply the name of a compiled

I think this is a bug; at least this behavior is not documented in plot or plot.default. plot.default resets xaxp, and leaves xaxp reset when it exits: par(xaxp=c(0,1,4)) print(par("xaxp")) plot(c(0,1),c(0.2,0.3)) print(par("xaxp")) R. Woodrow Setzer, Jr. Phone: (919) 541-0128 Experimental Toxicology Division

I found that predict.gnls failed with a wierd error message about a non-numeric argument to a binary vector in one of three nearly identical uses. Error in Inh/Ki : non-numeric argument to binary operator (Inh and Ki are arguments to the function used in the formula for the object whose predictions were requested). It turns out that the problem is in getCovariateFormula(). The final line in

Does the following indicate that I have a memory leak? > gc(TRUE) Garbage collection 22891 = 21012+1557+322 (level 2) ... 483257 cons cells free (59%) 11.0 Mbytes of heap free (82%) used (Mb) gc trigger (Mb) Ncells 334906 9.0 818163 21.9 Vcells 309342 2.4 1746173 13.4 > memory.size()/1024/1024 ## in MB [1] 643.8978 > sum(sapply(ls(all.names=TRUE),object.size)) [1] 1776

If a function needs to be passed as an argument to another function, default arguments to the function being passed are lost. Consider this example: fun1 <- function(x, A=c("power","constant")) { arg <- match.arg(A) cat(paste("A is:",paste(A,collapse=", "),"\narg is:",arg,"\n")) cat("formals:\n") print(formals())

I think qqline does not do exactly what it is advertised to do ("`qqline' adds a line to a normal quantile-quantile plot which passes through the first and third quartiles."). Consider the graph: tmp <- qnorm(ppoints(10)) qqnorm(tmp) qqline(tmp) The line (which I expected go through all the points), has a slightly shallower slope than does the points plotted by qqnorm. I think

A couple of minor points about the odesolve package (which I am otherwise enjoying very much): 1. "scalar" is misspelled as "scaler" in the definitions of the rtol and atol parameters 2. it is possible to crash R by doing something dumb, e.g failing to read the documentation carefully enough and (a) returning only a vector of derivatives and not a list of (derivatives,

I'm trying to use odesolve for integrating various series of coupled 1st order differential equations (derived from a system of enzymatic catalysis and copied below, apologies for the excessively long set of parameters). The thing that confuses me is that, whilst I can run the function rk4: out <- rk4(y=y,times=times,func=func, parms=parms) and the results look not unreasonable:

Probably I've just misread the documentation, but I don't understand the behavior of .Alias. Consider this (on R-1.1.1, both Windows and SGI): > tmp <- matrix(nrow=3,ncol=2) > new <- .Alias(tmp) > new[1,1] <- 1 > tmp [,1] [,2] [1,] NA NA [2,] NA NA [3,] NA NA > new [,1] [,2] [1,] 1 NA [2,] NA NA [3,] NA NA I expected tmp[1,1] to

Hello there, Suppose you want to solve the following system of ODE's (a simple Lotka-Volterra predator prey model) dP/dt = beta*P*V - mu*P dV/dt = r*V - beta*P*V where P and V are the numbers of predators and prey. Now, this is easy to do, but suppose you have a system of equations like this, dP1/dt = beta1*P1*V1 - mu1*P1 dP2/dt = beta2*P2*V2 - mu2*P2 dV1/dt = r1*V1 - beta1*P1*V1

The binary package (version of 6/17) of rw1051 for Windows is missing html documentation for dyn.load() The file is present in $R_HOME/library/base/html, but has zero bytes. dyn.load help IS accessible through the windows help system. platform i386-pc-mingw32 arch i386 os mingw32 system i386, mingw32 status major 1 minor 5.1 year 2002 month 06 day 17 language R

Back on January 16, a message on R-help from Ravi Varadhan described a problem with gnls using weights=varPower(). The problem was that the fit failed with error Error in eval(expr, envir, enclos) : Object "." not found I can reliably get this error in version 2.0.1-patched 2004-12-09 on Windows XP and 2.0.1-Patched 2005-01-26 on Linux. The key feature of that example is that the

I am beginning a much-delayed update of odesolve to include several ordinary differential equation solvers from the Livermore package ODEPACK. These are much-used and reliable Fortran codes, and I plan (as I did for lsoda in the current odesolve package) to make as few changes as possible to the Fortran 77 code. However, recently someone wanted to make nested calls to lsoda, which will not work,

I'm probably misunderstanding something in "Writing R Extensions" version 1.0.0. In the chapter on the R API, section 4.7, it is stated that the functions listed in R_ext/Linpack.h are available to users' Fortran code. I am developing a developing a library of ode solvers, based on lsoda and ddassl, and which in turn call some routines from linpack and double precision blas. I

This is in Windows 98. The problem exists in both rw1021 and rw1023, using ActiveState Perl builds 620 and 626: C:\>Rcmd REMOVE --library=C:/R/Rlibs RBMDS Unknown option: library Usage: Rcmd REMOVE [options] pkgs Remove the add-on packages specified by pkgs from the default library tree (/library) or the tree specified via `--library'. Options: -l, --library=LIB remove packages