similar to: Error "singular gradient matrix at initial parameter estimates" in nls

Hi R mailing list, probably a very basic problem here, I try to do the following: > Q<-c(1,2,3) > P<-c(4,5,6) > A<- data.frame(Q,P) > A Q P 1 1 4 2 2 5 3 3 6 this is my simplified data.frame (matrix) now I try to create following loop for subtraction of element within the data.frame: > for(i in length(A[,"P"]-1){ delta[i]<-

Dear mailing list, I have following array: X2 Y2 [1,] 422.7900 6.0 [2,] 469.8007 10.5 [3,] 483.9428 11.0 [4,] 532.4917 25.5 [5,] 596.1942 33.5 [6,] 630.8496 40.5 [7,] 733.2996 45.0 [8,] 946.4779 32.0 [9,] 996.8068 35.5 [10,] 1074.3310 23.0 I do afterwards the following: plot.new() plot.window(xlim=c(min(X1)-50,max(X1)+50),

I have a data set which contains two columns. The first column is a list of countries, and the second column contains their political risk ratings. I would like to create one large plot that contains 5 different sections, each with a scatter plot. To clarify, I have divided the countries into 5 groups. For each group (continent), I would like to have the name of the continent on the x-axis, and

Hi R-mailing list I would have the following set-up below with a simplified data-frame. Through a loop which includes certain criteria for the densities I would like to plot the different density-distributions in the same plot. Of course I hope I don't do any mistakes with all the indexes of the dataframe. All I would like to have is the different densities in the same plot with a general

Hi We have a Replica 2 + Arbiter Gluster setup with 3 nodes Server1, Server2 and Server3 where Server3 is the Arbiter node. There are several Gluster volumes ontop of that setup. They all look a bit like this: gluster volume info gv-tier1-vm-01 [...] Number of Bricks: 1 x (2 + 1) = 3 [...] Bricks: Brick1: Server1:/var/data/lv-vm-01 Brick2: Server2:/var/data/lv-vm-01 Brick3:

What's the logic behind the following, and where can I find any documentation about it? In particular, why are 2:9 - as characters - not regarded as being smaller than 10? # R-Code: a <- as.character(1:12) a < 10 # [1] TRUE FALSE FALSE FALSE FALSE FALSE FALSE FALSE FALSE FALSE FALSE FALSE Thanks in advance! Niklaus

Hi all, I ran into a problem in some of my code that could be traced back to 'mean' sometimes returning NA and sometimes NaN, depending on the value of na.rm: > mean(c()) [1] NA > mean(c(NA),na.rm=T) [1] NaN However, I don't understand the reasoning behind this and would appreciate and explanation. I understand that the mean of an empty vector is not definied, but I

Hi, well after some ever further investigation of some bugs, I ended up cleaning up the code, again. I also fixed some really heavy memory leaks in csvToList and made it in general more stable. Regards, Patrick "Marex" Niklaus -------------- next part -------------- A non-text attachment was scrubbed... Name: ini.c.patch Type: text/x-patch Size: 8275 bytes Desc: not available Url :

Dear all, I am not sure if this mail is for R-help or should be sent to R-devel instead, and therefore post to both. While using nlrob from package 'robustbase', I ran into the following problem: For psi-functions that can become zero (e.g. psi.bisquare), weights in the internal call to nls can become zero. Example: d <- data.frame(x=1:5,y=c(2,3,5,10,9)) d.nlrob <-

Hi all, I am trying to produce some panels with dots in an X/Y plane where the diameter of the dots indicates a Z value (like e.g. earthquake maps where dot sizes indicate magnitudes and X/Y the location). This works fine with xyplot, e.g.: xyplot(1:3~1:3,cex=1:3,pch=16) However, when I do this with a panel variable, e.g.: x<-rep(1:3,5) y <- rep(1:3,5) sz <- rep(1:5,each=3) grp

Hello! A couple of days ago I submitted the problem report shown below to the support forum for WinSCP. I got a reply (shown at the end of this e-mail) saying that this in part was a WinSCP problem, but also that there appeared to be something wrong with the replies from OpenSSH-3.9p1 under AIX 4.3.3. The full dialog including the reply can be seen at

Hi all, Is there a way to modify the optional arguments (...) passed to a function, so that these can be passed in modified form to a subsequent function call? I checked "Programming with Data" but could not find a solution there. What I'd like is something along these lines: test <- function(x,y,...) { if(!hasArg(xlab)) { ___add xlab to ...___ } if(hasArg(xlab)) {

Hi, For a n-level factor, I'd like to specify the first contrast and have the remaining n-2 constructed automatically so that the set is orthogonal. I then test the contrasts with summary.lm(anova-object). In S-Plus, the following works: >y.anova <- aov( y ~ C(CO2,c(1,0,-1)) ) >summary.lm(y.anova) In R, it fails with the following error: >levels(CO2) [1]

Dear all, I hope this is the right mailing list for my question -- I felt that this was too technical for R-help. I am trying to compile R-2.4.0 on a HP-UX system: ./configure MAKE=gmake --prefix=$HOME --without-x --without-tcltk --disable-R-profiling --without-readline --disable-multibyte R is now configured for ia64-hp-hpux11.23 Source directory: . Installation directory:

Our previous attempts to select a project name for CompComm have not gone well. We tried a poll but the names in the poll were not pre-qualified to determine if they were usable and the ?winning? name was also a strongly disliked name. Kristian and I each proposed a plan to select a name, but that effort got sidetracked by a massive flame war on the mailing list. While our community has some

pypy-0.7.0: first PyPy-generated Python Implementations ============================================================== What was once just an idea between a few people discussing on some nested mailing list thread and in a pub became reality ... the PyPy development team is happy to announce its first public release of a fully translatable self contained Python implementation. The 0.7 release

A new compiz release 0.7.2 is now available from: http://xorg.freedesktop.org/archive/individual/app/compiz-0.7.2.tar.gz which can be verified with: http://xorg.freedesktop.org/archive/individual/app/compiz-0.7.2.tar.gz.sha1 da3dbe8c5be19aa03529e8d3c7e32365c82a973b compiz-0.7.2.tar.gz

Hi, recently I took a look at ini to find a bug which made it crash at startup. While fixing the bug I realized I could do some improvements to the option reading code (espacally the action part). So I went on and here we are. This patch should make ini more robust, clean it up and fix also some more crashes I had here. Additionaly, I have a attached a second patch from Maniac which should fix

Hi all, Is there a function to check if a particular value is contained in a vector? I've looked at grep in the hope that I could use a Perl-like syntax, but obviously it's different... I'd like to do something like: y <- c("a","b","c") if("a" in y) { # "a" is not in y } Also, is there a way to

Hi all, Given a number of points (x,y) in a plane, I'd like to plot a map of polygons, so that 1) each polygon contains exactly one point 2) the polygon defines the area for which this specific point is closer than any other point. It's a bit like a map of areas "influenced" by that point, and it's obviously a matter of intersecting the perpendicular bisectors