similar to: CDF for pareto distribution

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

Hi, I'm a second year Master's student in Applied Statistics. I am doing a project using average weekly U.S. regular gasoline prices (in cents, per gallon) from an Excel file (from the years 1990- May 2010). I want to find the probability that the average weekly U.S. regular gasoline prices (in the long term) goes over 400 cents a gallon (or $4.00 a gallon). I am using the

Hi, I'm trying to evaluate a Monte Carlo p-value (using truncgof package) on a left truncated sample. >From an empirical sample I've estimated a generalized pareto distribution parameters (xi, beta, threshold) (I've used fExtremes pkg). I'm in doubt on what of the following command is the most appropriate: Let: x<-sample t<-threshold xt<-x[x>t] xihat<-gpdFit(x,

Dear all, I am trying to get the emperical bayes estimates together with their standard errors out of lme4. Up to now I have used MLwiN to get these estimates. I have fitted the following - very simple - model, just to find out how this works. test<-lmer(y~(1|subject),data,REML=F) ranef(test,postVar=T) str(ranef(test,postVar=T) If I use the formulation of the emperical bayes estimates and

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}

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 friends. In Carlin and Louis "Bayes and emperical Bayes methods.." 1996 the classical example of 12 independent tosses of a fair coin producing 9 heads and 3 tails is given. If the situation is seen as a fixed sample of 12, a binomial lieklihood is used, and Carlin et al reports a probability of 0.075. Using sum(dbinom(9:12,12,.5)) I obtain 0.073 Likewise, if the experiment is

Hi all, I would like to plot part of the emperical CDF. Suppose the variable is x, I just need the part when x>1,therefore, I am using the following codes. tail <- x>1 plot(ecdf(x[tail]), do.points=FALSE, verticals=TRUE) The "x" value starts from 1, but the yaxs still begins from 0, not the corresponding value when "x" is 1. How can I make it match? Could anyone

Hi all. Belove is the example from the cluster-help page wtih the output. I simply cannot figure out how to relate the estimate and robust Std. Err to the p-value. I am aware this a marginal model applying the sandwich estimator using (here I guess) an emperical (unstructered/exchangeable?) ICC. Shouldent it be, at least to some extend, comparable to the robust z-test, for rx :

Hello. I am not certain even how to search the archives for this particular question, so if there is an obvious answer, please smack me with a large halibut and send me to the URLs. I have been experimenting with fitting curves by using both maximum likelihood and maximum spacing estimation techniques. Originally, I have been writing distribution-specific functions in 'R' which work

Hello all, I am trying to estimate the cumulative distribution function for a single stock return time series. A piecewise estimation is composed of three parts: parametric generalized Pareto (GP) for the lower tail (10% of the observation), non-parametric kernel-smoothed interior (80% of the observations), and GP for the upper tail (10%). I wonder if anyone has clue about this in R. The

Dear R users, Could you please help me out. I am in trouble as I am unable to model graphs to explain the GARCH (1 , 1) model, the Hill estimator (of alpha), and the Pareto estimator. I just got introduce to R. I am working on a paper which must be worked from R. You look at the difficulty I had from the text below. [1] "DAX" "DAX_CAC" "DAX_CAC40"

Dear All, I am trying to fit some data to a Pareto distribution and would like to estimate the parameters with the fitting. I have come across some options so far. Unfortunately I haven't managed to get any of them to make the right fits (as is evident when I check with the goodness of fit). One such option is: library(VGAM) b1 <- read.table(file("FitPareto_Values.txt",

I would like to fit a Pareto Distribution and I am using the following codes. First, 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 the same value. fit=vglm(ycf1 ~ 1, pareto1(location=alpha), trace=TRUE, crit="c") coef(fit, matrix=TRUE) summary(fit) fitted(fit) Secondly, how can I plot the

Dear R-users, I am trying to fit pareto distribution to some data but i've one problem. Using optim to calculate the maximum of the likelihood function of the pareto I use as start parameters the moments method(using the distribution function in the package actuar): media=mean(x) var=mean(x^2)-media^2 scale=2*var/(var-media^2) shape=(scale-1)*media

Dear R users, the first 1.0 version of the rPref package is now on CRAN. rPref allows to select the Pareto-optimal tuples from a data set, also called Skylines in the database community. For example, optimal tuples from mtcars according to "high(mpg) * high(hp)" (where "*" is the Pareto operator) are those cars, for which no other dominating car exists. There, a car

Dear all: First, apologies for cross-posting multiplicities and for a query that is more analytically related than S-language related. The bottom-line wish is: Could you please provide and advice, references, etc on S software approaches for fitting a distribution with density: p*g(x) + (1-p)*f(x) where g(x) is the familiar lognormal 2-parameter density and f(x) is Pareto as defined below?

Ukech U. Kidi wrote: > dax<- diff(log(DAX_CAC$DAX[1:1865])) > m1<- garch(dax) > Error: couldn't find function "garch" > m1<- garch(dax[1:1865]) > Error: couldn't find function "garch" > m1<- garch(dax[1:1865]) I am sorry, but I forgot to change the addres to r-help in the reply. Well, I am not sure, wheere do you want to get

Hi, I'm fitting a generalized Pareto distribution to POT exceedances of a data set. The practical stuff works ok, but I have a question regarding theory. Is there an equation relating parameters of a gpd tail to its (first) moments? According to theory for certain parameters either the first moment does not exist or the distribution has an upper bound, but I haven't found the

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