similar to: inconsistency between gamma and choose functions

.. the round function seems to get round it (excuse the pun): > round(gamma(21)/(gamma(6)*gamma(16))) == round(choose(20,5)) [1] TRUE So I presume the problem is related to gamma being a real function. However, I'm still not sure why there's a break from TRUE to FALSE for certain values. Any comments would be appreciated! _____________________________________ Paul S.P. Cowpertwait

Hi all, I have this code, but the marginal distribution plot doesn?t appear aligned with the left plot. I think could be something about layout or par() mar. The code was programmed by me time ago. Can anyone help me to get the marginal distribution on the center (more higher centered) id.txt Could have this code: 05/01/2016;9335,200195 06/01/2016;9197,400391 07/01/2016;9059,299805

Hi Pedro, As a one-off, you just shove the coordinates around a bit: par(mar=c(11,0,6,6)) barplot(fhist$counts,axes=FALSE, space=0,horiz=TRUE,col="lightgray", ylim=c(0,24)) However, I don't think that this plot illustrates quite what you think it does. Jim On Wed, Jun 7, 2017 at 4:01 PM, Pedro p?ramo <percentil101 at gmail.com> wrote: > Please, I'm trying to put

Are there any functions in R for fitting distributions? In particular, I would like to fit the weibull and extreme-value distributions when: i) the data are given and the maximum likelihood estimates are required; ii) the sample moments are available but the complete data record is not (so moment estimates are required). Any suggestions will be gratefully received.

Many thanks Jim. What I,m trying to show with the fhist plot is the empirical distribution of the values of the left plot simulation. You say: However, I don't think that this plot illustrates quite what you think it does. Can you give me a clue to try to illustrate better if it is not showing what I believe it shows a better way to show it? Many thanks in advance. El 7 jun. 2017

Please can you send me some orientation? Many thanks in advance. Only if posible one book o similar example to understand why it is not what I try. El 8 jun. 2017 7:50 PM, "Pedro p?ramo" <percentil101 at gmail.com> escribi?: > Many thanks Jim. > > What I,m trying to show with the fhist plot is the empirical distribution > of the values of the left plot simulation.

Please, I'm trying to put the right plot higher and centered on the left values but I don't achive. I would appreciate so much your help El 6 jun. 2017 22:37, "Pedro p?ramo" <percentil101 at gmail.com> escribi?: > Hi all, > > I have this code, but the marginal distribution plot doesn?t appear > aligned with the left plot. > > > I think could be

I'd like to substitute the Subversion revision number for a string in a package man page every time the package is installed or built. I can assume it is being installed or built from a Subversion working copy. I tried putting a target that depends on FORCE into Makefile or Makevars in the src directory, but it is not being built. What sort of make magic do I need, and where do I put it

Hi Pedro, If you keep that same margins for the second plot: par(mar=c(10,0,6,6)) barplot(fhist$counts,axes=FALSE, space=0,horiz=TRUE,col="lightgray") it looks reasonably well aligned to me. Because you are plotting the counts of the values in Simulation, the ordinate (vertical axis) of the bar plot is in quite different units from that of the plot on the left side. Jim On Wed, Jun

Please can someone advise me how I can adjust correlations using bonferroni's correction? I am doing manny correlation tests as part of an investigation of the validity/reliability of a psychometric measure. Help would be so appreciated! Cheers, Louise

Dear list fellows, I want to study time series and use R to analyse time series of fishing data from several species (landings and cpue) investigating the correlation between them and with environmental factors (water temperature, wind, etc.). Searching at Amazon I found three books with examples in R: Time Series Analysis: With Applications in R by Jonathan D. Cryer and Jonathan D. Cryer

Hi there, I have subscribed to R-help but am not sure how to view or post questions? I think this is the right way. I am planning on doing a multivariate regression investigating the relationship between depression (a continuous variable) and social support variables (mostly continuous, some categorical) among older people. I have a number of demographic and health-related variables that I am

Can you suggest a plot on the right side so that the right plot has the same units of left plot and reflect the counts of the number of time the values repeat, something like the empirical distribution? 2017-06-14 12:30 GMT+02:00 Jim Lemon <drjimlemon at gmail.com>: > Hi Pedro, > If you keep that same margins for the second plot: > > par(mar=c(10,0,6,6)) >

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

Dear R users, I have to calculate gamma functions for negative numbers beyond -171.4. e.x. gamma(-500.4) I got following: > gamma(-170.4) [1] -5.824625e-308 > gamma(-171.4) [1] 0 Warning message: underflow occurred in 'gammafn' I have tried to use a recursion getting values a little futher -180. How could I solve this problem? Thank you beforehand. Chuse.

The following seems very strange: > lgamma(-0.04) [1] 3.243307 > gamma(-0.04) [1] -25.61830 Jim -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-devel 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:

I have to calculate some formula like: gamma(x)/(gamma(x+y) and I observed that for relatively big values of x, R returns infinity and so cannot compute the formula. Is it possible to force R to give the real value of gamma(x) instead of Inf ? thanks

I'm looking for goodness of fit tests for gamma distributions with large data sizes and for different data. I have a matrix with around 4.000 data values in it and i have fitted a gamma distribution with "fitdistr". You can see the example: > fitdistr(corpo,"gamma",lower=0.001) Errore in optim(x = c(5000, 5000, 5000, 5000, 5000, 5000, 5000, 5000, 5000, :

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

I need to know exactly how Wine tells the X server to change gamma settings... so gamma controls in games can work. Basically I need to know.. when a game has like a gamma adjust bar, and it moves... what specifically happens behind the scenes. I need to know this because right now Xquartz on Mac OS X does not support doing this at all, and I'd like to be able to add that ability into it...