similar to: Bessel function

Hello everybody I need some help. Currently, I'm translating a code from Splus to R. In Splus code there is the function nlminb (Nonlinear Minimization subject to Box Constraints) used. Does anybody know an equivalent function in R? Thanks for help Kind regards, Beat Huggler _____________________________________________ Beat Huggler Quantitative Analysis RMF Investment Management

Hi guys I've got an easy question but couldn't find any quick solution. I woulk like to change the following matrix good bad worse Blue 1 2 2 Yellow 2 1 3 Black 3 4 4 Into the following structure good 1 Blue Bad 2 Blue Worse 2 Blue Good 2 Yellow Bad 1 Yellow Worse 2 Yellow Good 2 Black Bad 4 Black Worse 4 Black Thanks a lot. Regards Beat Huggler --- Beat Huggler Quantitative

In order to fit a probability distribution proposed by Sichel [Journal of the Royal Statistical Society. Series A (General), Vol. 137, No. 1. (1974), pp. 25-34], I need a modified Bessel function of the 2nd kind. I notice that the base package of "R" only has modified Bessel functions of the 1st and 3rd kind. Does a modified Bessel function of the 2nd kind exist anywhere? Many

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

I would be very grateful for any help from members of this list for what might be a simple problem... We are trying to simulate the behaviour of a clinical measurement in a series of computer experiments. This is simple enough to do in R if we assume the measurements to be Gaussian, but their empirical distribution has a much higher peak at the mean and the distribution has much longer tails.

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

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.

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,

Dear R People: There is a function from the fBasics library to get the probability and quantiles for the hyperbolic probability function. Is there one that will estimate parms of the hyperbolic probability function from a data set, please? Thanks in advance! Sincerely, Erin Hodgess mailto: hodgess at gator.uhd.edu R Version 2.0.1 windows

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!

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), +

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

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

hi i'm looking for a gibbs sampling algorithm for R for the case of mixture of K normals, and in particular for the case of bivariate normals. i'd be grateful if anyone could send its own R-routine, at least for the univariate case. thank you in advance matteo

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

Hello All, I would like to optimize a (log-)likelihood function subject to a number of linear constraints between parameters. These constraints are equality constraints of the form A%*%theta=c, ie (1,1) %*% 0.8,0.2)^t = 1 meaning that these parameters should sum to one. Moreover, there are bounds on the individual parameters, in most cases that I am considering parameters are bound between zero

Hi guys Recently a wanted to look at the opotunities which are suported by the labrary of tcltk. Unfortunatly I couldn't load the library and it occured the following error message: Error in firstlib(which.lib.loc, package) : TCL_LIBRARY is not set Error in library(pkg, character.only = TRUE) : .First.lib failed I've tried to load it in sveral ways. I reinstalt the

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