search for: ibeta

...t;=med.y,1,0) > y <- cens * y.compl + (1-cens)*med.y > library(locfit) > m <- locfit(y~x,cens=cens,family="gamma") /usr/local/src/R-0.61.1/bin/R.binary: can't resolve symbol 'igamma' I took a look to the original code, not the R/S one. There igamma (and also ibeta which the r distribution misses too) are defined in the file random.c. I inserted them in the locfit/src-c/random.c (see below) and then things seem to go well (= the previous example works and give sensible results). Perhaps, a better possibility is to define igamma and ibeta using the R pgamma a...

I have a parameter vector, ibeta, and a corresponding vector of loglikelihoods, llbet. If llbet contains no NAs or Inf's then I can extract the best parameter by index <- order(-llbet)[1] beta <- ibeta[index] or similar. The argument na.last of order() allows me to fix this up even if llbet contains some NAs which I...

...d pasted directly from the web resource - the only change I have made is to fill in my values for T, n, low/high se/sp and the alpha beta for the distributions) prevalence.bayes<-function(theta,T,n,lowse=0.5,highse=1.0, sea=1,seb=1,lowsp=0.5,highsp=1.0,spa=1,spb=1,ngrid=20,coverage=0.95) { ibeta<-function(x,a,b) { pbeta(x,a,b)*beta(a,b) } ntheta<-length(theta) bin.width<-(theta[ntheta]-theta[1])/(ntheta-1) theta<-theta[1]+bin.width*(0:(ntheta-1)) integrand<-array(0,c(ntheta,ngrid,ngrid)) h1<-(highse-lowse)/ngrid h2<-(highsp-lowsp)/ngrid...

...with many contributors. Type "?contributors" for a list. Type "demo()" for some demos, "help()" for on-line help, or "help.start()" for a HTML browser interface to help. Program received signal SIGFPE, Arithmetic exception. 0xa0bac in machar (ibeta=0xefffef7c, it=0xefffef78, irnd=0xefffef74, ngrd=0x200, machep=0xefffef6c, negep=0x0, iexp=0xefffef64, minexp=0xefffef60, maxexp=0xefffef5c, eps=0xefffef50, epsneg=0xefffef48, xmin=0xefffef40, xmax=0xefffef38) at machar.c:204 204 z = y * y; (gdb) where #0 0xa0bac...

...ct -O2 -I/usr/local/R-0.61.1/include -c density.c -o density.o cc -ieee_with_inexact -O2 -I/usr/local/R-0.61.1/include -c kappa0.c -o kappa0.o cc -ieee_with_inexact -O2 -I/usr/local/R-0.61.1/include -c simul.c -o simul.o ld: Warning: Unresolved: pnorm exp log sqrt fabs printf igamma lgamma ibeta pow calloc floor __remq __divq sin cos pchisq atan2 runif pf dchisq df fprintf rnorm rpois rexp __exc_add_pc_range_table __exc_remove_pc_range_table __exc_add_gp_range __exc_remove_gp_range R data DONE Installing documentation of package `locfit' ... >>> Building help pages for packa...

...se /* not NLS */ > < #define _(String) (String) > < #define N_(String) String > < #endif > < > --- >> char *EncodeString(SEXP, int, int, int); > 950,957d874 > < /* structure for caching machine accuracy values */ > < typedef struct { > < int ibeta, it, irnd, ngrd, machep, negep, iexp, minexp, maxexp; > < double eps, epsneg, xmin, xmax; > < } AccuracyInfo; > < > < extern AccuracyInfo R_AccuracyInfo; > < > > ______________________________________________ > R-devel at r-project.org mailing list > ht...

...#define N_(String) gettext_noop (String) < #else /* not NLS */ < #define _(String) (String) < #define N_(String) String < #endif < --- > char *EncodeString(SEXP, int, int, int); 950,957d874 < /* structure for caching machine accuracy values */ < typedef struct { < int ibeta, it, irnd, ngrd, machep, negep, iexp, minexp, maxexp; < double eps, epsneg, xmin, xmax; < } AccuracyInfo; < < extern AccuracyInfo R_AccuracyInfo; <