[R] fractional balanced design conjoint analysis

S Ellison S.Ellison at LGCGroup.com
Thu Sep 13 18:45:56 CEST 2012


 

>> package BHH2 includes routines for 2-level fractional 
>> factorials with specified generators.
>> package FrF2 provides a different 2-level fractional 
>> factorial design toolkit.
> 
> Something tells me you won't get a balanced 2x2x2x2x3 
> fractional in 16 or fewer runs though. I think 18 is possible.

Apologies; should have noted that an 18 observation design without aliasing is possible _with additional levels_ in factors B-D.

But after re-engaging brain you can also manage a 16-run 2x2x2x2x3 that allows separation of main effects using Montgomery's trick: generate a 16-run 6-factor fractional factorial coded +-1, then add the first two columns to get a 7th, three-level, factor coded +1, 0, -1 and remove the first two columns to give a 2x2x2x2x3.

Coding for that using BHH2 would be

B<-ffDesMatrix(6, gen=list(c(5,1,2), c(6, 1,3)))
B <- cbind(B, rowSums(B[,1:2])/2)[,3:7]  #Generates a 3-level column and then cuts off the first two

#Alias check:
y <- rnorm(16)
alias(y~., data=as.data.frame(B))

This is free from main effect aliasing for a linear model, though whether the design (which is not ' balanced' in as much as the 3-level factor does not have all levels equally represented) will behave adequately for a conjoint analysis I'm afraid I can't say.

S Ellison

*******************************************************************
This email and any attachments are confidential. Any use...{{dropped:8}}




More information about the R-help mailing list