similar to: Non linear modeling

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

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

Hi, I'm exploring the ismev package and in particular the pp.fit function. Documentation says that this function performs Maximum-likelihood fitting for the point process model, including generalized linear modelling of each parameter. Actually I'm afraid not to understand what it means with "Point process model". Any suggestion on possible references for further infos? thanks

Hi all, I've only now noticed that on our Opteron boxes running Suse Enterprise Server 8 that R-2.0.1 did not pass make chek-all. It bombed in tests/lapack.R at: > eigenok <- function(A, E, Eps = 1000 * .Machine$double.eps) { V <- E$vect lam <- E$values stopifnot(abs(A %*% V - V %*% diag(lam)) < Eps, .... [TRUNCATED] > Ceigenok <- function(A, E, Eps =

R helpers, I am working with the R program performing multivariate regression trees (MRT). I have a matrix with species and environmental variables saved as a CSV file (sprot_matrix.csv), I have 42 species and 8 environmental variables (SECCHI+PH+TA+PTOT+NTOT+CHLA+AREA+ MEANDEP) for 104 samples Title SpA SpB SpC SpD Varible1 Variable2 Variable3 Sample1 Sample 2

Hi people, Does anybody know some Density-Based Method for clustering implemented in R? Thanks, Fernando Prass _______________________________________________________

first:I'd like to have the choice between breaking out a repeat loop or continue it then: i'd would like my code after the end of my repeat loop not to be read unless I type 2 (to continue executing the content of my repeat loop) THAT'S THE PROBLEM Is there any way to go straight out of the loop (like a "goto" associated to a "label") or a pause like command to

How is it possible to "split" a data.frame in order to get classes with variable size. actually I'd like to get classes of data with classe size (cs) so that cs becomes bigger with an other increasing value. Guillaume Storchi

I'm doing a non linear regression with 8 parameters to be fitted: J.Tl.nls<-nls(Gw~(a1/(1+exp(-a2*Tl+a3))+a4)*(b1/(1+exp(b2*Tl-b3))+b4),data=Enveloppe, start=list(a1=0.88957,a2=0.36298,a3=10.59241,a4=0.26308, b1=0.391268,b2=1.041856,b3=0.391268,b4=0.03439)) First, I fitted my curve on my data by guessing the parameters'

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, How can I define a new psi function in the rlm comand? In particular I would like to implement the case for \rho(u)=0.5*(u^2) and \psi(u)=u in order to assume normally distributed errors. Any help? Thanks -- ======================================================== Angelo Secchi PGP Key ID:EA280337

I've cc'd this to the package maintainer, Andy Liaw <andy_liaw at merck.com>. I'm not sure he reads this list. Duncan Murdoch On 2024-10-23 1:26 a.m., Stevie Pederson wrote: > Hi, > > It appears there is an OSX-specific bug in the function > `randomForest.default()` Going by the source code at >

Hi R user, I am a kind of an intermediate user of R. Now I am using GLM model (library MASS, VEGUS). I used a backward stepwise logistic regression, but i got a problem in removing those predictors which are above 0.05. I don't want to include those variables which were above 0.05 in final backward stepwise logetsic regression model. for example: first I run the model,

Hi, It appears there is an OSX-specific bug in the function `randomForest.default()` Going by the source code at https://github.com/cran/randomForest/blob/master/R/randomForest.default.R the bug is on line 103 If the vector `cutoff` is formed using `cutoff <- rep(1/9, 9)` (line #101) the test on line 103 will fail on OSX as the sum is greater than 1 due to machine precision errors. sum(rep(1

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