similar to: Install polynom pack

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

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>

The trailing backslash seems to be the culprit. After looking at rw1020\bin\INSTALL, I noticed $pkg =~ s/\/$//; which I assume is to strip the trailing slash from a package name before obtaining the package name and cleaning out its directory I added $pkg =~ s/\\$//; immediately after this line. Rerunning Rcmd INSTALL C:\rw1011\src\library\bqtl\ with that change seemed to work OK

Hello! I tried to build the RMySQL package by myself on an NT machine. I have a NT machine with cygwin compiler and activestate Perl installed on it. Perl details: D:\prog\rw1011\src\gnuwin32>perl -v This is perl, version 5.005_03 built for MSWin32-x86-object (with 1 registered patch, see perl -V for more detail) Copyright 1987-1999, Larry Wall Binary build 522 provided by ActiveState Tool

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

Dear R People: Here is an interesting question(I think) Suppose I want to calculate Centered Moving Averages; i.e. x[1] <- ( sum(y[1:12]) )/12 x[2] <- ( sum(y[2:13]) )/12 and so on. Of course, this is easily done through loops. However, I have been trying to do this more elegantly, but have failed. I have tried things like j1 <- 1:109 j2 <- 12:120 x[1:109] <- ( sum( y[j1:j2])

Dear R People: I have the source for R version 1.1.1 on a UNIX system alpha-dec-osf4. The ./configure works fine. When I do the make, it works fine until it tries to build the package tcklc. I get an error of Tk_Saf_Init not identified Then it kicks out. Any suggestions, please? Thanks! Sincerely, Erin M. Hodgess, Ph.D. Associate Professor Department of Computer and Mathematical Sciences

Dear R People: Is there a function for exponential smoothing, please? R version 1.4.1 for Windows. Thanks in advance! Have a great weekend! Sincerely, Erin Hodgess Associate Professor Department of Computer and Mathematical Sciences University of Houston - Downtown 1 Main Street Houston, TX 77002 mailto: hodgess at uhddx01.dt.uh.edu

Dear R people is there a module for quality control, please? X bar charts, and such? perhaps a contributed package? Thanks in advance! Erin Hodgess Associate Professor Department of Computer and Mathematical Sciences University of Houston - Downtown 1 Main Street Houston, TX 77002 mailto: hodgess at uhddx01.dt.uh.edu

Dear R People: Here is a very simple question, please: I have a data frame with the state name, state division, and population. Note: state division is something like Atlantic, Mountain, etc. I want to produce a table with the state division and the total of the populations for the states within that division. I've been messing with "table" but so far, no luck. This is R1.3.1,

Hello, R People: How do I change the size of a graph, please? I tried par(cin=c(8.5,11)) and it said that Warning message: parameter "cin" can't be set in: par(args) > par(fin=c(8.5,11)) > plot(xa,ya) Error in plot.new() : Figure region too large Also, can I set for portrait or landscape, please? Thanks in advance! Erin M. Hodgess, Ph.D. Associate Professor Department

Dear R People: I'm trying to use the predict command on an arima0 object. I do the following: xm.arma <- arima0(xm2,order=c(1,0,1)) predict(xm.arma,n.ahead=2) and I get the message: Error in round(x, digits) : Non-numeric argument to mathematical function Any ideas what the problem might be, please? R version 1 4 1 on Windows. Thanks in advance! Sincerely, Erin Hodgess Associate

Dear R People: Do any of the linear model or regression function calculate the Durbin-Watson test statistic for serial correlation, please? I found the hat matrix, studentized residuals, and so on, but no D-W. Thanks so much! Sincerely, Erin M. Hodgess, Ph.D. Associate Professor Department of Computer and Mathematical Sciences University of Houston - Downtown One Main Street Houston, TX 77002

Hello, Although your message is not very explicit about the structure of your data frame, I think you can get what you want by using a combination of the function "split" and the "apply" functions family: apply, lapply and sapply. If you provide a more detailed organization of your data frame, perhaps we can give to you a more wiser advice. Thanks, Carlos. P.S: By the way,

Dear R People: I have R Version 0.90.1 for UNIX (binary) and the Windows version from a binary distribution as well. I would like to look at the source code for the C program setup_starma in the ts library. How would I find that, please? Thanks in advance. Sincerely, Erin M. Hodgess, Ph.D. Assistant Professor Department of Computer and Mathematical Sciences University of Houston - Downtown

Dear R People: Suppose I put in the following: > w1 <- expression(2*x) > v1 <- expression(x^2) and I want to get the derivative of w1 * v1 with respect to x > D(w1*v1,"x") Error in w1 * v1 : non-numeric argument to binary operator Is there a way around this, please? Thanks in advance! Sincerely, Erin M. Hodgess, Ph.D. Associate Professor Department of Computer and

Dear R People: Does anyone have any code for Fractals, chaos, or anything like that, please? This is strictly for demo purposes...decorative only. This is R version 1.5.1 for Windows. Thank you in advance! Have a great weekend! Sincerely, Erin Hodgess Associate Professor Department of Computer and Mathematical Sciences University of Houston - Downtown 1 Main Street Houston, TX 77002

Dear R People: Here is a non-R probability question, please: There is a problem, originally posed from an American game show, called "Let's Make a Deal." A contestant makes a selection of one door, from Doors 1 through 3. He or she knows that there are 2 goats and 1 car behind those doors. Once the contestant makes the selection, say Door 2, then the host, Monty Hall, tells him

I recently wrote a function to compute functions over moving windows. (See below). It varies from your request by using the last N points rather than the last N/2 and next N/2 points, so you may have to modify it to get what you want. It will be in the next gregmisc package, which will go out in a month or so when the publication review process here finishes. (I hate red tape!) -Greg

Hello R People: Suppose I create a function, x.e, via the fix command: > x.e <- fix(x.e) When I type in x.e, the function code appears, along with the name/number of the environment. However, if I create a function, y.e, and source it in, when I type in y.e, the function code only displays(no environment). Why is this, please? (It's not a bug....I just wondered). I am using