similar to: R beginner

I'm trying to do a linear regression between the columns of matrices. In example below I want to regress column 1 of matrix xdat with column1 of ydat and do a separate regression between the column 2s of each matrix. But the output I get seems to give correct slopes but incorrect intercepts and another set of slopes with value NA. How do I do this correctly? I'm after the slope and

Full_Name: bernardo moises lagos alvarez Version: 2.4.0 OS: Windows XP professional Submission from: (NULL) (152.74.219.16) I need obtained the MLE of weibull parameters using the nlm with exact gradient an hessian. I am doing. bug report :Erro en log(b) : el argumento "b" est? ausente, sin default 1.Construction to objectiv functin with n=1 data

[This email is either empty or too large to be displayed at this time]

Hi guys, I'm trying to use lmomco package. first I did the manual calculation on what is the estimates scale and location parameter given L-CV=0.2, L1=1000 L-moments and k (shape parameter) =- 0.1. so what i get is: location: 821.0445 scale: 260.7590 shape: -0.1000 #I assign this as GEV vectors using vec2par GEVpara2<-vec2par(c( 821.0445 , 260.7590 ,-0.1),'gev') #then I

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,

Dear forum users, It's 3:35am and I am swamped with statistics homework lol I'm terrible with R and this time I have no idea what the prof wants. Here is the question: Consider the (two-??sample) Wilcoxon rank statistic T = ?rank(Xi). For n1=106 and n2=192, determine by simulation the ?=.05 critical point for testing H0: ?=0, H1:?<0. We can do this as follows: For m=10000 (no wimpy

  Hi r-users,   I would like to use gev and my data (annual rainfall ) is as follows:   > head(dat,20) A B C D E F G H I J 1 45.1 41.5 58.5 50.1 46.0 49.1 37.7 49.1 59.8 54.0 2 50.3 39.8 49.4 56.4 49.4 48.8 42.1 49.8 49.4 58.3 3 41.7 39.3 44.6 39.1 35.7 41.5 40.8 40.8 38.5 45.6 4 50.7 33.9 48.4 28.2 35.5 39.1 61.4 17.0 30.7 38.3 5 39.3 30.6 46.9 23.8 25.8

i need HELPPP!! how do i calculate the RMSE value for two GEV models?first GEV is where the three parameters are constant.2nd GEV model a 4 parameter model with the location parameter is allowed to vary linearly with respect to time while holding the other parameters at constant. is there any programming code for this? i really really need help. please reply to me as soon as possible. thanks in

Hi All, I am a newcomer to R. Could anyone explain me how to define link functions for either mu/sigma to allow for quadratic trends in the same, when fitting non-stationary GEV distributions? Thanks -- View this message in context: http://www.nabble.com/QUADRATIC-TREND-FOR-LINK-FUNCTIONS-ON-NON-STATIONARY-GEV-tp23360751p23360751.html Sent from the R help mailing list archive at Nabble.com.

Hi, I'm new to R, and am trying to embed R into another application. I'm calling gev.fit() from the ismev package, and it is crashing somewhere inside it. gdb is not catching it, and valgrind is not showing any memory corruption issues. I suspect it's memory corruption, because it doesn't crash in exactly the same spot each time. I'm running R 2.12.2 on a 64 bit linux (Ubuntu

Can anyone explain to me what this error message means, why I'm getting it, and how to fix it? lines(lowess(xdat, ydat, f=.5), col=3) Error: NAs in foreign function call (arg 1) ______________________________________________________________________ Stuart Luppescu -=-=- University of Chicago ºÍʸ ¤ÈÃÒÆàÈþ¤ÎÉã(EUC) -=-=- s-luppescu at uchicago.edu

Dear R-users, I am very enthusiastic about the capacities offered by the np package and pretty fond of it. Nevertheless, trying to estimate Ichimura and Klein and Spady models on my data, I would like to be able to provide the npindex function with my guesses for the bandwidth (eventually computed by other means). The fact that the option 'bandwidth.compute' in the function npindexbw

Hello, I am using the lmomco package (lmom.ub and pargev) to compute the GEV parameters (location, scale, and shape), which are used to estimate return values. I was wondering how/if I can calculate upper and lower confidence (CI_u, CI_l) intervals for each return frequency using the GEV parameters to fill-in the table below? Xi (location) = 35.396 Alpha (scale) = 1.726 Kappa (shape) =

Dear All- I am looking at the bivariate point process model for extreme values as described in Section 8.3.2 of Coles book. I am stuck at a very low level in being unable to reproduce the initial step of transforming the data to unit Frechet margins. In particular, I am unable to reproduce the wave-surge pairs distribution after transformation to Frechet scale (Figure 8.6 of Coles

Can anyone show me how to add a log Pearson type III plot using the evdistq() command to an extreme value plot using the lmom package? Attached sample code below... Thanks in advance, Dave library(lmom) # annual maximum daily streamflows Mackenzie River mackenzieRiver = c(26600, 30300, 34000, 32000, 29200, 28300, 28600, 26400, 28300, 28800, 29000, 22100, 32900, 31800, 21600, 32100, 27000,

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

Hola a todos, Soy nuevo en R y estoy intentando modelizar una serie de datos no estacionarios usand la distribucion Generalizada de Valores Extremos GEV. ¿Podriais indicarme como se modeliza una tendencia polinómica (cuadrática, por ejemplo) en alguno de los 3 parámetros (situación, escala o forma)? He encontrado documentación a cerca de modelización linear o exponencial, pero no acabo de

Hi, I have been having difficulties in finding packages/ codes that simplify plotting of a GEV fitted to dataset (by L-moments) that would print out graph comprising dataset versus gumbel reduced variate n return period at the same. Anyone can help me on this? Thanks. [[alternative HTML version deleted]]

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.