[R] Non-finite finite difference error

2008-12-23 Thread JS . Augustyn
Hello, I'm trying to use fitdistr() from the MASS package to fit a gamma distribution to a set of data. The data set is too large (1167 values) to reproduce in an email, but the summary statistics are: Min. 1st Qu. Median Mean 3rd Qu. Max. 116.7 266.7 666.7 1348.0 1642.0 16720.0 The call I'm

Re: [R] Standard error of mean for aov

2008-12-14 Thread JS . Augustyn
Thanks for all of your help, David, I finally got it. Here's some generic syntax in case it helps someone else down the road (using a 4-way ANOVA with repeated measures on all factors): # LOAD DATA data <- read.table("PATH\\datafile.txt") # RUN THE ANOVA data.aov <- aov(y ~ factor(x1)*factor(

Re: [R] Standard error of mean for aov

2008-12-13 Thread JS . Augustyn
> Does not this give you what you need? > model.tables(rawfixtimedata.aov,"means", se=TRUE) I tried that, but get an error: SEs for type 'means' are not yet implemented Maybe I'm not using the correct terminology to describe what I need to do. Using the main effect of Marking as an example, I h

[R] Standard error of mean for aov

2008-12-12 Thread JS . Augustyn
Hi all, I'm quite new to R and have a very basic question regarding how one gets the standard error of the mean for factor levels under aov. I was able to get the factor level means using: summary(print(model.tables(rawfixtimedata.aov,"means"),digits=3)), where rawfixtimedata.aov is my aov m