Displaying 20 results from an estimated 200 matches similar to: "Buglet/omission in nls package (PR#6901)"
2009 Nov 17
1
Define lm/glm object without evaluating them
For e.g. lm/glm type models I would like to separate model specification and model fitting and then only fit the models later 'when data arrives'. To be specific, I would like make a specification like
m1 <- lm(rate~conc)
m2 <- lm(rate~I(conc^2))
and then later I want to 'put data into' the objects and evaluate (fit the model), e.g. something like
update(m1, data=Puromycin)
2006 Nov 29
3
R2.4 xyplot + panel.number problem
Hi all;
I'm trying to display a 2 panel plot for the Puromycin data from R
with 2 different non-linear models fitted to each group. The problem
is that as far as I know panel.number doesn't work in the latest
version of R. Can anyone give a hint how to solve this?
Here is the code that I used before and now doesn't work
xyplot(rate ~conc| state,Puromycin,
2011 Apr 20
1
How can I 'predict' from an nls model with a fit specified for separate groups?
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,
2003 Oct 28
2
Confidence ellipse for correlation
Hello,
SAS' point and click interface has the option of produce a scatterplot with a
superimposed confidence ellipse for the correlation coefficient. Since I
generally like R so much better, I would like to reproduce this in R. I've
been playing with the ellipse package. In order to have the points and the
ellipse on the same graph I've done the following.
(Load ellipse
2003 Oct 30
2
'nls' and its arguments
Dear R experts!
I'd to fit data by 'nls' with me-supplied function 'fcn'.
1) I'd like 'fcn' to accept arbitrary arguments, i.e. I defined it
as f(...) {<body>}. (Ok, that's not actually impotant).
2) Second, I would NOT like to supply every parameter in the formula.
To illustrate this, let's look at the last example of 'nls' help
2007 Jun 07
2
Nonlinear Regression
Hello
I followed the example in page 59, chapter 11 of the 'Introduction to R'
manual. I entered my own x,y data. I used the least squares. My function has
5 parameters: p[1], p[2], p[3], p[4], p[5]. I plotted the x-y data. Then I
used lines(spline(xfit,yfit)) to overlay best curves on the data while
changing the parameters. My question is how do I calculate the residual sum
of squares.
2004 Feb 02
1
Robust nonlinear regression - sin(x)/x?
You reall have only one parameter in your model, c = a/b. You can't
identify both a and b from your model, therefore, you should fit the
linear model: lm(z ~ c* sin(x)/x)
Ravi.
----- Original Message -----
From: cstrato <cstrato at aon.at>
Date: Monday, February 2, 2004 2:28 pm
Subject: [R] Robust nonlinear regression - sin(x)/x?
> Dear all
>
> Since I did not receive any
2008 Nov 03
2
standard errors for predict.nls?
Dear all,
Is there a way to retrieve standard errors from nls models? The help page tells me that arguments
such as se.fit are ignored...
Many thanks and best wishes
Christoph
--
Dr. rer.nat. Christoph Scherber
University of Goettingen
DNPW, Agroecology
Waldweg 26
D-37073 Goettingen
Germany
phone +49 (0)551 39 8807
fax +49 (0)551 39 8806
Homepage http://www.gwdg.de/~cscherb1
2005 Feb 02
4
(no subject)
can you recommend a good manual for R that starts with a data set and gives
demonstrations on what can be done using R? I downloadedR Langauage
definition and An introduction to R but haven't found them overly useful.
I'd really like to be able to follow some tutorials using a dataset or many
datasets. The datasets I have available on R are
Data sets in package 'datasets':
2005 Jul 11
2
Weighted nls
Dear R Community,
I am attempting to perform a weighted non-linear least squares fit. It has already been noted that the weights option is not yet implemented for the nls function, but no one seems to offer any suggestions for getting around this problem. I am still curious if a) anyone has code they have written which includes a weight options for nls, or b) if there is another model which
2004 Nov 08
2
Nonlinear weighted least squares estimation
Hi there,
I'm trying to fit a growth curve to some data and need to use a weighted least squares estimator to account for heteroscedasticity in the data. A weights argument is available in nls that would appear to be appropriate for this purpose, but it is listed as 'not yet implemented'. Is there another package which could implement this procedure?
Regards,
Robert Brown
2008 Jan 22
2
extension to nlme self start SSmicmen?
Dear list,
Has anyone created a version of SSmicmen that allows testing for group
differences? The basic Michaelis-Menten equation is:
(Bmax * X) / (Kd + X).
The nlme package allows modeling of random effects for Bmax and Kd as
needed, but I curious how I can build in group differences? I have
receptor binding data for strains of mice, and following Pinheiro and
Bates' lead in their
2004 Feb 04
1
Fitting nonlinear (quantile) models to linear data.
Hello.
I am trying to fit an asymptotic relationship (nonlinear) to some
ecological data, and am having problems. I am interested in the upper
bound on the data (i.e. if there is an upper limit to 'y' across a range
of 'x'). As such, I am using the nonlinear quantile regression package
(nlrq) to fit a michaelis mention type model.
The errors I get (which are dependant on
2008 Apr 03
1
Omission in sshd_config man page
[Not subscribed to this list, so please respond directly if you need to speak to me]
In man5/sshd_config.5, a permissible keyword in a 'Match' block is missing. It currently lists only:
AllowTcpForwarding, Banner, ForceCommand, GatewayPorts, GSSApiAuthentication, KbdInteractiveAuthentication, KerberosAuthentication, PasswordAuthentication, PermitOpen, PermitRootLogin,
2009 Sep 06
0
[LLVMdev] omission: clang-cc -MD -MF -c -fPIC
On Sep 5, 2009, at 11:36 AM, Bruce Korb wrote:
> Hi,
>
> I do not really have much direct control over the compile lines any
> more.
> I guess I could code up a wrapper for your wrapper so that the
> libtool wrapper
> would be happy, but could you all consider libtool built projects and
> embed the necessary option stripping in your clang-cc thing?
> Thank you.
>
2001 Apr 20
1
[patch] one warning, one omission, and two requests; portable CVS
Requests first: after a build, "make" always remakes all the .out files;
maybe we could add a few dependencies there? And could we possibly have a
few more '@if' lines in the Makefile to quiet it down?
These should be obvious:
--- #readpass.c Fri Apr 20 13:17:39 2001
+++ readpass.c Fri Apr 20 13:37:26 2001
@@ -43,7 +43,7 @@
#include "ssh.h"
char *
2017 Sep 12
0
CyberPower UPS Omission
On Sep 11, 2017, at 10:09 PM, Jason Berg <jbswman at gmail.com> wrote:
> I am happy to provide you the output of "upsc" but I do not know what that is. "Does it mean UPS console"?
Basically - "UPS client" was the original author's intent, I think.
https://doc.freenas.org/9.3/freenas_services.html#ups says that you can run it from Shell, which I assume
2009 Dec 01
1
documentation omission in Extract.Rd?
The help file for "[<-" (and "[[<-" doesn't contain those functions in
it's usage section ( .../src/library/base/man/Extract.Rd ). Although
that file lists "[<-" and "[[<-" as aliases, it doesn't document them in
a formal manner (consequently, the "value" argument is not described in
the "arguments" section.)
2017 Sep 11
2
CyberPower UPS Omission
The Hardware Compatibility List page on our website said to report any
omissions in the list to this address. I bought a CyberPower UPS Model
*LX1500GU*. It looks similar to the *CP1500AVRLCD*. My intention is to
configure it in FreeNAS (FreeBSD). The FreeNAS user guide says it uses NUT
to provide UPS support an contains a link to your website for a list of
drivers of supported UPS devices. I
2011 Feb 06
1
anova() interpretation and error message
Hi there,
I have a data frame as listed below:
> Ca.P.Biomass.A
P Biomass
1 334.5567 0.2870000
2 737.5400 0.5713333
3 894.5300 0.6393333
4 782.3800 0.5836667
5 857.5900 0.6003333
6 829.2700 0.5883333
I have fit the data using logistic, Michaelis?Menten, and linear model,
they all give significance.
> fm1 <- nls(Biomass~SSlogis(P, phi1, phi2, phi3), data=Ca.P.Biomass.A)