[R] Weird (wrong?) result for fft

2021-03-22 Thread Filippo Monari
Hi all, I am plying around with fft function from the stats package. Running the example that is listed in the documentation: x <- 1:4 fft(x) #output 1 fft(fft(x), inverse = TRUE)/length(x) #output 2 I was expecting that output 1 and two were the same but I get: > fft(x) [1] 10+0i -2+2i -2+0i -

[R] Fwd: Splines

2016-09-05 Thread Filippo Monari
Hi, does anybody know what header file includes the bspline functions used in R? Regards, Filippo Forwarded Message Subject:Splines Date: Sun, 04 Sep 2016 09:12:37 +0100 From: Filippo Monari To: R-help@r-project.org Hi, I would like to use the C spline

[R] Splines

2016-09-04 Thread Filippo Monari
Hi, I would like to use the C spline functions if R for a FORTRAN subroutine. What header file should I refer to? Regards Filippo Monari __ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see https://stat.ethz.ch/mailman/listinfo/r-help

[R] .C and .Fortran

2014-09-03 Thread Filippo Monari
Hi, I'd like to know what is the difference between the functions .C() and .Fortran. I noticed that by enclosing my F90 subroutines in modules .Fortran() can't find them any more in the load table, while .C() still can. I also checked that the subroutine was loaded with the is.loaded() function

[R] garchFit updating

2013-10-23 Thread Filippo Monari
NULL, NULL)) Original Message Subject:garchFit Date: Wed, 23 Oct 2013 17:21:34 +0100 From: Filippo Monari To: r-help@r-project.org Hi, I'm using "garchFit" function from packages fGarch, to fit some nose. The weird thing is that in some cases the estimate f

[R] garchFit

2013-10-23 Thread Filippo Monari
Hi, I'm using "garchFit" function from packages fGarch, to fit some nose. The weird thing is that in some cases the estimate for the alpha1 coefficient is equal to 1, which can't be otherwise the model is not stable. Even more alpha1 = 1 makes impossible to simulate the model. My model is an ARM

[R] reference class error

2013-10-22 Thread Filippo Monari
hi, I'm new to reference classes but as I have experience with python I decided to port some of my code in this frame work. I have the following two reference classes in separate files: #1# uvRndVar = setRefClass( Class = 'uvRndVar', fields = list( desc = 'character', ty

[R] strange boolean results

2013-06-13 Thread Filippo Monari
Hi, anyone can explain to me the following? > rho [1] 0.9452398 0.8792735 0.9641829 0.9876954 0.9993560 0.9826084 1.000 [8] 1.000 0.7982916 1.000 0.3361956 > any(rho >= 1) [1] FALSE thanks, regards Filippo __ R-help@r-project.org mailin

[R] Fortran modules

2013-06-11 Thread Filippo Monari
Hi, I have some subroutines using function and subroutine as well from Fortran modules. In the f90 source code I used the statement: use mymodule and it compile well through the R CMD SHL command. Anyway when I call dyn.load('myF90.so') form R I get the following error: unable to load shared

[R] help pairs

2013-03-20 Thread Filippo Monari
Hi everybody, I'm exploring likelihood function having several parameters and I'm trying to work with pairs but I have the following problem. I need to colour my dots depending on the value of the likelihood of just two parameters, thus basically I need different colour scale for each plot of t