similar to: R-beta: problem with locfit

I was wondering if someone could help me with a problem getting locfit running. I am running Digital Unix 3.2D on an alphastation, with R 0.61.1 I downloaded the CRAN distribution of locfit. I had no problems installing it with R INSTALL (I installed it into my library rather than the default base library; the output is attached at the end of this email message). I then tried to test it with the

Dear R-list, I try to transform a mathematica script to R. #######relevant part of the Mathematica script (* p_sv *) dd = NN (DsD - DD^2); lownum = NN (L-DD)^2; upnum = NN (H-DD)^2; low = lownum/(2s^2); up = upnum/(2s^2); psv = NIntegrate[1/(s^NN) Exp[-dd/(2s^2)] (Gamma[1/2,0,up] + Gamma[1/2,0,low]),{s,sL,sH}, MinRecursion->3]; PSV = psv/Sqrt[2NN]; Print["------------- Results

Hello all, In the course of trying to implement the CDF of an InverseGammaDistribution, I have run across the need for an igamma() function. Several others have needed this function but the answers I have found so far are not totally clear to me. I'm writing for three reasons: 1) to present a small error in the docs 2) to clarify the approach we are expected to take 3) to request,for the

Hello I am trying to evaluate an Incomplete gamma function in R. Library Zipfr gives the Igamma function. From Mathematica, I have: "Gamma[a, z] is the incomplete gamma function." In[16]: Gamma[9,11.1] Out[16]: 9000.5 Trying the same in R, I get > Igamma(9,11.1) [1] 31319.5 OR > Igamma(11.1,9) [1] 1300998 I know I have to understand the theory and the math behind it rather

Hello, fellow Rdicts, I have the code for the program below. I need to integrate a function of "x" and "p". I use integrate to integrate over "x" and pass "p" as an additional argument. "p" is specified and given default value in the argument list. Still, integrate() cannot read "p", unless I explicitly insert a numeric value in the

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 wish to ignore. Unfortunately my llbet contains Inf's, which

Hi all! I'm new to this forum so please excuse me if I don't conform perfectly to the protocols on this board! I'm trying to get an estimate of true prevalence based upon results from an imperfect test. I have various estimates of se/sp which could inform my priors (at least upper and lower limits even if with a uniform distribution) and found the following code on this website..

Hi -- I'm trying to package R for inclusion on some upcoming Red Hat Powertools distribution. I have started with the 0.61 version packaged by Martyn Plummer <plummer at iarc.fr> (Thanks Martyn!) and have upgraded to 0.62-2. I have succeeded in running R demo's on i386-redhat-linux-gnu, and now I am also trying to get the alpha and sparc versions working. On sparc-redhat-linux,

Full_Name: Richard Boyce Version: 2.1.-1 OS: Debian testing/unstable Submission from: (NULL) (128.95.123.29) While building a custom package using a modified version of Duncan's REventLoop with R version 2.1 (Debian package r-base, r-base-dev) and R source from apt-get source 2.1.1 I get the following error: $ R CMD build vjREventLoop * checking for file

It is a pleasure for me to announce the new built for Rmetrics Version 191.0057. The source files and Windows binary packages can be downloaded from www.rmetrics.org . The new built has also been submitted to the CRAN server. Some new functions and example files have been added. Unfortunately the user guides and reference guides are not yet updated, they have still the status of Version

In what way is this a bug in R? It looks like a bug in the package, and as Defn.h is not part of R's API any packge using it is `at risk' (and cannot be installed in a binary-only installation, or even an installed version of R). In particular, Defn,.h depends on config.h, and it seems you installed a binary version of R and used separate sources. I would suggest building R from

Running R.app on Mac OS X 10.4 > version _ platform powerpc-apple-darwin8.6.0 arch powerpc os darwin8.6.0 system powerpc, darwin8.6.0 status major 2 minor 3.1 year 2006

Running R.app on Mac OS X 10.4 > version _ platform powerpc-apple-darwin8.6.0 arch powerpc os darwin8.6.0 system powerpc, darwin8.6.0 status major 2 minor 3.1 year 2006

This is a simple SSA code generation 101 question. If I follow the IR code generation techniques in the Dragon book the statement x = y + z would translate into something like this in SSA/LLVM %0 = add %y, %z %x = %0 Obviously "copy instructions" like %foo = %bar are senseless in SSA since %foo and %bar are immutably fixed to the same value and there is no need for two aliases

Dear useRs and wizaRds, I am currently developing a set of functions using S4 classes. On the way I encountered the problem exemplified with the code below. For some reason the 'validity' method does not seem to work, i.e. does not check for errors in the specification of the slots of the defined class. Any hints? My understanding of the whole S4 system was that validity checks are made

Dear R-devel, I'm working on Prof. Loader's new version of locfit to try to get it pass R CMD check. I'm almost there, but I have a problem with some Rd files that I hope some one can help me resolve. Here's an example: In the package there's a function called locfit.censor(). This function can be used in a few different ways: locfit.censor(x, y, cens, ...)

From: Daniel Kaschek <daniel.kaschek at physik.uni-freiburg.de> > ... When I evaluate this list of functions by > another lapply/sapply, I get an unexpected result: all values coincide. > However, when I uncomment the print(), it works as expected. Is this a > bug or a feature? > > conditions <- 1:4 > test <- lapply(conditions, function(mycondition){ >

Hi All, In another thread Andy Liaw, who CRAN lists as locfit maintainer; said: <quote> From: "Liaw, Andy" <andy_liaw at merck.com> To: "Guy Green" <guygreen at netvigator.com>; <r-help at r-project.org> Subject: Re: Alternatives to linear regression with multiple variables Date: 22 February 2010 17:50 You can try the locfit package, which I believe

Hello!!! I'm just a beginner of R. I work with locfit to get the density from an numeric vector. It looks good when I plot the graph, but I need the output data of locfit. For example I need somthing like that (only integer values for the input): numeric vector (input of locfit) density (output of locfit) 1 0.05.... 2 0.04.. 3 0.01 ... Thank you for your help and your time.

Would introducing the new frame, with the call to local(), cause problems when you use frame counting instead of <<- to modify variables outside the scope of lapply's FUN, I think the frame counts may have to change. E.g., here is code from actuar::simul() that might be affected: x <- unlist(lapply(nodes[[i]], seq)) lapply(nodes[(i + 1):(nlevels - 1)],