similar to: Troubleshooting with "gpd" (Fit generalized pareto model)

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 all R users, My simulation function works correctly, but I have problems with plot function. You will find the following code using it. Thank you for your help ##################################################" simulation <- function(k, n){ conc <- seq(0,10,by=0.5) #choixg <- seq(1, length(conc)) choixg <- rep(0,length(conc)) for (i in 1:length(conc)){ choixg[i]

Hi R-users, I have to define a noise level function L and its energy in the various moment of the day by: if time is between 18:00:00 and 23:59:59 then L[j] <- L[j]+5 and W <- 10^((L+5)/10) if time is between 22:00:00 and 05:59:59 ==> L <- L+10 and W <- 10^((L+10)/10) else L=L and W = W Could someone help me to realize this function please? You will find my following

Hi! R users I would like to ask you where could we find the VGAM Package. I don't find it in the list of packages. Thak you for your help Lassana KOITA Etudes de S?curit? et d'Exploitation a?roportuaires / Safety Study & Statistical analysis Service Technique de l'Aviation Civile (STAC) / Civil Aviation Technical Department Direction G?n?rale de l'Aviation Civile (DGAC) /

Hi all R users, I have problems with my second loop for drawing the three curves in the same graphic. I need help please Thank you in advance ######################################################################### simulation <- function(k, n){ conc <- seq(0,100,by=0.5) #choixg <- seq(1, length(conc)) choixg <- rep(0,length(conc)) for (i in 1:length(conc)){ choixg[i] <- (k

Hi, Dear R users I have problem with the following code. The matrix result must be a matrix (3x3). But I have obtained a matrix(3x1) and I don't know why. So, I need your help Best regards ##################################################################### taille <- function (delta, t, prob = 0.2) { niv.conf <- c(0.90, 0.95, 0.99) if(niv.conf <- 0.90) { t <- 1.645 }

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 have a problem in fitting GPD distribution. i generate random numbers from gpd distribution from specific parameters using pot packege then i used fitgpd function to estimete the parameters.The estimated parameters are not matched with the given parameters i.e.from which i generate random numbers.I think estimated parameters should be matched with the given parameters.Also suggest me

Using 'Samba version 2.2.1a' Please could someone offer advice in order to configure the following. I need domain level authentication for most shares, however 1 share needs to be made available to a none domained NT client, such that perhaps share or user level security would be appropriate. How can this be configured on a single samba server. Thanks in advance to anyone who responds

hi, we are using a voip h323 switch. the switch sends all caals to our Gatekeeper (gnugk). gnugk musst send all calls to asterisk and asterisk must do his choice (sip endpoint or out to PSTN) Making calls to our h323 switch works fine over asterisk. what must i configure to get inboung h323 calls from our gnugk to asterisk? any hints for me? thx -- Thomas K?pper 01063 Telecom GmbH &

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

Hi List, we are using openh323 gatekeeper for voip telefony. We also have a voip over ss7 TELES Switch for voip into POSTN Network. Know we want to use Asterisk for converting SIP to h323. Now my question. Is Asterisk an full h323 gatekeeper like openh323? Do we need openh323 GK for astrerisk, too?. And how can i tell asterisk to sent all none SIP-ip calls to the gatekeeper over h323? thx

Dear All, I have a question about package "ismev", its function "gpd.fit", and interpretation of the results. I used the package ismev to do an extreme value analysis on a fire dataset. Two variables are used in the analysis. The focal variable is acreage burned per fire, ranging from 1 to 5000 acres per fire. In total, there are 69,980 observations. The date covers

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

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

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

with a h323 client over my gatekepper a call comes over asrerisk to my sip endpoint: == Spawn extension (sip-phones, 01634255122, 1) exited non-zero on 'SIP/0699073201-528d' -- Executing Dial("H323/ip$10.0.0.124:49638/18690", "SIP/0699073201") in new stack -- Called 0699073201 -- SIP/0699073201-dc61 is ringing -- SIP/0699073201-dc61 answered

Hi ! I have a small problem. When switching a call (pstn -> sip user), I get the sip phone ringing - ie. everything is OK, but I do not get a ringtone in the handset on the pstn side. Can anyone help me out in how to make * play tones ? My setup: E1 IP pstn ------ Asterisk ------ sip phone Regards, Dave

Hi all, I'd like to export the results of GPD or GEV analysis generated with the Extremes Toolbox for plotting in Grapher or Excel (for manipulation by our publications group). Is it possible to access the data used to generate the plots in extRemes, or do I need to code the analysis from scratch and write to file? Thanks, Ben Stay connected to the people that matter most with a

Hi, I'd like to export the results of GPD or GEV analysis generated with the Extremes Toolbox for plotting in Grapher or Excel (for manipulation by our publications group). Is it possible to access the data used to generate the plots in extRemes, or do I need to code the analysis from scratch and write to a file? Stay connected to the people that matter most with a smarter inbox.