search for: hermite

The search at www.r-project.org mentioned a function "gauss.hermite{rmutil}". However, 'install.packages("rmutil")' produced, 'No package "rmutil" on CRAN.' How can I find the current status of "gauss.hermite" and "rmutil"? Thanks, Spencer Graves

Hi, I am trying to find a function that returns simply the weights and points of an n point gauss hermite integeration, so that I can use them to fit a non-standard likelihood. I have found some documentation for the function 'gauss.hermite' written by jim lindley, but can't find the actual binary on CRAN I'm aware there are lots of functions like glmm, glmmML etc to fit mixed models...

Which function implements the piecewise cubic Hermite interpolation? I am looking for equivalent of matlab's interp1 with the method = 'pchip' Here is the reference http://www.mathworks.com/access/helpdesk/help/techdoc/index.html?/access/helpdesk/help/techdoc/ref/interp1.html& -- View this message in context: http://www.nabble.com/...

Does anyone know how to do hermite or Bezier splines in R? I can find Matlab routines, but really need to implement them in R. Failing that I'd be interested in other conformal splines. I need to smooth and interpolate animal tracking data. Antarctic Wildlife Research Unit School of Zoology University of Tasmania PO Box...

Hi all, I am using the 'gaussHermite' function from the 'pracma' library ############ CODES ########### library(pracma) cc=gaussHermite(10) cc$x^2 cc$x^5 cc$x^4 ############ CODES ########### as far so good. However, it does NOT work for any NON integer values, say ############ CODES ########### cc$x^(2.5) cc$x^(-2.5) #...

...en limits are exceeded are not always meaningful. For example: > x <- rep(0,50000*50000) Error in rep(0, 50000 * 50000) : invalid 'times' argument In addition: Warning message: In as.vector(data) : NAs introduced by coercion Here's another example: > library(orthopolynom) > hermite <- hermite.h.polynomials(1001) > hermite[[1001]] # should display 1000-th degree polynomial Error in if (signs[1] == "- ") "-" else "" : missing value where TRUE/FALSE needed Thanks, Dominick

...o use the Gram-Charlier series expansion to model some data. To do that, I need functions to: 1) Calculate 'n' moments from given data 2) Transform 'n' moments to 'n' central moments, or 3) Transform 'n' moments to 'n' cumulants 4) Calculate a number of Hermite polynomials Are there R-functions to do any of the above? (mean, sd and cum3 are very limited) Thank you for your help, Augusto -------------------------------------------- Augusto Sanabria. MSc, PhD. Mathematical Modeller Risk Research Group Geospatial & Earth Monitoring Division Geosc...

...-infty,infty) with weighting function W(x)=exp(-(x-mu)^2/sigma) to estimate an integral. Is there a way to do it in R? Is there a function already implemented which uses such weighting function. I have been searching in the statmode package and I found the function "gauss.quad(100, kind="hermite")" which uses the weighting function W(x)=exp(-x^2). Is there a more general version of this weighting function (using mu and sigma)? Thanks for your help Iazaiez

Can anyone help me on how to get the nodes and weights of the adaptive quadrature using R. Best wishes Boikanyo. ----- The University of Glasgow, charity number SC004401

Does anyone see why my code does not integrate to 1? library(statmod) mu <- 0 s <- 1 Q <- 5 qq <- gauss.quad(Q, kind='hermite') sum((1/(s*sqrt(2*pi))) * exp(-((qq$nodes-mu)^2/(2*s^2))) * qq$weights) ### This does what's it is supposed to myNorm <- function(theta) (1/(s*sqrt(2*pi))) * exp(-((theta-mu)^2/(2*s^2))) integrate(myNorm, -Inf, Inf)

I haven't found who does it for RHEL/CentOS yet, and there's a wrong attribution.... --------------------- Fortune Begin ------------------------ I'm N-ary the tree, I am, N-ary the tree, I am, I am. I'm getting traversed by the parser next door, She's traversed me seven times before. And ev'ry time it was an N-ary (N-ary!) Never wouldn't ever do a binary. (No

Hi, While I browsing the source of the Oggenc utility, I saw a comment says: /* A common error is to have an 18-byte format chunk with the last two * bytes 0. This is incorrect, but sufficiently common that we only warn * about it instead of refusing it. * Please, if you have a program that's creating these 18 byte chunks, send * a bug report to whoever makes it */ And it was

Hi all, We know that Hermite polynomial is for Gaussian, Laguerre polynomial for Exponential distribution, Legendre polynomial for uniform distribution, Jacobi polynomial for Beta distribution. Does anyone know which kind of polynomial deals with the log-normal, Student抯 t, Inverse gamma and Fisher抯 F distribution? Thank you...

Hi, I'm a student studing for Math. & Infor. Ing. in Tirana,Albania and a have a final project in R to finish in a week, so I badly need your help... The topic of the project is "Hermite Interpolation" and a I have allready done a script that finds the approximation in a certan point x, H(x) ~f(x) but i want to find the value of his derivative H'(x)~f'(x) as well and to do this I need the polynomial form of H(x)=a1*x^(2n+1) + a2*x^(2n) + ... + an but R requires x to be...

Hi all, Does anybody know any function that performs gaussian adapative quadrature integration of univariate functions? Thanks in advance, Regards, Caio __________________________________________________ [[alternative HTML version deleted]]

Dear All, I am trying to compile R-3.0 on Cray xe6 (HLRS) HERMIT, no success so far. Here is my experience: I use this to configure and make R: CC="cc" \ CXX="CC" \ F77="ftn" \ FC="ftn" \ CPPFLAGS="-I$PREFIX/include" \ LDFLAGS="-L$PREFIX/lib${LIBDIRSUFFIX}" \ ./configure --prefix=$PREFIX \ --exec-prefix=$PREFIX \ --bindir=$PREFIX/bin \

...in get an error because the function has the elements "0.d0" and "2.d0". How can I change it to run in R? The function can be found in page 230 from http://www.stat.wisc.edu/~mchung/teaching/stat471/stat_computing.pdf Function is as follows: gauher <- function(n) {# Gauss-Hermite: returns x,w so that #\int_-\infty^\infty exp(-x^2) f(x) dx \doteq \sum w_i f(x_i) EPS <- 3.e-14 PIM4 <- .7511255444649425D0 MAXIT <- 10 m <- trunc((n+1)/2) x <- w <- rep(-1,n) for (i in 1:m) { if (i==1) { z <- sqrt(2*n+1)-1.85575*(2*n+1)^(-.16667) } else if(i==2) { z <- z-1...

A new version, 0.65-1, of glmmML is now on CRAN. It is a major rewrite of the inner structures, so frequent updates (bug fixes) may be expected for some time. News: * The Laplace and adaptive Gauss-Hermite approximations to the log likelihood function are fully implemented. The Laplace method is made the default. It should give results you can compare to the results from 'lmer' (for the models that glmmML can handle). * Binomial responses can now be represented as a two-column matrix with No...

A new version, 0.65-1, of glmmML is now on CRAN. It is a major rewrite of the inner structures, so frequent updates (bug fixes) may be expected for some time. News: * The Laplace and adaptive Gauss-Hermite approximations to the log likelihood function are fully implemented. The Laplace method is made the default. It should give results you can compare to the results from 'lmer' (for the models that glmmML can handle). * Binomial responses can now be represented as a two-column matrix with No...

...estimators. I thing that I could use any package that allows to estimate generalized linear mixed model, because above models is equivalent to: x | exp(lambda) ~ Poisson (exp(lambda)) exp(lambda) ~ Normal(c,d) Does it exist any package that can estimate it? May be you know a package that do Gauss-Hermite quadrature for estimation or simply do estimation for the first model? Robert