similar to: can't get names from vector in nlm calls

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(

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've been trying to figure out how to build a list of terms from a nonlinear model (terms() returns a error). I need to compute and evaluate the partial derivatives (Jacobian) for each equaiton in a set of equations. For example: > eqn <- q ~ s0 + s1 * p + s2 * f + s3 * a > sv2 <- c(d0=3,d1=4.234,d2=4,s0=-2.123,s1=0.234,s2=2.123,s3=4.234) > names( sv2 ) [1] "d0"

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

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'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:

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'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

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'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

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 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

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

> 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 am new to SAMBA and I have what I'm not even sure is an issue. I am aware of the difference in time stamps between *nix and Windows. What I don't understand is this: I used touch to modify the time stamps of a large number of files on the file server from the server side to match the time in the file name. They were video files from my Digital Video Cam and the import program used

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 got a c function in a dll and have been able to call it from R. The function is : void __stdcall testfunc2( unsigned long *a, double *b, unsigned long *c, double *result ) { *result = 0; *result = (*a) * (*b) * (*c); } and the R code to call the function is as follows: dyn.load("c:/openfvs/biometrics/debug/biometrics.dll" ) tf2 <- function(