spss分段回归的问题

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

这是试验过程中的一个结果，但是不会用spss实现，所以找了个sas程序做的，大家看看在spss里能不能实现？

以下是sas的程序：

1. DATA A;
   
2. INPUT aa ry;
   
3. cards;
   
4. 
   
5. 0        49.9
   
6. 75        52.1
   
7. 150        53.5
   
8. 225        53.9
   
9. 300        54.2
   
10. 375        53.7
    
11. 
    
12. 
    
13. 
    
14. 
    
15. 
    
16. 
    
17. 
    
18. 
    
19. 
    
20. ;
    
21. TITLE1 'Yield Response Data';
    
22. TITLE2 'LRP Model';
    
23. PROC NLIN DATA=A BEST=3;
    
24. PARMS A=49.9 to 54.2 by 1 B=0 to 1 by 0.1 K=0 to 375 by 15;
    
25. IF AA<K THEN DO;
    
26. MODEL RY = A + B*AA;
    
27. DER.A = 1;
    
28. DER.B = AA;
    
29. DER.K = 0;
    
30. END;
    
31. ELSE DO;
    
32. MODEL RY = A + B*K;
    
33. DER.A = 1;
    
34. DER.B = K;
    
35. DER.K = B;
    
36. END;
    
37. FILE PRINT;
    
38. IF _OBS_ = 1 AND _MODEL_= 0 THEN DO;
    
39. PLATEAU = A + B*K;
    
40. PUT PLATEAU=;
    
41. END;
    
42. PLATEAU = A + B*K;
    
43. ID PLATEAU;
    
44. OUTPUT OUT=NEW P=PRY PARMS=A B K SSE=SSE;
    
45. RUN;
    
46. PROC PLOT;
    
47. PLOT RY*AA='+' PRY*AA='*'/OVERLAY;
    
48. RUN;
    
49. /*   For you information:
    
50.      */
    
51. /*   This part of the program outputs Rsquares, significance,
    
52. plateau,  */
    
53. /*   and soil-test value at the joining point into a nice dataset.
    
54.     */
    
55. PROC MEANS NOPRINT;
    
56. VAR RY SSE A B K PLATEAU;
    
57. OUTPUT OUT=NEW2 N=TDF
    
58.                 MEAN=RY SSE A B K PLATEAU
    
59.                 CSS=CSST;
    
60. DATA NEW3; SET NEW2;
    
61. INTERCPT=A; SLOPE=B; JOINT=K;
    
62. RSQ=(CSST-SSE)/CSST;
    
63. EDF=TDF-3;
    
64. SSR=CSST-SSE;
    
65. MSR=SSR/2;
    
66. MSE=SSE/EDF;
    
67. F=MSR/MSE;
    
68. PROBF=1-(PROBF(F,2,EDF));
    
69. KEEP INTERCPT SLOPE JOINT PLATEAU RSQ F PROBF;
    
70. PROC PRINT;
    
71. RUN;

复制代码

扫码加我 拉你入群

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

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

分享0 收藏0 回帖

关键词：SPSS PSS response Title parms 程序

记得“蓝石”有一本书里提到利用SPSS做阶跃函数的回归。按理说，SPSS应该也能做楼主这类回归吧。

非常感谢，我已经做出来了。

应该可以的

学习下下

IF AA<K THEN DO;这个语句似乎有问题啊

sxdtlili 发表于 2012-4-19 06:05
非常感谢，我已经做出来了。

您好，我也想用spss做分段线性回归，您能指导一下吗？但是我的数据和您的不同，我的数据量大，而且是连续变量，拐点不好确定，请问用spss进行分段回归，拐点需要自己设定，还是会自动检测拐点？

您是怎么做到分段回归的呀，好厉害