similar to: R-alpha: polynom add-on

Hello, I have 3 columns : a, b and a*b I would like to find the pair (a,b) so that a*b is the minimum but not from the points I measured but from the fit of the curve (I have more points that the ones given below but I fit only on this part because I know that the minimum a*b is in this interval). I thought doing it this way : - to fit a*b=f(a) abfit<-lm(ab ~ poly(a,8,raw=T)) - to use the

Dear R people: I am trying to download the package polynom to R version 1.0.1.1. on Windows. Here are my commands and output: > install.packages("polynom",lib="c:\rw1011\library",CRAN="http://www.r-project.org/src/contrib/") Error in start[k]:(start[k + 1] - 1) : NA/NaN argument In addition: Warning message: Download had nonzero exit status in:

Hi, I can not comprehend the linear fitting results of polynoms. For example, given the following data (representing y = x^2): > x <- 1:3 > y <- c(1, 4, 9) performing a linear fit > f <- lm(y ~ poly(x, 2)) gives weird coefficients: > coefficients(f) (Intercept) poly(x, 2)1 poly(x, 2)2 4.6666667 5.6568542 0.8164966 However the fitted() result makes sense: >

I've put all ``current'' add-on packages into CRAN's source/contrib tree and created an INDEX file (attached below). As you can see, currently we have acepack bootstrap ctest date e1071 fracdiff gee jpn snns splines survival4 (Yes, e1071 and jpn are new ... more on the latter in a later mail.) In the near future, I am hoping for the following: oz (Bill

I am using the 0.50 alpha version of R packaged (in 3 parts) by Kurt for Debian Linux. Package: r-base Status: install ok installed Priority: optional Section: local Maintainer: Kurt Hornik <Kurt.Hornik at ci.tuwien.ac.at> Version: 0.50b7-1 Depends: libc5, xlib6, libreadline2 Description: R, a language not entirely unlike the language S. Package: r-contrib Status: install ok

Hello, I have a question regarding fitting a polynomial to a data set, then constructing a polynom from the coefficients so that I can integrate it. I first use lm to fit the polynomial setting the coefficients to raw=TRUE - this appears to work fine. I plot the model and it is a true representation of the data. I then take the coefficients vector and construct a polynom from the

Hello, I'm trying to fit some points with a 8-degrees polynom (result of lm is stored in pfit). In most of the case, it is ok but for some others, some coefficients are "NA". I don't really understand the meaning of these "NA". And the problem is that I can't perform a derivation (pderiv<-as.function((deriv(polynomial(pfit$coefficients))))) on pfit due to the

Dear all, I am trying to transform polynomial coefficients from orthogonal form to the standard power basis. There's poly.transform in S-plus. Does anybody know how to do that in R ? I've found question about that in the archives of R-help but no real answer. Example : I'm doing polynomial regression of percentage of one insect in a community on altitude, precipitations,

Merci beaucoup, Jean, for the bug report -- which I'm no "completeing" to R-bugs >>>>> "Jean" == Jean Coursol <coursol@cristal.math.u-psud.fr> >>>>> on Thu, 24 Jun 2004 15:22:37 +0200 (CEST) writes: Jean> I was exploring the polynom library with students: <and found a segmentation fault from parsing a long expression>

Attached is a snapshot of the new version of the FAQ. What is still missing is something on eval and .Options versus options(). As always, feedback is greatly appreciated. Best, -k *********************************************************************** R FAQ Kurt Hornik v0.2-0, 1997/09/01 This document contains answers to some of the most frequently asked questions about R. Feedback

Please, I'm a beginner with the R language. I'm looking for a function to compute a Polynomfit for simple x-y Data. Who can help? Many Greetings E.A. Welge -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send "info", "help", or "[un]subscribe" (in the

I found the following 2 problems (no debugging, sorry). * qt(0.975, 3) never returns. Can someone please check that? Happens for me on Debian GNU/Linux/ix86 compiled with CFLAGS="-O2 -g". * Yet another stupid way to cause a segfault: R> x Error: Object "x" not found R> debug(t.test) R> t.test(rnorm(10)) debug: choices <- c("two.sided",

Try something like: install.packages("polynom") library(polynom) predict(polynomial(rev(p)), x) HTH, Andy > From: McClatchie, Sam (PIRSA-SARDI) > > > Background: > > OS: Linux Mandrake 9.1 > > release: R 1.9.0 > > editor: Xemacs 21.4 > > frontend: ESS 5.1.23 > > --------------------------------- > > > > Colleagues > >

Something strange is happening with update.packages on my machine (Win2000, R-1.5.0 installed from binaries). Some of the packages (see below) keep getting updated every time. When I issue update.packages() I get this every time for the last few days (I think since I installed 1.5 and the packages but I am not sure). This happens with CRAN shown below and also with CRAN set to

Full_Name: Jens Keienburg Version: 2.3.0 OS: Windows XP Submission from: (NULL) (193.174.53.122) I used the function lm() to calculate the coefficients of a polynome. If I used the function poly(t,2) to denote a polynome of form 1 + x + x^2, the coefficients are wrong. I appended an excerpt below: > t=1:100 > p=-20 - 10 * t + 2 * t^2 > p [1] -28 -32 -32 -28 -20 -8 8

Hi, I was just looking at the "CRAN Daily Package Check Results" [http://cran.r-project.org/src/contrib/checkSummary.html], and realized there are 1000 packages on CRAN as of today (look at row 3 in the table below). Yet another quite extraordinary milestone in R history. Last updated on 2007-04-12 11:48:32 Results for installing and checking packages using the three current flavors

Sorry to bring up this old topic again but ... The 2.0 release of Debian Linux is about to happen and the release team requires copyright statements on any code included in the release. Presently there are 4 packages that I have created from R source code and libraries and contributed to the release. I think this will help with the visibility of R when there is a well-regarded release of an

"Sanjay Patel" <spatel at rotateright.com> wrote: > IIUC, you want to use x86-specific bit-hacks (sbb masking) in cases like > this: > unsigned int foo(unsigned int crc) { > if (crc & 0x80000000) > crc <<= 1, crc ^= 0xEDB88320; > else > crc <<= 1; > return crc; > } To document this for x86 too: rewrite the function

Attached is an updated version of the CRAN announcement, also to be sent out as soon as 0.50 is out. Again, please have a look if possible. -k ************************************************************************ This is the first announcement of the Comprehensive R Archive Network (CRAN) CRAN is a collection of sites which carry identical material, consisting of the

I've already fixed the bugs, but as with the last one, this is not critical enough to allow breaking current R-devel's code freeze. I hope I will have corrected it for 1.5.1.. ## Here is code reproducing the problems; ## I use try(.) whenever I know the current versions of R would ## give an error: library(splines) x <- c(1:3,5:6) y <- c(3:1,5:6) (isP <- interpSpline(x,y))#