similar to: Modified Bessel Function - 2nd kind

Dear R-Group: Recently, I ran into a problem. I was using a function called "I.1", which evaluates the first-order modified Bessel function of the first kind, in the package "CircStats". This function is not vectorized, since it uses a couple of "if" conditions. However, when I called this function with a vector argument, I got no error/warning messages in

#bug 1: besselI() for nu<0 and expon.scaled=TRUE #tested with R-devel (2007-06-17 r41981) x <- 2.3 nu <- -0.4 print(paste(besselI(x, nu, TRUE), "=", exp(-x)*besselI(x, nu, FALSE))) #fix: #$ diff bessel_i_old.c bessel_i_new.c #57c57 #< bessel_k(x, -alpha, expo) * ((ize == 1)? 2. : 2.*exp(-x))/M_PI #--- #> bessel_k(x, -alpha, expo) * ((ize == 1)? 2. :

The besselK() function knows about overflows/underflows internally; there is a constant xmax_BESS_K in src/nmath/bessel.h (and referred to only in bessel_k.c), equal to 705.342, which is checked if expon.scaled is FALSE. (The equivalent number for bessel_i.c is 709, defined as exparg_BESS in bessel.h.) However, besselK(x) silently returns +Inf if x>705.342. This behavior is reasonable for

I have a small handheld pc having ARM process as a CPU. I installed debian and installed R using apt-get command. Everything worked great except for drawing even simple graphs x <- 1:10 plot(x) I got error messages 1: Nonfinite axis limits [GScale(nan,nan,1, .); log=0] 2: relative range of values = 9.0072e+15 * EPS, is small (axis 1). 3: Nonfinite axis limits

Dear all Currently, I'm implementing the generalized hyperbolic distribution into Splus. Unfortunately the Bessel function is not implemented in Splus. In R the Bessel function does exist but it is an internal function and I'm not able to look at the code. Is there any possibility to see the code of the Bessel function in R or does anybody has an implementation of the Bessel function in

Hello, I am searching for code in C++ or fortran for Modified Bessel function of third kind (fractional or real order). Can someone help me? Thank you --------------------------------- Découvrez un nouveau moyen de poser toutes vos questions quelque soit le sujet ! Yahoo! Questions/Réponses pour partager vos connaissances, vos opinions et vos expériences. Cliquez ici.

I am not yet a R user but I will be soon. I am looking for the R command and syntax to compute the roots of Bessel function i.e. computing the z values that lead to Jnu(z)=0 where J is a Bessel function or order nu. May You help me ? thanks in advance. Dr Catherine COUTAND Institut National de la Recherche Agronomique (INRA) umr Physiologie Int?grative de l'Arbre Fruitier et Forestier

Hi: There is a Modified Bessel function (third kind, real order) besselK in R. Is there a C version of this function? I can only find a Fortran function rkbesel on http://www.netlib.org/specfun/. However, this rkbesel fortran function is the second kind. Thank you! liu _______________________________ Declare Yourself - Register online to vote today!

I am looking for a method of dealing with the modified Bessel function K_\nu(x) for large \nu. The besselK function implementation of this allows for dealing with large values of x by allowing for exponential scaling, but there is no facility for dealing with large \nu. What would work for me would be an lbesselK function in the manner of lgamma which returned the log of K_\nu(x) for large

Debian tries to build its packages on a variety of platforms. The arm platform compiled 0.90.1 (the last Debian release before the Debian package required an Atlas library, something we no longer require) failed in 'make check'. The log snippet follows; I traced this to the example(Bessel) code. > matplot(nu, t(outer(xx,nu, besselI)), type = 'l', ylim = c(-50,200), +

Hi R-users, Does R has a topic on newton's method? Thank you for the info.

Dear all, I'm writing a code that requires Bessel functions with complex argument. Searching the list, I found the continuation of a thread I initiated a few months ago: http://tolstoy.newcastle.edu.au/R/e4/devel/08/03/0746.html As I understand, the most promising option would be to use the fortran or C implementation of Amos,

Here is a copy of the last few lines in base-Ex.Rout.fail: > x <- seq(3,500);yl <- c(-.3, .2) > plot(x,x, ylim = yl, ylab="",type='n', main = "Bessel Functions Y_nu(x)") > for(nu in nus){xx <- x[x > .6*nu]; lines(xx,besselY(xx,nu=nu), col = nu+2)} > legend(300,-.08, leg=paste("nu=",nus), col = nus+2, lwd=1) > > x <-

Hi, R guys: I'm using L-BFGS-B method of optim for minimization problem. My function called besselI function which need non-negative parameter and the besselI will overflow if the parameter is too large. So I set the constraint box which is reasonable for my problem. But the point outside the box was test, and I got error. My program and the error follows. This program depends on CircStats

> Dears users, > I'm trying to run R on a computer with OS X. I installed the R folder > and the dylibs, where they must be installed. I updated my path, and > R runs but tells me : > > dyld: /usr/local/R-1.4.0/lib/R/bin/R.bin version mismatch for > library: /usr/lib/libz.1.dylib (compatibility version of user: 1.1.3 > greater than library's version: 1.0.0)

Dear R users, I'm porting a piece of Matlab code to R, but I'm now stuck with the following: I need an equivalent of besselJ(x, nu) that can handle a complex argument x. I couldn't find any R implementation. I did find a possible fortran solution in SLATEC (< http://www.netlib.org/slatec/ > , CBESJ-C), however I've never tried to use external C or Fortran code

Dear R users, I'm porting a piece of Matlab code to R, but I'm now stuck with the following: I need an equivalent of besselJ(x, nu) that can handle a complex argument x. I couldn't find any R implementation. I did find a possible fortran solution in SLATEC (< http://www.netlib.org/slatec/ > , CBESJ-C), however I've never tried to use external C or Fortran code

Dear list, Could someone enlighten me a bit on the exact format used for text help pages as I see them on Windows in the help/ folders of (compiled) package roots. When opening an example of these files in a text editor (?Bessel in GNU emacs), the file is displayed as follows: _^HB_^He_^Hs_^Hs_^He_^Hl_^H _^HF_^Hu_^Hn_^Hc_^Ht_^Hi_^Ho_^Hn_^Hs Bessel package:base

Hi, I have written out the log-likelihood function to fit some data I have (called ONES20) to the non-central chi-squared distribution. >library(stats4) >ll<-function(lambda,k){x<-ONES20; 25573*0.5*lambda-25573*log(2)-sum(-x/2)-log((x/lambda)^(0.25*k-0.5))-log(besselI(sqrt(lambda*x),0.5*k-1,expon.scaled=FALSE))} > est<-mle(minuslog=ll,start=list(lambda=0.05,k=0.006))

Previously I had managed to get R-0.99.0 to compile and pass the tests on a Debian Linux/Alpha 2.2 system by recompiling src/main/optim.c by hand omitting the -O2 flag. I just tried R-1.0.0 on a similarly configured machine, although not the same machine as before, and it compiled ok but failed the tests. This time it fails on bessel> x <- seq(0, 4, len = 501) bessel> plot(x, x, ylim