An alternative is to use the ffDesMatrix(..) procedure in the BHH2 
library (from CRAN mirrors) - it's drawn from _Statistics for 
Experimenters_ by Box, Hunter, and Hunter (hence the name). It will 
general a variety of design matrices - best for factorial and fractional 
factorial. It has a bunch of other useful experimental design goodies.

David Stevens

On 1/31/2012 12:34 PM, Sarah Goslee wrote:
> Which part? This might help:
> rep(1:5, each=9)
>
>
> On Tue, Jan 31, 2012 at 1:22 PM, Daniel Negusse
> <daniel.negusse@my.mcphs.edu>  wrote:
>>   can someone tell me if there is an easier way to do this in R - create a design matrix? thanks.
>> --->  design<- model.matrix(~ -1+factor(c(1,1,1,1,1,1,1,1,1,2,2,2,2,2,2,2,2,2,3,3,3,3,3,3,3,3,3,4,4,4,4,4,4,4,4,4,5,5,5,5,5,5,5,5,5,5,5,5,5))) # Creates an appropriate design matrix
>> can someone tell me if there is an easier way to do this matrix?
>> thanks.

-- 
David K Stevens, P.E., Ph.D., Professor
Civil and Environmental Engineering
Utah Water Research Laboratory
8200 Old Main Hill
Logan, UT  84322-8200
435 797 3229 - voice
435 797 1363 - fax
david.stevens@usu.edu




	[[alternative HTML version deleted]]

