similar to: how to fit in R

Dear All, I run the goodness of fit test using goodfit() in vcd package. The result is as follow: Goodness-of-fit test for poisson distribution X^2 df P(> X^2) Pearson 1.053348 2 0.5905661 Warning message: In summary.goodfit(gf) : Chi-squared approximation may be incorrect I want to save the the test statistics(X^2), df, and p-value. How can I save the result.

Dear R-Users, How can I use chisq.test() as a goodness of fit test? Reading man-page I?ve some doubts that kind of test is available with this statement. Am I wrong? X2=sum((O-E)^2)/E) O=empirical frequencies E=expected freq. calculated with the model (such as normal distribution) See: http://www.itl.nist.gov/div898/handbook/eda/section3/eda35f.htm for X2 used as a goodness of fit test. Any

Hi everybody, somehow i dont get the shapiro wilk test for normality. i just can?t find what the H0 is . i tried : shapiro.test(rnorm(5000)) Shapiro-Wilk normality test data: rnorm(5000) W = 0.9997, p-value = 0.6205 If normality is the H0, the test says it?s probably not normal, doesn ?t it ? 5000 is the biggest n allowed by the test... are there any other test ? ( i know qqnorm

Hi there, I am trying to perform different normality tests in R. I have no problem when I use Shapiro.test( ). However, when I try to use any normality tests from the Nortest package, I kept getting this error " Error in read.table(source, header = T, fill = T) : 'file' must be a character string or connection". I have no idea why because there was no problem with

Hi All, I have a dataset consisting of abundance counts of a fish and I want to test if my data are poisson in distribution or normal. My first question is whether it is more appropriate to model my data according to a poisson distribution (if my test says it conforms) or use transformed data to normalise the data distribution? I have been using the vcd package gf<-goodfit(Y,type=

Hello, Sample code: library("vcd") dummy <- rnbinom(200, size=1.5, prob=0.8) gf <- goodfit(dummy, type="nbinomial", method="MinChisq") plot(gf) I would like to: 1. add a lgened stating the bars show the actual counts and the red dots - the fit. 2. show the goodness-of-fit values calculated somewhere on an empty white space ob the plot. But... the legend

Hi, is there any method for goodness of fit testing of an (as general as possible) univariate distribution with parameters estimated, for normal, exponential, gamma distributions, say (e.g. the corrected p-values for the Kolmogorov-Smirnov or Chi-squared with corresponding ML estimation method)? It seems that neither ks.test nor chisq.test handle estimated parameters. I am aware of function

Dear r list,   I am using glm.nb in the MASS package to fit negative binomial models to data on manta ray abundance, and AICctab in the bbmle package to compare model IC.  However, I need to test for the goodness of fit of the full model, and have not been able to find a Pearson's Chi Squared statistic in any of the output.  Am I missing it somewhere?  Is there a way to run the test using

Dear R-Help: As you can see, Prof. Friendly refers me to your site for an executable version of vcd. I don't mean to be obtuse, but 15 minutes spent exploring your site failed to locate a downloadable version of the vcd package to which he referred. I know plainly what this application can do. What I need to know is how to obtain the application itself. My thanks in advance for any

I have written a function to emulate minitab's QQ plotting output (with SW test and AD test results on the graph): mtab.norm<-function(x) { library(nortest) library(lattice) x<-as.numeric(x) x<-as.vector(x) plot.ht<-4.6 plot.wd<-4.6 pt.ht=plot.ht/5 txt.sz<-(plot.ht/7.5) X11(width=plot.wd, height=plot.ht, bg='gray96') qqplot(x, pch=16, cex=pt.ht,

Try changing your method to "ML" and try again. I tried the run the first example from the documentation and it failed with the same error. Changing the estimation method to ML worked. @List: Can anyone else verify the error I got? I literally ran the following two lines interactively from the example for goodfit: dummy <- rnbinom(200, size = 1.5, prob = 0.8) gf <- goodfit(dummy,

Good morning, i'm drawing a barplot (data are stored in a matrix,nrow=3,ncol=10). So i get 10 groups each containing 3 sub-bars. I need to change the range of the y-axis when i plot respectively the first and the other two sub-bars in each group. I can't manage. Could you help me? Thank you. Angela

I have troubles understanding how goodfit() function in the vcd package computes the Pearson coefficient. Can anybody provide more information on the computation? In particular, for HorseKicks data in vcd package, goodfit() yields > oo <- goodfit(HorseKicks,type="poisson",method="MinChisq") > summary(oo) Goodness-of-fit test for poisson distribution

I have a dataset which I want to model using a Poisson distribution, with a given parameter. I would like to know what is the proper way to do a ''goodness of fit'' test using R. I know the steps I''d take if I were to do it ''manually'': grouping the numbers into classes, calculating the expected frequencies using ''ppois'', then

Hi, from what you're writing: "The logaritmic transformation "shapiro.test(log10(y))" says: W=0.9773, p-value= 2.512e-05." it seems the log-values are not distributed normally and so original data are not distributed like a log-normal: the p-value is extremally small! Other tests for normality are available in package: nortest compare the log-transformation of your ecdf

I believe there to be a bug in the cvm.test module of the nortest package authored by Juergen Gross. I do not know how to contact the author directly. I've been running some normality tests using the nortest package. For some of my datasets the Cramer-von Mises normality test generates an extremely high probability (e.g., 1.637e+31) and indicates normality when the other tests do

Hello all: Is there a way to detect in R if a dataset is normally distributed or skewed without graphically seeing it? The reason I want to be able to do this is because I have developed and application with Visual Basic where Word,Access and Excel "talk" to each other and I want to integrate R to this application to estimate confidence intervals on fish sizes (mm). I basically want to

Hello, Do you know a R implemented normality test like the shapiro test but more suitable for large data set ? Thanks, _________________________________________________________________ Découvrez toutes les possibilités de communication avec vos proches [[alternative HTML version deleted]]

Hello all, Would somebody be kind enough to show me how to do a KS test in R for a lognormal distribution with ESTIMATED parameters. The R function ks.test()says "the parameters specified must be prespecified and not estimated from the data" Is there a way to correct this when one uses estimated data? Regards, Kwabena. -------------------------------------------- Kwabena Adusei-Poku

I simulate a uniform data and then transformed into Burr(1,3,1) data, which is of pdf: f(x)=[3*(x^2)] / [(1+x^3)^2], x>0 How can I perform a goodness-of-fit test (k-s, anderson-darling,chisq,cramer-von mises,...) on it (should highly accept) to get test-statistics & p-values? Thanks! Sincerely, Shelton Jin -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-