Displaying 20 results from an estimated 1000 matches similar to: "nls making R "not responding""
2011 Feb 28
3
nls not solving
I am running the following nls equation. I tried it with data that excel was
fitting and got the error:
singular gradient matrix at initial parameter estimates
I thought it was due to a low number of points (6), but when I create a
dataset, I get the same problem. If I remove the parameter "a," then it can
find a solution. Does anyone know what I can do to fit this model?
2010 Sep 02
1
How using the weights argument in nls2?
Good morning gentlemen!
How using a weighted model in nls2? Values with the nls are logical since
values with nls2 are not. I believe that this discrepancy is due to I did
not include the weights argument in nls2.
Here's an example:
MOISTURE <- c(28.41640, 28.47340, 29.05821, 28.52201, 30.92055,
31.07901, 31.35840, 31.69617, 32.07168, 31.87296, 31.35525, 32.66118,
33.23385,
2005 Dec 21
0
Help with Krige.conv using linear models
A majority of my data makes a kriged map perfectly using an exponential
model for the semivariogram to fit my data and then going through the
commands variofit() to define the model and then krige.conv() to use the
model to predict values in a grid. But?one set of my data appears to be
linearly correlated for the first 5000 meters and not correlated beyond
that. I have been having problems
2010 May 11
1
nls() and nls2() behavior?
first, apologies for so many posts yesterday and today. I am
wrestling with nls() and nls2(). I have tried to whittle it down to a
simple example that still has my problem, yet can be cut-and-pasted
into R. here it is:
library(nls2)
options(digits=12);
y= c(0.4334,0.3200,0.5848,0.6214,0.3890,0.5233,0.4753,0.2104,0.3240,0.2827,0.3847,0.5571,0.5432,0.1326,0.3481)
x=
2006 Mar 27
0
not a problem: a submission (PR#8714)
Full_Name: TOBY MARTHEWS
Version: 2.1.1
OS: Windows XP
Submission from: (NULL) (139.133.7.38)
I think you should have better examples on the ?plot man page, if you don't mind
me saying. Can I suggest something like this, which would probably stop so many
emails to the R help about how to put on error bars
Cheers,
Toby M
******************************
xvalsT1=1:3
yvalsT1=c(10,20,30)
2010 Nov 24
1
The nls2 function automatically prints the object!
Good morning gentlemen!
When I use the function nls2, and store it in an object, that object is
automatically printed, without the summary or to draw the object. For
example.
model <- nls2 (...)
Number of iterations to convergence: ...
Achieved convergence tolerance: ...
Nonlinear regression model
model: ... ~ ...
Date: NULL
The B k
... ... ...
residual sum-of-squares: ...
Number
2005 Jul 28
2
problem with an IF statement?
Can somebody please take a look at this and tell me whats going wrong?
It seems to be parsing wronly around the 'if' statement and gives me a
directory listing.
Thanks in advance
Tom
N.B. datan is an invented dataset
xvals<-c(1,0.4,0.2)
datan<-data.frame(s1=c(3,4,5),s2=c(5,5,5),s3=c(21,55,34),s4=c(5,3,2))
datan$sint<-NA
datan$sgrad<-NA
for(icount in 1:dim(datan)[1]) {
2013 Jan 02
1
Need help with self-defined function to perform nonlinear regression and get prediction interval
Dear All,
I was trying to call a self-defined function that performs nonlinear
regression and gets the corresponding prediction upper limit using nls2
package. However, weird thing happened. When I called the function in the
main program, an error message "fitted(nlsmodel): object 'nlsmodel' not
found" came up. But when I directly ran the codes inside the function, no
error came
2006 Sep 11
3
Extracting overdispersion estimates from lmer amd glm objects
Dear list,
I am needing to extract the estimate of overdispersion (deviance / residual degrees of freedom or c-hat) from multiple model objects - so they can then be used to compare the extent of overdispersion among alternative models as well as calculate qausi-AIC values. I have been unable to do this, despite consulting a number of manuals and searching the R-help. I am imaging that in
2007 Aug 01
1
Problem to remove loops in a routine
Dear R-users,
I have written the following code to generate some trellis plots. It
works perfectly fine except that it is quite slow when it is apply to my
typical datasets (over several thousands of lines). I believe the
problem comes from the loops I am using to subset my data.frame. I read
in the archives that the tapply function is often more efficient than a
loop in R. Unfortunately ,
2011 May 18
1
How to make array of regression objects
Dear all,
I have made couple logistic regressions, what making a distribution of some
event.
Currently, i store it like this:
o1 <- lrm(...)
o2 <- lrm(...)
o3 <- lrm(...)
...
Then, i have made a function to peak required regression object from this
variables by it number:
get_object <- function(obj_name, nModel) {
eval (parse(text=paste("o <- ", obj_name, nModel,
2007 Jun 15
0
No subject
using Asterisk.
=20
Is this all you want Asterisk to do? (eg as an application service
rather than provide telephony for the rest of the library as well), or
are you looking to have it replace your existing telephony equipment?
=20
As a suggestion if you google Trixbox and Nerd Vittles you will find a
fairly detailed explanation of how to set your Trixbox server (a version
of Asterisk) up to
2012 Jan 27
0
[LLVMdev] Double spills with Greedy regalloc
On Jan 27, 2012, at 10:20 AM, Borja Ferrer wrote:
> Hello,
>
> I noticed the following interesting code sequence while compiling a piece of code with the backend I'm developing. Probably this issue is for Jakob, but anyways this is what I'm getting:
>
> STDWPtrQRr <fi#12>, 0, %R25R24; mem:ST2[FixedStack12](align=1)
> STDWPtrQRr <fi#12>, 0, %R25R24;
2012 Jan 27
2
[LLVMdev] Double spills with Greedy regalloc
Hello,
I noticed the following interesting code sequence while compiling a piece
of code with the backend I'm developing. Probably this issue is for Jakob,
but anyways this is what I'm getting:
STDWPtrQRr <fi#12>, 0, %R25R24; mem:ST2[FixedStack12](align=1)
STDWPtrQRr <fi#12>, 0, %R25R24; mem:ST2[FixedStack12](align=1)
STDWPtrQRr <fi#13>, 0, %R23R22;
2003 Nov 16
2
prevent conversion to factors in aggregate?
I've been trying to figure out how to prevent a column that is the result of
an aggregate function call so that I can use it in further calculations. For
example, I would like to aggregate the expf for the data.frame by sp
(character) and dbh (double d=rounded to integer) using the command:
> st2 <- aggregate( ntrs$expf, by=list(sp=ntrs$sp,dbh=ntrs$dbh), sum )
> st2$expf <- st2$x
2006 Oct 21
0
tcltk: multiple listboxes, selection
Dear list,
I have multiple (BWidget) listboxes in the same toplevel window.
The problem is, if I select (by left clicking) on one of those
listbox elements, the current selection in the *other* listboxes is
cleared!
Anybody knows how I can prevent this?
Here's my code (sorry not complete):
(E.g. If I select an X value, I'd lose the Y value I selected before)
gui.create.tab.general <-
2011 Oct 15
0
shoulda installation troubles
I cannot manage to have shoulda working.
I''m using rvm with rails 3.1 and ruby 1.9.2
I added
group :test do
gem "shoulda"
end
to the Gemfile, and then run bundle install, thus the gem is installed,
and it is featured in the list running gem list --local.
but whenever I add something like
class ModelTest < Test::Unit::TestCase
should_not
2000 Dec 15
1
resolution of windows metafiles
I am having trouble getting smooth-appearing curves in figures produced as
windows metafiles using R (R-1.1.1 on Windows 98). When I import them into
a word processor (either MS Word 97, or WordPerfect 9.0), the figures
appear only slightly bumpy on screen, but when printed, there is a clear
jaggedness reminiscent of aliasing. Is there a fix for this (if the answer
is "upgrade to 1.2",
2004 Aug 20
0
Proposed (minor) change to plot.TukeyHSD
Attached follows a patch to a minor change in the plot method of
the TukeyHSD class (package stats). Basically it defines main= and xlab=
as formal arguments of the plot function, with reasonable default
values, passing them to title(), instead of using hard-coded only values
for main= and xlab=. This way it's possible to change the title and xlab
of the plots created using plot.TukeyHSD().
2020 Jun 19
2
Aliasing and forwarding optimization
----Snip--
struct st1{
int a;
};
struct st2{
int b;
};
struct st {
struct st1 obj1;
struct st2 obj2;
}Obj;
int test1(struct st1 * ptr1 , struct st2 * ptr2, struct st2 *ptr3) {
ptr1->a = 10;
*ptr3 = *ptr2;
return ptr1->a;
}
--Snip---
For the above case GCC is able to store forward the value 10 to the return
place.
LLVM is not doing this.
GCC
https://godbolt.org/z/FCjCXy
LLVM