search for: draper

Dear all, Can any one tell me how can i perform Orthogonal Polynomial Regression parameter estimation in R? -------------------------------------------- Here is an "Orthogonal Polynomial" Regression problem collected from Draper, Smith(1981), page 269. Note that only value of alpha0 (intercept term) and signs of each estimate match with the result obtained from coef(orth.fit). What went wrong? -------------------------------------------- Data: > Z<-1957:1964 > Y<-c(0.93,0.99,1.11,1.33,1.52,1.60,1.47,1.33) &g...

...> > On Tue, Aug 19, 2014 at 1:40 PM, Jun Koi <junkoi2004 at gmail.com > > <mailto:junkoi2004 at gmail.com>> wrote: > > > > > > > > > > > > On Mon, Aug 18, 2014 at 9:18 PM, Carback, Richard T., III > > <rcarback at draper.com <mailto:rcarback at draper.com>> wrote: > > > > Hi All, > > > > > > > > 2 of my summer interns (Aimee Dipietro and Greg Simpson) used > > their time over the summer to resurrect the LLVM C Backend: > > > > &g...

I'm sure this must be a nls() newbie question, but I'm stumped. I'm trying to do the example from Draper and Yang (1997). They give this snippet of S-Plus code: Specify the weight function: weight < - function(y,x1,x2,b0,b1,b2) { pred <- b0+b1*x1 + b2*x2 parms <- abs(b1*b2)^(1/3) (y-pred)/parms } Fit the model gmfit < -nls(~weight(y,x1,x2,b0,b1,b2), observe,list("starting value&quo...

...ld be in modern C, C++, or any other high level language or assembly language with a translator to LLVM. On Tue, Aug 19, 2014 at 1:40 PM, Jun Koi <junkoi2004 at gmail.com<mailto:junkoi2004 at gmail.com>> wrote: On Mon, Aug 18, 2014 at 9:18 PM, Carback, Richard T., III <rcarback at draper.com<mailto:rcarback at draper.com>> wrote: Hi All, 2 of my summer interns (Aimee Dipietro and Greg Simpson) used their time over the summer to resurrect the LLVM C Backend: https://github.com/draperlaboratory/llvm-cbe Improvements include recovery of simple for/while loops (instead of g...

...program could be in modern C, C++, or any other high level language or assembly language with a translator to LLVM. On Tue, Aug 19, 2014 at 1:40 PM, Jun Koi <junkoi2004 at gmail.com> wrote: > > > > On Mon, Aug 18, 2014 at 9:18 PM, Carback, Richard T., III < > rcarback at draper.com> wrote: > >> Hi All, >> >> >> >> 2 of my summer interns (Aimee Dipietro and Greg Simpson) used their time >> over the summer to resurrect the LLVM C Backend: >> >> >> >> https://github.com/draperlaboratory/llvm-cbe >> >...

Hi All, 2 of my summer interns (Aimee Dipietro and Greg Simpson) used their time over the summer to resurrect the LLVM C Backend: https://github.com/draperlaboratory/llvm-cbe Improvements include recovery of simple for/while loops (instead of goto), better variable naming, inline asm support, and making it work on a more recent version of llvm. I believe they used the repository here as a starting point: https://github.com/glycerine/llvm/tree/cbe_re...

...ying to use R that someone installed in his own home directory (which is located on a common filesystem). Now, he is able to use it with no problems, but when I try to run it I get the following messages (the first line is just the command line, and the way I call the program): an2 1> /projects/draper/R-1.2.1/bin/R cannot find system Renviron R : Copyright 2001, The R Development Core Team Version 1.2.1 (2001-01-15) R is free software and comes with ABSOLUTELY NO WARRANTY. You are welcome to redistribute it under certain conditions. Type `license()' or `licence()' for distribution deta...

Dear R-experts! I have y(x) data, dim(y) == dim(x) == c(2000, 2) I'd like to fit them with nls: fit.result <- nls ( y ~ f(x, p1, p2, p3), start = list(p1 = ... , p2 = .. , p3 = ..) ) Actually I want to fit y[,1] ~ x[,1] and y[,2] ~ x[,2] *simulaneously*, with the same parameters set {p1, p2, p3}. I tried to feed R tha above formula, R errors with:

Hello all, How do I actually use the output of predict.lm(..., type="terms") to predict new term values from new response values? I'm a chromatographer trying to use R (2.15.1) for one of the most common calculations in that business: - Given several chromatographic peak areas measured for control samples containing a molecule at known (increasing) concentrations, first

I'm trying to fit a model to solve a biological problem. There are multiple independent variables, and also there are multiple responses. Each response is a function of all the independent variables, plus a set of parameters. All the responses depend on the same variables and parameters - just the form of the function changes to define each seperate response. Any ideas how I can fit

...how this situation has come about, why it is important to the future of research, and what can be done about it. Robert McGreevy will draw on his extensive experience at the STFC ISIS Facility, and explain how these issues are being addressed for the benefit of research science globally. Nicholas Draper, consultant at Tessella, will then expand on this, using the example of the Mantid project at ISIS. " Tessella (www.tessella.com) is a technology and consultancy firm, based in Oxford. ISIS (International Species Information System) (www.isis.org) has as its mission the facilitation of &quo...

I''ve seen a few examples of dummy Rails apps (for testing, so they live under test or spec dirs, typically) for use with the Appraisals gem that supposedly work with both Rails 3.x and Rails 4, but they seem hackish and not fully functional. It is somewhat expected, as it is a stripped down Frankenstein monster that is trying to be compatible with various versions of Rails 3 as well

Hi R users! One thing I cannot understand with R is the frequently appearing " Correlation of coefficients". I do not find a through going explanation of this concept either on the help pages of R or on the internet. A Google search has shown many examples of R/S-Plus outputs but no explanations of meaning. I am afraid this turns out to be a silly and unworthy question. Would anyone be

What is the difference between the two R-squareds returned for a linear regression by summary(lm...)? When might one report multiple vs. adjusted R-squared? Thank you, Ben Osborne -- Botany Department University of Vermont 109 Carrigan Drive Burlington, VT 05405 benjamin.osborne at uvm.edu phone: 802-656-0297 fax: 802-656-0440

...IX: <cy@FreeBSD.org> Web: http://www.FreeBSD.org BC Government: <Cy.Schubert@gov.bc.ca> "Lift long enough and I believe arrogance is replaced by humility and fear by courage and selfishness by generosity and rudeness by compassion and caring." -- Dave Draper

...IX: <cy@FreeBSD.org> Web: http://www.FreeBSD.org BC Government: <Cy.Schubert@gov.bc.ca> "Lift long enough and I believe arrogance is replaced by humility and fear by courage and selfishness by generosity and rudeness by compassion and caring." -- Dave Draper ------- Forwarded Message Date: Wed, 16 Nov 2005 20:24:21 -0500 From: Tom Yu <tlyu@mit.edu> To: kerberos-announce@mit.edu Subject: krb5-1.4.3 is released - -----BEGIN PGP SIGNED MESSAGE----- The MIT Kerberos Team announces the availability of MIT Kerberos 5 Release 1.4.3. Thi...

...IX: <cy@FreeBSD.org> Web: http://www.FreeBSD.org BC Government: <Cy.Schubert@gov.bc.ca> "Lift long enough and I believe arrogance is replaced by humility and fear by courage and selfishness by generosity and rudeness by compassion and caring." -- Dave Draper

Dear all, Does any one know if any R package or function can do Ordinary Least Products regression? Many thanks! Bill

...> lack of fit. The error term will be partitioned into pure error (error within replicates) > , and a lack of fit error and the F-test can be used to test if you have chosen an > adequate regression model. See Neter, Kutener, Nachtsheim, and Wasserman fourth edition > page 115, or Draper and Smith. Bill Venables wrote "...It makes it > impossible to write code to do it automatically, but if you know > what you are doing, the procedure is simple with the software you > have. As with so many things in statistics, it is not a matter > of good software so much as of...

Dear friends. On the Bayesian averaging homepage http://www.research.att.com/~volinsky/bma.html I found some S code some of which perhaps may run in R. There was a call to an algorithm possibly within S but not supported by R 64.1: "leaps and bounds". I guess it is a minimization step. Can anyone clarify the algorithm and perhaps even give a pointer to some code ? I guess this may be