Displaying 20 results from an estimated 900 matches similar to: "weighted 2 or 3 parameter weibull estimation?"
2003 Sep 30
1
can't get names from vector in nlm calls
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 ) )
2001 Nov 13
2
how do I install add on packages?
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:
2002 Jan 04
1
plotting maps?
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
2002 May 22
5
plotting a 1000 plots?
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
2003 Oct 17
2
nlm, hessian, and derivatives in obj function?
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
2002 Mar 20
3
tex/latex output?
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
2002 May 15
1
can't get rcmd to work on windows2k....
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
2002 Feb 28
1
trying to build matrix for tests...
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
2003 Feb 16
1
__stdcall funcitons called using .C() on win32 chokes
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
2002 Jan 16
1
getting the response (dependent) from formula?
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 ),
2001 Nov 12
8
help uses only netscape?
> 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
2002 Jan 18
2
easier way to update packages i'm developing?
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
2003 Oct 06
1
getting names of p vector in nlm function...
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(
2001 Dec 08
1
plotting values "by"
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
2001 Dec 08
5
no _?
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
2003 Aug 15
1
rterm not shutting down from ESS on Win32
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
2002 May 04
1
lattice help....
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
2002 Oct 08
1
dyn.load and c-function
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(
2001 Dec 11
4
journal form?
Are these posts available in journal form so I could get only one message a
day? I really enjoy getting the lowdown on R and I'm really putting it to
work and have almost stopped using that other package whose initials are
"SAS", but I' like to unclutter the inbox. Or are these archived at the end
of the day, week, or month? I wouldn't want to go for a week without
learning
2001 Nov 27
2
overlaying qqnorm plots...
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 )