similar to: nls, restrict parameter values

Hi, I am trying to fit a sigmoid curve to some data with nls but I am getting into some trouble. Seems that the optimization method is getting down to some parameter estimates that make the equation unsolvable. This is an example: >growth<-data.frame(Time=c(5,7,9,11,13,15,17,19,21,23,25,27),BodyMass=c(45,85,125,210,300,485,570,700,830,940,1030,1120))

Hi all I try to make a weibull survival analysis on R. I know make this on GLIM, and now I try to make the GLIM exercice GLEX8 on R to learning and compare the test. The variables are: time censor group bodymass In GLIM I make: $calc %s=1 $ to fit weibull rather than exponential $input %pcl weibull $ $macro model group*bodymass $endmac$ $use weibull t w %s $ Then, GLIM estimate an alpha for the

Hi The best is definitively (X)Emacs (my opinion). For Windows have a look at (GNU Emacs): http://www.math.auc.dk/~dethlef/Tips/introduction.html And for XEmacs with Cygwin (XEmacs, Latex, ..., for Ess simple save (require 'ess-site) in your init.el): http://home.imf.au.dk/hellmund/indhold.html For native Xemacs for Windows: www.xemacs.org Xemacs is easy to install, use and configure

Which is the best editor for .R files? I currently use kate on my linux as it has R highlighting and allows me to split the window into two: in one I edit the .R file and in the other I have a shell so I run R and can easily copy and paste the code. There are some features that I don't like and I am having a look on some alternatives. I've heard wonders of emacs with ess but I am a

Why do x<-b%*%ginv(A) and x<-solve(A,b) give different results?. It seems that I am missing some basic feature of matrix indexing. e.g.: A<-matrix(c(0,-4,4,0),nrow=2,ncol=2) b<-c(-16,0) x<-b%*%ginv(A);x x<-solve(A,b);x Thanks in advance, Angel

I want to 1) Select for each of the n points in a matrix A, those of the m points in B that lay within a given radius. 2) Of those points within the radius, select the k nearest ones. What I now do is 1) Create an n*m matrix C were I put the distances from all the points in B to the points in A and make NA those cells were the distance is larger than the radius. (The points are geographical

I needed a symbolic and algebraic computation software. The best solution I've found so far is maxima computer algebra system. I was wondering if there is any R package that has similar features as maxima's or somebody could point me to reference manuals on how to use R for this purpose. Thanks, Angel

I have a WLS regression with 1 dependent variable and 3 independent variables. I wish to constrain the predicted values (the fitted values) so that they are greater than zero (i.e. they are positive). I do not know how to impose this constraint in R. Please respond if you have any suggestions. There are some previous postings about constraining the coefficients, but this won't accomplish

[Apologies for the slightly-off-topic posting.] BizConf is a one-of-a-kind event that I am hosting about 3 months from now here in northern Florida. It''s a 2-day conference specifically crafted for independent consultants and owners/managers of small to mid-sized web consultancies. Because of who I am and the audience I can reach, I expect the majority of the attendees to be involved

I am using R on a Windows XP professional platform. The following code is part of a bigger one CODE press=function(y,x){ library(qpcR) models.press=numeric(0) cat("\n") dep=y print(dep) indep=log(x) print(indep) yfit=dep-PRESS(lm(dep~indep))[[2]] cat("\n yfit\n") print(yfit) yfit.orig=yfit presid=y-yfit.orig press=sum(presid^2)

hello short description: ssh client sends a wrong all-zero terminal info when requesting forced pseudo-terminal allocation while stdin is not a terminal. sshd then sets the terminals ospeed to 0 because it receives this information from the client. on solaris this means that the pseudo-terminal is closed and output of any remote command is dropped. longer description: what we're doing is

Hi. I have a question about activerecord, computed columns and to_xml. I have a simple model called person that has a first_name and a last_name field (and thus a first_name and last_name field in the Person table). When I :find a model and call to_xml, I would like the result to not only send back first name and last name as xml but also add an additional attribute called full_name that is

I am looking into fitting a so-called double von Bertalanffy function to fish length-at-age data. Attempting to simplify the situation, the model looks like this ... Y ~ f(X; a,b,c) if x < Z Y ~ g(X; a,d,e) if x >= Z where * f and g are non-linear functions (the "traditional" "single" von Bertalanffy growth function), * Y (length) and X (age) are observed variables,

http://bugzilla.mindrot.org/show_bug.cgi?id=1199 Summary: ssh sends invalid ttymodes when stdin is not a tty Product: Portable OpenSSH Version: -current Platform: All OS/Version: All Status: NEW Severity: normal Priority: P2 Component: ssh AssignedTo: bitbucket at mindrot.org ReportedBy:

Hi, the following R lines work fine in R 2.4.0 alpha (and older R versions), but not in R 2.4.0 beta (details below): library(drc) # to load the dataset 'PestSci' library(nlme) ## Starting values sv <- c(0.328919, 1.956121, 0.097547, 1.642436, 0.208924) ## No error m1 <- nlme(SLOPE ~ c + (d-c)/(1+exp(b*(log(DOSE)-log(e)))), fixed =

Following an example on p 111 in 'Nonlinear Regression with R' by Ritz & Streibig, I have been fitting nls models using square brackets with the grouping variable inside. In their book is this example, in which 'state' is a factor indicating whether a treatment has been used or not: > Puromycin.m1 <- nls(rate ~ Vm[state] * + conc/(K[state] + conc), data = Puromycin,

http://bugs.freedesktop.org/show_bug.cgi?id=16148 Summary: primary display is borked with RandR12 enabled Product: xorg Version: unspecified Platform: Other OS/Version: All Status: NEW Severity: normal Priority: medium Component: Driver/nouveau AssignedTo: nouveau at lists.freedesktop.org

> Date: Friday, August 28, 2015 16:47:43 +0000 > From: Tony Mountifield <tony at softins.co.uk> > > In article > <CAOZy0enqddiPvpd+M-Ltwih9dPmA7b_ro4-_5bQ=u1GAALDebA at mail.gmail.co > m>, Tim Dunphy <bluethundr at gmail.com> wrote: >> Hey guys, >> >> Sorry for the failed attempts at obscuring the company I work >> for. My boss

Hello, I'm a beginner with R and it's the first time I'm using the R-help list... I hope I'm in the right place, if not: Sorry!! I need to do non linear regressions on a data set which columns are: "river.name" "Portata" "PTG.P" "PO4.P" "NT.N" "NH4.N" "NO3.N" "BOD5" "SiO2"

Assuming I have a firewire card with a supported chipset, is there anything special I need to do in order for the kernel to recognize the card? Shouldn't kudzu "wake up" and configure the card when I boot up after installing it? This is a brand spanking new 4.4 system and a supported Agere chipset firewire card. I'm having trouble tickling devices connected to the card.