[交流] 【求助】SAS问题

> options(contrasts=c("contr.SAS", "contr.poly"))
> summary(lm(SOC ~ landform, soil))

Call:
lm(formula = SOC ~ landform, data = soil)

Residuals:
    Min      1Q  Median      3Q     Max
-0.6867 -0.6267 -0.1533  0.3608  1.2967

Coefficients:
                  Estimate Std. Error t value Pr(>|t|)   
(Intercept)         7.8967     0.3165  24.948 2.45e-10 ***
landformslopeland   2.0167     0.4476   4.505  0.00113 **
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

Residual standard error: 0.7753 on 10 degrees of freedom
Multiple R-squared: 0.6699,     Adjusted R-squared: 0.6369
F-statistic:  20.3 on 1 and 10 DF,  p-value: 0.001134

> summary(lm(SOC ~ landuse, soil))

Call:
lm(formula = SOC ~ landuse, data = soil)

Residuals:
   Min     1Q Median     3Q    Max
-1.685 -0.905  0.250  0.630  2.315

Coefficients:
                 Estimate Std. Error t value Pr(>|t|)   
(Intercept)        8.8950     0.5509  16.146 1.72e-08 ***
landusegrassland   0.0200     0.7791   0.026     0.98   
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

Residual standard error: 1.349 on 10 degrees of freedom
Multiple R-squared: 6.589e-05,  Adjusted R-squared: -0.09993
F-statistic: 0.0006589 on 1 and 10 DF,  p-value: 0.98

> anova(aov(SOC ~ landform, soil))
Analysis of Variance Table

Response: SOC
          Df  Sum Sq Mean Sq F value   Pr(>F)   
landform   1 12.2008 12.2008  20.297 0.001134 **
Residuals 10  6.0113  0.6011                    
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
> anova(aov(SOC ~ landuse, soil))
Analysis of Variance Table

Response: SOC
          Df  Sum Sq Mean Sq F value Pr(>F)
landuse    1  0.0012  0.0012   7e-04   0.98
Residuals 10 18.2109  1.8211               
> anova(aov(SOC ~ landform + landuse, soil))
Analysis of Variance Table

Response: SOC
          Df  Sum Sq Mean Sq F value   Pr(>F)   
landform   1 12.2008 12.2008 18.2706 0.002067 **
landuse    1  0.0012  0.0012  0.0018 0.967113   
Residuals  9  6.0101  0.6678                    
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

                                          The SAS System          20:33 Tuesday, June 20, 2009  21

                                        The REG Procedure
                                          Model: MODEL1
                                     Dependent Variable: SOC

                                       Analysis of Variance

                                              Sum of           Mean
          Source                   DF        Squares         Square    F Value    Pr > F

          Model                     1       12.20083       12.20083      20.30    0.0011
          Error                    10        6.01127        0.60113
          Corrected Total          11       18.21210


                       Root MSE              0.77532    R-Square     0.6699
                       Dependent Mean        8.90500    Adj R-Sq     0.6369
                       Coeff Var             8.70661


                                       Parameter Estimates

                                    Parameter       Standard
               Variable     DF       Estimate          Error    t Value    Pr > |t|

               Intercept     1        9.91333        0.31652      31.32      <.0001
               landform1     1       -2.01667        0.44763      -4.51      0.0011

                                          The SAS System          20:33 Tuesday, June 20, 2009  22

                                        The REG Procedure
                                          Model: MODEL1
                                     Dependent Variable: SOC

                                       Analysis of Variance

                                              Sum of           Mean
          Source                   DF        Squares         Square    F Value    Pr > F

          Model                     1        0.00120        0.00120       0.00    0.9800
          Error                    10       18.21090        1.82109
          Corrected Total          11       18.21210


                       Root MSE              1.34948    R-Square     0.0001
                       Dependent Mean        8.90500    Adj R-Sq    -0.0999
                       Coeff Var            15.15416


                                       Parameter Estimates

                                    Parameter       Standard
               Variable     DF       Estimate          Error    t Value    Pr > |t|

               Intercept     1        8.89500        0.55092      16.15      <.0001
               landuse1      1        0.02000        0.77912       0.03      0.9800

data ex;
input landform $        landuse $        SOC;
if landform eq 'tablelan' then landform1 = 1; else landform1 = 0;
if landform eq 'slopelan' then landform2 = 1; else landform2 = 0;
if landuse eq 'grasslan' then landuse1 = 1; else landuse1 = 0;
if landuse eq 'orchard' then landuse2 = 1; else landuse2 = 0;
datalines;
tableland        grassland        8.28
tableland        grassland        9.10
tableland        grassland        8.25
tableland        orchard        7.21
tableland        orchard        7.27
tableland        orchard        7.27
slopeland        orchard        10.20
slopeland        orchard        11.21
slopeland        orchard        10.21
slopeland  grassland      9.23
slopeland  grassland      9.31
slopeland  grassland    9.32
;
proc reg data = ex;
model SOC = landform1;
/* model SOC = landform2; */
run;
proc reg data = ex;
model SOC = landuse1;
/* model SOC = landuse2; */
run;
quit;