[R] DOE in R

Spencer Graves spencer.graves at pdf.com
Thu Aug 3 19:50:45 CEST 2006


<see in line>

Petr Pikal wrote:
> See
> 
> ?aov
> ?lm
> 
> something like
>  
> lm(result~yourfactor1+yourfactor2+yourfactor3+yourfactor5, 
> data=yourdataframe)
>
or

lm(result~(yourfactor1+yourfactor2+yourfactor3+yourfactor5)^2,
data=yourdataframe)

if you want all interactions.  Make sure all "yourfactors" are of class 
"factor".

	  Hope this helps.
	  Spencer Graves

  or
> aov(result~yourfactor1+yourfactor2+yourfactor3+yourfactor5, 
> data=yourdataframe)
> 
> but the exact structure of lm or aov construction depends on what you 
> want to test.
> 
> HTH
> Petr
> 
> 
> On 28 Jul 2006 at 21:43, dvrecko at sfu.ca wrote:
> 
> Date sent:      	Fri, 28 Jul 2006 21:43:38 -0700
> From:           	dvrecko at sfu.ca
> To:             	r-help at stat.math.ethz.ch
> Subject:        	[R] DOE in R
> Send reply to:  	dvrecko at sfu.ca
> 	<mailto:r-help-request at stat.math.ethz.ch?subject=unsubscribe>
> 	<mailto:r-help-request at stat.math.ethz.ch?subject=subscribe>
> 
>> Hi.
>>
>> I'm a student in a graduate program at Simon Fraser University in
>> Canada.
>>
>> I am trying to run a simple screening experiment with some simulated
>> data.
>>
>> I simply want to do an ANOVA of an experiemnt with 5 factors (4 have 2
>> levels, the last has 3 levels) and 48 runs (ie, full factorial).
>>
>> The thing is that I have multiple observations for each level
>> combination (run).
>>
>> So,
>>
>> 1) How do I do the anova based on the setup above?
>>
>> and 
>>
>> 2) More importantly, because of convergence issues for my simulations,
>> I will likely have an unequal number of observations for the 48 runs.
>> How can I do this?
>>
>> Seems like a straightforward enough situation.
>>
>> I am trying to avaoid writing my own C code to do the analysis since I
>> am working under some pretty tight time constraints.
>>
>> ANy help would be appreciated.
>>
>> Thanks very much.
>>
>> Dean Vrecko
>>
>> ______________________________________________
>> R-help at stat.math.ethz.ch mailing list
>> https://stat.ethz.ch/mailman/listinfo/r-help
>> PLEASE do read the posting guide
>> http://www.R-project.org/posting-guide.html and provide commented,
>> minimal, self-contained, reproducible code.
> 
> Petr Pikal
> petr.pikal at precheza.cz
> 
> ______________________________________________
> R-help at stat.math.ethz.ch mailing list
> https://stat.ethz.ch/mailman/listinfo/r-help
> PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
> and provide commented, minimal, self-contained, reproducible code.



More information about the R-help mailing list