Displaying 20 results from an estimated 1000 matches similar to: "nlrob problem"
2009 Aug 12
1
psi not functioning in nlrob?
Hi all,
I'm trying to fit a nonlinear regression by "nlrob":
model3=nlrob(y~a1*x^a2,data=transient,psi=psi.bisquare,
start=list(a1=0.02,a2=0.7),maxit=1000)
However an error message keeps popping up saying that the function
psi.bisquare doesn't exist.
I also tried psi.huber, which is supposed to be the default for nlrob:
model3=nlrob(y~a1*x^a2,data=transient,psi=psi.huber,
2013 Mar 15
2
nlrob and robust nonlinear regression with upper and/or lower bounds on parameters
I have a question regarding robust nonlinear regression with nlrob. I
would like to place lower bounds on the parameters, but when I call
nlrob with limits it returns the following error:
"Error in psi(resid/Scale, ...) : unused argument(s) (lower = list(Asym
= 1, mid = 1, scal = 1))"
After consulting the documentation I noticed that upper and lower are
not listed as parameter in
2011 May 16
1
nlrob(...) returns error message
Dear all,
?
I implemented a non-linear model using nls(...) and it works just fine. I now tried to run the same model using nlrob(...) which basically does the same but uses a more robust estimation procedure. My problem: I cannot seem to get nlrob(...) running. Irrespective of how I try to call the function, I always get the error message "Error in is.null(x) : 'x' is missing".
2008 Jan 19
1
How do we get two-tailed p-values for rlm?
How do we get 2-tailed p-values for the rlm summary?
I'm using the following:
> fit <- rlm(oatRT ~ oatoacData$erp, psi=psi.bisquare, maxit=100,
na.action='na.omit')
> fitsum <- summary(fit, cor=F)
> print(fitsum)
Call: rlm(formula = oatRT ~ oatoacData$erp, psi = psi.bisquare, maxit = 100,
na.action = "na.omit")
Residuals:
Min 1Q Median
2018 Apr 06
1
Fast tau-estimator line does not appear on the plot
R-experts,
I have fitted many different lines. The fast-tau estimator (yellow line) seems strange to me?because this yellow line is not at all in agreement with the other lines (reverse slope, I mean the yellow line has a positive slope and the other ones have negative slope).
Is there something wrong in my R code ? Is it because the Y variable is 1 vector and should be a matrix ?
Here is the
2018 Apr 07
0
Fast tau-estimator line does not appear on the plot
You need to pay attention to the documentation more closely. If you don't
know what something means, that is usually a signal that you need to study
more... in this case about the difference between an input variable and a
design (model) matrix. This is a concept from the standard linear algebra
formulation for regression equations. (Note that I have never used RobPer,
nor do I regularly
1999 Sep 17
1
Tukey's biweight
I want to estimate the center of a distribution with lots of outliers in one
tail, and thought I would use a function such as S-plus's location.m() with
psi.fun=bisquare (as per MASS 3 p. 131). However, R seems not have such a
function, so my questions are:
1) Is there an R equivalent to location.m()?
2) Would huber() give me results that are similar (i.e., close enough)?
Thanks.
2018 Mar 31
0
Fast tau-estimator line does ot appear on the plot
On 31/03/2018 11:57 AM, varin sacha via R-help wrote:
> Dear R-experts,
>
> Here below my reproducible R code. I want to add many straight lines to a plot using "abline"
> The last fit (fast Tau-estimator, color yellow) will not appear on the plot. What is going wrong ?
> Many thanks for your reply.
>
It's not quite reproducible: you forgot the line to create
2018 Mar 31
2
Fast tau-estimator line does ot appear on the plot
Dear R-experts,
Here below my reproducible R code. I want to add many straight lines to a plot using "abline"
The last fit (fast Tau-estimator, color yellow) will not appear on the plot. What is going wrong ?
Many thanks for your reply.
##########
Y=c(2,4,5,4,3,4,2,3,56,5,4,3,4,5,6,5,4,5,34,21,12,13,12,8,9,7,43,12,19,21)
2005 Mar 24
1
Robust multivariate regression with rlm
Dear Group,
I am having trouble with using rlm on multivariate data sets. When I
call rlm I get
Error in lm.wfit(x, y, w, method = "qr") :
incompatible dimensions
lm on the same data sets seem to work well (see code example). Am I
doing something wrong?
I have already browsed through the forums and google but could not find
any related discussions.
I use Windows XP and R
2010 Jun 10
1
nls model fitting errors
What am I failing to understand here?
The script below works fine if the dataset being used is
DNase1 <- DNase[ DNase$Run == 1, ] per the example given in
help(nlrob).
Obviously, I am trying to understand how to use nls and nlrob to fit
curves to data using R.
#package=DAAG
attach(codling)
plot(pobs~dose)
#next command returns 'step factor reduced below min factor
2012 Nov 22
1
help in M-estimator by R
hi guys and gals ... How are you all ...
i have to do something in robust regression by R programm , and i have some
problems as following:
*the first :*
suppose
w(r) =1/(1 r^2) and r <- c(7.01,2.07,7.061,5.607,8.502,54.909,12.222)
and i want to exclude some values from r so that (abs(r)>4.9 )...
after ,i want to used (w) to get on coefficients beta0 and beta1 (B1 <-
2007 Nov 21
1
equivalent of Matlab robustfit?
Hi,
I've been using the Matlab robustfit function for linear regressions
where I suspect some data points are outliers. Is there an equivalent
function in R?
Take care, Darren
PS, This is the Matlab help on robustfit:
>> help robustfit
ROBUSTFIT Robust linear regression
B = ROBUSTFIT(X,Y) returns the vector B of regression coefficients,
obtained by performing robust
2004 Oct 11
3
split and rlm
Hello, I'm trying to do a little rlm of some data that looks like this:
UNIT COHORT perdo adjodds
1010 96 0.39890 1.06894
1010 97 0.48113 1.57500
1010 98 0.36328 1.21498
1010 99 0.44391 1.38608
It works fine like this: rlm(perdo ~ COHORT, psi=psisquare)
But the problem is that I have about 100 UNITs, and I want to do a
2018 Apr 13
0
cvTools for 2 models not working
Dear R-experts,
I am trying to do cross-validation for different models using the cvTools package.
I can't get the CV for the "FastTau" and "hbrfit". I guess I have to write my own functions at least for hbrfit. What is going wrong with FastTau ?
Here below the reproducible example. It is a simple toy example (not my real dataset) with many warnings, what is important to
2002 Mar 19
3
plot with the axes at xlim and ylim
I'd like to make a plot with the axes drawn at the values of xlim and
ylim. The default plot draws the axes slightly outside these values. I
have been experimenting with the par(mgp) setting, but specifying
par(mgp=c(2,0,-1)) gives me
"invalid value specified for graphics parameter "mgp".".
A more complicated mathod that almost seems to work is doing the plot
with
2008 Dec 22
0
nlsrob fails with puzzling error message on input accepted by nls
I have a nonlinear model estimation problem with ~50,000 data records
and a simple 3 parameter model (logistic type - please don't tell me
that there are linear methods for such a problem). I run nls with
constraints once to get a good initial parameter guess, then try to run
nlrob to get improved estimates. The model is well-behaved for the
parameters that come from nls - no huge values, NAs
2005 Aug 23
1
Robust M-Estimator Comparison
Hello,
I'm learning about robust M-estimators right now and had settled on the
"Huber Proposal 2" as implemented in MASS, but further reading made clear,
that at least 2 further weighting functions (Hampel, Tukey bisquare) exist.
In a post from B.D. Ripley going back to 1999 I found the following quote:
>> 2) Would huber() give me results that are similar (i.e., close
2010 Oct 11
2
Nonlinear Regression Parameter Shared Across Multiple Data Sets
I'm working with 3 different data sets and applying this non-linear
regression formula to each of them.
nls(Y ~ (upper)/(1+10^(X-LOGEC50)), data=std_no_outliers,
start=list(upper=max(std_no_outliers$Y),LOGEC50=-8.5))
Previously, all of the regressions were calculated in Prism, but I'd like to
be able to automate the calculation process in a script, which is why I'm
trying to move to
2018 Apr 25
0
Zero errors : Bug in my R code ?
Dear R-experts,
I guess I have a problem with my fast function (fast tau estimator) here below. Indeed, zero errors look highly suspicious. I guess there is a bug in my R code. How could I correct my R code ?
# install.packages( "robustbase" )
# install.packages( "MASS" )
# install.packages( "quantreg" )
# install.packages( "RobPer" )
#