similar to: simultaneous equation estimation?

I'm trying to build a flexible OLS/SUR/2SLS/3SLS package and I'm having trouble getting some information out of a formula. For example.... # set up the system of equations demand <- q ~ p + d supply <- q ~ p + f + a inst <- ~ d + f + a systemeq <- list( demand, supply ) ...blah, blah, blah... # get the number if instruments... # or build the "bigX"... #I

I know this topic has had plenty of discussion in the last couple of days, but.... I've been trying to compare the effects of different fitted methods for systems of equations (OLS, SUR, 2SLS, 3SLS ) and would like to compare the residual plots (easy) and the qqnorm/qqplot of the fits for the different fitted methdos. For example, qqnorm( residuals( lm( q ~ p + f + a ) ) ) par( new = TRUE )

I've figured out how to use optim (barely) to estimate 2 parameter = weibull distributions. I can't get over how easy this is. What I need to = do is use a weight in the observations..... For example,=20 the tree diameters and weights are are=20 4.70 , 100 6.00, 98 7.10, 75.0 8.10, 86.3 8.60, 80.456 8.90, 20.5 9.50, 16.6 11.40, 12.657 11.80, 12.47 14.50,

I've been trying to figure out how to get the names of the parameter vector variables when inside the function that nlm calls to return the objective function value: knls <- function( theta, eqns, data, fitmethod="OLS", instr=NULL, S=NULL ) { ## print( names( theta ) ) # returns NULL ## get the values of the parameters for( i in 1:length( theta ) )

I thought I had plotted a map of oregon and washington and overlayed some points using R. Is that possibe since I can't seem to find the source file I used or commands for plotting maps and points? If not, sorry for the bother... Thanks, Jeff. Jeff D. Hamann Hamann, Donald & Associates, Inc. PO Box 1421 Corvallis, Oregon USA 97339-1421 Bus. 541-753-7333 Cell. 541-740-5988 jeff_hamann

I'm having some trouble installing the add on packages. I would like to install the RMySQL_0.4-3.tar.gz add on package. I've downloaded it and unzipped the package. The INSTALL file says R CMD INSTALL [-l<alternate-dir>] RMySQL_<version>.tar.gz does this mean in the R system enter the following? > INSTALL -l /root/RMySQL_0.4-3.tar.gz (?) when I do this I get... Error:

Dear R programming folks: I'm trying to finish off a package for non-linear simultaneous system estimation and I've been trying to figure out how to get the names of the parameter vector variables when inside the function that nlm calls to return the objective function value: knls <- function( theta, eqns, data, fitmethod="OLS", instr=NULL, S=NULL ) { ## print(

Is it possible to write the output/results (redirect) to a latex file? Jeff. Jeff D. Hamann Hamann, Donald & Associates, Inc. PO Box 1421 Corvallis, Oregon USA 97339-1421 Bus. 541-753-7333 Cell. 541-740-5988 jeff_hamann at hamanndonald.com www.hamanndonald.com -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read

I'm trying to make sure my R package will install properly (build help files, etc) and when I try to run the rcmd check (or any rcmd option) I get the following error. Does dcf.pm stand for perl module? If not, why do I get this message? C:\systemfit\r_code>rcmd check --help Can't locate R/Dcf.pm in @INC (@INC contains: c \PROGRA~1\r\rw1041\share\perl; /usr/lib/perl5/5.6.1/cygwin-multi

Okay, I stumped. I'm trying to build a matrix or data.frame that consists of values by group. For example, resid, group 574 0.536196370873122 E 575 0.93191093696265 E 576 -1.27116094894117 A 577 -1.05159666055027 A 578 0.398273458301184 D 579 0.586373126102573 E 580 0.806575304513802 E 581 0.897258393118338 E 582 0.342393324256838 D 583 0.781245511773227 E I

I'm trying to build a vector (fitting systems of equations) and I'm having a little trouble getting the response variables from a list of formulas... bdgmodel <- lbdg ~ d1sqr + ld1 + lhg hgmodel <- lhg ~ ht2 + lht + whc inst <- ~ d1sqr + ld1 + ht2 + lht + whc systemeq <- list( bdgmodel, hgmodel ) # manually generate the y matrix... y <- cbind( matrix( lbdg ),

I would like to plot lots of plots (2x2 plots) of a set of data for many sets. I can do it just fine for one set of data, but I would like to generate a small mountain of plots for others to examine. Basically, I would like to generate the same par( mcol=c(2,2) ) plot for about 1000 datasets. I've been iterating through them in the plot window, but I would like to send all the plots to a eps

I would like to produce a set of values, on the same chart, from an sql table that is structured like... species dbh ht expf DF 1.2 8.9 10.0 DF 2.4 17.3 12.4342 DF 3.1 20.9 56.76 PP 2.3 16.9 100.0 PP 12.8 97.3 40.3 PP 8.2 63.0 98.34 . . . SS blah, blah, blah... is it possible to, using a single command in the plot command to plot the different groups on the same plot or will I have to iterate

> I don't have netscape installed on my machine (FreeBSD 4.4) and when I type > the command help.start(), I get > > > netscape: not found > > error. I do have mozilla installed. Is this a defect? can I run any web > browser? > > jeff. > > > ----- Original Message ----- > From: <Majordomo at stat.math.ethz.ch> > To: <jeff_hamann at

I see the uderscore "_" is not allowed in R. This make R a real drag when trying to use with SQL packages and c code. Why is the underscore not allowed and will it be allowed in a future release? Jeff. Jeff D. Hamann Hamann, Donald and Associates PO Box 1421 Corvallis, Oregon USA 97339-1421 Bus. 541-753-7333 Cell. 541-740-5988 jeff_hamann at hamanndonald.com www.hamanndonald.com

I've been working on a new package and I have a few questions regarding the behaviour of the nlm function. I've been (for better or worse) using the nlm function to fit a linear model without suppling the hessian or gradient attributes in the objective function. I'm curious as to why the nlm requires 31 iterations (for the linear model), and then it doesn't work when I try to add

I'm frightfully new at developing r packages. I'm developing an econometrics package, and each time I update my beta package, I close R, delete the package dir in the library dir, rebuild the zip file, rerun R, reinstall the package and run the functions. Is there a more efficient way? Jeff. Jeff D. Hamann Hamann, Donald & Associates, Inc. PO Box 1421 Corvallis, Oregon USA 97339-1421

I've been having problems with Rterm.exe not shutting down when I exit an R (1.7.0 and 1.7.1) session from within emacs when using ESS. I've just upgraded to 5.1.24 and still have the same problems. I'm running ntemacs and winxp. I don't recall having these troubles with version 1.6.2 of R. Jeff. --- Jeff D. Hamann Hamann, Donald and Associates, Inc. PO Box 1421 Corvallis, Oregon

I'm trying to create some trellis plots (like the last plot in the lattice demo). When I run the demo( "lattice" ), the demo runs just fine, but when I cut and paste some or all of the demo code into a r file of mine (and nothing else), I can't see anything on a device and it doesn't look like any of the demo is being run. Has anyone else seen this? Jeff. Jeff D. Hamann

I've developed a dll (win32) of a bunch of functions prototypes as: /* for building on *nix */ #ifndef WIN32 #define __stdcall /*nothing*/ #endif void __stdcall dbh_2_height( const unsigned long *func_idx, const unsigned long *metric, const unsigned long *species, const double *dbh, double *pred_height ); and have created a R interface