Daten Tschernobyl


Daten:
------

             ort regen     dist      bq
1        Pripjet     0     3.16 5300.00
2  Chistogalovka     0     5.62 2000.00
3          Lelev     0     7.94 2100.00
4    Tschernobyl     0    14.13 2000.00
5          Rudki     0    15.85  800.00
6       Orevichi     0    28.18 2000.00
7           Kiew     0    89.13   21.00
8    Tschernikow     0   125.89   55.00
9    Tscherkassy     0   281.84   12.00
10         Minsk     0   316.23   20.00
11        Donezk     0   707.95    6.00
12          Wien     0  1000.00    3.00
13   Oesterreich     1  1000.00   53.00
14     Stockholm     0  1122.02    1.50
15        Gaevle     1  1258.93   31.00
16    SuedBayern     1  1258.93   81.00
17      Konstanz     1  1584.89   31.00
18        Irland     1  1584.89   16.00
19     Stuttgart     0  1590.54    1.50
20       Chilton     0  1995.26    1.50
21    Schottland     1  1995.26   17.00
22         Japan     0 12589.25    0.15
23        Japan2     1 12589.25    0.85


Radioaktivitaet mit Distanz allein erklaeren (einfache Regression):
-------------------------------------------------------------------

> NurDistanz<-lm(log(bq)~log(dist))
> aov(NurDistanz)
Call:
   aov(formula = NurDistanz)

Terms:
                log(dist) Residuals
Sum of Squares  152.87486  34.39141
Deg. of Freedom         1        21

Residual standard error: 1.279721 
Estimated effects may be unbalanced
> summary(NurDistanz)

Call:
lm(formula = log(bq) ~ log(dist))

Residuals:
    Min      1Q  Median      3Q     Max 
-1.8592 -1.1348 -0.1033  1.1650  2.3797 

Coefficients:
            Estimate Std. Error t value Pr(>|t|)    
(Intercept)   9.8027     0.7022  13.960 4.26e-12 ***
log(dist)    -1.0911     0.1129  -9.662 3.53e-09 ***
---
Signif. codes:  0 `***' 0.001 `**' 0.01 `*' 0.05 `.' 0.1 ` ' 1 

Residual standard error: 1.28 on 21 degrees of freedom
Multiple R-Squared: 0.8164,     Adjusted R-squared: 0.8076 
F-statistic: 93.35 on 1 and 21 DF,  p-value: 3.527e-09 


Radioaktivitaet mit Regen allein erklaeren (einfache Regression):
-----------------------------------------------------------------

> NurRegen<-lm(log(bq)~ regen)
> aov(NurRegen)
Call:
   aov(formula = NurRegen)

Terms:
                    regen Residuals
Sum of Squares    3.31193 183.95434
Deg. of Freedom         1        21

Residual standard error: 2.959684 
Estimated effects may be unbalanced
> summary(NurRegen)

Call:
lm(formula = log(bq) ~ regen)

Residuals:
    Min      1Q  Median      3Q     Max 
-5.6755 -2.3332 -0.1205  2.1735  4.7970 

Coefficients:
            Estimate Std. Error t value Pr(>|t|)    
(Intercept)   3.7784     0.7399   5.107 4.66e-05 ***
regen        -0.8247     1.3412  -0.615    0.545    
---
Signif. codes:  0 `***' 0.001 `**' 0.01 `*' 0.05 `.' 0.1 ` ' 1 

Residual standard error: 2.96 on 21 degrees of freedom
Multiple R-Squared: 0.01769,    Adjusted R-squared: -0.02909 
F-statistic: 0.3781 on 1 and 21 DF,  p-value: 0.5452 


Volles Modell (zweifache Regression):
-------------------------------------

> FullModel<-lm(log(bq)~log(dist)+ regen)
> aov(FullModel)
Call:
   aov(formula = FullModel)

Terms:
                log(dist)     regen Residuals
Sum of Squares  152.87486  26.78701   7.60440
Deg. of Freedom         1         1        20

Residual standard error: 0.6166197 
Estimated effects may be unbalanced
> summary(FullModel)

Call:
lm(formula = log(bq) ~ log(dist) + regen)

Residuals:
     Min       1Q   Median       3Q      Max 
-1.37021 -0.37287 -0.05441  0.21542  1.61984 

Coefficients:
            Estimate Std. Error t value Pr(>|t|)    
(Intercept) 10.52246    0.34904  30.147  < 2e-16 ***
log(dist)   -1.36027    0.06316 -21.536 2.62e-15 ***
regen        2.72254    0.32436   8.394 5.51e-08 ***
---
Signif. codes:  0 `***' 0.001 `**' 0.01 `*' 0.05 `.' 0.1 ` ' 1 

Residual standard error: 0.6166 on 20 degrees of freedom
Multiple R-Squared: 0.9594,     Adjusted R-squared: 0.9553 
F-statistic: 236.3 on 2 and 20 DF,  p-value: 1.219e-14