spss语法程序如何编写？

libname a ‘e:\data\’

data a.student

infile ‘e:\data\student.txt’

input name height weight

以上程序将目录“e:\data\”下的文本文件“student.txt”中的数据输入数据集student中，该数据集存放于目录“e:\data\”下。

（2）直接输入方式

数据量较少或 *** 作者意志力坚强的情况下采用此种输入方式，在data语句之后写入如下语句：

input变量名1变量名2 …变量名n

datalines（在以前的版本下为cards，新版本下两者可通用）

… … … …（数据行）

… … … …（数据行）

… … … …（数据行）

datalines语句用于直接输入数据，标志着数据块的开始。

注意：这里的数据行中数据之间以空格分隔，当然也可以其它东东如逗号等来分隔，这里大家先以空格来分隔好了。因为不同的分隔方式下input语句要采取相应的控制选项，这些我们以后再讨论，这里我们还是省省力气吧。另外数据行输完后不能像其它语句那样直接在后面加上分号，而要另起一行输入分号，这样SAS才认为这是在输入原始数据而不是在搞别的什么。

SPSS编程

通过菜单选择：File-New-Syntax，打开语句编辑窗口(SyntaxEditor)编写程序(也可以在word或其他文本编辑软件中编写，通过复制粘贴方式将程序粘贴到此窗口)，若已建立程序，可通过File--+Open--+Syntax直接打开。

程序一（7组）：

data list free/Hc r1 r2 r3 r4 r5 r6 r7 N n1 n2 n3 n4 n5 n6 n7.

begin data

* 请于此行输入Kruskal-Wallis H 检验的结果（即Hc r1 r2 r3 r4 r5 r6 r7 N n1 n2 n3 n4 n5 n6 n7，尽量精确）

end data.

compute H=(12*((r1*n1)**2/n1+(r2*n2)**2/n2+(r3*n3)**2/n3+(r4*n4)**2/n4+(r5*n5)**2/n5+(r6*n6)**2/n6+(r7*n7)**2/n7))/(N*(N+1))-3*(N+1).

compute c=H/Hc.

compute x12=(r1-r2)**2/((N*(N+1)/12)*(1/n1+1/n2)*c).

compute x13=(r1-r3)**2/((N*(N+1)/12)*(1/n1+1/n3)*c).

compute x14=(r1-r4)**2/((N*(N+1)/12)*(1/n1+1/n4)*c).

compute x15=(r1-r5)**2/((N*(N+1)/12)*(1/n1+1/n5)*c).

compute x16=(r1-r6)**2/((N*(N+1)/12)*(1/n1+1/n6)*c).

compute x17=(r1-r7)**2/((N*(N+1)/12)*(1/n1+1/n7)*c).

compute x23=(r2-r3)**2/((N*(N+1)/12)*(1/n2+1/n3)*c).

compute x24=(r2-r4)**2/((N*(N+1)/12)*(1/n2+1/n4)*c).

compute x25=(r2-r5)**2/((N*(N+1)/12)*(1/n2+1/n5)*c).

compute x26=(r2-r6)**2/((N*(N+1)/12)*(1/n2+1/n6)*c).

compute x27=(r2-r7)**2/((N*(N+1)/12)*(1/n2+1/n7)*c).

compute x34=(r3-r4)**2/((N*(N+1)/12)*(1/n3+1/n4)*c).

compute x35=(r3-r5)**2/((N*(N+1)/12)*(1/n3+1/n5)*c).

compute x36=(r3-r6)**2/((N*(N+1)/12)*(1/n3+1/n6)*c).

compute x37=(r3-r7)**2/((N*(N+1)/12)*(1/n3+1/n7)*c).

compute x45=(r4-r5)**2/((N*(N+1)/12)*(1/n4+1/n5)*c).

compute x46=(r4-r6)**2/((N*(N+1)/12)*(1/n4+1/n6)*c).

compute x47=(r4-r7)**2/((N*(N+1)/12)*(1/n4+1/n7)*c).

compute x56=(r5-r6)**2/((N*(N+1)/12)*(1/n5+1/n6)*c).

compute x57=(r5-r7)**2/((N*(N+1)/12)*(1/n5+1/n7)*c).

compute x67=(r6-r7)**2/((N*(N+1)/12)*(1/n6+1/n7)*c).

compute p12=1-cdf.chisq(x12,2).

compute p13=1-cdf.chisq(x13,2).

compute p14=1-cdf.chisq(x14,2).

compute p15=1-cdf.chisq(x15,2).

compute p16=1-cdf.chisq(x16,2).

compute p17=1-cdf.chisq(x17,2).

compute p23=1-cdf.chisq(x23,2).

compute p24=1-cdf.chisq(x24,2).

compute p25=1-cdf.chisq(x25,2).

compute p26=1-cdf.chisq(x26,2).

compute p27=1-cdf.chisq(x27,2).

compute p34=1-cdf.chisq(x34,2).

compute p35=1-cdf.chisq(x35,2).

compute p36=1-cdf.chisq(x36,2).

compute p37=1-cdf.chisq(x37,2).

compute p45=1-cdf.chisq(x45,2).

compute p46=1-cdf.chisq(x46,2).

compute p47=1-cdf.chisq(x47,2).

compute p56=1-cdf.chisq(x56,2).

compute p57=1-cdf.chisq(x57,2).

compute p67=1-cdf.chisq(x67,2).

execute.

复制代码

注：少于7组的，删掉相应的行与变量即可。

******

举例3组：

data list free

/Hc r1 r2 r3 N n1 n2 n3.

begin data

9.94 8.40 18.78 19.27 30 10 9 11

end data.

COMPUTE H = (12*((r1*n1)**2/n1+(r2*n2)**2/n2+(r3*n3)**2/n3))/(N*(N+1))-3*(N+1).

compute C = H/Hc.

compute x12 = (r1-r2)**2/((N*(N+1)/12)*(1/n1+1/n2)*c).

compute x13 = (r1-r3)**2/((N*(N+1)/12)*(1/n1+1/n3)*c).

compute x23 = (r2-r3)**2/((N*(N+1)/12)*(1/n2+1/n3)*c).

compute p12 = 1-cdf.chisq(x12,2).

compute p13 = 1-cdf.chisq(x13,2).

compute p23 = 1-cdf.chisq(x23,2).

execute.

结果：

举例4组：

data list free

/Hc r1 r2 r3 r4 N n1 n2 n3 n4.

begin data

27.625 44 33.46 23.29 13.25 56 14 14 14 14

end data.

COMPUTE H = (12*((r1*n1)**2/n1+(r2*n2)**2/n2+(r3*n3)**2/n3)+(r4*n4)**2/n4)/(N*(N+1))-3*(N+1).

compute C = H/Hc.

compute x12 = (r1-r2)**2/((N*(N+1)/12)*(1/n1+1/n2)*c).

compute x13 = (r1-r3)**2/((N*(N+1)/12)*(1/n1+1/n3)*c).

compute x14 = (r1-r4)**2/((N*(N+1)/12)*(1/n1+1/n4)*c).

compute x23 = (r2-r3)**2/((N*(N+1)/12)*(1/n2+1/n3)*c).

compute x24 = (r2-r4)**2/((N*(N+1)/12)*(1/n2+1/n4)*c).

compute x34 = (r3-r4)**2/((N*(N+1)/12)*(1/n3+1/n4)*c).

compute p12 = 1-cdf.chisq(x12,2).

compute p13 = 1-cdf.chisq(x13,2).

compute p14 = 1-cdf.chisq(x14,2).

compute p23 = 1-cdf.chisq(x23,2).

compute p24 = 1-cdf.chisq(x24,2).

compute p34 = 1-cdf.chisq(x34,2).

execute.

结果：

---