similar to: Legendre polynomials

Hi All, I need to evaluate a series expansion using Legendre polynomials. Using the 'orthopolinom' package I can get a list of the first n Legendre polynomials as character strings. > library(orthopolynom) > l<-legendre.polynomials(4) > l [[1]] 1 [[2]] x [[3]] -0.5 + 1.5*x^2 [[4]] -1.5*x + 2.5*x^3 [[5]] 0.375 - 3.75*x^2 + 4.375*x^4 But I can't figure out how to

Hi, Anybody can tell me why I don't have colors in my picture? library(lattice) my.pch<-c(17,16,15,18,20,9,10,12) my.fill<-c('yellow','black','green','red','blue','cyan','magenta','pink') wzorzec1<- c(1,3,2,4,6,4,2,3) wzorzec2<- c(5,3,6,6,2,3,4,1) a<-factor(wzorzec1) b<-factor(wzorzec2)

Dear list, This is quite a specific question requiring the package orthopolynom. This package provides a nice implementation of the Legendre polynomials, however I need the associated Legendre polynomial which can be readily expressed in terms of the mth order derivative of the corresponding Legendre polynomial. (For the curious, I'm trying to calculate spherical harmonics [*]).

Dear list, If have some repeated measurement data which looks something like: time <- rep(1:5 , each=2*4) groups <- rep(c("Case", "Control"), each=4) subjects <- factor(rep(1:(2*4), 5)) responses <- time + rnorm(5*2*4) + as.integer(factor(groups)) data <- data.frame(responses, time, groups, subjects) Now I want to plot each subject in a separate panel:

## Not run: from <- sprintf("<sendmailR@ to <- "<olafm at datensplitter.net>" subject <- "Hello from R" msg <- "It works!" sendmail(from, to, subject, msg, control=list(smtpServer="ASPMX.L.GOOGLE.COM")) ## End(Not run) the above commands are provided in this document ie http://cran.r-project.org/web/packages/sendmailR/sendmailR.pdf

Dear All, I have found in the poly help this sentence: The orthogonal polynomial is summarized by the coefficients, which can be used to evaluate it via the three-term recursion given in Kennedy & Gentle (1980, pp. 343–4), and used in the predict part of the code. My question: which type of orthogonal polynomials are used by this function? Hrmite, legendre.. TIA Giovanni [[alternative HTML

Have you ever worked in R with bivariate polynomials? How did you implement simple operators like addition/multiplication? I found a package called multipol that seems to support these kinds of operators but I do keep receiving error. Check for example the following snippet of code (you can copy & paste) require('orthopolynom') require('polynom') require('multipol')

Hola a todos, Por casualidad alguno que tenga el libro de Legendre de Numerical Ecology with R sabe dónde o cómo bajar los data sets que se mencionan para seguir los ejemplos? Muchas gracias por su ayuda! ---- *Angela Andrea Camargo Sanabria* Estudiante Doctorado en Ciencias Biológicas Laboratorio de Ecología de poblaciones y comunidades tropicales Centro de Investigaciones en Ecosistemas

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]]

There are still places available on the following courses. For more detailed information and a registration form please contact training@mango-solutions.com, or visit our website at www.mango-solutions.com <http://www.mango-solutions.com/> . R for Finance **************************************************** Dates: 21st to 23rd October 2008. Duration: 3 days. Location: London.

BaselR - The new R meeting We are pleased to announce the new R meeting to be held in Basel, Switzerland. BaselR will be held from 6:30-9:30pm on Tues, Apr 27 at TransBARent: http://transbarent.business.sv-group.ch Doors open at 6:30,pm with the presentations starting at 7:00pm Introduction: What is Basel R? Andreas Krause:... Graphing Pharma Data Yann Abraham: Graphics Charles Roosen: Web

Dear all, Lattice provides automatic coloring for subgroups on each panel by the simple use of a groups statement. For an application I want to change these colors to a predifined set. This works well using a panel function in stead of the default as long as there are only points in the graphs. When I set type="b" things get messed up. Any idea why? I include sample code for

Good Evening, I would like to work with multivariate polynomials (x and y variables). I know that there is a package called multipol but I am not sure that supports my needs. I use a function (in reality legendre.polynomials) which creates me the polynomials I want. For example the following returns > legendre.polynomials(2)[[2]] x (first order polynomial) I would like to calculate the

Greetings to all, and my apologies for a question that is mostly about statistics and secondarily about R. I have just started a new job that (this week, apparently) requires statistical knowledge beyond my training (as an epidemiologist). The problem: - We have 57 food production facilities in three categories - Samples of 4-6 different foods were tested for listeria at each facility - I need

Pierre Legendre has developed a beta version of a new redundancy analysis package called RdaTest that is available on his web page at the Universit® de Montréal. The test example that is included with the package is based on the example provided in his book (Numerical Ecology, Chapter 11 (Legendre & Legendre 1998)) I have downloaded the package and am attempting to run it so that I might

Hi, I would like to solve a double integral of the form \int_0^1 \int_0^1 x*y dx dy using Gauss Quadrature. I know that I can use R's integrate function to calculate it: integrate(function(y) { sapply(y, function(y) { integrate(function(x) x*y, 0, 1)$value }) }, 0, 1) but I would like to use Gauss Quadrature to do it. I have written the following code (using R's statmod package)

I am trying to use the cca/rda/capscale functions in vegan to analyse genetic distance data ( provided as a dist object calculated using dist.genpop in package adegenet) with geographic distance partialled out ( provided as a distance object using dist function in veganthis method is attempting to follow the method used by Geffen et al 2004 as suggested by Legendre and . FORTIN (2010). I

Dear R-Help Does "The R Package for Multivariate and Spatial Analysis Version 4.0 (Casgrain and Legendre, 2001)" exist on CRAN and under what name? It supposedly has a chronological clustering program ,CHRONO, that I would like to use. Alternatively, I would ask if there is a R based program that performs chronological clustering? Thanks Alex Alex Hanke Department of Fisheries and

Hi every body: I need some help with kendall.global. The example in the manual seems not working well, and cannot used with my data, always the same error. data(mite) > mite.hel <- decostand(mite, "hel") > > # Reproduce the results shown in Table 2 of Legendre (2005), a single group > mite.small <- mite.hel[c(4,9,14,22,31,34,45,53,61,69),c(13:15,23)] >

Hi, I hope someone knows the answer to this or has a real good reference about it (I am using Legendre & Legendre, Numerical Ecology, 1998).... My data is a data.frame with locations as rows and vegetation assemblages / species as columns. I've done a PCA, a correspondance analysis (CA) using ca in ca package and a detrended correspondance analysis (DCA) using decorana from vegan package.