search for: wle

Hi all -- I'm seeing the following commands reliably produce a crash in RGui, version 2.0.1, for both my home and office machine: > rm(list = ls(all = TRUE)); > load("dataset.R"); > library("wle"); > data.wle = wle.lm(abortion ~ year * lib.con + age + gender + + urbanism + census + income + church.att + children + educ + + religion.imp, data = data.set); dataset.R is moderately sized (about 200k compressed), and the command works just fine with "lm" rather than &qu...

-------------- next part -------------- An embedded message was scrubbed... From: Hella Heikki <Heikki.Hella at bof.fi> Subject: R code and robust regression Date: Tue, 12 Dec 2000 12:12:14 +0200 Size: 2187 Url: https://stat.ethz.ch/pipermail/r-help/attachments/20001212/c3361c7d/attachment-0001.mht

Dear [R]-list, I'm looking for a classic equivalent of the wle.gamma function (library wle) that estimate robustly the shape and the scale parameters of gamma data. I have a vector of iid gamma rv : >data=rgamma(100,shape=10,scale=3) and a vector of their weights: >weights=c(rep(.5/70,70),rep(.25/20,20),rep(.25/10,10)) and want to estimate the scale...

Subject: Eliminate repeated components from a vector X-Mailer: VM 7.00 under 21.4 (patch 6) "Common Lisp" XEmacs Lucid Reply-To: fjmolina at lbl.gov FCC: /home/f/.xemacs/mail/sent Does anyone know how I can eliminate repeated elements from a vector?

...) and I installed a few packages using install.packages() . Now, I am trying to upgrade these packages using update.packages() but I believe I am running into some kind of error. See below. Any tips will be greatly welcome. Thanks in advance. Francisco PS: update.packages() only tries to update wle from version 0.3 to 0.4, gives me an error message at the end, and installed.packages() continues to list version 0.3 as the installed version. I ran install.packages("package-name") as root, and am running update.packages() as root. I must be missing something simple. ------------...

Hi, all When I am using mle.cv(wle), I find a interesting problem: I can't do leave-one-out cross-validation with mle.cv(wle). I will illustrate the problem as following: > xx=matrix(rnorm(20*3),ncol=3) > bb=c(1,2,0) > yy=xx%*%bb+rnorm(20,0,0.001)+0 > summary(mle.cv(yy~xx,split=nrow(xx)-1,monte.carlo=2*nrow(xx),verb...

Hello, I have a question pertaining to the stepwise regression which I am trying to perform. I have a data set in which I have 14 predictor variables accompanying my response variable. I am not sure what the difference is between the function "mle.stepwise" found in the wle package and the functions "step" or "stepAIC"? When would one use "mle.stepwise" versus "step/stepAIC"? Other than the references in the R-software, is there anything which discusses the use of "mle.stepwise"? Thank you kindly in advance for a...

R users   Thanks in advance:   I need to use the function "mle.cp" in R.  When I run this function I recive   result <- mle.cp(y.hald~x.hald) Error: could not find function "mle.cp"  could some one  give me a help?   thanks [[alternative HTML version deleted]]

..._idlc_ -- fracdiff /usr/bin/ld: multiple definitions of symbol _gammfd_ (and others) -- odesolve -- princurve lamix: Error on line 116: Declaration error for v: adjustable dimension on non-argument -- splancs /usr/bin/ld: multiple definitions of symbol _bounds_ -- wle mleaic: Error on line 229: Declaration error for work: adjustable dimension on non-argument -- subselect configure: error: R was configured without a FORTRAN compiler ERROR: configuration failed for package `subselect' The multiple defined symbols are a well known OS X problem. There ar...

Buenas tardes, Estoy intentando ajustar distribuciones utilizando un vector de ponderación en los datos (Weighted Likelihood). ¿Existen paquetes en R que resuelven esto? He mirado ya el paquete "wle" pero no me permite introducir los pesos mediante los cuales ponderar los datos. En un primer momento, se me ha ocurrido realizar lo siguiente: repetir cada elemento del vector datos tantas veces como indique el vector pesos ?modificado? (multiplicado por una potencia de 10, para que sean to...

...ust') Warning message: In update.views("Robust") : The following packages are not available: covRobust, distr, FRB, MASS, mblm, multinomRob, mvoutlier, quantreg, RandVar, rgam, RobAStBase, robfilter, RobLox, RobRex, robust, RobustAFT, robustbase, ROptEst, ROptRegTS, rrcov, sandwich, wle > library(covRobust) > It is puzzling because — as the second command shows — the package covRobust is installed. Here is my > sessionInfo() R version 2.14.2 (2012-02-29) Platform: x86_64-apple-darwin9.8.0/x86_64 (64-bit) locale: [1] en_US.UTF-8/en_US.UTF-8/en_US.UTF-8/C/en_US.UTF-8/en...

...dows/contrib are now for rw1020. If you are upgrading to rw1020, please replace all the compiled packages you have downloaded by ones from this set. There a few packages available for the first time: RMySQL mlbench netCDF and a few that are not yet running on 1.2.0 on any platform hdf5 locfit wle (crashes with an array overrun). Beware of the following (from the current rw-FAQ) 3.5 Loading a package fails. ============================ Is the package compiled for this version of R? Many of the packages need to be compiled for a fairly recent version, and in particular a...

...dows/contrib are now for rw1020. If you are upgrading to rw1020, please replace all the compiled packages you have downloaded by ones from this set. There a few packages available for the first time: RMySQL mlbench netCDF and a few that are not yet running on 1.2.0 on any platform hdf5 locfit wle (crashes with an array overrun). Beware of the following (from the current rw-FAQ) 3.5 Loading a package fails. ============================ Is the package compiled for this version of R? Many of the packages need to be compiled for a fairly recent version, and in particular a...

Is there a painless way to find the names of all packages on CRAN that "Depend" on a specified package? url: www.econ.uiuc.edu/~roger Roger Koenker email rkoenker at uiuc.edu Department of Economics vox: 217-333-4558 University of Illinois fax: 217-244-6678 Champaign, IL 61820

...tiv mvnmle mvtnorm netCDF nlme nlrq norm odesolve oz panel pcurve permax pinktoe pixmap polymars polynom princurve pspline quadprog quantreg rmeta rpart rpvm scatterplot3d sem sgeostat sm sma sn sna splancs sptests spweights strucchange survival tensor tree tripack tseries twostage vegan wavethresh wle xgobi xtable zmatrix -------------- next part -------------- An HTML attachment was scrubbed... URL: https://stat.ethz.ch/pipermail/r-help/attachments/20011012/5b0887c1/attachment.html

Hello, does R have functions for setting and testing bit values? I want to conserve memory for storing presence/absence data for large multiple arrays within a single array, using element values like present[x,y] <- ntharray[x,y]*(2^n) where presence is 1, non-presence is 0 and n is the nth array e.g. 1*(2^0) + 0*(2^1) + 0*(2^2) + 1*(2^3) + 0(2^4) for storing the value 9 for presence in

Dear R-helpers, I was wondering if anyone has or knows someone who might have an implementation of algorithm for estimating SETAR models including the lag-order. For some reason my code gives me a bit wrong results. I am fighting with it for a week and cannot bring it down. Thanks a million in advance, Sincerely, Evgueni McGill University Department of Economics

Hola! Anybody knows if there exists somewhere in R some implementation of robust Poisson regression, where robust is taken in the sense as usen in rlm(MASS). I found something in the package wle, but only for the Poisson distribution, not for regression. For the moment I try to use linear models with the square-root transformation, and rlm. Kjetil -- Kjetil Halvorsen. Peace is the most effective weapon of mass construction. -- Mahdi Elmandjra

...e regression (forward as well as the backward regression). >From other software (i.e. minitab), I know only X1 and X2 are significant so my regression equation should be Y = constant + (b1) X1 + (b2) X2. I am not sure whether I am using the correct R commands, which are as follows - library (wle) ONS <- mle.stepwise(Y ~ X1+X2+X3+X4, data=ONS) summary(ONS) The output looks something like this OUTPUT Forward selection procedure F.in: 4 Last 3 iterations: (Intercept) X1 X2 X3 X4 [1,] 1 0 0 0 1 15.57 [2,] 1 1 0 0 1 14.80 [3,] 1...

Here's a more general question following up on the specific question I asked earlier: Can anybody recommend an R command other than mle.aic() (from the wle package) that will give back a ranked list of submodels? It seems like a pretty basic piece of functionality, but the closest I've been able to find is stepAIC(), which as far as I can tell only gives back the best submodel, not a ranking of all submodels. Thanks in advance, Alexandra