similar to: Fitting Distributions Directly From a Histogram

We obtained some disturbing results from convolve() (inaccuracies and negative probabilities). We'll try to make the context clear in as few lines as possible... Our function panjer() (code below) basically computes recursively the probability mass function of a compound Poisson distribution. When the Poisson parameter lambda is very large, the starting value of the recursive scheme ---

Hi just wondering if there is a package that can get the maximum likelihood or method of moments estimator for distributions with censored data? The distributions I'm interested in are: Exponential, pareto, beta, gamma and lognormal. -- View this message in context: http://www.nabble.com/MLE-for-censored-distributions-in-R-tp15022863p15022863.html Sent from the R help mailing list archive at

Given a numeric vector of observations, does R have any generic way to estimate the parameters of commonly used distributions (exponential, gamma, etc.) without numerically optimizing the likelihood function? Thanks, David _______________________________________ David R. Bickel http://davidbickel.com Research Scientist Pioneer Hi-Bred International Bioinformatics & Exploratory Research 7250

The hist() command produces this message on my machine... Error in title(main = main, sub = sub, xlab = xlab, ylab = ylab, ...) : X11 font at size 14 could not be loaded How can I either (a) determine what font is required, or (b) specify one of the fonts I have available? This problem is specific to hist(), plot() works fine. I am using R on SuSE Linux v9.3, from the KDE desktop.

I have grouped data in this format Size -- Count 0-10 -- 15 10-20 -- 25 20-50 -- 10 50-100 -- 5 I've been trying to find a way to set this up with the proper histogram heights, but can't seem to figure it out. So any help would be much appreciated! -- View this message in context: http://www.nabble.com/Histogram-for-grouped-data-in-R-tp21624806p21624806.html Sent from the R help

Is there any neat way to add a curve (frequency function or the like) to a histogram or other plots? I haven't found one yet... Robert

Hi, how can I treat data organised in classes and frequencies? Ex. class frequency 20-23 9 23-25 7 26-28 5 29-31 5 32-34 3 Thanks Angelo Scozzarella

Dear All, I am trying to fit some data both as a negative binomial and a zero inflated binomial. For the first case, I have no particular problems, see the small snippet below library(MASS) #a basic R library set.seed(123) #to have reproducible results x4 <- rnegbin(500, mu = 5, theta = 4) #Now fit and check that we get the right parameters fd <- fitdistr(x4, "Negative

Is there a canonical way to check if all elements of a vector or matrix are the same? Solutions below work, but look hackish to me. > x <- rep(1, 10) > all(x == x[1]) # == operator does not provide for small differences [1] TRUE > isTRUE(all.equal(x, rep(x[1], length(x)))) # ugly [1] TRUE Best, Vincent -- Vincent Goulet, Associate Professor ??cole d'actuariat

UseRs, actuar is a package for Actuarial Science. A rather preliminary version (0.1-3) of the package has been available on CRAN since February 2006. We now announce the immediate availability of version 0.9-2 sporting a large number of new features. Non actuaries behold! There can be some features of interest for you, especially those related to new probability distribution and to the

UseRs, actuar is a package for Actuarial Science. A rather preliminary version (0.1-3) of the package has been available on CRAN since February 2006. We now announce the immediate availability of version 0.9-2 sporting a large number of new features. Non actuaries behold! There can be some features of interest for you, especially those related to new probability distribution and to the

Dear R developers, I'm a former APL programmer. In that language, the "take" (up arrow) and "drop" (down arrow) operators were extensively used to, well, take and drop elements of vectors. Functions head() and tail() are equivalents in R for the "take" operator, but nothing seems to mimic the "drop" operator. I think it would be useful. For

[Version fran??aise plus bas] To the R community, A quick word to announce the publication of my document "Introduction ?? la programmation en S". It is available in the French section of the Contributed Documentation page of CRAN. Many of the documents or books currently available on S-Plus and/or R present the software in a statistical analysis context. My document rather focuses

Dear All, I need to fit a custom probability density (based on the symmetric beta distribution B(shape, shape), where the two parameters shape1 and shape2 are identical) to my data. The trouble is that I experience some problems also when dealing with the plain vanilla symmetric beta distribution. Please consider the code at the end of the email. In the code, dbeta1 is the density of the beta

Dear R People: I have finally seen the error of my ways and have decided to use ESS for R and S + stuff. However, I have a question right from the beginning. I'm somewhat confused by the installation instructions. Do I install XEMACS or ESS first, please? Windows XP R Version 2.0.1 (S + 6.2) Thanks so much! Sincerely, Laura Holt mailto: lauraholt_983 at hotmail.com

I created a type of object similar to a data frame. In some circumstances, It needs special methods for "[" and "[<-" and rbind() (but not cbind()). Then I found this in the cbind()/rbind() man page: The method dispatching is _not_ done via 'UseMethod()', but by C-internal dispatching. Therefore, there is no need for, e.g., 'rbind.default'.

Hy all, Is there a direct way to build the complete function call of an arbitrary function? Here's what I want to do. A function will build a function which will itself call a probability density function for some law given in argument to the first function: > f("gamma", 1000) will return, say, function(x, shape, rate, scale = 1/rate) dgamma(x + 1000, shape, rate,

Full_Name: marco zamboni Version: no R version OS: ubuntu gusty Submission from: (NULL) (87.9.174.188) I have just finish to install gusty ubuntu on my AMD sempron. I would like to: $ sudo apt-get install r-base but ... I seguenti pacchetti hanno dipendenze non soddisfatte: r-base: Dipende: r-base-core (>= 2.6.1-1gutsy0) ma 2.5.1-1 sta per essere installato Dipende:

Full_Name: marco zamboni Version: no R version OS: ubuntu gusty Submission from: (NULL) (87.9.174.188) I have just finish to install gusty ubuntu on my AMD sempron. I would like to: $ sudo apt-get install r-base but ... I seguenti pacchetti hanno dipendenze non soddisfatte: r-base: Dipende: r-base-core (>= 2.6.1-1gutsy0) ma 2.5.1-1 sta per essere installato Dipende:

Please guide me through to resolve the error message that I get this is what i have done. >x1<- rnorm(100,2,1) >x1fitbeta<-fitdistr(x1,"beta") Error in fitdistr(x1, "beta") : 'start' must be a named list Yes, I do understand that sometime for the distribution to converge to the given set of data, it requires initial parameters of the distribution, to