Displaying 15 results from an estimated 15 matches for "sabaric".
Did you mean:
sabaris
2008 Mar 30
2
convert weekly time series data to monthly
...lass(ts).
I can perform the conversions outside of R and then import but I would
rather perform all conversions within R. I have looked through the zoo and
Rmetrics packages but without success.
Any help will be appreciated.
Thanks,
Richard Saba
Department of Economics
Auburn University
Email: sabaric at auburn.edu
Phone: 334 844-2922
2008 Jan 03
1
R procedure similar to STATA heckprob?
Is anyone aware of an R procedure similar to STATA's "heckprob" procedure?
"Heckprob" fits maximum likelihood probit models correcting for sample
selection bias.
Thanks,
Richard Saba
Department of Economics
Auburn University
Email: sabaric@auburn.edu
[[alternative HTML version deleted]]
2008 Mar 21
1
tseries(arma) vs. stats(arima)
...stats" or in other time
series packages like fArima, forecast, or FinTS? They all take a "lag"
argument. I would like to have the ability to estimate models like the one
above while utilizing the "xreg" argument available in the other arima
functions .
Thanks,
Richard Saba
sabaric at auburn.edu
2007 Nov 28
2
extracting year an month from ts data set
I have an ascii data set of monthly observation starting in Jan 1946 with a
header.
hstarts
57
65
95
103
103
97
94
.
.
.
Which I read with the following code
tab6.1<-ts(read.table(fname, header=TRUE),frequency=12,start=c(1946,1))
I would like to run a time series model with dummy variables for each month.
If I had a variable which take values from 1 to 12 indicating the month I
could use
2007 Dec 05
1
Working with "ts" objects
...I found the Shumway "Time series analysis
and its applications with R Examples" website very helpful but many
practical questions involving manipulation of time series data still remain.
Any help will be appreciated.
Thanks,
Richard Saba
Department of Economics
Auburn University
Email: sabaric at auburn.edu
Phone: 334 844-2922
anscombe<-read.table(fname, header=TRUE)
names(anscombe)<-c("x1","y1","x2","y2","x3","y3","x4","y4")
reg1<-lm(y1~1 + x1, data=anscombe)
reg2<-lm(y2~1 + x2, da...
2008 Feb 12
2
Formulae for R functions
...e function help files without
success. For example I have attempted to replicate by hand the se.fit
calculation from a lm object calculated by a call to the predict function
and have not been able to reproduce the results.
Thanks,
Richard Saba
Department of Economics
Auburn University
Email: sabaric@auburn.edu
[[alternative HTML version deleted]]
2008 Apr 30
2
Bug? in summary( ) function base package
There seems to be an error in the summary() function when applied to "ts"
class objects. The results of a call to summary( ), on the R "ts" data set
USAccDeaths , reports the wrong value for Max. The value reported by the
summary function is 11320. The max( ) function returns the correct value
11317, the July 1993 value. Coercing the data to a data.frame and calling
summary
2010 Feb 07
2
Reading hierarchical data
I would like to read the following hierarchical data set. There is a family
record followed by one or more personal records.
If col. 7 is "1" it is a family record. If it is "2" it is a personal
record.
The family record is formatted as follows:
col. 1-5 family id
col. 7 "1"
col. 9 dwelling type code
The personal record is formatted as follows:
col.
2007 Oct 26
1
Newey-West and SUR regression models
Is anyone aware of a procedure to apply Newey-West corrections for
autocorrelation to a SUR regression model? The SANDWICH package seems to be
applicable only to LM or GLM models.
Thanks,
Richard Saba
Department of Economics
Auburn University
Email: sabaric at auburn.edu
2008 Jan 11
1
question about xreg of arima
Hi,
I am trying to understand exactly what xreg does in arima. The documentation
for xreg says:"xreg Optionally, a vector or matrix of external regressors,
which must have the same number of rows as x." What does this mean with
regard to the action of xreg in arima?
Apparently somehow xreg made the following two arima fit equivalent in R:
arima(x, order=c(1,1,1), xreg=1:length(x))
is
2010 Apr 08
1
reshape panel data
I have a data set with observations on 549 cities spanning an 18 year
period. However, some of cities did not report in one or more of the 18
years. I would like to implement the procedure suggested by Wooldridge
section 17.1.3 in his "Econometric analysis of cross section and panel data"
to correct for attrition. For example the table below indicates that the 3rd
and the 7th cities in
2011 Nov 07
0
vars impulse responce function output
Does anyone know if the bootstrap CI intervals generated by the irf() function (impulse response function) in the " vars" package are bias corrected?
Thanks,
Richard Saba
[[alternative HTML version deleted]]
2011 Aug 15
1
Read variable column width data
Reading data with variable column widths.
Here are several lines of a txt data set I would like to read.
The number of variables is fixed at 13 . The problem is how to read the
first variable when it can contain blank space-- for example " Alabama
(Seasonally Adjusted)" , "St. Clair", etc.
Alabama (Seasonally Adjusted) 2,168,870 2,162,604 2,122,787 1,954,895
1,956,026
2013 Feb 26
1
problem with nested loops
Each of the data sets contains monthly observations on price indices for 7 countries. I use the fitted values from reg1 in the reg2 model. The interior loop executes without error as long as I explicitly specify the data set, i.e. data=dat70. However the code fails to execute if I specify the model in the form of the commented line, i. e reg1 <-dynlm(form1,data=Dnames[j])
I get the following
2009 Mar 05
2
Overriding contributed package functions
The "tsdiag" function in the TSA package overrides the "tsdiag" function in
the "stats" package. There are a few annoying bugs in the TSA's version of
the function so I would like to use the "stats" function but still have
access to other TSA functions. I have tried using stats::tsdiag( ) but as
long as the TSA package is attached the function from