similar to: numericDeriv and ecdf

> On only ten points, what did you expect ? Even with 1000 > observations, estimating a density is difficult, and has > been the subject of a century of research. Kernel density > estimates are among the most successful. For your immediate > application, try plot(density(rnorm(10)), type="l"), etc. wait, you misunderstood me! I'd like to see 10 or 9 points with

Hi all, is there any package in R capable of smooth approximation of CDF basing on given sample? (Thus, I am not speaking about ecdf) In particular, I expect very much that the approximation should subject to the property: f(x0)<=f(x1) for x0<x1, where x0 and x1 belong to range of the sample given. Polynomial approximation could be OK for me as well. P.S.

Hi all, once more about the old subj :-) My data has too much various distribution families and for every particular experiment I need just to decide whether the data is "quite homogeneous" or it has two or more clusters. I've revisited the following libraries: amap, clust, cclust, mclust, multiv, normix, survey. And I didn't find any ready-to-use general

Dear Prof. Brian Ripley, first of all thank you for your answer, I do appreciate how do you manage to keep successfully all your activities and answer posts in this forum! > An empirical CDF is a step function: it does not have a > derivative at the jump points, and has a zero > derivative everywhere else. of course! Let me add few words concerning my simple motivation. 1.

Hi All, We'd like to use functions provided in R in our application. Our application is written in C/C++ and currently runs on win32, Linux and Mac. We'd be happy to attach the whole R ( i.e. not just transfer some function by hand). It is important that we deal with big amount of data, so "command line"-like invocations won't be very interesting. We'd

Hi All, We'd like to use functions provided in R in our application. Our application is written in C/C++ and currently runs on win32, Linux and Mac. We'd be happy to attach the whole R ( i.e. not just transfer some function by hand). It is important that we deal with big amount of data, so "command line"-like invocations won't be very interesting. We'd

> Almost any method of fitting a density estimate would work on > integrating (numerically) the result. it is a nice idea concerning the monotony property, which will be obtained automatically, but I am going to use results of approximation analytically > In particular, look at package polspline, where > p(old)logspline does the integration for you. thank you, I am going to

Hi All is there "homals" package prepared for win32? kind regards, Valery A.Khamenya --------------------------------------------------------------------------- Bioinformatics Department BioVisioN AG, Hannover

Hi All, while looking why the cclust(cclust) doesn't work for 1-dimensional data, I've found unpleasant behavior in semantics of R. Indeed: is.matrix(matrix(cbind(c(1,2,3,4)),ncol=2)[1:2,]) == TRUE but: is.matrix(matrix(c(1,2))[1:2,]) == FALSE kind regards, Valery A.Khamenya --------------------------------------------------------------------------- Bioinformatics

> Well, that is in all good texts on R, together with the > solution: drop=FALSE. See ?"[" for the on-line details. OK. Thank you a lot. Now patched cclust and clustIndex work fine for 1D case. BTW, why not to apply the "drop=F" to these functions? I guess other users need 1D case as well. kind regards, Valery A.Khamenya

Hi devels, i did not find at this page: http://developer.r-project.org/ideas.txt any ideas concerning incorporating documentation possibilities (say, Sweave-based) into R-scripts. Was it discussed already? (If discussed, then what is the decision/conclusion then?) thanks, Valery

Hi all, Situation: there is a function `f' already defined by someone and provided in package. `f' looks like that: f <- function() { x+1 } i.e. `f' is not closed i.r.t. term `x' now I have my own function `g', where I'd like to override variable `x' while calling `f': x <- "dummy gloabal value" g <- function() { x

Hi all, I took a look at last 2 months post in R-help maillist and surfed through the R-project.org . Unfortunately, I can't find some page with roadmap/statements about major changes coming in R-2.0 in comparison to R-1.9 Could anyone point me to the right URL? Thank you in advance. -- Valery

Hi all, I took a look at last 2 months post in R-help maillist and surfed through the R-project.org . Unfortunately, I can't find some page with roadmap/statements about major changes coming in R-2.0 in comparison to R-1.9 Could anyone point me to the right URL? Thank you in advance. -- Valery

Dear R developers, I've run into a weird behavior of the numericDeriv function (from the stats package) which I also posted on StackOverflow (question has same title as this email, except for the version of R). Running the code bellow we can see that the numericDeriv function gives an error as the derivative of x^a wrt a is x^a * log(x) and log is not defined for negative numbers. However,

I have to compute some standard errors using the delta method and so have to use the command "numericDeriv" to get the desired gradient. Befor using it on my complicated function, I've done a try with a simple exemple : x <- 1:5 numericDeriv(quote(x^2),"x") and i get : [1] 1 8 27 64 125 216 attr(,"gradient") [,1] [,2] [,3] [,4] [,5] [,6] [1,] Inf

Hi all, What is the right portable way to save graphics in batch mode? Remarks: 1. Problem is STFWed and RTFMed. In particular a short note about png() is found in R-FAQ. In fact, there were stated that png() is not reliable under Linux in batch mode. 2. savePlot under windows is quite convenient, but not supplied under Linux. 3. pdf() + postscript() < savePlot()

Hi, I am stuck on something for a couple days, I am almost about to give up. This looks simple, but I can't figure out. I hope I can get some help here. I am trying to do some symbolic and numerical derivations. Let me explain the problem. Let's say, I have a matrix as follows: > load <- matrix(c(3,0,1,4,1,3),nrow=3,ncol=2,byrow=TRUE) > > load [,1] [,2] [1,] 3 0

R Help List -- I have defined two time-series-vector-valued-functions, let them be f and g, and want to find the numeric derivative of f with respect to the variable x where f depends on x through g: (d/dx)(f (g(x) ) Moreover, x is a vector I tried this out the long way (naming every element of the x vector and then making the 'theta' argument in numericDeriv() the character vector of

Dear all As far as I could trace, looking at the function C function numeric_deriv, this unwanted behavior comes from the inner most loop in, at the very end of the function, for(i = 0, start = 0; i < LENGTH(theta); i++) { for(j = 0; j < LENGTH(VECTOR_ELT(pars, i)); j++, start += LENGTH(ans)) { SEXP ans_del; double origPar, xx, delta; origPar = REAL(VECTOR_ELT(pars, i))[j];