search for: hermits

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

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:

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 252-05 Hobart TAS, 7001 Australia Phone: (0)3 6226 2645 Fax:

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) ############ CODES ########### But just think about it

Is there documentation on R limits? That is, max matrix size, etc.? Diagnostics when 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 <-

Good day everyone, I want to 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?

Hi All, I am trying to use A Gaussian quadrature over the interval (-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

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)

...t 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 sir!) I'm 'er eighth tree that was N-ary. N-ary the tree I am, I am, N-ary the tree I am. -- Stolen from Paul Revere and the Raiders Hermans' Hermits. Paul Revere and the Raiders never did it.... mark "somewhere, somebody on the 'Net said something wrong!"

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 in advance! David [[alternative HTML version deleted]]

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

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 \

I have a function written for Splus, when I run it in R I obtain 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

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

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

Hi R Users, I'm going to estimate via. ML the parameters in Poisson Lognormal model. The model is: x | lambda ~ Poisson(lambda) lambda ~ Lognormal(a,b) Unfortunately, I haven't found a useful package allowing for such estimation. I tried to use "poilog" package, but there is no equations and it's hard to understand what exactly this package really does. Using it I get the