similar to: MLE for censored distributions in R

Hi, I would like to use the following codes to plot the CDF for pareto distribution. Before doing this, I have plot the emperical one. x <- seq(1.6, 3, 0.1) lines(x,pgpd(x, 1.544,0.4477557,), col="red") Could anyone give me some advice whether the above codes are correct? Many thanks. -- View this message in context:

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

Hello: Could anyone point me to a nice example where someone has created methods for "[" on a user defined Class? I looked at the package Matrix but that was a little daunting. I'm looking for someone a little more introductory. I've tried to search the help section and the web but its difficult since "[" isn't searchable. Thanks in advance! Greg [[alternative

Dear developeRs, Compilation of the latest version (0.9-5) of my actuar package fails with r-release MacOS_X ix86 on CRAN; see http://www.R-project.org/nosvn/R.check/r-release-macosx-ix86/actuar-00check.html All errors come from accented letters in comments in latin-1 encoded files (except hierarc.R which is in UTF-8, my bad). Encoding is declared as latin-1 in DESCRIPTION. The package

Dear R Users, I am trying to obtain a best-fit analytic distribution for a dataset with 11535459 entries. The data range in value from 1 to 300000000. I use: fitdistr(data, "gamma") to obtain mle's for the parameters. I get the following error: Error in optim(start, mylogfn, x = x, hessian = TRUE, ...) : non-finite finite-difference value [1] And the following warnings:

Dear All, A simple question: packages like fitdistr should be ideal to analyze samples of data taken from a univariate distribution, but what if rather than the raw data of the observations you are given directly and only a histogram? I was thinking about generating artificially a set of data corresponding to the counts binned in the histogram, but this sounds too cumbersome. Another question is

Package grDevices has a function cm() to convert from inches to centimeters. Its definition is: cm <- function(x) 2.54*x As far as I can tell, the function is not used anywhere in the R sources (I grepped for "cm(", "cm)" and ", cm".) I did not check for all packages on CRAN, though. May I propose to remove this apparently useless function from the sources? [I

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

Hi R users: I will like to know if somebody works on insurance statistics (actuarial problems) and had use TRICAST, and can tell me if with all the R tools it can be build a solution like TRICAST or similar. In a word: Do you think that R has all the statistical tools (I mean modeling tools) to make a job similar to TRICAST? Does TRICAST has modeling tools that are not implemented on R yet?

Dear All I am trying to replicate a numerical application (not computed on R) from an article. Using, ks.test() I computed the exact D value shown in the article but the p-values I obtain are quite different from the one shown in the article. The tests are performed on a sample of 37 values (please see "[0] DATA" below) for truncated Exponential, Pareto and truncated LogNormal

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

Dear List, As someone who is in the process of trying to migrate from Excel, I'd appreciate any help on this question: I have a data set and want to fit, say, three distributions to it. I would like to create a plot that shows my data points against all three fitted curves (estimated d.f.). Basically, I lookint to creat a plot that looks like the one presented in the attached paper (Figure

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

Dear List Members, I need a bit help about fitting some theoretical distributions (such as geometric, exponential, lognormal or weibull distribution) to the following *dry spell*, *wet spell*, *cycles (Wet-Dry or Dry-Wet)* from my meteorological (daily rainfall) data http://www.angelfire.com/ab5/get5/R.rainfall.txt only for rainy seasen (july - september) of 14 years only:

[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 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: