Re: [R] Aggregation of data frame with calculations of proportions

Dear Jim, On 26/06/07, jim holtman [EMAIL PROTECTED] wrote: I think something like this will do it for you. set.seed(1) n - 10 x - data.frame(a=sample(1:100,n), b=sample(1:100,n),d=sample(1:100,n)) x$a[c(1,5)] - NA x$b[c(7,6,4)] - NA x$d[c(5,8)] - NA x$total - rowSums(x, na.rm=TRUE)

[R] Is Landau-H available

Does anybody knows whether Landau-H index is available in ? I did not found anything about with r Regards Knut __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide

Re: [R] A really simple data manipulation example

For merging, selecting and aggregating, R is not too bad: I believe the following code is more or less equivalent to your Vilno example, isn't? # to input data, replace the following two lline by 'read.table' labresults - data.frame(patient.id=gl(10,5), visit.num=gl(5,10), sodium=rnorm(50)) demo

Re: [R] Looking for parallel functionality between Matlab and R

El-ad David Amir wrote: I'm slowly moving my statistical analysis from Matlab to R, and find myself missing two features: a) How do I mimic Matlab's 'hold on'? (I want to show several plots together, when I type two plots one after the other the second overwrites the first) b) How do I

Re: [R] Looking for parallel functionality between Matlab and R

This zooming function on the R-Wiki page was very neat: http://wiki.r-project.org/rwiki/doku.php?id=tips:graphics-misc:interactive_zooming Also, to answer question (a), maybe these examples might help? ## add elements to plot plot(1:10,1:10) lines(1:10,(1:10)/2) points(1:10,(1:10)/1.5) ## add

Re: [R] Strange RODBC problem

On 6/26/07, Richard M. Heiberger [EMAIL PROTECTED] wrote: I am using RODBC to collect data from an ODBC connection to an MS Access Database. Everything seems to be working well except datetimes between March 12, 2006 02:00 and 02:59 get moved one hour forward. This does not seem to be

[R] Another loop avoidance question.

Hi I want to sum over one of the dimensions of a n x k1 x k2 array in which each column is the product of the corresponding columns from two matrices with dimensions n x k1 and n x k2. I can see two approaches: a loop on k1 and a loop on k2. But I cannot figure a solution that avoids the loop? Is

Re: [R] Aggregation of data frame with calculations of proportions

I can drink a virtual beer in a bar in Second Life and put it on your tab. On 6/27/07, Mark Wardle [EMAIL PROTECTED] wrote: Dear Jim, On 26/06/07, jim holtman [EMAIL PROTECTED] wrote: I think something like this will do it for you. set.seed(1) n - 10 x -

[R] Import from excel

Hello, I have imported data from Excel using the command cours=read.delim(w:/apprentissage/cours_2.txt) after transforming my initial file with tab delimiters It seemed to work It is 2-dimensionnal. When I type cours[5,5], I get this type of message : [1] 0,9760942761824 Levels:

[R] cmdscale(eurodist)

Dear all, I have a question regarding the 'y - -loc[,2]' in ?cmdscale. Although I see that the plot is more sensible when using the '-loc' instead of just 'y - loc[,2]', I don't understand if there is a statistical reason to do '-loc[,2]'. So is this just to make the graph look better, or should

Re: [R] Import from excel

Bruce Willy wrote: Hello, I have imported data from Excel using the command cours=read.delim(w:/apprentissage/cours_2.txt) after transforming my initial file with tab delimiters It seemed to work It is 2-dimensionnal. When I type cours[5,5], I get this type of message : [1]

Re: [R] Import from excel

Thank you, but it still doesn't work completely Thanks to you and the dec=, option, I can now do cours[5,5]*5 and get the exact value But I still cannot do matrix operations like cours[2,]%*%cours[5,] It says the arguments are neither matrix nor vectors. :( Date: Wed, 27 Jun 2007

Re: [R] Import from excel

Hi [EMAIL PROTECTED] napsal dne 27.06.2007 15:43:53: Thank you, but it still doesn't work completely Thanks to you and the dec=, option, I can now do cours[5,5]*5 and get the exact value But I still cannot do matrix operations like cours[2,]%*%cours[5,] As only you have cours we can

Re: [R] Import from excel

Thank you, it works now I'm sorry, I have worked all year with R but with simulated data (except one time actually) Subject: RE: [R] Import from excel Date: Wed, 27 Jun 2007 10:01:51 -0400 From: [EMAIL PROTECTED] To: [EMAIL PROTECTED] cours is a dataframe not a matrix. Try:

[R] lme correlation structures

Hi all, I've been using SAS proc mixed to fit linear mixed models and would like to be able to fit the same models in R. Two things in particular: 1) I have longitudinal data and wish to allow for different repeated measures covariance parameter estimates for different groups (men and women),

Re: [R] lme correlation structures

Please read ?lme carefully -- the info you seek is there. In particular, the weights argument for changing variance weighting by covariates and the correlation argument for specifying correlation structures. Pinheiro and Bates's MIXED EFFECT MODELS IN S... is the canonical reference (which you

[R] [R-pkgs] New package tradeCosts

We would like to announce the availability of the 'tradeCosts' package in R for analysing transaction costs of trading. Version 0.1-0 is now available on CRAN. To take a look, you can: install.packages(tradeCosts) ... vignette(tradeCosts) and play around. This is the first and very basic

[R] moving-window (neighborhood) analysis

Hello all I was wondering what would be the best way to do a moving-window analysis of a matrix? By moving-window I mean that kind of analysis common in GIS, where each pixel (matrix element) of the resulting map is a function of it neighbors, and the neighborhood is a square matrix. I was hoping

[R] how to use chi-square to test correlation question

Hi There, There are 300 boy students and 100 girl students in a class. One interesting question is whether boy is smarter than girl or not. first given the exam with a difficulty level 1, the number of the student who got A is below 31 for boy, 10 for girl. Then we increase the difficulty

Re: [R] moving-window (neighborhood) analysis

See the Spatial section under CRAN's Task views Bert Gunter Genentech Nonclinical Statistics -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Carlos Guâno Grohmann Sent: Wednesday, June 27, 2007 8:27 AM To: r-help@stat.math.ethz.ch Subject: [R]

[R] xyplot with par

All, Is there are a simple way to plot multiple xyplots on the same page in the code below (it currently overwrites the first plot w/ the second). I searched the archives and saw a similar question but the answer didn't seem to work. thanks dave x1 = rnorm(10) x2 = rnorm (10) y1 = rnorm(10) y2 =

[R] CORRECTION: Re: BioC2007, August 6-7 ** early registration ends July 1 **

BioC2007: Where Software and Biology Connect August 6-7 in Seattle, WA, USA The BioC2007 developer-focused meeting is on August 8, the day *after* the main conference. In addition, I neglected to mention the poster session on Monday evening. Best, Martin Morgan -- Bioconductor / Computational

Re: [R] xyplot with par

On 6/27/07, Afshartous, David [EMAIL PROTECTED] wrote: All, Is there are a simple way to plot multiple xyplots on the same page in the code below (it currently overwrites the first plot w/ the second). I searched the archives and saw a similar question but the answer didn't seem to work.

[R] Solving a system of nonlinear equations

Hi All, Here is a modification of nlsolve() that I had submitted before. The new version of nlsolve() has the option to generate multiple random starting values in order to improve the chances of locating the zeros of the nonlinear system. The last test example in the attached file, the

[R] no applicable method

I'm getting started in R, and I'm trying to use one of the gradient boosting packages, mboost. I'm already installed the package with install.packages(mboost) and loaded it with library(mboost). My problem is that when I attempt to call glmboost, I get a message that Error in glmboost() :

[R] stepAIC on lm() where response is a matrix..

dear R users, I have fit the lm() on a mtrix of responses. i.e M1 = lm(cbind(R1,R2)~ X+Y+0). When i use summary(M1), it shows details for R1 and R2 separately. Now i want to use stepAIC on these models. But when i use stepAIC(M1) an error message comes saying that dropterm.mlm is not

Re: [R] how to use chi-square to test correlation question

someone else might have another viewpoint but I think your categories needs to be score categories ( 60-65, 65-70 etc ) and the data needs to be the number of girls and boys that fall into a category. I've never seen this below a certain value Methodology but maybe it's popular and I just

[R] Matlab end operator

Dear list members, I use both R and Matlab and find that each has its own strengths. Matlab definitely has the edge when it comes to the interactivity of its graphs. In addition I find the little operator end extremely useful in indexing arrays. (as in x(1:end,) ) The notation is MUCH more

[R] Timer

Hello, This might be a very basic question but I was not sure how to go about it. I just wanted to calcluate the time it takes to run my program. Basically I was to put a timer at the start and the end of the program in order to see how much time it takes to give the output (similar to tic-tac in

Re: [R] Timer

see ?system.time -- Henrique Dallazuanna Curitiba-Paraná-Brasil 25° 25' 40 S 49° 16' 22 O On 27/06/07, suman Duvvuru [EMAIL PROTECTED] wrote: Hello, This might be a very basic question but I was not sure how to go about it. I just wanted to calcluate the time it takes to run my program.

Re: [R] Strange RODBC problem

On 6/26/07, Richard M. Heiberger [EMAIL PROTECTED] wrote: I am using RODBC to collect data from an ODBC connection to an MS Access Database. Everything seems to be working well except datetimes between March 12, 2006 02:00 and 02:59 get moved one hour forward. This does not seem to be

[R] Re : Help please for graphics!

To whom it may concern, I want to plot two or more graphs in the same window by the means of the plot command. I have tried the option add=TRUE but this last does not work! Have you an hit for me please? Thank you very much for your attention. Bernard Colin [[alternative HTML

Re: [R] Matlab end operator

Hello Markus, On 6/27/07, Markus Loecher [EMAIL PROTECTED] wrote: Dear list members, I use both R and Matlab and find that each has its own strengths. Matlab definitely has the edge when it comes to the interactivity of its graphs. I also use both. R definitely has the edge when it comes to

Re: [R] Re : Help please for graphics!

par(mfrow=c(2,1)) plot(...plot1...) plot(plot2...) On 6/2/07, Bernard Colin [EMAIL PROTECTED] wrote: To whom it may concern, I want to plot two or more graphs in the same window by the means of the plot command. I have tried the option add=TRUE but this last does not work! Have you an

[R] read.xls problem

Hi to all, I have a strange problem There was a Excel file with 5 sheets I deleted sheet 2 and 4. as all the times before I tried to read the sheets with data1- read.xls(excelfile_2,sheet=1, as.is = TRUE ,verbose=FALSE, perl=C:/perl/bin/perl.exe) data2- read.xls(excelfile_2,sheet=2, as.is = TRUE

[R] Condensed PCA Results

Hello all, I'm currently using R to do PCA Analysis, and was wondering if anyone knew the specific R Code that could limit the output of the PCA Analysis so that you only get the Principal Component features as your output and none of the extraneous words or numbers that you don't want. If that

Re: [R] Condensed PCA Results

if x is your matrix, then prcomp(x)$rotation gets you what you want Katharina Wayne Aldo Gavioli schrieb: Hello all, I'm currently using R to do PCA Analysis, and was wondering if anyone knew the specific R Code that could limit the output of the PCA Analysis so that you only get the

[R] skeleton for C code?

Dear R experts---I would like to write a replacement for the read.csv function that is less general, but also more efficient. could someone please provide me with a skeleton function that shows me how to read the arguments and return a data frame for a call to a C function that handles

[R] running saved scripts

Hi All, I have a rather naive question. I need to run some simple calculations such as read.table, resid, and quantile on a data file. How can I save these commands in a text file and ask R to run this text file? Thanks. Jiong Zhang The email message (and any attachments) is for the sole use

Re: [R] running saved scripts

put into a file called abc.R or whatever, source(abc.R) check ?source On 6/27/07, Jiong Zhang, PhD [EMAIL PROTECTED] wrote: Hi All, I have a rather naive question. I need to run some simple calculations such as read.table, resid, and quantile on a data file. How can I save these commands

Re: [R] running saved scripts

Jiong Zhang, PhD schrieb: Hi All, I have a rather naive question. I need to run some simple calculations such as read.table, resid, and quantile on a data file. How can I save these commands in a text file and ask R to run this text file? Thanks. There is a very comfortable editor

[R] SEM model fit

I wonder if someone could explain why, when I perform confirmatory factor-analysis model using polychoric correlations why I do not get an estimated confidence interval for the RMSEA. My experience with these type models is that I would obtain a confidence interval estimate. I did not get any

[R] exaustive subgrouping or combination

Dear Colleagues, I am looking for a package or previous implemented R to subgroup and exaustively divide a vector of squence into 2 groups. For example: 1, 2, 3, 4 I want to have a group of 1, (2,3,4) (1,2), (3,4) (1,3), (2,4) (1,4), (2,3) (1,2,3), 4 (2,3), (1,4) ... Can someone help me as

Re: [R] running Rcmdr

Yes, of course, i know it, and i have installed all packages i need, but this is not what i am searching. I want to do that only with a command line without have to type only first R and after library(Rcmdr); i only want to type something like R library(Rcmdr) but i know it´s doesn´t work

[R] error message survreg.fit

Dear All, I am doing a parametric survival analysis with: fit - survreg(Surv(xyz$start, xyz$stop, xyz$event, type=interval) ~ 1, dist='loglogistic') At this point I do not want to look into covariates, hence the '~1' as model formulation. As event types I have exact, interval, and right

Re: [R] Matlab end operator

Hi Markus, Christophe I also use both matlab and R. I agree with Christophe: you can define the 'end' functionality by nrow or length also have a look at the following link that may be useful. http://mathesaurus.sourceforge.net/octave-r.html good luck AA. - Original Message - From:

Re: [R] running Rcmdr

Either one of these works under Linux. 'r' is our own littler (see Google for it), Rscript comes with R 2.5.0 or greater $ r -e 'suppressMessages(library(Rcmdr)); while (TRUE) Sys.sleep(1)' $ Rscript -e 'suppressMessages(library(Rcmdr)); while (TRUE) Sys.sleep(1)' Hth, Dirk -- Hell, there

Re: [R] running Rcmdr

I am very much a novice R user, and this is my first post to the List, but given that disclaimer, perhaps you could put library(Rcmdr) as a line in your Rprofile.site file? I think this contains commands that are run every time R is opened. [At least, I noticed after I installed my Tinn-R and

[R] Correlation ratio

Hi, I wanted to know how to compute the correlation ratio (eta) between two variables using R. Is there any function to compute the correlation ratio. Any help will be very much appreciated. Thanks, Suman [[alternative HTML version deleted]]

Re: [R] exaustive subgrouping or combination

Do you realize that for n items, there are 2^(n-1) such groups -- since you essentially want all possible subsets divided by 2: all possible subsets and their complements repeats each split twice, backwards and forwards. So this will quickly become ummm... rather large. If you really want to do

[R] Loading problem with XML_1.9

Hello all, I have loading problem with XML_1.9 under 64 bit R2.3.1, which I got from http://R.research.att.com/. XML_1.9 works fine under 32 bit R2.5.0. I thought that could be installation problem, and I tried install.packages or biocLite, every time the package installed fine, except some

[R] Gaussian elimination - singular matrix

Hello, I hope it is not a too stupid question. I have a singular matrix A (so not invertible). I want to find a nontrivial nonnegative solution to AX=0 (kernel of A) It is a special matrix A (so maybe this nonnegative solution is unique) The authors of the article suggest a Gaussian

Re: [R] Gaussian elimination - singular matrix

I am sorry, there is just a mistake : the solution cannot be unique (because it is a vectorial space) (but then I might normalize it) can R find one anyway ? This is equivalent to finding an eigenvector in fact From: [EMAIL PROTECTED] To: r-help@stat.math.ethz.ch Date: Wed, 27 Jun 2007

[R] How to discard data out of polr object

Hi, I would like to discard in-sample data that are stored in the polr object out after estimation. All I need are only the coefficients for estimation. I'm wondering if there is anyway to do this? Thank you. - adschai __ R-help@stat.math.ethz.ch

Re: [R] Correlation ratio

On Wed, 27 Jun 2007, suman Duvvuru wrote: Hi, I wanted to know how to compute the correlation ratio (eta) between two variables using R. Is there any function to compute the correlation ratio. Any help will be very much appreciated. help.search(correlation) ?cancor ?lm

[R] levelplot in lattice

Hi, I'm new to lattice. So please kindly be patient with me. I'm trying to arrange groups of levelplots into 3 rows as follows: Row1: Probabilities as functions of x and y, and conditioned on an event factor vector factor(a,b,c) Row2: Number of days as functions of x and y, and conditioned

Re: [R] exaustive subgrouping or combination

On Wed, 27 Jun 2007, Waverley wrote: Dear Colleagues, I am looking for a package or previous implemented R to subgroup and exaustively divide a vector of squence into 2 groups. For example: 1, 2, 3, 4 I want to have a group of 1, (2,3,4) (1,2), (3,4) (1,3), (2,4) (1,4), (2,3)

Re: [R] levelplot in lattice

On 6/27/07, [EMAIL PROTECTED] [EMAIL PROTECTED] wrote: Hi, I'm new to lattice. So please kindly be patient with me. I'm trying to arrange groups of levelplots into 3 rows as follows: Row1: Probabilities as functions of x and y, and conditioned on an event factor vector factor(a,b,c)

Re: [R] Gaussian elimination - singular matrix

All the nontrivial solutions to AX = 0 are the eigenvectors of A corresponding to eigenvalue 0 (try eigen function). The non-negative solution may or may not exist. For example, if A is a 2x2 matrix Aij = 1 for 1 =i,j =2 then the only non-trivial solution to AX = 0 is X = a*(1,-1), where a is any

[R] restructuring matrix

Hi all, let's say I have matrix PeopleDescValue Mary Height50 Mary Weight 100 FannyHeight 60 Fanny Height200 Is there a quick way to form the following matrix? People HeightWeight Mary 50 100 Fanny 60200 (Assuming I

Re: [R] restructuring matrix

If your original matrix is A then unique(A$People) and unique(A$Desc) will produce a vector of different people and a vector of different descriptions. --- yoo [EMAIL PROTECTED] wrote: Hi all, let's say I have matrix PeopleDescValue Mary Height50 Mary

[R] Sweave bug? when writing figures / deleting variable in chunk

I have found a quite strange (to me) behaviour in Sweave. It only occurs in the following situation: 1. define a variable in one chunk 2. use it within a subsequent figure-generating chunk 3. delete it at the end of that same chunk Then the Sweave driver chokes, not finding the variable name

Re: [R] restructuring matrix

Yea... let's say I constructed a matrix with rownames/colnames be those unique elements.. then what should I do? I don't want to do mapply, etc to find the field.. I'm wondering if there's a smarter way using row/col.. etc... Thanks! Moshe Olshansky-2 wrote: If your original matrix is A

Re: [R] Sweave bug? when writing figures / deleting variable in chunk

On 6/27/07, D G Rossiter [EMAIL PROTECTED] wrote: I have found a quite strange (to me) behaviour in Sweave. It only occurs in the following situation: 1. define a variable in one chunk 2. use it within a subsequent figure-generating chunk 3. delete it at the end of that same chunk Then the

Re: [R] Sweave bug? when writing figures / deleting variable in chunk

On 27/06/2007 9:06 PM, D G Rossiter wrote: I have found a quite strange (to me) behaviour in Sweave. It only occurs in the following situation: 1. define a variable in one chunk 2. use it within a subsequent figure-generating chunk 3. delete it at the end of that same chunk Then the

Re: [R] Sweave bug? when writing figures / deleting variable in chunk

On Wed, 27 Jun 2007 21:06:04 -0400 D G Rossiter wrote: I have found a quite strange (to me) behaviour in Sweave. It only occurs in the following situation: 1. define a variable in one chunk 2. use it within a subsequent figure-generating chunk 3. delete it at the end of that same chunk

Re: [R] Sweave bug? when writing figures / deleting variable in chunk

I have found a quite strange (to me) behaviour in Sweave. It only occurs in the following situation: You need to understand what Sweave does when it creates pictures: = sel - 1:5 @ fig=T= plot(trees[sel,]) rm(sel) @ By default, a eps and pdf version of the graphic is made. That is, this

[R] unequal variance assumption for lme (mixed effect model)

Dear Douglas and R-help, Does lme assume normal distribution AND equal variance among groups like anova() does? If it does, is there any method like unequal variance T-test (Welch T) in lme when each group has unequal variance in my data? Thanks, Shirley

Re: [R] restructuring matrix

Is this what you want? x - PeopleDescValue + Mary Height50 + Mary Weight 100 + FannyHeight 60 + Fanny Weight200 x - read.table(textConnection(x), header=TRUE, as.is=TRUE) reshape(x, direction='wide', idvar=People, timevar=Desc) People Value.Height

Re: [R] stepAIC on lm() where response is a matrix..

I see several options for you: 1. Write a function 'dropterm.mlm', copying 'dropterm.lm' and modifying it as you think appropriate. The function 'dropterm.lm' is hidden in a namespace, which you can see from 'methods(dropterm)'. To get it, either use getAnywhere(dropterm.lm) or

Re: [R] skeleton for C code?

On 6/27/07, ivo welch [EMAIL PROTECTED] wrote: Dear R experts---I would like to write a replacement for the read.csv function that is less general, but also more efficient. could someone please provide me with a skeleton function that shows me how to read the arguments and return a data frame

Re: [R] Gaussian elimination - singular matrix

RSiteSearch(generalized inverse, fun) produced 194 hits for me just now, including references to the following: Ginv {haplo.stats} ginv {MASS} invgen {far} Ginv {haplo.score} At least some of these should provide a solution to a singular system. You could further provide side constraints as

Re: [R] levelplot in lattice

Thank you Deepayan. Let's do the following exercise to reproduce the problem I'm facing. In the following code, I tried to have plot.rw1 and plot.rw2 shows in two separate rows on my plot windows. However, after the call to the 2nd levelplot, plot of plot.rw2 replace my first plot and the

Re: [R] levelplot in lattice

Sorry. My email editor from my ISP always screws up the text after sending out. Below is my response to you (in plain text). Thank you Deepayan. Let's do the following exercise to reproduce the problem I'm facing. In the following code, I tried to have plot.rw1 and plot.rw2 shows in two

Re: [R] epitools and R 2.5

--- Prof Brian Ripley [EMAIL PROTECTED] wrote: On Mon, 11 Jun 2007, Pietro Bulian wrote: At work after updating to R 2.5 I get an error using epitab from package epitools, when at home (R 2.4) I get no error. Could someone help me? The maintainer: this is a long-standing bug in the

Re: [R] unequal variance assumption for lme (mixed effect model)

The default settings for lme do assume equal variances within groups. You can change that by using the various varClasses. see ?varClasses. A simple example would be to allow unequal variances across groups. So if your call to lme was: lme(...,random=~1|group,...) then to allow each group to

Re: [R] levelplot in lattice

On 6/27/07, [EMAIL PROTECTED] [EMAIL PROTECTED] wrote: Sorry. My email editor from my ISP always screws up the text after sending out. Below is my response to you (in plain text). Thank you Deepayan. Let's do the following exercise to reproduce the problem I'm facing. In the following code, I

Re: [R] Correlation ratio

Hi Bruce, correlation ratio (eta) is different from correlation coefficient (rho). While correlation coefficient captures only a linear relationship btw variables, correlation ratio captures both linear and non-linear relationships. It is the defined as the ratio of the variance between arrays to

Re: [R] unequal variance assumption for lme (mixed effect model)

Hi Simon, Thanks for your reply. Your reply reminds me that book. I've read it long time ago, but haven't try the weights option in my projects yet:) Is the heteroscedastic test always less powerful because we have to estimate the within group variance from the given data? Should we check

Re: [R] stepAIC on lm() where response is a matrix..

On Wed, 27 Jun 2007, Spencer Graves wrote: I see several options for you: 1. Write a function 'dropterm.mlm', copying 'dropterm.lm' and modifying it as you think appropriate. The function 'dropterm.lm' is hidden in a namespace, which you can see from 'methods(dropterm)'. To get

Re: [R] Loading problem with XML_1.9

Please don't post to multiple lists: I have removed the BioC-devel list. This is about MacOS X, and the appropriate list is R-sig-mac. There is no intrinsic 64-bit problem: package XML 1.9-0 (sic) works fine on 64-bit versions of Solaris and Linux. Most likely there was an installation

[R] WEIBULL FRAILTY MODEL HELP

Dear R users, I try to write an rcode for a paranetric weibull model with unobserved heterogeneity. The data I have are giuven below. ( The data comes from Keiding and Klein artricle published in the Statistics in Medicine 1996) . While the simple weibull model runs perfectly ( I get

Re: [R] Correlation ratio

Suman, Try this: # some data example(aov) # the summary table anova(npk.aov) # extract sums of squares SS - anova(npk.aov)$Sum Sq SS #[1] 343.295 189.2816667 8.4016667 95.2016667 21.2816667 33.135 0.4816667 185.287 # eta-squared for factor N #SS factor N/SS Total