similar to: Re: gregmisc version 0.7.3 now available

Version 0.8.0 of the gregmisc package is now (or will shortly be) available on CRAN. New in this release: - Enhanced and bug-fixed 'CrossTable' function (contributed by Marc Schwartz) - Augmented the 'barplot2' function with an 'add' argument to allow for the addition of a barplot to an existing graphic. (contributed by Marc Schwartz) - Added the

Version 0.6.0 of the gregmisc package has just been released on CRAN. Description: =========== Package: gregmisc Description: Misc Functions written/maintained by Gregory R. Warnes Title: Greg's Miscellaneous Functions Version: 0.6.0 Date: 2002/09/23 Maintainer: Gregory R. Warnes <Gregory_R_Warnes at groton.pfizer.com> Author: Gregory R. Warnes. Includes code provided by Ben

Version 0.6.0 of the gregmisc package has just been released on CRAN. Description: =========== Package: gregmisc Description: Misc Functions written/maintained by Gregory R. Warnes Title: Greg's Miscellaneous Functions Version: 0.6.0 Date: 2002/09/23 Maintainer: Gregory R. Warnes <Gregory_R_Warnes at groton.pfizer.com> Author: Gregory R. Warnes. Includes code provided by Ben

> From: kjetil halvorsen [mailto:kjetilh at umsanet.edu.bo] > > Martin Hoyle wrote: > > > > Dear Greg, > > Thanks for the new release. The decomposition of the SSQ is > just what I need! > > > I thought decomposition of the SSQ can be obtained with the split= > argument to summary.aov? That argument seems tro be new --- it didnt > exist two years

Hello All, Is there a way of producing an ANOVA table split into contrasts, thus showing the contrasts sums of squares and associated p-values? Thanks, Martin. Martin Hoyle, School of Life and Environmental Sciences, University of Nottingham, University Park, Nottingham, NG7 2RD, UK Webpage: http://myprofile.cos.com/martinhoyle

Error: could not find function "hier.part" things I have tried: 1. reinstall R (lastest version, on windows XP) 2. install package "gtools" 3. include: library(gtools) require(gtools) 4. how I call this function: hier.part(c$Y, xdata, fam = "gaussian", gof = "Rsqu") 5.when I try to check what's in the package "gtools", I get (hier.part is

Hi all, I have a problem using rdirichlet{gtools}. For Dir( a1, a2, ..., a_n), its mode can be found at $( a_i -1)/ ( \sum_{i}a_i - n)$; The means are $a_i / (\sum_{i} a_i ) $; I tried to study the above properties using rdirichlet from gtools. The code are: ############## library(gtools) alpha = c(1,3,9) #totoal=13 mean.expect = c(1/13, 3/13, 9/13) mode.expect = c(0, 2/10, 8/10) #

I've thrown a couple of useful functions into a new package "gregmisc" and uploaded them to ftp://ftp.ci.tuwien.ac.at/incoming/gregmisc_0.1.1.tgz Here's DESCRIPTION: Package: gregmisc Description: Misc Functions written/maintained by Gregory R. Warnes Title: Greg's Miscellaneous Functions Version: 0.1 Date: 2001/05/24 Depends: R Maintainer: Gregory R. Warnes

Hi How do I install "gregmisc" packages? I did- % sudo R > install.packages("gregmisc") . . > barplot2() but, Error: couldn't find function "barplot2" -- atuya Mac OSX 10.2.6 R 1.7.1

Hi everyone, Can someone explain me how to calculate SAS type III sums of squares in R? Not that I would like to use them, I know they are problematic. I would like to know how to calculate them in order to demonstrate that strange things happen when you use them (for a course for example). I know you can use drop1(lm(), test="F") but for an lm(y~A+B+A:B), type III SSQs are only

Dear All, I'm not sure of the interpretation of interactions with contrasts. Can anyone help? I do an ANCOVA, dryweight is covariate, block and treatment are factors, c4 the response variable. model<-aov(log(c4+1)~dryweight+treatment+block+treatment:block) summary(model); Df Sum Sq Mean Sq F value Pr(>F) dryweight 1 3.947 3.947 6.6268 0.01076 *

Dear R Users, I suppose this is a school boy question, but here it is anyway. I'm trying to re-create the residuals for a poisson GLM with simulated data; x<-rpois(1000,5) model<-glm(x~1,poisson) my.resids<-(log(x)- summary(model)$coefficients[1]) plot(my.resids,residuals(model)) This shows that my calculated residuals (my.resids) are not the same as residuals(model). p 65 of

Hello R Users, I am investigating the basic use of the LME function, using the following example; Response is Weight, covariate is Age, random factor is Genotype model.lme <- lme (Weight~Age, random=~ 1|Genotype) After summary(model.lme), I find that the estimate of Age is 0.098 with p=0.758. I am comparing the above model with the AOV function; model.aov <- aov (Weight~Age + Genotype)

Hola comunidad, Les consulto mi duda. Quisiera generar N números aleatorios tal que su suma esté predeterminada a un valor k. Sé que con la distribución Dirichlet, por ejemplo, se podrían obtener números aleatorios tal que su suma es la unidad. Por ejemplo MCMCpack::rdirichlet(1,rep(10,7)) genera 7 números tal que su suma es 1 (el parámetro rep(10,7) es cualquier cosa). Pero me pregunto si se

I am a pretty new R user, I am running the latest linux version on xandros, updated with some extra debian packages, and I also run the latest windows version, but prefer linux. I am having trouble downloading "bayesm", it won't do it all from any of the sites on the web, I resorted to this one, http://packages.debian.org/unstable/math/r-cran-bayesm. and got slightly further, but I

Looking over the contents of various packages, including my own, it is clear that lots of things end up 'hidden away' in packages where they don't belong. My gregmisc package is a particularly egregious example, containing something from almost every functional category. I propose that from time to time the R community go through the complete set of packages and 'refactor'

Full_Name: John Hendrickx Version: 1.6.0 OS: Windows 98 Submission from: (NULL) (137.224.174.216) CrossTable in the "gregmisc" package fails when the fisher.exact test produces an error (I suspect this is because the number of cases is too large). This can be fixed using "FTt <- try(fisher.test(t, alternative = "two.sided"))" or by making the test optional.

What I am doing : create a simple correlation matrix on 41 variables, then plot an heatmap with this program : library(gplots) mat<-cor(temp.alln,use="pairwise.complete.obs") hm<-heatmap.2(mat,symm=T) HM<-format(round(mat[hm[[1]],hm[[2]]],2)) library(RColorBrewer) brewer.pal(10,"Spectral")->mp heatmap.2(mat, symm = TRUE, col = mp,

I'm seeing an odd error message when running R CMD check gregmisc with R-1.9.0beta_2004-03-22.tar.gz * checking package dependencies ... WARNING Error in switch(type, code = c("R", "r", "S", "s", "q"), data = c("R", : Argument "type" is missing, with no default Execution halted See the information on DESCRIPTION

Hi, I'm trying to use the command "glm.nb" in library(MASS) to test for a significant difference in the aggregation parameter "theta" between the three levels of a factor. Any help gratefully received! Martin. Martin Hoyle, School of Life and Environmental Sciences, University of Nottingham, University Park, Nottingham, NG7 2RD, UK Webpage: