similar to: Incomplete Gamma function

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, 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 installed the locfit package under Linux (gcc 2.7.2). Installation was ok but > x <- runif(200) > y.compl <- 10*x*x*rgamma(200,3) > med.y <- median(y.compl) > cens <- ifelse(y.compl<=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

Could anyone help me with the following issue. Using the GSL library in R I define the following code: ######### library(gsl); S<-function(n) { r<-0:n; ans<-sum(gsl_sf_choose(n,r)*(-1)^r*2^(2*r)*gamma_inc(6-2*r,2)) ans } ######### >SS(10) yields 34.91868 >SS(40) yields 5.340422 >SS(60) yields 180.3162 Doing the same computations in maple I get

Dear all, I would be very grateful if you could help me with: Given the regularized gamma function Reg=int_0^r (x^(k-1)e^(-x))dx/int_0^Inf (x^(k-1)e^(-x))dx ; 0<r<Inf (which is eventually the ratio of the Incomplete gamma function by the gamma function), does anyone know of a package in R that would evaluate the derivative of the inverse of Reg with respect to k? I am aware that the

Can someone please enlighten me as to why the following happens? > -2.7^8.6 [1] -5125.407 > p<- -2.7 > q<- 8.6 > p^q [1] NaN R seems perfectly able to calculate -2.7^8.6, but fails when the exact same values are assigned to variables and then the computation is repeated. Thanks in advance for any suggetsions. Kris.

Hi, I am new to using R for Time series analysis. I was wondering if there are any functions that can convert ARMA or ARIMA time series into their corresponding AR or MA time series representations (by calculating the corresponding AR or MA coefficients). Thanks a lot Kris.

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

Hello, I am trying to call a FORTRAN subroutine from R. The Fortran code is @: http://lib.stat.cmu.edu/apstat/206 It performs a bivariate isotonic regression on a rectangular grid (m X n) matrix. I used the g77 compiler and successfully created a dll file and it also loads successfully from R. But somehow the programs fails to run properly. (I do get the correct result when I compile the

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

Hi all, I am a biostatistician and I have developed my own ranking system for clinical data. I would like to test the efficiency of it w.r.t. to other ranking systems. I would like to simulate the data and after assigning ranks to my observed scores(after neglecting dropouts), observe the type I error. If I want to do a Kruuskal Wallis type of test, what test statistic should I use to test for a

I'm testing a FC5 machine for use in a student lab. R 2.3.1 is installed and seems to work fine. There is one peculiarity - the logins are authenticating to a server, and a "verbose" flag is set somewhere, leading to lots of spurious messages like this request done: ld 0xa227598 msgid 1 which may be confusing R. However, R CMD check seems to fail for packages with no

Hi, I have a question on generating random variables based on zipf-mandelbrot distribution. So when I execute the following lines: ZM = lnre ("zm", alpha = 2/3, B=0.1) zmsample = rlnre (ZM, n =100) zmsample It generates 100 random values based on a zipf-mandelbrot distribution as below. But how do I make sure the generated random number is within the range of 1 - 6000 only? Can I

Hi, I''m trying to program custom type for managing Couchbase buckets. Problem is that every RW operation to Couchbase cluster (creating, deleting, resizing bucket) requires admin privileges (username/password). Do you have any suggestion how should I pass the username/password to provider? If I make it a param, then every resource of type ''couchbucket'' will need to

So far we have a similar situation, for each different server one fabric and one puppet file, where the fabric file simply applies it in a brutal way. with settings(user=''root''): put(''qa.pp'', ''qa.pp'') put(''puppet apply qa.pp'') And puppet files don''t use anything like classes or modules, but

Dear R-users, I am new to R and would like to use it for fitting the zipf distribution to some numeric data that I have. Here's the snippet that I use: library(VGAM) X <- read.table(file("~\\mydata.txt", encoding="latin1")) w <- as.vector(t((X[2]))) w <- w/sum(w) y <- (1:length(w)) fit = vglm (y ~ 1, zipf, tra=TRUE, weight=w) zipf(N=NULL,

Using R, I plotted a log-log plot of the frequencies in the Brown Corpus using plot(sort(file.tfl$f, decreasing=TRUE), xlab="rank", ylab="frequency", log="x,y") However, I would also like to add lines showing the curves for a Zipfian distribution and for Zipf-Mandelbrot. I have seen these in many articles that used R in creating graphs. Thank you! [[alternative HTML

Hey all, Just a sanity check, but should the shorecap script in shorewall6-lite be sourcing /usr/share/shorewall6-lite/lib.base rather than /usr/share/shorewall-lite/lib.base like it does currently? In fact shouldn''t there be a general s/shorewall-lite/shorewall6-lite/ in shorecap in shorewall6-lite? Maybe there is more of that lurking about as well. Also, the first line of the

Hello R, To analyze multi-level data, I started learning and using lmer. So far so wonderful. I then found some useful functions in package languageR. But then the following problem ocurred: Whenever I load and use the languageR package, then save the workspace - or quit R with saving the workspace - I am unable to reload that workspace in a later session. That is, R doesn't start at all