similar to: Goodness of fit for gamma distributions

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 with this question again: I would like to test few distributional assumptions for some behavioral response data. There are few theories about true distribution of those data, like: normal, lognormal, gamma, ex-Gaussian (exponential-Gaussian), Wald (inverse Gaussian) etc. The best way would be via qq-plot, to show to students differences. First two are trivial: qqnorm(dat$X)

Hi there, I'm a newbie, plesae bear with me. I have a dataset with about 10000 ~ 30000 data points. Would like fit to both Gamma and Normal distribution to see which one fits better. How do I do this in R? Or I could do a normality test of the data, if it's normal, I then will do a normal fit, otherwise, a gamma fit. But again, I don't know how to do this either. Please help! David

Hi, I have the following vector which is created from 3 distinct distribution (three components) of gamma: x=c(rgamma(30,shape=.2,scale=14),rgamma(30,shape=12,scale=10),rgamma(30,shape=5,scale=6)) I want to plot the density curve of X, in a way that it shows a distinct 3 curves that represent each component. How can I do that? I tried this but doesn't work: lines(density(x)) Please

Hi, I want to create upper and lower 95% confidence intervals for a p-p plot of an empirical distribution with a theoretical gamma distribution. This is my code: x<-rgamma(100,shape=2, rate=1) # empirical data fitdistr(x,"gamma") # fit a gamma distribution dist<-pgamma(x,shape=1.9884256 ,rate=0.8765314 ) # fitted distribution, using the loglikelihood estimated parameters

Previously the code relied on device->pri to be NULL and to fail probing later. We really should just return an error inside nvkm_device_ctor for unsupported GPUs. Fixes: 24d5ff40a732 ("drm/nouveau/device: rework mmio mapping code to get rid of second map") Signed-off-by: Karol Herbst <kherbst at redhat.com> Cc: dann frazier <dann.frazier at canonical.com> ---

By any chance is any one aware of a website, book, paper, etc. or combinations of those sources that show plots of different distributions? After reading a pretty good whitepaper I became aware of the benefit of I the benefit of doing Q-Q plots and histograms to help assess a distribution.   The whitepaper is called: "Univariate Analysis and Normality Test Using SAS, Stata, and SPSS*" ,

On Mon, Oct 05, 2020 at 03:52:30PM +0200, Karol Herbst wrote: > Previously the code relied on device->pri to be NULL and to fail probing > later. We really should just return an error inside nvkm_device_ctor for > unsupported GPUs. > > Fixes: 24d5ff40a732 ("drm/nouveau/device: rework mmio mapping code to get rid of second map") > Signed-off-by: Karol Herbst

On Tue, Oct 06, 2020 at 06:46:35PM +0200, Karol Herbst wrote: > On Tue, Oct 6, 2020 at 6:43 PM dann frazier <dann.frazier at canonical.com> wrote: > > > > On Mon, Oct 05, 2020 at 03:52:30PM +0200, Karol Herbst wrote: > > > Previously the code relied on device->pri to be NULL and to fail probing > > > later. We really should just return an error inside

Hi Everyone, I want to make 3 lines on the same graph (not as subplots, all in the same window, one on top of each other) and I want them to be quantile-quantile plots (qqplot). Essentially, I am looking for the equivalent of Matlab's "hold on" command to use with qqplot. I know I can use 'points' or 'lines', but these do not give me a qqplot (only appear to work

Dear People, In a qqplot I am doing, I get lines/points that are very thick. I've tried setting the lwd variable to 0.1, but it doesn't seem to have any effect. Also, I have set the value of lty to dashed, but I still get dots. The command looks like qqplot(cdf.inv(seq(0,1,length=size),theta,pos,len),empmargdistvec(len,theta,pos,size), xlim=c(-theta,theta), ylim=c(-theta,theta),

Hi, We could use qqplot to see how two distributions are different from each other. To show better how they are different (departs from the straight line), how is it possible to plot the straight line that goes through them? I am looking for some thing like qqline for qqnorm. I thought of abline but how to determine the slope and intercept? Best wishes, Carol

I am a beginner in R, so please don't step on me if this is too simple. I have two data sets datax and datay for which I created a qqplot qqplot(datax,datay) but now I want a line that indicates the perfect match so that I can see how much the plot diverts from the ideal. This ideal however is not normal, so I think qqnorm and qqline cannot be applied. Perhaps you can help? Ralf

I am trying to do something in R and would appreciate a push into the right direction. I hope some of you experts can help. I have two distributions obtrained from 10000 datapoints each (about 10000 datapoints each, non-normal with multi-model shape (when eye-balling densities) but other then that I know little about its distribution). When plotting the two distributions together I can see that

Hello, I would like to use qqplot() to compare two experimental distributions. But I do not understand how qqplot() compute quantiles. In fact, quantile() do not return the same results. Thank you for your help. Vincent.

>>>>> On Tue, 31 Aug 1999 14:57, Werner Stahel <stahel@stat.math.ethz.ch> said: WSt> Here is a suggestion. It seems that qqplots, comparing a sample WSt> to a distribution other than the normal, are not explicitly WSt> available in S or R. I found (in S-plus / Trellis it is, see below) WSt> qqplot(y, rt(300, df = 5)) WSt> as an

Hi All, I am trying to plot 4 graphs on to 1 page using layout(...), or par(mfcol = c(...)); with the function QQplot from the package PerformanceAnalytics. The problem is that, no matter what order I use, it only plots 3 graphs on to 1 page and the last QQplot is shunted to the next page. Also, this only happens to the QQplot, i.e. there is no problem with 4 Histograms. set.seed(1033) data

Hi Group, Here is some data. p <- runif(1000) # sample data groups <- rep(c(1,2),each=500) #conditioning variable mydata <- cbind(p,groups) n <- length(p) u <- (1:n)/(n + 1) # uniform distribution reference for qqplot logp <- -log(p,base=10) logu <- -log(u,base=10) qqplot(logp,logu) How can I make the above qqplot in lattice and/or ggplot2. The sample is uniform, and I take

Full_Name: Matthias von Davier Version: 1.2.0 OS: linux Submission from: (NULL) (144.81.31.148) pairs(NxK_Matrix,panel=qqplot) produces a message Error in pairs.default(NxK_Matrix,panel=qqplot) : The panel function made a new plot best regards and a happy new year -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-devel mailing list -- Read

hi everyone! season's greetings! is there any way that i can create a normal qqplot showing, aside from the qqline, the 95% confidence limits? thank you very much.. happy holidays! -- View this message in context: http://n4.nabble.com/how-to-create-normal-qqplot-with-the-95-confidence-interval-tp977727p977727.html Sent from the R help mailing list archive at Nabble.com.