similar to: nls() newbie convergence problem

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

I can't see why you'd want to do this, no. -eric On Tue, Aug 19, 2014 at 8:02 PM, Isaac Dupree < ml at isaac.cedarswampstudios.org> wrote: > Is the C backend at all suitable to be adapted to emit OpenCL code? Or > do the target-dependence, and/or things that C can do but OpenCL can't, > make that hopeless? > -Isaac > > On 08/19/2014 03:08 PM, Carback,

It provides a useful starting point, but I agree with Jim that it is not a complete solution and requires rework of the results in a lot cases. I think we could improve it further to address these issues but that work is nontrivial. If you are deciding between a quick and dirty implementation of a custom backend vs. the C backend, then the C backend is sometimes preferable in my experience

An embedded and charset-unspecified text was scrubbed... Name: n?o dispon?vel URL: <https://stat.ethz.ch/pipermail/r-help/attachments/20080829/57df9fc7/attachment.pl>

I don't know how good it is, but the applications seem obvious. e.g. compiling programs in any of a number of original formats to run natively on CPUs that have a working simple C compiler (maybe only K&R or C89) but don't have an LLVM back end. The source program could be in modern C, C++, or any other high level language or assembly language with a translator to LLVM. On Tue, Aug

Hi. I'm trying 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

Hi, I have a data.frame in wide format which I would like to reshape into a long format: example (nonsense) data: > dput(perus2) structure(list(id = c(30L, 38L, 21L, 12L, 22L, 28L, 31L, 44L, 8L, 47L, 23L, 20L, 41L, 42L, 29L, 50L, 5L, 33L, 4L, 17L, 11L, 1L, 18L, 6L, 9L, 32L, 16L, 14L, 39L, 48L, 37L, 43L, 25L, 27L, 35L, 40L, 45L, 49L, 2L, 34L, 10L, 15L, 24L, 26L, 46L, 7L, 13L, 3L, 19L, 36L),

Hi: How do I find average across three columns such as weight1, weight2, weight3 (20 people, each person has three different measures)? They have NA's also. Thanks Tom, Univ of Pittsburgh, ttrut at yahoo.com

I came across this notice of an upcoming webinar. The issues identified in the first paragraph below seem to me exactly those that the R project is designed to address. The claim that "most research software is barely fit for purpose compared to equivalent systems in the commercial world" seems to me not quite accurate! Comments! WEBINAR SERIES A Crack in the Code: Why software

Hi, I am an R newbie and one thing I am having trouble with binding variables that I have created within one data frame into a new data frame when using na.omit(). To illustrate this problem I will give the example I am working on and the approah I have been using:- data.frame1<-filepath.... attach(data.frame1) #create a new variable using a function new.variable<-rep(1,length(weight3))

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

I'm trying to plot a series of densities using/comparing differing weights. I see the reference to weights and subscripts, but I don't understand how to implement that. My data are of the form: I, J, Actual, Distance, Subset, Weight1, Weight2, ... I'm trying to see the effect of the distance distribution (Actual by Distance) compared to the various weighted distributions

Hi all, I have a data set with missing value. I would like to estimate those missing value by using normal ratio method. Below is part of my data: AS BL Serdang Jhr Phg Target station 0 0.0 12.8 0.0 23.7 0.0 6 0.0 81.7 0.2 0.0 NA 0 1.5 60.9 0.0 0.0 15.5 1 13.0 56.8 17.5 32.8 6.4 4 3.0 66.4

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

I have an updated tripwire port which I'd like to release for testing prior to replacing the existing and currently broken tripwire port. I've tested it under the upcoming 6.1, which it seems to work well in my test environment and will commence testing under 4.11-STABLE, the upcoming 5.5, and 7.0-CURRENT. The things that are on my todo list are: - Fully test under 4.11-STABLE. -

I will be updating the KRB5 port shortly. The updated port, as it stands right now, builds and installs cleanly on the i386 version of 6.0. I need to verify it builds and installs on 5.4, 4.11, and 7.0-CURRENT (my "checkout" of 7.0 as of three days ago currently doesn't build here). I should have the upgraded port committed sometime during the weekend. Cheers, Cy Schubert

While the ports freeze is still in effect, I have a new MIT KRB5-1.4 port ready for commit. Anyone willing to use this preview version of the port or test it is welcome to use it. I've put it at http://komquats.com/ports/krb5- 1.4.tar.bz2. If you have any problems or issues with it, just send me an email. Cheers, Cy Schubert <Cy.Schubert@komquats.com> Web: http://www.komquats.com

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