search for: kaschek

...en I uncomment the print(), it works as expected. Is this a bug or a feature? conditions <- 1:4 test <- lapply(conditions, function(mycondition){ #print(mycondition) myfn <- function(i) mycondition*i return(myfn) }) sapply(test, function(myfn) myfn(2)) Cheers, Daniel -- Daniel Kaschek Institute of Physics Freiburg University Room 210 Phone: +49 761 2038531

From: Daniel Kaschek <daniel.kaschek at physik.uni-freiburg.de> > ... When I evaluate this list of functions by > another lapply/sapply, I get an unexpected result: all values coincide. > However, when I uncomment the print(), it works as expected. Is this a > bug or a feature? > > conditions &...

...defining the derivatives. No problem so far. However, now I need to implement an ODE that contains spline functions, i.e. the derivatives at given time points depend on the value of a spline function at this time point. What is an efficient way to implement this in precompiled C code? -- Daniel Kaschek <daniel.kaschek at physik.uni-freiburg.de> -------------- next part -------------- A non-text attachment was scrubbed... Name: not available Type: application/pgp-signature Size: 490 bytes Desc: This is a digitally signed message part URL: <https://stat.ethz.ch/pipermail/r-help/attachments...

Hi, When I fit the mixed-effects model by lme, I got the convergence error code = 1 message = iteration limit reached without convergence (9). Even after I increase the iteration numbers, it still did not work. The only way I found to avoid the error message is to increase the number of iterations for the EM algorithm, which need to be at least 1100. But the weird thing is that the correlation

Dear Martyn, On Fr, Jan 15, 2016 at 4:01 , Martyn Plummer <plummerm at iarc.fr> wrote: > > Alternatively, you may be able to control the maximum number of > threads > by setting and exporting an appropriate environment variable depending > on what backend you are using, e.g. OPENBLAS_NUM_THREADS or > MKL_NUM_THREADS. Thanks a lot. Running export OPENBLAS_NUM_THREADS

...code optimization might reduce the overhead in some instances (e.g. if FUN is a BUILTIN), but it's not clear that would be worth while. Variants of this issue arise in a couple of places so it may be worth looking into. Best, luke On Tue, 24 Feb 2015, Radford Neal wrote: > From: Daniel Kaschek <daniel.kaschek at physik.uni-freiburg.de> >> ... When I evaluate this list of functions by >> another lapply/sapply, I get an unexpected result: all values coincide. >> However, when I uncomment the print(), it works as expected. Is this a >> bug or a feature? >&gt...

...tances (e.g. if FUN is a BUILTIN), but it's > not clear that would be worth while. > > Variants of this issue arise in a couple of places so it may be worth > looking into. > > Best, > > luke > > > On Tue, 24 Feb 2015, Radford Neal wrote: > > From: Daniel Kaschek <daniel.kaschek at physik.uni-freiburg.de> >> >>> ... When I evaluate this list of functions by >>> another lapply/sapply, I get an unexpected result: all values coincide. >>> However, when I uncomment the print(), it works as expected. Is this a >>>...

On 23/02/2015 3:57 PM, Daniel Kaschek wrote: > Hi everybody, > > with the following code I generate a list of functions. Each function > reflects a "condition". When I evaluate this list of functions by > another lapply/sapply, I get an unexpected result: all values coincide. > However, when I uncomment the...

On Fri, 2016-01-15 at 15:03 +0100, Daniel Kaschek wrote: > Dear all, > > I run different R versions (3.2.1, 3.2.2 and 3.2.3) on different > platforms (Arch, Ubuntu, Debian) with a different number of available > cores (24, 4, 24). The following line produces very different behavior > on the three machines: > > for(i in...

Dear all, I have a function "fn" with its own environment, i.e. env <- environment(fn) and env is not .GlobalEnv. And another function getValue <- function(x) environment(x)$mylist which returns the list object "mylist" which is in "env". If I want to modify "mylist", I could write 'getValue<-' <- function(x, value) {

I am trying to calculate monthly means by year of phosphates and nitrates from a multi year data set. Can anybody help me out with the most effective way to do this? My data looks like this: Collection_Date                   Test.Name              Value 2000-01-24 17:00:00            Phosphate               0.108 2000-01-24 17:00:00            Nitrate                     0.037 2001-11-12

Dear all, I run different R versions (3.2.1, 3.2.2 and 3.2.3) on different platforms (Arch, Ubuntu, Debian) with a different number of available cores (24, 4, 24). The following line produces very different behavior on the three machines: for(i in 1:1e6) {n <- 100; M <- matrix(rnorm(n^2), n, n); M %*% M} On the Ubuntu and Arch machine one core is used, but on the Debian machine ALL

Hi all, I need to test the significnce of difference between slopes of two regression lines and regression line with theoretical line. I try to use Slope.test function from emu package, but an error occured... library(emu) d1<-data.frame(P1=c(1,2,3,5,7,8,9,13,14,15), P2=c(1,2,5,8,11,13,15,15,18,24), R=c(2,7,8,9,16,21,27,31,33,36)) # First data set m1<-lm(R~P1+P2+P1*P2,data=d1) # Regr.