153 questions
2
votes
0
answers
89
views
Quantify sum of squared differences in a robust multilinear regression in R
I would like to obtain estimates of the variance explained by each predictor in multiple regression using robust linear regression (for instance with the R function lmrob from robustbase R package or ...
0
votes
0
answers
160
views
Calculation of robust standard errors with gam model
I would like to run a generalised additive model (GAM) regression with some clustered survey data, where I fit a spline for one of the variables. I would like to estimate robust standard errors, ...
0
votes
1
answer
121
views
Error using mcp2atm function for robust pairwise comparisons - undefined columns selected
I used a robust two-way ANOVA using the WRS2 package (R v.4.4.1), using the following code as suggested by Mair and Wilcox (2020) - the full dataset is on the bottom of this message:
library(WRS2)
...
0
votes
1
answer
36
views
How can I store list of t1waybt ANOVAs (WRS2 package) into a data frame (R)
I have computed 12 ANOVAs using WRS2 package, with the t1waybt function. I would like to be able to put all of these results into a data frame so I can view/print them all.
I have tried the following ...
1
vote
1
answer
2k
views
Difference between rlm() and lm_robust
Due to outliers, I would like to use a robust regression method instead of lm().
I can't decide whether to use lm_robust (estimatr package) or rlm (MASS package).
Are there mathematical differences ...
0
votes
0
answers
528
views
Replacing Standard Errors in a GLM Model in R [Follow-up]
My post is a follow-up question to this one: Replacing Standard Errors in a Reg Model in R
I am in search of a way to directly replace the standard errors in a regression model with my own standard ...
0
votes
0
answers
71
views
How can I compare the results of two lmRob() models?
I am looking for a way to compare the results from two lmRob() functions statistically.
Here is an explanation of what I am trying to do: My professor wants me to compare the results of two ANCOVAs (...
1
vote
2
answers
60
views
apply a function to a robust data frame with similar variables
Perhaps someone could help me to elaborate better on the title of this post.
I am trying to apply the exact same formula to several sites (Australia, Mexico and France). Each site has the same ...
1
vote
0
answers
19
views
How to flag Outlyingness using Stanhel-Donoho Estimator in R
Please, I need help to successfully run the Stanhel and Donoho Estimator on my data to flag outlyingness. I have a four column data with the summary presented below:
Data Summary
I did some study on ...
1
vote
2
answers
485
views
What the purpose of 'bufp += nread;' in the Robust reading rio_read() function?
I'm reading the code below in Chapter 10 in the book called Computer System: A Programmer's Perspective. I understand all the code below except the line: bufp += nread; I found this line is pointless, ...
1
vote
0
answers
26
views
Lmrob returning infinity for t-value of the parameter incorrectly
I have time series with outliers, and I want to get a robust estimate of the trend and its t-value. For that I use lmrob, which works fine in most cases. However, in some cases, it returns -inf for ...
0
votes
2
answers
365
views
Error in round(x[[y]], digits = digits) : non-numeric argument to mathematical function
I was trying to fit a linear regression model using the lm function on a dataset (in R):
model <-lm(DLP~gender+Kvp120+mAs+length+gantry+device_age+detectors+scanner2,data = data)
and here were the ...
1
vote
1
answer
51
views
Variable in a function is not used as the value
I am using the package robust.arima in R, which works fine when I call it in a script. However, I want to organize my files and therefore call robust arima in a function. Here all of a sudden the ...
1
vote
0
answers
203
views
Manually calculate robust standard errors of 2SLS regression
I used the "systemfit" function in R to estimate a 2SLS model as it allows to specify the first and second stage separately which is important for my estimation. I need robust standard ...
1
vote
0
answers
111
views
automatic selection of lag order in robust arima R always returns maximum number of AR lags
I am using library robustarima, which has an option auto.ar=TRUE. According to the description and the book (https://onlinelibrary-wiley-com.eur.idm.oclc.org/doi/pdf/10.1002/9781119214656 chapter 8 ...