Re: [R] BBH2 and FrF2 packages

2010-06-30 Thread Andrea Bernasconi DG
; function is in both package BsMD > and FrF2. Try > > library(FrF2) > > and then run your code again; it worked for me... > > If you check the list of functions in BHH2 under HTML help, you'll find that > none of the plot functions you used below are found in that

[R] BBH2 and FrF2 packages

2010-06-24 Thread Andrea Bernasconi DG
Hi R HELP, I consider the 2^3 factorial experiment described at page 177 of the book Statistics for Experimenters: Design, Innovation, and Discovery by George E. P. Box, J. Stuart Hunter, William G. Hunter (BHH2). This example use the following data in file BHH2-Data/tab0502.dat at ftp://ftp.wile

Re: [R] Which is the easiest (most elegant) way to force "aov" to treat numerical variables as categorical ?

2010-06-14 Thread Andrea Bernasconi DG
I think I found the solution ! > cc<-factor(cars) > dd<-factor(driver) > MODEL<-y~cc+dd+additive > summary(aov(MODEL,data=DATA)) On 14 Jun, 2010, at 2:52 PM, Andrea Bernasconi DG wrote: > Hi R help, > > Hi R help, > > Which is the easiest (most elegant) w

[R] Which is the easiest (most elegant) way to force "aov" to treat numerical variables as categorical ?

2010-06-14 Thread Andrea Bernasconi DG
Hi R help, Hi R help, Which is the easiest (most elegant) way to force "aov" to treat numerical variables as categorical ? Sincerely, Andrea Bernasconi DG PROBLEM EXAMPLE I consider the latin squares example described at page 157 of the book: Statistics for Experimenters: Design,

[R] What is the best way to plots surfaces in 3 dimensions?

2010-04-28 Thread Andrea Bernasconi DG
Hi R help, What is the best way to plots surfaces in 3 dimensions? I also have the following availability problem with "plot3d" and "scatterplot3d", and "wireframe": install.packages("scatterplot3d") Warning: unable to access index for repository http://cran.ch.r-project.org/bin/macosx/leopard

Re: [R] How to obtain the coefficients from a summary of aov ?

2010-04-21 Thread Andrea Bernasconi DG
, David Winsemius wrote: > > On Apr 21, 2010, at 8:37 AM, Andrea Bernasconi DG wrote: > >> Dear Madame, Dear Sir, >> >> I am able to obtain the coefficients from a 'summary' of 'lm', but NOT from >> a 'summary' of 'aov'. &

[R] How to obtain the coefficients from a summary of aov ?

2010-04-21 Thread Andrea Bernasconi DG
Dear Madame, Dear Sir, I am able to obtain the coefficients from a 'summary' of 'lm', but NOT from a 'summary' of 'aov'. The following example shows my steps. ## Initialize rm(list = ls()) # remove (almost) everything in the working environment utils::data(npk, package="MASS") # get data model <