similar to: fExtreme packages

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

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,

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

Good day good people, I have troubble with printing. I have samba 3.0.14a-Debian, on Debian Stable Sarge. There is my smb.conf # Global parameters [global] dos charset = CP852 unix charset = ISO8859-2 workgroup = IRIXAS server string = Samba Irix a.s. null passwords = Yes obey pam restrictions = Yes pam password change = Yes passwd

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

Prof Brian Ripley <ripley at stats.ox.ac.uk> wrote: mkdir fBasics unzip fBasics.zip -d fBasics rm fBasics/src/*.o R CMD check fBasics and that took me about 3 minutes. Now me, I just did unzip -a fBasics_190.10051.zip R CMD INSTALL fBasics rm -rf fBasics in a naive and trusting manner. It took me considerably longer than 3 minutes to learn that this was what I should do, and

Hello, and please advise regarding any errors/omissions on my part. However, the documentation for R's tcplot function in the evd package appears to contain an error. I am using evd version 2.2-4 in R 2.14.0 with Windows XP. > data(portpirie) > mrlplot(portpirie) ## No Error > tlim <- c(3.6, 4.2) > tcplot(portpirie, tlim) ## Line from documentation Error in fpot(data, u[1],

I am estimating a large model by groups. How do you save the results and?returns the associated quantiles? For this example I need a data frame n?? ?xi??????? mu????????beta 1?? 0.1033614? 2.5389580 0.9092611 2? ?0.3401922? 0.5192882 1.5290615 3?? 0.5130798? 0.5668308 1.2105666 I also want to apply gevrlevelPlot() for each "n" or group. ? #Example n <- c(1, 1, 1, 1, 1, 1, 2, 2, 2,

rpcclient localhost -d 4 -U'root%Passwd' -c 'adddriver "Windows NT x86"

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

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

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

Hi all, Samba 3.0.7 and using CUPS printing. I am logged in as root to an XP Pro system, I have used the APW to install the XP drivers to all of my printers. I went to Shared and clicked "Additional Drivers" and checkmarked both 9x and NT/2K boxes. Using drivers straight from the HP website, the 9x driver installs fine, but it seems like I can't get any NT/2K drivers to

Dear all, I'm trying to produce a log-likelihood profile for a GPD estimate with the ismev package. When following the examples, everything is working fine, i.e. I get a nice parabolic log-likelihood curve. When using my own data (intensity of hurricanes at landfall), the log-likelihood curve looks like a step function. Am I doing something wrong here? Is my data not suited for the GPD fit?

I've read in the documentation that Samba does not use ports like Win2k. So how does Samba know where to find a network printer if it has it's own IP address (stand-alone printer)? I have created a printer in my /etc/printcap file: HP4000:\ :sd=/var/spool/lpd/hp4000:\ :mx#0: :sh:\ :rm:192.168.2.13:rp=raw: :rp=:\ added this to my smb.conf [printers] comment = All Printers printable =

-----BEGIN PGP SIGNED MESSAGE----- Hash: SHA1 | [root@enigma 0 3]$ rpcclient enigma -U 'root%XXXXXXXXXX' -c 'enumdrivers' | | [Windows NT x86] | Printer Driver Info 1: | Driver Name: [Lexmark Z53 Color Jetprinter] | | | [Windows NT x86] | Printer Driver Info 1: | Driver Name: [Lexmark Z53 Color Jetprinter] | | [root@enigma 0 3]$ rpcclient enigma -U

In the process of learning R, with a specific interest on financial time series. While I continue to get through the documents I am more a fan of learning by example and then looking up how each function is used. Any websites which post sample code for R? Thanks in advance. -- Jonathan jqm475 at gmail.com

Hi, I am new to bioconductor, trying to install KEGGSOAP package, but got warnings() when installing and error message when trying to load the package, can anyone suggest what went wrong? many thanks John > source("http://bioconductor.org/biocLite.R") Bioconductor version 2.11 (BiocInstaller 1.8.3), ?biocLite for help > biocLite("KEGGSOAP") BioC_mirror: