R-beta: problem with locfit

Clive Loader clive at bell-labs.com
Thu Mar 26 16:27:17 CET 1998


guido at sirio.stat.unipd.it wrote:
> 
> I installed the locfit package under Linux (gcc 2.7.2). Installation
> was ok but
...
> > m <- locfit(y~x,cens=cens,family="gamma")
> /usr/local/src/R-0.61.1/bin/R.binary: can't resolve symbol 'igamma'

> Perhaps, a better possibility is to define igamma and ibeta using
> the R pgamma and pbeta function. However, I have not tried it.

placing (in the locfit/src/family.c file)

extern double pgamma(), pbeta();
#define igamma(x,p) pgamma(x,p,1.0)
#define ibeta(x,a,b) pbeta(x,a,b)

seems to solve the problem. (Big plus of R over S: pgamma() etc. can be
called directly,
instead of cribbing source code!).


A second problem is that the R (0.61.1) model.extract() function differs
from the
S version, and the cens argument won't be interpreted properly. Adding
the following
modified model.extract() to the locfit/R/lfrsup.r file should work...

model.extract <-
function (frame, component)
{
        component <- as.character(substitute(component))
        rval <- switch(component, response = model.response(frame),
                offset = model.offset(frame), weights =
frame$"(weights)",
                start = frame$"(start)", cens=frame$"(cens)",
                base=frame$"(base)")
        if (length(rval) == nrow(frame))
                names(rval) <- attr(frame, "row.names")
        else if (is.matrix(rval) && nrow(rval) == nrow(frame)) {
                t1 <- dimnames(rval)
                dimnames(rval) <- list(attr(frame, "row.names"),
                        t1[[2]])
        }
        return(rval)
}


Regards,
Clive.
-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-
r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html
Send "info", "help", or "[un]subscribe"
(in the "body", not the subject !)  To: r-help-request at stat.math.ethz.ch
_._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._



More information about the R-help mailing list