求助！广义相加模型！（简单sas程序老是报错） [推广有奖]

11楼

zyq7547182

发表于 2013-5-26 11:33:43

                                                                                                               LINP_
Obs year month day rain pressure wind  T  Tmax Tmin rh sun Case  P_Case    P_rain P_pressure P_wind P_sun    Case

  2536 2011 12 11 0 1028.2 0.9 1.5  9.8 -4.5 78 8.6  11  11.0816 .003314639  -0.004458 -0.004497  0.001095 2.40528
  2537 2011 12 12 0 1025.2 0.9 4.6 11.7  0.1 79 6.4  19  11.4616 .003314639  -0.040043 -0.004497  0.007407 2.43900
  2538 2011 12 13 0 1023.1 0.9 4.9 13.3 -0.1 79 7.4 7  12.0103 .003314639  -0.054170 -0.004497 -0.000246 2.48576
  2539 2011 12 14 0 1024.2 1.1 5.7 12.3 -0.1 85 1.4  15  11.4731 .003314639  -0.047978  0.003446  0.028692 2.44000
  2540 2011 12 15 0 1029.0 3.3 5.6  9.4  3.5 75 0.0  10 9.8483 .003314639 0.007345 -0.005554 -0.019852 2.28730
  2541 2011 12 16 0 1032.7 2.0 2.1  9.3 -2.9 76 8.0  11  10.4939 .003314639 0.069248  0.002346 -0.000752 2.35079
  2542 2011 12 17 0 1030.5 1.4 1.8  8.3 -3.4 75 3.6 9  10.7963 .003314639 0.031331  0.007100  0.040819 2.37920
  2543 2011 12 18 0 1028.5 1.0 4.9 11.1  1.0 68 4.4 8  10.9823 .003314639  -0.000125  0.000193  0.032029 2.39629
  2544 2011 12 19 0 1024.9 1.3 2.3 11.2 -3.2 81 7.3  15  11.6676 .003314639  -0.042652  0.006575  0.000105 2.45682
  2545 2011 12 20 0 1024.0 1.0 6.7 11.3  2.4 70 0.4  11  11.0164 .003314639  -0.049304  0.000193 -0.003435 2.39939
  2546 2011 12 21 0 1025.4 2.1 6.4  8.9  4.1 69 0.0 5  10.4857 .003314639  -0.038195  0.000421 -0.019852 2.35001
  2547 2011 12 22 0 1027.7 2.7 5.7 10.5  2.4 58 4.6  11  11.0402 .003314639  -0.011408 -0.003432  0.029634 2.40155
  2548 2011 12 23 0 1027.5 1.1 1.0  8.9 -3.6 76 5.9 9  11.1270 .003314639  -0.014088  0.003446  0.013438 2.40938
  2549 2011 12 24 0 1028.2 1.1 2.0 11.5 -4.3 72 7.9  12  11.0712 .003314639  -0.004458  0.003446 -0.000839 2.40434
  2550 2011 12 25 0 1028.7 1.0 3.1 10.7 -2.6 75 1.3 7  10.5701 .003314639 0.002826  0.000193  0.026269 2.35803
  2551 2011 12 26 0 1028.5 1.4 4.6  7.8  0.1 65 0.0 6  10.0596 .003314639  -0.000125  0.007100 -0.019852 2.30853
  2552 2011 12 27 0 1026.1 1.3 5.5  7.0  4.4 71 0.0 8  10.4282 .003314639  -0.031060  0.006575 -0.019852 2.34451
  2553 2011 12 28 0 1025.6 1.3 4.9 11.5 -1.0 76 6.4 7  11.4984 .003314639  -0.036261  0.006575  0.007407 2.44221
  2554 2011 12 29 0 1026.6 3.1 5.8  8.5  3.1 77 0.0  10  10.2162 .003314639  -0.025378 -0.003705 -0.019852 2.32397
  2555 2011 12 30 0 1029.3 1.5 4.8  9.4 -0.2 86 0.8 5  10.3520 .003314639 0.011967  0.007244  0.011407 2.33718
  2556 2011 12 31 0 1026.2 0.6 5.7 10.8  3.1 83 0.0 3  10.0900 .003314639  -0.029961 -0.025432 -0.019852 2.31154

书山有路勤为径，学海无涯苦作舟！

12楼

jingju11 发表于 2013-5-26 12:00:01

Not sure exactly what you need.
But
P_** is predicted ...;
ods output ParameterEstimates = XXX gives a table 'XXX' containing parameters estimates. By the way, that should be almost globally working on all procedures about regression models in SAS.

Jingju

13楼

zyq7547182

发表于 2013-5-26 12:06:04

类似于图片所示的结果

1.jpg (86.32 KB)

书山有路勤为径，学海无涯苦作舟！

14楼

zyq7547182

发表于 2013-5-26 12:14:27

jingju11 发表于 2013-5-26 12:00
Not sure exactly what you need.
But
P_** is predicted ...;

变量参数估计值标准误t值P值
Intercept 2.90476 0.03391 85.66 <0.0001
Linear(RH) -0.00178 0.00043 -4.17 <0.0001
Linear(T) 0.00585 0.00055 10.70 <0.0001
Linear(WIND) 0.00161 0.00678 0.24 0.8127
Linear(RAIN) 0.00008 0.00044 0.19 0.8487
Linear(DAY) 0.00000 0.00001 0.53 0.5987

书山有路勤为径，学海无涯苦作舟！

15楼

jingju11 发表于 2013-5-26 12:15:36

zyq7547182 发表于 2013-5-26 12:06
类似于图片所示的结果

all of those are SAS default outputs in listing window.
since SAS 9.3, a html ouput with htmlblue style was default. that is, you can get those tables automatically but with a little different style (if you did not reset your SAS output system).
If sas version was 9.1 and later, using ods like
ODS HTML ....;
proc ...;
ods html close;
Please refer to SAS doc for more details.
Jingju