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