similar to: ismev package

Hi, I'm using the ineq package to calculate some concentration measures (Gini, Herfindal, ...) and I was wondering if there's around also a function to calculate standard error on these measures. If not, is anybody aware of where I can find a reference on this point? Thanks. -- ======================================================== Angelo Secchi PGP Key ID:EA280337

Hi, after looking around I was not able to get info about these two questions: 1. Is there a function to have a "jackknifed corrected var/cov estimate" (as described in MacKinnon and White 1985) in a standard OLS regression? 2. Does R possess a LAD (Least Absolute Deviation) regression function? Any help? Thanks -- ======================================================== Angelo

AFAIK most model fitting techniques will only deal with additive errors, not multiplicative ones. You might want to try fitting: log(y-x) = a*x + e which is linear. Andy > From: Angelo Secchi > > Hi, > is there a way in R to fit a non linear model like > > y=x+exp(a*x)*eps > > where a is the parameter and eps is the error term? > Thanks > Angelo > >

Hi, I'm doing robust regression with the following command rlm(dip~ind1+ind2-1,method="M",psi=psi,maxit=1000,acc=1e-15) now when I ask for a summary summary(rlm(dip~ind1+ind2-1,method="M",psi=psi,maxit=1000,acc=1e-15)) I get Coefficients: Value Std. Error t value ind1 -0.0377 0.0000 -24203.1415 ind2 1.0370 0.0000 668735.7195 taht is

Hi, I'm trying to write a simple function like case1 <- function (m, cov, Q, R) { theta <- (acos(R/sqrt(Q^3))) beta <- (-2)*sqrt(Q)*cos(theta/3)+m[1]/3 rho1 <- (-2)*sqrt(Q)*cos((theta+2*pi)/3)+m[1]/3 rho2 <- (-2)*sqrt(Q)*cos((theta-2*pi)/3)+m[1]/3 stderrb <- deltamethod( ~(-2)*sqrt(Q)*cos(theta/3)+x1/3,m,cov) stderrr1 <- deltamethod(

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

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

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?

Hello there! Perhaps this is a little off-topic, but I notice this only on the Centos box. I'm running Centos 4 on an AMD64 which has the following entries in the fstab to connect to NFS shares on a Fedora3 box: 192.168.1.12:/home/angelo/ /home/angelo/NFS_share1 nfs rw,addr=192.168.1.12 0 0 192.168.1.12:/home/angelo/data /home/angelo/NFS_share2 nfs rw,addr=192.168.1.12 0 0

Hi, I am trying the following: library(ismev) library(evd) fit <- gev.fit(x,show=FALSE) ks.test(x,pgev,fit$mle[1],fit$mle[2],fit$mle[3]) but I am getting: Warning message: cannot compute correct p-values with ties in: ks.test(x, pgev, fit$mle[1], fit$mle[2], fit$mle[3]) where x is: [1] 239 38 1 43 22 1 5 9 15 6 1 9 156 25 3 100 6 [18] 5 100

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

Friends, Sorry my english, i'am speaking from brazil. I have a problem with an delphi application running on wine 0.93. The application open exactly, but the buttons not it appears correctly. I send the link of screenshot: http://200.195.11.195/~prolinx/screen01.jpg any help is coming well. Thanks!!! Angelo Braga

---------- Mensagem reenviada ---------- Subject: Re: [Wine] Delphi Application Date: Dom 01 Jan 2006 14:32 From: Bruno <cbruno@ferreiracosta.com.br> To: wine-users@winehq.org Cc: truls@asheim.ne Olah Angelo Braga, tambem estou no Brasil, e meu ingles eh pessimo, Acho que sei qual o problema com seus icones! Jah tive um problema parecido, e resolvi instalando fontes TrueType Originais,

One more follow on, The combination of kernel.panic=60 and kernel.printk=7 4 1 7 seems to have netted the culrptit: E01-netconsole.log:Jan 18 09:45:10 E01 (10,0):o2hb_write_timeout:137 ERROR: Heartbeat write timeout to device dm-12 after 60000 milliseconds E01-netconsole.log:Jan 18 09:45:10 E01 (10,0):o2hb_stop_all_regions:1517 ERROR: stopping heartbeat on all active regions.

Dear all, first of all I apologize for not having changed the object. I just re-used an old email I sent some time ago. The let us go into the question. Our architecture is the following: A (set of) java programs running under a jvm that passes data and instructions to an R instance via RNI and at the end of the process extracts calculated data from R. The question is: can I

Hi, I'm new to R and have a question which seems like it should be easy. I loaded the R port of my boot library and installed it. All of the functions and help files are available but I can't seem to access the data. The data directory was created and contains files with the extension .rda. I tried to use the data() function but I get > data(city) Error: restore file corrupted -- no

Hi, Suppose that I have a list where each component is a list of two matrices. I also have a vector of weights. How can I collapse my list of lists into a single list of two matrices where each matrix in the result is the weighted sum of the corresponding matrices. I could use a loop but this is a nested calculation so I was hoping there is a more efficient way to do this. To help clarify,

Hi, I want to sort a matrix by a specific variable without changing the row binding between variables. Ex. Name Sex Age Fred M 24 John M 18 Mary F 21 ordered by Age Name Sex Age John M 18 Mary F 21 Fred M 24 Thanks Angelo

Hi, how can I treat data organised in classes and frequencies? Ex. class frequency 20-23 9 23-25 7 26-28 5 29-31 5 32-34 3 Thanks Angelo Scozzarella

Hi, Why this function doesn't work? function (x) { if (is.factor(x)) { if (!is.ordered(x)) { warning("La mediana non si puo' calcolare!!!") return(NA) } me <- median(unclass(x)) if (me - floor(me) != 0) { warning("Mediana indeterminata") return(NA) } else