在SPSS中能不能做多样本的趋势秩检验

相似文件 换一批

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

像这题：
一般来说，年龄越大的人的β脂蛋白的含量越大。现观察三组人，他们都是男性。第一组人的年龄在20到30岁之间，第二组人的年龄在30到40岁之间，第三组人的年龄在40到50岁之间。他们的β脂蛋白的测量值如下表，问这三组人的测量值是否符合人们的经验：年龄越大的人的β脂蛋白的含量越大？

第  1 组	260	200	240	170	270	205	190	200	250	200
第  2 组	310	310	190	225	170	210	280	210	280	240
第  3 组	320	260	360	310	270	380	240	295	260	250

扫码加我 拉你入群

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

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

分享0 收藏0 回帖

关键词：SPSS PSS 脂蛋白 测量值 检验 样本

自己顶下，在SPSS中有没有这个功能？老师让做下这个，实验要用，自己又忙着看考研书籍，心思都不在上面，SPSS不能的话，有没有办法SAS实现？

坐等高手回答，谢谢！

自己后面找资料弄了下，先在mathematica中算出一些必要的值，之后化为标准正态分布：

1. a = {};
   
2. Subscript[m, 1] = {260, 200, 240, 170, 270, 205, 190, 200, 250, 200};
   
3. Subscript[m, 2] = {310, 310, 190, 225, 170, 210, 280, 210, 280, 240};
   
4. Subscript[m, 3] = {320, 260, 360, 310, 270, 380, 240, 295, 260, 250};
   
5. t = Sort[Join[Subscript[m, 1], Subscript[m, 2], Subscript[m, 3]]];
   
6. Do[Subscript[b, i] = {}, {i, 1, 3}];
   
7. Do[Do[AppendTo[Subscript[b, i],
   
8.     Position[t, Subscript[m, i][[j]]]], {j, 1,
   
9.     Length[Subscript[m, i]]}], {i, 1, 3}];
   
10. Do[Subscript[t, i] =
    
11.    Table[Apply[Plus, Subscript[b, i][[j]]]/
    
12.     Length[Subscript[b, i][[j]]], {j, 1,
    
13.      Length[Subscript[m, i]]}], {i, 1, 3}];
    
14. m = Split[t];
    
15. a = {};
    
16. For[i = 1, i <= Length[m], i++,
    
17.   If[Length[m[[i]]] > 1, AppendTo[a, Length[m[[i]]]]
    
18.    ]];
    
19. Print[a];
    
20. Do[Subscript[n, i] = Length[Subscript[t, i]], {i, 1, 3}];
    
21. Do[Subscript[R, i] = Apply[Plus, Subscript[t, i]][[1]], {i, 1, 3}];
    
22. Do[Subscript[w, i] = 2 \!\(
    
23. \*UnderoverscriptBox[\(\[Sum]\), \(t = 1\), \(i\)]
    
24. \*SubscriptBox[\(n\), \(t\)]\) - Subscript[n, i], {i, 1, 3}];
    
25. Do[Print[{Subscript[n, i], Subscript[R, i], Subscript[w, i]}], {i, 1,
    
26.   3}]
    
27. T = \!\(
    
28. \*UnderoverscriptBox[\(\[Sum]\), \(i = 1\), \(3\)]\((
    
29. \*SubscriptBox[\(w\), \(i\)]*
    
30. \*SubscriptBox[\(R\), \(i\)])\)\)
    
31. n = \!\(
    
32. \*UnderoverscriptBox[\(\[Sum]\), \(i = 1\), \(3\)]
    
33. \*SubscriptBox[\(n\), \(i\)]\);
    
34. ET = n^2  (n + 1)/2
    
35. Do[Subscript[k, i] = (a[[i]])^3 - a[[i]], {i, 1, Length[a]}];
    
36. DT = N[(n*(n^2 - 1) - \!\(
    
37. \*UnderoverscriptBox[\(\[Sum]\), \(i = 1\), \(Length[a]\)]
    
38. \*SubscriptBox[\(k\), \(i\)]\))*(\!\(
    
39. \*UnderoverscriptBox[\(\[Sum]\), \(i = 1\), \(3\)]\(
    
40. \*SubscriptBox[\(n\), \(i\)]*
    
41. \*SubscriptBox[\(w\), \(i\)]*
    
42. \*SubscriptBox[\(w\), \(i\)]\)\) - n^3)/(12 (n - 1)), 9]
    
43. y = -Abs[(T - (ET))/Sqrt[(DT)]]

复制代码

然后将算出的y值代入SAS程序中（mathematica中不懂怎么算p值，话说回来，求p值和求概率一样吗？（有点迷糊，基础不好啊））：

1. data;
   
2. p=probnorm(-1.768463653);
   
3. put p=;
   
4. run;

复制代码

最后算出p值为0.038491711
（PS：做完这个后发现教科书上解答有错）
趋势的秩检验不是很主流，多样本一般都是用W-H检验，所以也没简单的菜单操作，只能通过编程实现了

希望后面有人重新改进下代码，弄成用SAS或mathematica一次性就能解决并且输出的。

在spss中怎么做，期待高手临场指导！！！！！！！！！

用spss做一下，描述看看