Displaying 20 results from an estimated 200 matches similar to: "Deviance and AIC in weighted NLS"
2012 Sep 19
0
Discrepancies in weighted nonlinear least squares
Dear all,
I encounter some discrepancies when comparing the deviance of a weighted and
unweigthed model with the AIC values.
A general example (from 'nls'):
DNase1 <- subset(DNase, Run == 1)
fm1DNase1 <- nls(density ~ SSlogis(log(conc), Asym, xmid, scal), DNase1)
This is the unweighted fit, in the code of 'nls' one can see that 'nls'
generates a vector
2004 Jul 16
1
Does AIC() applied to a nls() object use the correct number of estimated parameters?
I'm wondering whether AIC scores extracted from nls() objects using
AIC() are based on the correct number of estimated parameters.
Using the example under nls() documentation:
> data( DNase )
> DNase1 <- DNase[ DNase$Run == 1, ]
> ## using a selfStart model
> fm1DNase1 <- nls( density ~ SSlogis( log(conc), Asym, xmid, scal ),
DNase1 )
Using AIC() function:
>
2004 Jul 16
0
Does AIC() applied to a nls() object use the correctnumber of estimated parameters?
Thanks Adaikalavan, however the problem remains.
Considering AIC() as applied to the linear model in AIC() help
documentation:
> data(swiss)
> lm1 <- lm(Fertility ~ . , data = swiss)
> AIC(lm1)
[1] 326.0716
Clearly this includes the estimation of the residual standard error as
an estimated parameter, as this gives the correct score:
> -2*logLik(lm1) + 2*(length(coef(lm1))+1)
2007 May 31
1
predict.nls - gives error but only on some nls objects
Dear list,
I have encountered a problem with predict.nls (Windows XP, R.2.5.0), but I am not sure if it is a bug...
On the nls man page, an example is:
DNase1 <- subset(DNase, Run == 1)
fm2DNase1 <- nls(density ~ 1/(1 + exp((xmid - log(conc))/scal)),
data = DNase1,
start = list(xmid = 0, scal = 1))
alg = "plinear", trace =
2009 Nov 09
1
Parameter info from nls object
Hi!
When checking validity of a model for a large number
of experimental data I thought it to be interesting
to check the information provided by
the summary method programmatically.
Still I could not find out which method to
use to get to those data.
Example (not my real world data, but to show the point):
[BEGIN]
> DNase1 <- subset(DNase, Run == 1)
> fm1DNase1 <- nls(density ~
2006 Sep 11
4
syntax of nlme
Hello,
How do I specify the formula and random effects without a startup object
? I thought it would be a mixture of nls and lme.
after trying very hard, I ask for help on using nlme.
Can someone hint me to some examples?
I constructed a try using the example from nls:
#variables are density, conc and Run
#all works fine with nls
DNase1 <- subset(DNase, Run == 1 )
fm2DNase1 <- nls(
2017 Apr 01
6
Intervalos de confianza de la varianza de los residuos en un modelo no lineal.-
Hola amigos,
Supongamos que se quiere ejecutar un modelo no lineal con nls. Pensemos en
el ejemplo de la ayuda:
DNase1 <- subset(DNase, Run == 1)
fm1DNase1 <- nls(density ~ SSlogis(log(conc), Asym, xmid, scal), DNase1)
summary(fm1DNase1)
Aquí se está modelando la densidad óptica de un ensayo relacionada de forma
no lineal (logística) con (el logaritmo) de la concentración de una
proteína.
2009 Jul 23
1
Network from package functions
Dear R-helpers,
does anyone know of some package/function that can build a network from
the functions that are implemented in a package,
i.e. visualize the cross-references from one function to another in the
same or some dependent package?
An example would be a function like 'nls' on top of the hierarchy and
then a network of nodes from the functions that are called within 'nls'
2018 May 05
0
Bug in profile.nls with algorithm = "plinear"
Dear sirs
It seems like there is a bug in `profile.nls` with `algorithm =
"plinear"` when a matrix is supplied on the right hand side. Here is
the bug and a potential fix
#####
# example where profile.nls does not work with `plinear` but does with
# `default`
require(graphics)
set.seed(1)
DNase1 <- subset(DNase, Run == 1)
x <- rnorm(nrow(DNase1))
f1 <- nls(density ~ b1/(1 +
2012 Jan 20
1
nobs() and logLik()
Dear all,
I am studying a bit the various support functions that exist for
extracting information from fitted model objects.
From the help files it is not completely clear to me whether the number
returned by nobs() should be the same as the "nobs" attribute of the
object returned by logLik().
If so, then there is a slight inconsistency in the methods for 'nls'
objects with
2008 Mar 25
1
Error propagation
Dear R-helpers,
I´m in the context of writing a general function for error propagation
in R.
There are somehow a few questions I would like to ask (discuss), as my
statistical knowledge is somewhat restricted.
Below is the function I wrote, the questions are marked.
Many thanks in advance.
propagate <- function(expr, varList, type = c("stat", "raw"), cov =
TRUE)
{
2006 Apr 18
1
Nonlinear Regression model: Diagnostics
Hi,
I am trying to run the following nonlinear regression model.
> nreg <- nls(y ~ exp(-b*x), data = mydf, start = list(b = 0), alg = "default", trace = TRUE)
OUTPUT:
24619327 : 0
24593178 : 0.0001166910
24555219 : 0.0005019005
24521810 : 0.001341571
24500774 : 0.002705402
24490713 : 0.004401078
24486658 : 0.00607728
24485115 : 0.007484372
2006 Jul 18
4
How can I extract information from list which class is nls
Hello!
I work with :
R : Copyright 2006, The R Foundation for
Statistical Computing
Version 2.3.1 (2006-06-01)
On Windows XP Professional (Version 2002) SP2.
At this moment I use the function "nls" combined
with a selfStar model (SSmicmen, related to
Michaelis-Menten equation, and provided by the
"stats" package).
When I realise the following operation (cf. p 59
of the
2008 Sep 10
3
writing simple function through script
Hi all,
I try to write a simple function in a script. The script is as follows
yo<-function(Xdata)
{
n<-length(Xdata[,1])
Lgm<-nls(formula=LgmFormula,
data=Xdata,
start=list(a=1500,b=0.1),weights=Xdata$Qe)
return(Lgm)
}
After the execution of the script, when I call the function yo on data
called NC60.DATA I get an error.
#yo(NC60.DATA)
Erreur dans eval(expr, envir, enclos)
2010 Jun 10
1
nls model fitting errors
What am I failing to understand here?
The script below works fine if the dataset being used is
DNase1 <- DNase[ DNase$Run == 1, ] per the example given in
help(nlrob).
Obviously, I am trying to understand how to use nls and nlrob to fit
curves to data using R.
#package=DAAG
attach(codling)
plot(pobs~dose)
#next command returns 'step factor reduced below min factor
2000 Feb 14
2
Error in the inverse of a diagonal matrix?
I?m new to R so maybe this issue has been asked before and I still could not read the complete set of past messages sent to the list.
I found a weird behabiour that I will explain with a simple example. Lets consider the following block of commands:
> x <- diag(c(1,4,10))
> x
[,1] [,2] [,3]
[1,] 1 0 0
[2,] 0 4 0
[3,] 0 0 10
> invx <- x^-1
> invx
2006 Apr 23
1
lme: null deviance, deviance due to the random effects, residual deviance
A maybe trivial and stupid question:
In the case of a lm or glm fit, it is quite informative (to me) to have
a look to the null deviance and the residual deviance of a model. This
is generally provided in the print method or the summary, eg:
Null Deviance: 658.8
Residual Deviance: 507.3
and (a bit simpled minded) I like to think that the proportion of
deviance 'explained' by the
2011 Nov 10
1
Sum of the deviance explained by each term in a gam model does not equal to the deviance explained by the full model.
Dear R users,
I read your methods of extracting the variance explained by each
predictor in different places. My question is: using the method you
suggested, the sum of the deviance explained by all terms is not equal to
the deviance explained by the full model. Could you tell me what caused
such problem?
> set.seed(0)
> n<-400
> x1 <- runif(n, 0, 1)
> ## to see problem
2010 Jul 30
2
svydesign syntax and deviance!
Un texte encapsul? et encod? dans un jeu de caract?res inconnu a ?t? nettoy?...
Nom : non disponible
URL : <https://stat.ethz.ch/pipermail/r-help/attachments/20100731/ac3b9e43/attachment.pl>
2012 Jan 10
0
rpart vs. tree and deviance calculations
Hi Everyone,
I'm working on building some classification trees, and up to this point
I've been using rpart. However, I recently discovered the tree package,
and found that it had some useful functions (in particular deviance(),
which I would really like to use for my project). I can't seem to find
an equivalent function for rpart. I've considered using tree() in place
of