similar to: help with memory problem in SystemFit

Hi - I encounter two problems with SystemFit. I have a matrix of 20 variables (380000 observations each). I am trying to fit using "2SLS" in Systemfit because I want to be able to simulate new observations with the resulting model. The first problem I found is that it ran out of memory given that I have 2GB of RAM. How can I circumvent this? Can I do some mixing or bootstrapping with

Hello !   I’m trying to estimate a system of equation (demand and supply) using the systemfit package.  My program is:   library(systemfit) demand <- tsyud ~ tsyud1 + tsucp + tspo + tssn supply <- tscn ~ tsyn + tsqn + tsksn + tsucp system <- list(demand=eqdemand, learning = eqsupply) labels <- list(demand="eqdemand", learning="eqsupply") inst <- ~ tsupp1 + tsupp2

Hi,I have a fundamental questions that I'm a bit confused. If any guru from this circle could help me out, I would really appreciate.I have a system of equations in which some of the endogs appear on right hand sides of some equations. To solve this, one needs a technique like 2SLS or FIML to circumvent inconsistency of the estimated coefficients. My question is that if I apply the nonlinear

Hi, I am trying to obtain the heteroskedasticity consitent standard errors (HCSE) after apply 2SLS. I obtain 2SLS through tsls from package sem or systemfit: #### tsls #### library (sem) Reg2SLS <-tsls(LnP~Sc+Ag+Ag2+Var+R+D,~I2+Ag+Ag2+Var+R+D) summary (Reg2SLS) #### systemfit #### library (systemfit) RS <- LnP~Sc+Ag+Ag2+Var+R+D Inst <- ~I2+Ag+Ag2+Var+R+D labels

The 'systemfit' function in systemfit 0.7-3 CRAN package seems to have a problem with formulas that contain transformed (eg. log) variables. If I have my data in a data frame, apparently systemfit doesn't "pass" the information of where the variables should be taken to the transforming function. I'm not entirely sure if this is a bug or just a limitation, I was just

Hi, I am using tsls function from sem package to estimate a model which includes large number of data. Among its predictors, it has a nominal data which has about 10 possible values. So I expand this parameter into 9-binary-value predictors with the coefficient of base value equals 0. I also have another continuous predictor. The problem is that, whenever I run the tsls, I will get 'System

Dear R users, the systemfit package contains functions for fitting systems of simultaneous equations by various estimation methods (e.g. OLS, SUR, 2SLS, 3SLS). Currently version 0.8 of systemfit is available on CRAN. However, shortly we will upload version 1.0, which is NOT BACKWARD COMPATIBLE. The changes that broke backward compatibility were necessary to make systemfit() more similar to

Dear R users, the systemfit package contains functions for fitting systems of simultaneous equations by various estimation methods (e.g. OLS, SUR, 2SLS, 3SLS). Currently version 0.8 of systemfit is available on CRAN. However, shortly we will upload version 1.0, which is NOT BACKWARD COMPATIBLE. The changes that broke backward compatibility were necessary to make systemfit() more similar to

Dear R users, The authors of the systemfit package have released a new version of this package with substantial enhancements. The systemfit package contains functions for fitting simultaneous systems of linear equations using Ordinary Least Squares (OLS), Weighted Least Squares (WLS), Seemingly Unrelated Regressions (SUR), Two-Stage Least Squares (2SLS), Weighted Two-Stage Least Squares

Dear R users, The authors of the systemfit package have released a new version of this package with substantial enhancements. The systemfit package contains functions for fitting simultaneous systems of linear equations using Ordinary Least Squares (OLS), Weighted Least Squares (WLS), Seemingly Unrelated Regressions (SUR), Two-Stage Least Squares (2SLS), Weighted Two-Stage Least Squares

Hi, can anybody tell me how the Hausman test for endogenty works? I have a simulated model with three correlated predictors (X1-X3). I also have an instrument W for X1 Now I want to test for endogeneity of X1 (i.e., when I omit X2 and X3 from the equation). My current approach: library(systemfit) fit2sls <- systemfit(Y~X1,data=data,method="2SLS",inst=~W) fitOLS <-

A copy of this question can be found on Cross Validated: https://stats.stackexchange.com/questions/645610 I am estimating a system of seemingly unrelated regressions (SUR) with `gmm::sysGmm` in R. Each of the equations has one unique regressor and one common regressor. The common regressor is a dummy variable indicating the last observation (n-1 zeros followed by 1). I impose a restriction that

Hi - I use stepAIC to automatically select the model. The stepAIC was applied on polr as follow:objPolr <- polr(formula=myformula, data=dat, method=METHOD);objPolr.step <- stepAIC(objPolr, trace=T);Then R complaints that it doesn't know about 'dat' when it executes the second line. Below is the exact error that I got when executing the stepAIC line above:Error in eval(expr,

I might be sorry for asking this question :-) I have two equations and I tried to estimate them individually with "lm" and "gls", and then in a system (using systemfit) with "OLS", "WLS" and "SUR". Quite surprisingly (for myself at least) the results are identical to the last digit. Could someone (please!) give a hint as to what am I

Hi, I have a model which tries to fit a set of data with 10-level ordered responses. Somehow, in my data, the majority of the observations are from level 6-10 and leave only about 1-5% of total observations contributed to level 1-10. As a result, my model tends to perform badly on points that have lower level than 6. I would like to ask if there's any way to circumvent this problem or not.

Hi all - I'm on R 2.5.1 for XP. in the systemfit package, the summary is set to print the McElroy's measure of fit unless it's NULL. When the option saveMemory = TRUE, the McElroy isn't included, instead it defaults to NA. Thus I am unable to use summary.systemfit. > library(systemfit) > example(systemfit) > surfit2 <-

I am having problem using output of lm() function for further analysing using systemfit package. Basicaly, the problem s following - I generate several formulas using lm() > fo1 <- lm(r98[,2] ~ f98[,1] + f98[,2] + ... + f98[,43]) > fo2 <- lm(r98[,1] ~ f98[,1] + f98[,2] + ... + f98[,43]) and than I want to estimate a general model using package systemfit. > fitsur <-

Hello. I am still a newbie in R. Excuse me if I am asking something obvious. My efforts to get an answer through browsing the mailing archives failed. I want to perform an augmented Dickey-Fuller test and to obtain AIC and BIC and to be able to impose some linear restrictions on the ADF regression so as to decide the correct order of autoregression. However I could find no obvious way to impose

Hi, I added a new function "predict.systemfit" to our package "systemfit" to make it closer to other packages (e.g. lm). Now "R CMD check" complains that the generic function "predict" has only the argument "object", while our function "predict.systemfit" has more arguments. However, the function "predict.lm" has also more

A friend of mine sent me below so I am posting below. If it is not enough information, please just let me know and I will tell him. Thanks. -----Original Message----- Sent: Friday, April 13, 2007 3:58 PM To: Leeds, Mark (IED) Subject: R question I am using the "predict" function after I have done a simultaneous estimation of a system using "systemfit". fitsur is the output