similar to: Non-homogeneity of variance - decreasing variance

Hi, All! To calculate Levene's Test for Homogeneity of Variance I use R Commander, and this is the output: > levene.test(Dataset$age, Dataset$sex) Levene's Test for Homogeneity of Variance Df F value Pr(>F) group 1 0.8739 0.3567 33 I am not sure what means "Pr(>F)"? Can anyone explain/translate this? Regards, Iurie Malai Department of Psychology and

Has the Levene test of homogeneity of variance been implemented in any library in R? Thanks, Maurice Haynes National Institute of Child Health and Human Development Child and Family Research Section 6705 Rockledge Drive Bethesda, MD 20892 Voice: 301-496-8180 Fax: 301-496-2766 E-Mail: mh192j at nih.gov

Dear all, I'm a beginner of R and I need to carry out some three-way mixed ANOVAs. Following examples at http://personality-project.org/r/r.anova.html, I managed to get the ANOVA part, but I don't know how can I check data normality and homogeneity of variance in R (since they're the required assumptions of ANOVA analysis). Are there any special command/packages? Could anyone give me

Dear all a couple of months ago i've found threads regard test that verify AnOVa assumption on homogeneity of variances. Prof. Ripley advice LDA / QDA procedures, many books (and many proprietary programs) advice Hartley's F_max, Cochran's minimum/maximum variance ratio (only balanced experiments), K^2 Bartlett's test, Levene's test. Morton B. Brown and Alan B. Forsythe in a

Dear R users I have been looking for functions that can deal with overdispersed poisson models. Some (one) of the observations are negative. According to actuarial literature (England & Verall, Stochastic Claims Reserving in General Insurance , Institute of Actiuaries 2002) this can be handled through the use of quasi likelihoods instead of normal likelihoods. The presence of negatives is not

I am testing the homogeneity of variances via bartlett.test and fligner.test. Using the following example, how should I interpret the p-value in order to accept or reject the null hypothesis ? set.seed(5) x <- rnorm(20) bartlett.test(x, rep(1:5, each=4)) Bartlett test of homogeneity of variances data: x and rep(1:5, each = 4) Bartlett's K-squared = 1.7709, df = 4, p-value =

En innebygd og tegnsett-uspesifisert tekst ble skilt ut... Navn: ikke tilgjengelig Nettadresse: https://stat.ethz.ch/pipermail/r-help/attachments/20080214/5a98c799/attachment.pl

Dear r-helpers, Prior to analysis of variance, I ran the Boxcox function (MASS library) to find the best power transformation of my data. However, reading the Boxcox help file, I cannot figure out if this function (through its associated log-likelihood function) corrects for * normality only * or if it also induces * homogeneity of variances *. I found in Biometry (Sokal and Rohlf, p. 419)

Dear all, sorry if this is a bit on the sidetrack for R-help. As a regular R user I have developed quite a lot of custom R functions, to the point of not always remembering what I have already programmed, where the file is and so on. I was wondering what other people do in this regards. A basic file with all your functions, or a custom R package, or directly integrated into a profile file

Dear all, I have been trying to find a function to calculate the Breslow-Day test for homogeneity of the odds ratio in R. I know the test can be preformed in SAS but i was wondering if anyone could help me to perform this in r. In addition i have the fullrefman file to search for functions in the basic R packages, does anyone have any suggestions of an efficient way of searching for

Un texte encapsul? et encod? dans un jeu de caract?res inconnu a ?t? nettoy?... Nom : non disponible Url : https://stat.ethz.ch/pipermail/r-help/attachments/20060609/ab46b028/attachment.pl

Hi, I am trying to run a generalized linear regression using a negative binomial error distribution. However, I want to use an overdispersion parameter that varies (dependent on the length of a stretch of road) so glm.nb will not do. >From what I've read I should be able to do this using GLM by specifying my own quasi family and describing the variance function using varfun, validmu,

Hello, We've got a dataset with several variables, one of which we're using to split the data into 3 smaller subsets. (as the variable takes 1 of 3 possible values). There are several more variables too, many of which we're using to fit regression models using lm. So I have 3 models fitted (one for each subset of course), each having slope estimates for the predictor variables.

I have a few questions relating to overdispersion in a sex ratio data set that I am working with (note that I already have an analysis with GLMMs for fixed effects, this is just to estimate dispersion). The response variable is binomial because nestlings can only be male or female. I have samples of 1-5 nestlings from each nest (individuals within a nest are not independent, so the response

Hello, I have a question regarding estimation of the dispersion parameter (theta) for generalized linear models with the negative binomial error structure. As I understand, there are two main methods to fit glm's using the nb error structure in R: glm.nb() or glm() with the negative.binomial(theta) family. Both functions are implemented through the MASS library. Fitting the model using these

Dear R users I have been looking for functions that can deal with overdispersed poisson models. According to actuarial literature (England & Verall, Stochastic Claims Reserving in General Insurance , Institute of Actiuaries 2002) this can be handled through the use of quasi likelihoods instead of normal likelihoods. However, we see them frequently in this type of data, and we would like to

Dear R-helpers, a S-code for piecewise regressions was provided by Toms & Lesperance (2003) Ecology, 84, 2034-2041 (paper can be found on the web).The code is quite complete with different types of transitions around breakpoints and model selection fonctions. It doesn't work directly under R due to some "translation" problems I guess. However I reckon that it would be a

 Hello everyone, Apologies in advance, as this is partially a stats question and partially an R question.  I have been using a GAM to model the activity level of bats going into and coming out from a forested edge.  I had eight microphones set up in a line transect at each of eight sites, and I am hoping to construct a model for each of 7 species.  My count data has a reverse J-shaped skew and

Hello! I am trying to model data on species abundance (count data) with a poisson error distribution. I have a fixed and a random variables and thus needs a mixed model. I strongly doubt that my model is overdispersed but I don't know how to get the overdispersion parameter in a mixed model. Maybe someone can help me on this point. Secondly, it seems that quasi models cannot be implemented

I would like to try a simple parametric bootstrap, but unfortunately (stupidly?) my models are "overdispersed" gams & glms. I'm hoping for a function that generates overdispersed poisson or negative binomial data with a given mean, scale (& shape parameter). The loose definition I'm using is overdispersed poisson produces integer values with variance=const*mean &