Displaying 20 results from an estimated 2000 matches similar to: "Non-linear Regression : Error in eval(expr, envir, enclos)"
2012 May 15
2
Renaming names in R matrix
I have the following matrix:
> dat
[,1] [,2] [,3] [,4]
foo 0.7574657 0.2104075 0.02922241 0.002705617
foo 0.0000000 0.0000000 0.00000000 0.000000000
foo 0.0000000 0.0000000 0.00000000 0.000000000
foo 0.0000000 0.0000000 0.00000000 0.000000000
foo 0.0000000 0.0000000 0.00000000 0.000000000
foo 0.0000000 0.0000000 0.00000000 0.000000000
and given this:
2011 Jul 27
2
fitting sine wave
Dear R-helpers
?
I have 7 data points that I want to fit a continuous curve to, that should look similar to a sine wave
My data points would mark the local minima and maxima respectively.
This is what I?ve got so far. And I would keep doing so, but sadly nls() then says that it has reached the maximum number of Iterations?
?
2009 May 06
1
'RG' looks like a pre-2.4.0 S4 object: please recreate it
I would like to load ApoAI.RData. During the operation of reading this
data an error occurs. There is also a problem with STF file.
> library (limma)
> load("ApoAI.RData")
Warning message:
'RG' looks like a pre-2.4.0 S4 object: please recreate it
> objects()
[1] "RG"
> names(RG)
[1] "R" "G" "Rb"
1997 Apr 30
2
R-alpha: New Incomplete Beta Function
Here is a drop-in replacement for the R incomplete beta function.
src/math/pbeta.c
It is a slightly modified version of the cephes library one from
Netlib. In the few cases I tried it seems to give at least 14
digit agreement with the one in S-PLUS (its hard to get more).
I'm not sure what performance is like. I'd like to know if it
helps with some of the problems which have been
2008 Feb 01
2
the "union" of several data frame rows
Hi,
I have a question about how to obtain the union of several data frame
rows. I'm trying to create a common key for several tests composed of
different items. Here is a small scale version of the problem. These
are keys for 4 different tests, not all mutually exclusive:
id q1 q2 q3 q4 q5 q6
1 A C
2 B D
3 A D B
4 C D B D
I would like
2008 Mar 06
2
Help with parsing a data file
Hi All,
I need to parse data from a file, example shown below. The first two lines
can be skipped, the third line contains the column names. The next 13 lines
can be skipped. The next line "1991" is a year value, with the following 13
values data for that year. The file then repeats this format with (year, 13
lines of data for that year). I would ideally like to end up with an
2004 Oct 04
3
Beginners problem
Hi,
I'm new to R and have a problem with a little test program (see below).
Why doesn't <<- in function rk4
assign the new value to y so that it is seen in rktest. I thought that
<<- does exactly this. But it seems that I
didn't get it right. I would be very appreciative for an explanation of
that behaviour of <<-. I know how to
write the whole thing so that it
2009 Oct 19
1
updating columns using other column as reference
Dear R-gurus,
Just supose I have a dara.frame that looks like
myDF<-read.table(stdin(),head=T,sep=",")
codID,namesp,k1,k2,k3,k4
1,spA,2,5,6,3
2,spB,4,5,4,6
3,spC,2,1,5,6
4,spC,5,4,3,2
5,spD,1,2,3,4
6,spE,2,4,3,1
I need to update the columns k1-k4 with the namesp, but
considering the math between Kx and codID.
My desired output must looks like:
codID,namesp,k1,k2,k3,k4
2006 Nov 29
2
How to solve differential equations with a delay (time lag)?
Hi,
I would like to solve a system of coupled ordinary differential equations,
where there is a delay (time lag) term. I would like to use the "lsoda"
function "odesolve" package. However, I am not sure how to specify the
delay term using the syntax allowed by odesolve.
Here is an example of the kind of problem that I am trying to solve:
> library(odesolve)
2010 Nov 21
1
solve nonlinear equation using BBsolve
Hi r-users,
I would like to solve system of nonlinear equation using BBsolve function and
below is my code. I have 4 parameters and I have 4 eqns.
mgf_gammasum <- function(p)
{
t <- rep(NA, length(p))
mn <- 142.36
vr <- 9335.69
sk <- 0.8139635
kur <- 3.252591
rh <- 0.896
# cumulants
k1 <- p[1]*(p[2]+p[3])
k2 <- p[1]*(2*p[2]*p[3]*p[4] +p[2]^2+p[3]^2)
k3 <-
2001 Jun 06
0
snk authentication
Here is a little patch against 2.9p1 that performs the SNK (also known
as TIS authserv) challenge-response automaticly instead of asking the
user. hope you find it useful.
--larry
-------------- next part --------------
diff -NuBw openssh-2.9p1/Makefile.in openssh/Makefile.in
--- openssh-2.9p1/Makefile.in Thu Apr 26 20:31:08 2001
+++ openssh/Makefile.in Wed Jun 6 16:15:56 2001
@@ -43,9 +43,9
2005 Feb 05
0
Rép : 2 small problems: integer division and the nature of NA
Thanks to the many R users who convinced me that the sum of NAs should
be zero and gave me a solution if I did not want it to be zero.
Thank you also for the explanations of rounding errors with floating
point arithmetics. I did not expect it. This small error was a real
problem for me as I was trying to find a way to recode numeric values
into intervals. Because I wanted to retain numeric
2013 Apr 25
2
Vectorized code for generating the Kac (Clement) matrix
Hi,
I am generating large Kac matrices (also known as Clement matrix). This a tridiagonal matrix. I was wondering whether there is a vectorized solution that avoids the `for' loops to the following code:
n <- 1000
Kacmat <- matrix(0, n+1, n+1)
for (i in 1:n) Kacmat[i, i+1] <- n - i + 1
for (i in 2:(n+1)) Kacmat[i, i-1] <- i-1
The above code is fast, but I am curious about
2008 Mar 28
1
Device problems in a loop.
I have tried the following code to plot and save many box plots in files. The code works by itself, but does not run in a loop. The error indicates that there is no device. The error is printed below the code.
Code:
for(i in 1:len)
{
A = as.numeric(delta[Delta[i,1]])
B = as.numeric(delta[Delta[i,2]])
C = as.numeric(delta[Delta[i,3]])
D = as.numeric(delta[Delta[i,4]])
data
2005 May 30
3
sapply following using by with a list of factors
Background:
OS: Linux Mandrake 10.1
release: R 2.0.0
editor: GNU Emacs 21.3.2
front-end: ESS 5.2.3
---------------------------------
Colleagues
I am having some trouble extracting results from the function by, used to
average variables in a data.frame first by one factor (depth) and then by a
second factor (station). The real data.frame is quite large
> dim(data.2001)
[1] 32049 11
Here is a
2012 Jun 14
0
fixed trimmed mean for j-group
Hello...i want to find the empirical rate for type 1 error using fixed
trimmed mean. To make it easy, i'm referring to journal given by this
website
http://www.academicjournals.org/ajmcsr/PDF/pdf2011/Yusof%20et%20al.pdf.
I already run the programme and there is no error in it but i got zero for
the empirical rate of type 1 error. The empirical rate for the type 1 error
given in the journal
2011 Nov 15
2
How can rpm "%{SUMMARY}" not be consistent?
I have been seeing something for quite some time which has confused me
considerably for over a year, perhaps one of you can help me understand.
Assumed: rpm queries are against _a_ database.
Assumed: database queries against the same database, without changes to
the data in the database, will return the same data.
Confusion: then why are some of the summaries reported by rpm different?
Each day
2012 Jul 07
0
fixed trimmed mean for group
Hello,
I haven't found errors in your code. I implemented the test in the paper
(the first, fixed symetric mean) and it also gives me zero Type I
errors, when alpha = 0.05. Try to see the value of min(pv) or to plot
the histogram of 'pv', hist(pv) and you'll see that there are no
significant p-values, at that level.
Anyway I'll continue to look at it, but my first
2000 Dec 23
1
Look what I found under the Xmas tree!
Hello people,
Looks like Santa Claus thinks I've been a good boy this year.
Here's the third in my performance patch series: d.m.l
Apply after applying d.o.n and d.n.m; I don't know how much of
those got applied to the CVS tree.
What's inside:
Request for help! Look in os.h if you're using a compiler or
processor I don't use (I use gcc on K5, K7, G3).
New MDCT! Now we
2001 Sep 11
2
Differential Equations Using R?
To whom it may concern,
I am a student at Macaleste College, and next semester Macalester
is going to offer a course for CellBio that is mainly statistically based.
For the most part the students will be using R for analysis. The problem is
there will be some simple differential equations for the students to solve.
The committee that in charge of the classes corriculam would like only to