Displaying 20 results from an estimated 700 matches similar to: "Exception in NeweyWest - Pre-Whitening necessary?"
2010 Jun 27
1
NeweyWest
I want to calculate Newey West robust standard error using NeweyWest. Comparing the results to what I get in STATA, in order to get the same results in I need to specify "prewhite=0". Can someone explain what this prewhite command means?
Thanks
[[alternative HTML version deleted]]
2011 Aug 27
3
Exception while using NeweyWest function with doMC
Dear R users,
I am using R right now for a simulation of a model that needs a lot of
memory. Therefore I use the *bigmemory* package and - to make it faster -
the *doMC* package. See my code posted on http://pastebin.com/dFRGdNrG
Now, if I use the foreach loop with the addon %do% (for sequential run) I
have no problems at all - only here and there some singularities in
regressor matrices which
2011 Jul 25
1
biglm() and NeweyWest()
Dear all,
I am working on a large dataset and need to use biglm() to perform OLS
regressions. I have detected significant ARCH effects which I try to account
for using the Newey-West correction.
So far, I have worked with NeweyWest() in the sandwich package. NeweyWest()
however seems to be unable to handle an object of class "biglm".
Looking into the code, I figured out that
2011 Sep 28
1
Robust covariance matrix with NeweyWest()
Dear R-users,
I would like to compute a robust covariance matrix of two series of realizations of random variables:
###Begin Example###
data <- cbind(rnorm(100), rnorm(100))
model <- lm(data ~ 1)
vcov(model)
library(sandwich)
NeweyWest(model) #produces an error
###End Example###
NeweyWest() produces an error but sandwich(), vcovHAC(), kernHAC, weave(),... do not produce any errors. It
2010 Sep 23
1
Newey West and Singular Matrix + library(sandwich)
thank you, achim. I will try chol2inv.
sandwich is a very nice package, but let me make some short
suggestions. I am not a good econometrician, so I do not know what
prewhitening is, and the vignette did not explain it. "?coeftest" did
not work after I loaded the library. automatic bandwidth selection
can be a good thing, but is not always.
as to my own little function, I like the
2010 Sep 22
1
Newey West and Singular Matrix
dear R experts: ?I am writing my own little newey-west standard error
function, with heteroskedasticity and arbitrary x period
autocorrelation corrections. ?including my function in this post here
may help others searching for something similar. it is working quite
well, except on occasion, it complains that
Error in solve.default(crossprod(x.na.omitted, x.na.omitted)) :
system is
2009 Feb 16
2
Whitening Time Series
Hi R users,
I am doing cross correlation analysis on 2 time series (call them y-series
and x-series) where I need the use the model developed on the x-series to
prewhiten the yseries.. Can someone point me to a function/filter in R that
would allow me to do that?
Thanks in advance for any help!
--
View this message in context:
2002 Nov 26
1
floor curve question (whitening filter)
<!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 4.0 Transitional//EN">
<HEAD><META http-equiv=Content-Type content="text/html; charset=iso-8859-1"></HEAD>
<BODY bgColor=#ffffff><FONT face="verdana,arial" size="2">
Hi there,</P>
one step in the encoding process (if I got it right), would becomputing a "floor
2007 May 14
1
Free Colgate Max Fresh Whitening Toothpaste
http://www.colgate.toothpaste-sample.com
Get New Colgate Max Fresh Whitening Toothpaste. Cool Mint Flaor with
Breath Strips.
--
Posted via http://www.ruby-forum.com/.
--~--~---------~--~----~------------~-------~--~----~
You received this message because you are subscribed to the Google Groups "Ruby on Rails: Talk" group.
To post to this group, send email to
2009 Dec 10
2
Problem with coeftest using Newey West estimator
Hi,
I want to calculate the t- and p-values for a linear model using the Newey West estimator.
I tried this Code and it usually worked just fine:
> oberlm <- lm(DYH ~ BIP + Infl + EOil, data=HU_H)
> coeftest(oberlm, NeweyWest(oberlm, lag=2))
t test of coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) 0.1509950 0.0743832 2.0300 0.179486
BIP
2007 Feb 26
1
Partial whitening of time series?
I have a time series with a one year lag, ar=0.5. The series has some
interesting events that disappear when the series is whitened (i.e.,
fitting an AR process and looking at the residuals). I'd like to remove
the autocorrelation in stages to see the effect on the time series. Is
there a way to specify the autocorrelation term while fitting an AR
process?
For instance, given the following:
2006 Mar 11
1
Analyze Fmri package
Hi,
Can you please help me?
Whenever I try to run f.ica.fmri.gui(), R crashes.
In the GUI, I selected an *.img file and I checked the 'Create mask' box and
I typed a name for the object and then I press 'Start' button. Then in the
R Console window it said 'Reading...Whitening', but after the Whitening
message an error window popped up and then R crashed.
Did I do something
2010 Oct 13
1
robust standard errors for panel data
Hi,
I would like to estimate a panel model (small N large T, fixed effects),
but would need "robust" standard errors for that. In particular, I am
worried about potential serial correlation for a given individual (not so
much about correlation in the cross section).
>From the documentation, it looks as if the vcovHC that comes with plm
does not seem to do autocorrelation, and the
2011 Apr 07
1
transform() on selective names. Is it possible?
Hi all,
I am whitening my data:
# code begins
N <- 300
M <- 2
x <- matrix(data=rnorm(N*M, 0, 3)-10, ncol=M, nrow=N)
y <- matrix(c(1,-2,-2,1), ncol=M, nrow=M)
z <- data.frame(x %*% y)
colnames(z) <- c('x','y')
par(mfrow=c(1,3))
plot(z, pch=5, col="blue")
whiten <- function(x) { (x-mean(x))/sd(x) }
zz <- transform(z, x=whiten(x), y=whiten(y))
2008 Nov 20
0
A Problem while Calculating Newey-West HAC
Hi,
Does anyone read Verbeek's "A Guide to Modern Econometrics"? In its Section
4.11, how does the last two equations' HAC calculate? I've tried several
groups of parameters in sandwich::NeweyWest, but I still cannot get the same
result. I've tried lag=2 and lag=3, as long as prewhite=FALSE and
prewhite=TRUE yet, but...
Sincerely
Hsiao-nan Cheung
2010 May 02
1
question about 2SLS
Hi All,
I am using R 2.11.0 on a Ubuntu machine. I estimated a model using "tsls"
from the package "sem". Is there a way to get Newey West standard errors for
the parameter estimates?
When estimating the model by OLS, I used "NeweyWest" from the package
"sandwich" to get HAC standard errors. But, I am not able to use the same
method with the results of the
2009 Mar 10
1
HAC corrected standard errors
Hi,
I have a simple linear regression for which I want to obtain HAC corrected
standard errors, since I have significant serial/auto correlation in my
residuals, and also potential heteroskedasticity.
Would anyone be able to direct me to the function that implements this in R?
It's a basic question and I'm sure I'm missing something obvious here. I
looked up this post:
2014 Aug 11
1
[PATCH 1/2] rngd: add udev rule to source from hwrng if virtio-rng present
On Mon, Aug 11, 2014 at 12:45:27PM +0530, Amit Shah wrote:
>
> What's the suggested value for rng->quality, though, for virtio-rng
> that I can use to ensure the kthread starts?
>
> Should I use the 700 (70%) as proposed in the original patchset? I'm
> not exactly sure how that value will be used as well..
There is no such thing as a suggested value, every number
2014 Aug 11
1
[PATCH 1/2] rngd: add udev rule to source from hwrng if virtio-rng present
On Mon, Aug 11, 2014 at 12:45:27PM +0530, Amit Shah wrote:
>
> What's the suggested value for rng->quality, though, for virtio-rng
> that I can use to ensure the kthread starts?
>
> Should I use the 700 (70%) as proposed in the original patchset? I'm
> not exactly sure how that value will be used as well..
There is no such thing as a suggested value, every number
2008 May 22
1
How to account for autoregressive terms?
Hi,
how to estimate a the following model in R:
y(t)=beta0+beta1*x1(t)+beta2*x2(t)+...+beta5*x5(t)+beta6*y(t-1)+beta7*y(t-2)+beta8*y(t-3)
1) using "lm" :
dates <- as.Date(data.df[,1])
selection<-which(dates>=as.Date("1986-1-1") & dates<=as.Date("2007-12-31"))
dep <- ts(data.df[selection,c("dep")])
indep.ret1