楼主: patac_chong
10065 40

SAS可以做复杂试验设计吗? [推广有奖]

  • 0关注
  • 0粉丝

已卖:29份资源

高中生

95%

还不是VIP/贵宾

-

威望
0
论坛币
49 个
通用积分
0
学术水平
0 点
热心指数
2 点
信用等级
0 点
经验
158 点
帖子
39
精华
0
在线时间
35 小时
注册时间
2010-10-28
最后登录
2014-5-5

楼主
patac_chong 发表于 2010-10-28 15:40:03 |AI写论文

+2 论坛币
k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群
赵安豆老师微信:zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

感谢您参与论坛问题回答

经管之家送您两个论坛币!

+2 论坛币
想通过SPSS生成较大规模的试验设计表,总试验次数在200-300之间,变量数较多,水平数不同但都是在9以下。不知道是否能通过SPSS实现。
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

关键词:试验设计 SPSS PSS 大规模 不知道 规模

沙发
情迷仲夏夜 发表于 2010-10-30 04:46:07
做正交设计和正交分析其实有很多方法:
1。用下面的这个SAS程序可以生成3x2^2正交表:
proc plan;
   factors A=3 b=3 c=3;
   output out=fact A nvals=(0 1 2)
   b nvals=(0 1 2)
   c nvals=(0 1 2);
data rep1; set fact; rep=1; block=mod(A+b+c,3)+1;
data rep2; set fact; rep=2; block=mod(A+2*b+c,3)+1;
data rep3; set fact; rep=3; block=mod(A+b+2*c,3)+1;
data rep4; set fact; rep=4; block=mod(A+2*b+2*c,3)+1;
data bdesign; set rep1 rep2 rep3 rep4;
   if b<2 and c<2;
proc sort data=bdesign; by rep block;
proc print data=bdesign; run;
所以,我想你当然也可以根据这个程序写出你想要写出的更复杂的程序来的!

2。用SPSS也可以进行正交设计及正交分析:
做法如下: 设要做二因素的正交设计,A因素有三个水平,B因素有两个水平。则选择Data-->Orthogonal Design-->generate,弹出的就是正交设计窗口: Factor name框:输入A:单击ADD钮:单击Define value钮:分别在Value列的头三行输入1、2和3,单击continue钮,这样就定义好了变量A。
按类似的方法定义好变量B的2个水平。单击OK,系统就输出一个新定义的数据集,前两个变量就是要分析的A和B,各个水平已经按正交设计的要求排列好了。后面的status_和card_变量是系统产生的LOG变量,可以不管它。现在你再建立一个结果变量,输入实验结果,就可以进行正交设计的分析了。 正交设计的分析用GLM模块进行。具体操作如下: Analyze->General Linear Model->Univariate... dependent中选入应变量,fixed factor中选入自变量。然后进入model钮进行模型设置,这一步非常重要!设置模型为custom,然后选择需要分析的主效应和交互作用。然后确认,就可以得到所需要的结果。 请注意,如果model钮进行模型设置时选择错误,则得到的结果肯定是不正确的。

3。使用Excel也可以进行有交互作用的正交设计方差分析:
在Excel中进行方差分析方法如下:
1) 新建Excel空白工作薄,输入试验所得数据,其中A、B、C、D表示不同因素,A×B表示A、B因素间互作,1、2表示各因素中不同的水平,如图1。
2) K1、K2表示某因素中不同水平试验结果之和,即得率之和。使用公式“=SUMIF(range, criteria, [sum_range])”;如A因素中K1使用公式“=SUMIF(B$3:B$10,1,$I$3:$I$10)”,其中1表示水平为1。T2/n为校正值,可用:“=SUM(I3:I10)^2/8”。
3) A因素的Sj使用公式“=SUMSQ(B$11:B$12)/4-$I$12”;MSj使用公式“=B13/B14”;F使用公式“=B$15/$I$15”,分别将其复制到其他各因素的Sj、MSj、F中自动得到各因素的结果。
4) MSe使用公式“=SUM(F13:G13)/2”;F值显示因素A、B、A×B、D均不显著,将其与误差项e合并形成新的误差项均方MSe' “=(SUM(B13:D13) SUM(F13:H13))/6”。
5) F(α)为临界值用公式FINV(α,dfj,dfe)。
方差分析结果如图1所示。最后得到结果因素C即配比达到显著水平,应使用水平2即2:1的比例。
3 将此处理方法保存为模板
模板是Excel的一个重要方法,图1的工作表界面可以保存为Excel模板,方便以后的正交试验和方差分析。

藤椅
patac_chong 发表于 2010-11-1 09:32:56
很感谢您的帮助,还请继续赐教。
SPSS我是试过的,但是它生成的正交实验表我是怀疑的,比如需要一个3^1,5^2的3因素分表有3水平和5水平的实验设计,SPSS得到的是一个25行的表,25 不能 被3整除,所以肯定不是正交的,这个如何理解呢?

我想应用到的领域可能涉及到的影像因素很多(比如某一性能可能有10几个,需要同时考虑多种性能)
我想通过软件实现设计,通过输入变量数,及不同水平数,由软件生成一个最小的试验次数及正交表。

板凳
patac_chong 发表于 2010-11-1 09:49:22
弱弱的问一下:3x2^2是表示3行2变量2水平的正交试验吧?
上面提到的SAS程序中,A,b,c代表什么呢?
proc plan;
   factors A=3 b=3 c=3;
   output out=fact A nvals=(0 1 2)
   b nvals=(0 1 2)
   c nvals=(0 1 2);

报纸
patac_chong 发表于 2010-11-1 11:15:32
还请高手关注一下这个疑难问题啊

地板
patac_chong 发表于 2010-11-1 13:55:59
没有知道的同学路过吗?55~~~~~~~

7
yxj166 发表于 2010-11-1 17:07:00
有点难度!

8
情迷仲夏夜 发表于 2010-11-2 00:39:06
当然,这个问题不是一个简单的小问题!已经有很多人都在研究!就我知道:这个问题现在就可以分成有简单固定的(比如:2^3)和混合的之分(3x2^3),平衡的或不平衡的等等从简单到更复杂的。但话又说回来:问题是到底你是想要解决你现在的实际问题,是否真正需要一个非常复杂的正交设计表来完成你的试验和课题吗?还是想综合来探讨产生的方法。如果有谁想更深入研究这个议题,我现在上传一个附件用三种方法(SAS Macro, R 和MatLab)来探讨这个问题。请各位网友多多支持呦!

如果有人还想了解更多:能否用SAS系统界面来直接获取正交设计表,R Package Buildin正交列表,我还会继续陆续上传有关议题,方式,程序或附件。还有,如果有些人想要知道是否有现成软件可以作出正交表,就我现在知道网上有个软件叫“DOE++”(DOE就是Design of Experiments)就可以用简单和混合的方式来解决这个问题,原价可要$795美元,花这么高价,是否能做出非常复杂的正交表不得而知。不知道哪位高手能得到破译版为我们大家来分享!当然,还有另一个国产软件,叫“正交设计助手”,如果谁感兴趣,也不妨一试。

1.zip
下载链接: https://bbs.pinggu.org/a-785502.html

233 KB

需要: 2 个论坛币  [购买]

本附件包括:

  • v25i06-data.zip
  • paper.pdf
  • selc.sas
  • selc.R
  • selc.m

9
情迷仲夏夜 发表于 2010-11-2 03:52:37
上传一个好工具:正交设计表生成器!这个工具可以根据试验的具体需要从12种不同的模式中产生下列复杂的正交设计表!

L12(3^1 2^4)
1  1  0  1  1  
1  0  1  0  1  
0  0  0  1  0  
0  1  1  0  1  
1  1  1  1  0  
0  1  0  0  0  
2  1  1  1  0  
2  0  1  0  0  
2  1  0  0  1  
1  0  0  0  0  
0  0  1  1  1  
2  0  0  1  1  

L18(2^1 3^7)
2  0  2  2  1  1  1  0  
0  2  0  1  1  1  2  1  
0  1  2  2  0  0  2  0  
2  1  2  1  2  2  0  1  
1  2  1  2  2  2  2  0  
2  2  1  2  0  1  0  1  
0  2  2  0  0  2  1  1  
1  1  0  0  0  1  0  0  
2  1  1  0  1  0  2  1  
2  0  0  1  0  2  2  0  
0  0  0  2  2  0  0  1  
0  1  1  1  2  1  1  0  
0  0  1  0  1  2  0  0  
1  2  2  1  1  0  0  0  
2  2  0  0  2  0  1  0  
1  0  2  0  2  1  2  1  
1  1  0  2  1  2  1  1  
1  0  1  1  0  0  1  1  

L20(5^1 2^8)
2  1  1  1  0  0  0  0  1  
3  0  0  0  0  0  1  1  0  
3  1  1  0  1  1  1  0  1  
0  1  1  0  0  0  1  0  1  
2  1  0  0  1  0  0  1  0  
0  1  0  1  1  1  0  1  1  
3  1  1  1  1  0  0  1  0  
4  1  0  1  0  0  1  0  0  
4  0  1  0  0  0  0  1  1  
3  0  0  1  0  1  0  0  1  
2  0  0  1  0  1  1  1  1  
4  0  1  1  1  1  0  0  0  
1  0  0  0  1  0  0  0  1  
1  1  0  0  0  1  0  0  0  
0  0  0  1  1  0  1  0  0  
4  1  0  0  1  1  1  1  1  
0  0  1  0  0  1  0  1  0  
2  0  1  0  1  1  1  0  0  
1  0  1  1  1  0  1  1  1  
1  1  1  1  0  1  1  1  0  

L25(5^6)
1  0  2  0  2  0  
3  3  3  0  3  3  
2  3  4  1  2  2  
3  4  4  3  4  0  
0  2  3  1  0  0  
1  2  4  2  3  4  
3  0  0  1  1  4  
0  1  4  0  1  1  
4  0  4  4  0  3  
0  4  0  2  2  3  
4  4  2  1  3  1  
2  0  3  2  4  1  
2  4  1  0  0  4  
1  4  3  4  1  2  
0  3  2  4  4  4  
0  0  1  3  3  2  
1  1  1  1  4  3  
2  1  0  4  3  0  
3  1  2  2  0  2  
4  1  3  3  2  4  
4  3  1  2  1  0  
2  2  2  3  1  3  
3  2  1  4  2  1  
1  3  0  3  0  1  
4  2  0  0  4  2  

此外,还有下列其他正交设计表,不再一一举例,如果网友对此生成器有兴趣,不妨一试(附件中有3种具体安装说明):
L60(2^15 6^1 10^1)
L84(2^14 6^1 14^1)
L'6(3^1 2^3)
L'10(5^1 2^5)
L'12(4^1 3^4)
L'12(2^3 3^4)
L'12(6^1 2^5)
L'12(6^1 2^6)
L'12(3^1 2^9)
L'12(2^1 3^5)
L'12(3^2 2^7)
L'12(3^2 2^7)
L'12(2^5 3^3)
L'15(5^1 3^5)
L'18(2^1 3^8)
L'18(3^7 2^3)
L'18(9^1 2^8)
L'20(5^1 2^15)
L'24(8^1 3^8)
L'24(3^1 2^21)
L'24(6^1 2^15)
L'24(6^1 2^15)
L'24(6^1 2^18)
L'24(2^1 3^11)
L'24(3^1 4^7)
L'24(4^4 2^4)
L'84(2^8 6^1 14^1 3^2)
L'100(10^4 2^4 3^2)
S16(2^15)
S11(2^11)
S8(2^14)
S12(2^18)
S12(6^5 4^4)
S20(5^5 4^1)
S16(2^16 8^8)
S20(2^19 4^1 5^5)

G.zip

111.19 KB

需要: 2 个论坛币  [购买]

本附件包括:

  • Array.class
  • Bound.class
  • Cyclic.class
  • Design.class
  • Dialog.class
  • Gendex.class
  • Layout.class
  • Matrix.class
  • OMEP.class
  • Tools.class
  • Util.class
  • Utils.class
  • Vector.class
  • alpha.class
  • cibd.class
  • cut.class
  • feado.class
  • ibd.class
  • lhd.class
  • menu.class
  • miga.class
  • noa.class
  • rat.class
  • rcd.class
  • rrcd.class
  • sod.class
  • sudoku.class
  • install.txt

10
patac_chong 发表于 2010-11-3 12:42:41
非常感谢你的答复,看来高人出现了
我的实际问题是:手头有一些混合水平正交表,
正交试验助手我用过,还有Minitab,他们做不了这么复杂的表的。
SPSS据我这段时间的了解,好像也不可以。

您需要登录后才可以回帖 登录 | 我要注册

本版微信群
加好友,备注cda
拉您进交流群
GMT+8, 2025-12-20 22:56