similar to: How to create a 'fit' plot

Hello all, I would like to plot the emperical CDF, normal CDF and pareto CDF in the same graph and I amusing the following codes. "z" is a vector and I just need the part when z between 1.6 and 3. plot(ecdf(z), do.points=FALSE, verticals=TRUE, xlim=c(1.6,3),ylim=c(1-sum(z>1.6)/length(z), 1)) x <- seq(1.6, 3, 0.1) lines(x,pgpd(x, 1.544,0.4373,-0.2398), col="red") y

Up to now, I have recognized problems with "gpd(..)", the function from the package "evir" I think that all these functions that estimate the parameters xi, beta for the GPD by given threshold mu use the function "optim(..)" ( gpd, fitgpd, ...) "Error" example: data1 <- rgpd(1000, xi= -1.5, mu=1000, beta=100) so the created poinnts take place in about

Dear R helpers, I need to use KS and AD test for Generalized Pareto and Generalized extreme value. E.g. if I need to use KS for Weibull, I have teh syntax ks.test(x.wei,"pweibull", shape=2,scale=1) Similarly, for AD I use ad.test(x, distr.fun, ...) My problem is fir given data, I have estimated the parameters of GPD and GEV using lmom. But I am not able to find out the distribution

Hi: I am going through some of the xtable examples and I can't make the one below work. I need to create a longtable on the fly keeping the column headers for all the pages and I thought this example could give some ideas on how to do it. I am using Sweave and xtable to create my tables and graphics. I wonder if someone could tell me what's wrong. Thanks ## Not run: \begin{small}

I would like to fit a Pareto Distribution and I am using the following codes. I thought the fitted (fit1) should be the fitted value for the data, is it correct? As the result of the "fitted" turns out to be a single value for all. fit=vglm(ycf1 ~ 1, pareto1(location=alpha), trace=TRUE, crit="c") fitted(fit) The result is fitted(fit) [,1] [1,] 0.07752694

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

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:

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

I'm looking for a company that offers Guatemala DID's. I saw that Lingo does, but Lingo isn't easily compatible w/ Asterisk, so they're a last resort. Thanks in advanced, Phil Astin.

Hello, First of all, let me apologize that my statistics background is modest at best. I am doing some extreme value analysis on model output (WRF) which have the following dimensions: speed(time,lat,lon) I am trying to fit the GPD (gpd.fit) to each point (time,lat,lon) to get a return level plot with values at each grid point. (Map with return level by location.) Here is some code I

Dear List, I have, say, a 2000x1 numeric vector and would like to split it into, say, a 200x10 matrix. Any help is appreciated.

Hello fellow developers, I wanted to follow a more up-to-date llvm version and saw the google/stable branch that seems to be updated every two weeks or so. But it seems that the last update of the google/stable branch is from 2016-08-15 [1]. Did it move to somwhere else and I am lookint at the wrong spot, or maybe there is something terribly wrong with the current trunk? Regards. [1]

Hi, I am trying to migrate print queues from a windows server to a samba share. 1. I ran printmig.exe on the windows server. 2. Backed up the printers on the windows servers to a cab file. 3. Tries to restore the printers to the sambashare by specifying the target as //sambaShare. But I get the following error: 2006:01:25 15:35:38 Access Granted to: \\sambaShare 2006:01:25 15:35:38

Hello, I am a new user of R. I am trying to use the packages fBasics and fExtremes when i am running the examples I get few error. Could someone tell me what is happenig? Thank you beforehand. from Fbasics packages: xmpfBasics() Error in file(file, "r") : unable to open connection In addition: Warning message: cannot open file '/usr/lib/R/library/fBasics/demoIndex'

Hello dear R-helpers, I'm working with R-2.15.2 on Windows 7 OS. I'm stucked with a merge of two data frames by characters. In each data frame I got two different list of names, that is my main-key to be merged. To figure out what I'm saying, I build up a modified "?merge" example, with errors by purpose: # Data for authors: authors <- data.frame( surname =

Dear Specialists in R, May be somebody has experiment in using pakage for the book Quantitative Risk Management by McNeil? This package is writen in R. I have run this package for fitting the data to Nornal Inverse Gaussian distribution and fased with following problem. > Return<-read.csv("data.csv") > Transpose<-t(Return) > fit.NH(Transpose, case="NIG",

Dear All, I am trying to generate some Pareto random variates using the inverse method. This is really straightfoward and my R function looks as below : pareto <- function(c, a, cnt=1000) { u <- runif(cnt) x <- (c / ((u ^ (1 / a)))) mean.theo <- ((c * a) / (a - 1)) mean.gen <- mean(x) cat('Pareto mean : theoritical', mean.theo, 'generated', mean.gen,

Hi, guys, I need to simulate Pareto distribution. But I found 'rpareto' didn't exist in R. And it seems that Pareto distribution don't have mathematical relationships with other distributions. What can I do? Thanks a lot. Ni --------------------------------- [[alternative HTML version deleted]]

Dear all, I am trying to simulate from truncated Pareto distribution. I know there is a package called PtProcess for Pareto distribution...but it is not for truncated one. Can anyone please help me with this? Thanks in advance. Cassie [[alternative HTML version deleted]]

I've gone ahead and updated both stable and testing for the branch. Sorry that it hasn't been updated in a while - I'll make sure it happens more regularly now. Thanks! -eric On Tue, Nov 22, 2016 at 8:25 AM David Blaikie via llvm-dev < llvm-dev at lists.llvm.org> wrote: > I think the process has bitrotted a bit and not sure it's going to > continue to be updated