Re: [R] Modifying Internal C Files

2005-11-08 Thread Ken Kelley
, Ken Kelley wrote: Hi Duncan and others. Thanks for your insight. Actually, I did change the function name in the pnbeta2.c code (I should have said so). When I check the package I get no errors. I build the package and all seems well. When I install and then load the package in R, I get

[R] Modifying Internal C Files

2005-11-07 Thread Ken Kelley
and loaded it it, the function doesn't work (it isn't even recognized). I have no idea why this is failing. Any information to help me figure out what I need to do to modify internal C code generally, or specifically as it applies to this scenario would be most helpful. Thanks, Ken -- Ken Kelley

Re: [R] Modifying Internal C Files

2005-11-07 Thread Ken Kelley
pnbeta2 internal function doesn't seem to be there. What I want to do is simple (modify an interal C file), but it is proving to be quite difficult for me to implement. Thanks for any thoughts, Ken Duncan Murdoch wrote: On 11/7/2005 5:17 PM, Ken Kelley wrote: Hi All. I want to tweak a few

Re: [R] Doubly Non-Central F-Distribution

2005-11-01 Thread Ken Kelley
-project.org/posting-guide.html -- Ken Kelley, Ph.D. Inquiry Methodology Program Indiana University 201 North Rose Avenue, Room 4004 Bloomington, Indiana 47405 http://www.indiana.edu/~kenkel __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch

[R] Problems with pf() with certain noncentral values/degrees of freedom combinations

2005-10-24 Thread Ken Kelley
until then? Thanks, Ken Kelley # Begin example code X - seq(10, 600, 10) # Gets stuck at .99135 round(pf(X, 10, 1000, 225), 5) round(pf(X, 10, 200, 225), 5) round(pf(X, 5, 1000, 225), 5) round(pf(X, 5, 200, 225), 5) round(pf(X, 1, 1000, 225), 5) round(pf(X, 1, 200

[R] Sometimes having problems finding a minimum using optim(), optimize(), and nlm() (while searching for noncentral F parameters)

2005-10-10 Thread Ken Kelley
, or even a reasonable, way? Are there other functions I'm missing that would be more appropriate given what I'm trying to do? Any help would most certainly be appreciated. Thanks, Ken -- Ken Kelley, Ph.D. Inquiry Methodology Program Indiana University 201 North Rose Avenue, Room 4004 Bloomington

[R] Returning singular nlme objects.

2004-08-02 Thread Ken Kelley
Hi everyone. I'm working with nlme and I have a question regarding nlme fits that fail because of singularity issues. Specifically, there a way to return an nlme object when the estimation process runs into a singular matrix? For example, can the results up to the point of an error such as Error

[R] Individual log likelihoods of nlsList objects.

2004-07-01 Thread Ken Kelley
Hello all. I was wondering if the logLike.nls() and logLike.nlme() functions are still being used. Neither function seems to be available in the most recent release of R (1.9.1). The following is contained in the help file for logLik(): classes which already have methods for this function

[R] Ignoring Errors in Simulations

2003-11-05 Thread Ken Kelley
Hello all. I'm doing a simulation study and every so often I get an error that stops the simulation. I would like to ignore the errors *and* identify the particular iterations where they occurred. I have tried: options(error = expression(NULL)) which I thought would ignore the error, but the

[R] Generating Data Sets -- Each with a Different Name

2003-10-23 Thread Ken Kelley
Hello all. I was wondering if anyone had insight into how I might generate a large number of data sets/vectors, where each of the data sets/vectors have a different name? For example, suppose I wanted to generate N data sets/vectors, each with a different name such as: Random.Data.1 -

[R] Denominator Degrees of Freedom in lme() -- Adjusting and Understanding Them

2003-10-21 Thread Ken Kelley
Hello all. I was wondering if there is any way to adjust the denominator degrees of freedom in lme(). It seems to me that there is only one method that can be used. As has been pointed out previously on the list, the denominator degrees of freedom given by lme() do not match those given by SAS