similar to: choose incorrect for fractional and some negative integer values (PR#10766)

choose(-5,-7) uses integer arguments (as specified in Help) and returns a numeric value that is incorrect. Either the function or the documentation should be fixed. If the function is not fixed, a warning or an error would be helpful. The fact that choose(n,k) usually returns choose(n,round(k,0)) is not obvious from either the output or the documentation. I suggest issuing a warning when

Full_Name: Jerry W. Lewis Version: 2.7.0 OS: Windows XP Professional Submission from: (NULL) (198.180.131.16) If you are saying that there is no need to solve for the noncentrality parameter, please justify this amazing assertion. If you are saying that this need is already adequately addressed in R, then please enlighten me.

Full_Name: Jerry W. Lewis Version: 2.7.0 (2008-03-23 r44847) OS: Windows XP Professional Submission from: (NULL) (71.184.230.48) choose(n,k) = choose(n,n-k) is not satisfied if either 1. n is a negative integer with k a positive integer (due to automatically returning 0 for n-k<0) 2. n is not an integer with k a positive integer (due to rounding n-k to an integer, compounded by

Hi Firstly , happy Christmas to R-Help! Secondly, I wonder if anyone can help me with the following query: I am trying to reproduce some explicit probability calculations performed in APPL (a Maple extension for computational probability). For instance, in APPL, to compute the probability that the sum of 10 iid uniform variables [0,1] will be between 4 and 6, (i..e Pr( 4 < \sum_{i=1}^{10}X_i

While I agree that the reported results from Mathematica have only 10-13 correct digits, that does not mean that pt() in R is any better for these calculations. For instance the following three calculations are mathematically equivalent, but pt() disagrees at the 13th figure in R v2.6.2 pt(1e-4,13) pf(1e-8,1,13)/2+0.5 pbeta(1/(1+13/1e-8),.5,6.5)/2+0.5 Using 1/(1+n/x^2) and reversing the

Full_Name: Jerry W. Lewis Version: 2.6.1 OS: Windows XP Professional Submission from: (NULL) (24.147.191.250) pchisq(0,0,ncp=lambda) returns 0 instead of exp(-lambda/2) pchisq(x,0,ncp=lambda) returns NaN instead of exp(-lambda/2)*(1 + SUM_{r=0}^infty ((lambda/2)^r / r!) pchisq(x, df + 2r)) qchisq(.7,0,ncp=1) returns 1.712252 instead of 0.701297103 qchisq(exp(-1/2),0,ncp=1) returns 1.238938

Full_Name: Jerry W. Lewis Version: 2.10.1 OS: Windows XP Professional Submission from: (NULL) (198.180.131.21) In the line "The quantile is right continuous: qpois(q, lambda) is the smallest integer x such that P(X <= x) >= q." "q" is used as a probability when the Arguments section defines it to be a quantile. Also there are some representation problems where the

The data is of two forms, ie one numeric and another ordinal ie 1.5 and <.5. The X Value is a factor. I'm trying to create a lattice conditioned plot with the following characteristics: 1) The plot is conditioned using the form (Conductivity~Day|Valve) 2) The plot should use a barplot for the ordinal (<.5) and dots for the numeric (1.5) 3) A line should be created specifying a

I'd like to fit the (1-displaced) negative hypergeometric distribution to data samples such as the following:- x <- c(370, 311, 299, 266, 265, 232, 197, 198, 170, 154, 133, 123, 120, 103, 80, 72, 69, 67, 67, 50, 36, 35, 26, 23, 15, 11, 9, 6, 5, 3, 3, 2, 2, 2) i.e., I want to estimate the parameter values of K and M (with my data, n would usually be the same as the length of the data

I have a question. I can find a function about "Gaussian hypergeometric function" in sofeware Maple V. I need this function to calculate. Is there a similar function in R ? Or I have to design such a function. Thanks! Jimmy Cheng

Hi, I need Beta binomial and Beta negative binomial functions but in R there is only SuppDists package which provide this distributions using a limited parameter space of the generalized hypergeometric distribution (dghyper & Co.) which provide a limited parameter space for Beta binomial and Beta negative binomial functions (e.g. alpha + beta <1 in the Beta negative binomial). I've

Dear R Helpers, I have searched for any R package or code for simulating multivariate fractional Brownian motion (mFBM) or multivariate fractional Gaussian noise (mFGN) when a covariance matrix are given. Unfortunately, I could not find such a package or code. Can you suggest any solution for multivariate FBM and FGN simulation? Thank you for your help. Best Regards, Ryan ----- Wonsang You

Dear R Users, There exists a non-central hypergeometric distribution function in the (MCMCpack) package, and a hypergeometric distribution function in the (stats) package. Is there a function for sampling from an extended hypergeometric distribution? Thanks, Narcyz This message is intended for the addressee named and may con...{{dropped}}

[This email is either empty or too large to be displayed at this time]

Hi, I hope somebody can help me on how to use the hypergeometric function. I did read through the R documentation on hypergeometric but not really sure what it means. I would like to evaluate the hypergeometric function as follows: F((2*alpha+1)/2, (2*alpha+2)/2 , alpha+1/2, betasq/etasq). I'm not sure which function should be used- either phyper or qhyper or dhyper Where

Dear folk, I am stempting to estimate a vector error correction model using a seemingly fractionally integrated multivariate time series. The *fracdiff *package provides tools to estimate degree of fractional integration. But *fracdiff *can't help me to: 1. test equality of two degrees of fractional integration, say d1=d2? 2. estimate a multivariate cointegrating error correction model,

folks, does anyone know if the SEMIFAR model has been implemented in R? i see that there's a S-FinMetrics function SEMIFAR() that does the job, but I have no access to that software. essentially, this semiparametric fractional autoregressive model introduces a deterministic trend to the FARIMA(p,d,0) model (which, as i understand it, takes care of the random trend and short and long memory).

Hello, I am searching for code in C++ or fortran for Modified Bessel function of third kind (fractional or real order). Can someone help me? Thank you --------------------------------- Découvrez un nouveau moyen de poser toutes vos questions quelque soit le sujet ! Yahoo! Questions/Réponses pour partager vos connaissances, vos opinions et vos expériences. Cliquez ici.

Dear All, I have analysed time to event data for continuous variables by considering the multivariable fractional polynomial (MFP) model and comparing this to the untransformed and log transformed model to determine which transformation, if any, is best. This was possible as the Cox model was the underlying model. However, I am now at the situation where the assumption that the competing risks

Hi, I created a 4 Meg Byte primanry DOS patition onto a hard disk. This partition was formatted with format /s. The hard drive was able to boot DOS. After that I connected the drive to a Linux machine and created and image with the command dd if=/dev/had1 of=hard.img When I send this image to my PXE client I see on the screen: MEMLINUX: Image seems to have fractional end cylinder MEMLINUX: