H:\public_html\sta6208\rpd>C:\R-2.9.0\bin\Rterm --vanilla
R version 2.9.0 (2009-04-17)
Copyright (C) 2009 The R Foundation for Statistical Computing
ISBN 3-900051-07-0
R is free software and comes with ABSOLUTELY NO WARRANTY.
You are welcome to redistribute it under certain conditions.
Type 'license()' or 'licence()' for distribution details.
Natural language support but running in an English locale
R is a collaborative project with many contributors.
Type 'contributors()' for more information and
'citation()' on how to cite R or R packages in publications.
Type 'demo()' for some demos, 'help()' for on-line help, or
'help.start()' for an HTML browser interface to help.
Type 'q()' to quit R.
> pdf("rpd12_1r.pdf")
>
> xy <- matrix( c(
+ 1,0,0.36,
+ 1,0,0.47,
+ 1,0,0.3,
+ 1,72,3.54,
+ 1,72,4.35,
+ 1,72,4.88,
+ 1,144,6.13,
+ 1,144,6.49,
+ 1,144,7.04,
+ 2,0,0.32,
+ 2,0,0.63,
+ 2,0,0.51,
+ 2,72,4.95,
+ 2,72,4.45,
+ 2,72,1.5,
+ 2,144,6.46,
+ 2,144,4.35,
+ 2,144,2.18,
+ 3,0,0,
+ 3,0,0.43,
+ 3,0,0.72,
+ 3,72,4.26,
+ 3,72,3.89,
+ 3,72,6.54,
+ 3,144,5.93,
+ 3,144,6.29,
+ 3,144,9.62,
+ 4,0,0.54,
+ 4,0,0.7,
+ 4,0,0.49,
+ 4,72,3.69,
+ 4,72,2.81,
+ 4,72,4.08,
+ 4,144,5.68,
+ 4,144,4.68,
+ 4,144,5.79,
+ 5,0,0.44,
+ 5,0,0.42,
+ 5,0,0.39,
+ 5,72,3.01,
+ 5,72,4.08,
+ 5,72,4.54,
+ 5,144,6.06,
+ 5,144,6.05,
+ 5,144,6.97,
+ 6,0,0.55,
+ 6,0,0.57,
+ 6,0,0.45,
+ 6,72,2.32,
+ 6,72,3.57,
+ 6,72,3.59,
+ 6,144,4.32,
+ 6,144,6.11,
+ 6,144,6.49,
+ 7,0,0.2,
+ 7,0,0.51,
+ 7,0,0.27,
+ 7,72,3.16,
+ 7,72,3.17,
+ 7,72,3.75,
+ 7,144,4.79,
+ 7,144,5.74,
+ 7,144,5.95,
+ 8,0,0.31,
+ 8,0,0.44,
+ 8,0,0.84,
+ 8,72,2.8,
+ 8,72,2.96,
+ 8,72,2.04,
+ 8,144,10.58,
+ 8,144,4.44,
+ 8,144,1.7), byrow=T,ncol=3)
>
> nobs <- nrow(xy)
> y <- xy[,3]
> x <- xy[,2]
> family <- factor(xy[,1])
>
> for (i in 1:nobs) {
+ if (y[i]==0) y[i] <- 0.01
+ if (x[i]==0) x[i] <- 1
+ }
>
> btreg.1 <- lm(y ~ family + x -1)
>
> summary(btreg.1)
Call:
lm(formula = y ~ family + x - 1)
Residuals:
Min 1Q Median 3Q Max
-3.89298 -0.54133 0.03935 0.49320 4.98702
Coefficients:
Estimate Std. Error t value Pr(>|t|)
family1 1.011978 0.454650 2.226 0.02961 *
family2 0.099756 0.454650 0.219 0.82704
family3 1.470867 0.454650 3.235 0.00194 **
family4 0.445312 0.454650 0.979 0.33110
family5 0.834200 0.454650 1.835 0.07125 .
family6 0.390867 0.454650 0.860 0.39321
family7 0.343089 0.454650 0.755 0.45329
family8 0.184200 0.454650 0.405 0.68674
x 0.037561 0.002522 14.893 < 2e-16 ***
---
Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
Residual standard error: 1.249 on 63 degrees of freedom
Multiple R-squared: 0.9213, Adjusted R-squared: 0.91
F-statistic: 81.94 on 9 and 63 DF, p-value: < 2.2e-16
>
> z <- 0.03756*x*log(x)
>
> btreg.2 <- lm(y ~ family + x +z -1)
>
> summary(btreg.2)
Call:
lm(formula = y ~ family + x + z - 1)
Residuals:
Min 1Q Median 3Q Max
-3.71333 -0.54549 0.09222 0.38316 5.16667
Coefficients:
Estimate Std. Error t value Pr(>|t|)
family1 0.77248 0.46728 1.653 0.10336
family2 -0.13974 0.46728 -0.299 0.76590
family3 1.23137 0.46728 2.635 0.01061 *
family4 0.20582 0.46728 0.440 0.66114
family5 0.59471 0.46728 1.273 0.20788
family6 0.15137 0.46728 0.324 0.74707
family7 0.10359 0.46728 0.222 0.82528
family8 -0.05529 0.46728 -0.118 0.90619
x 0.09488 0.03254 2.916 0.00493 **
z -0.30483 0.17254 -1.767 0.08220 .
---
Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
Residual standard error: 1.229 on 62 degrees of freedom
Multiple R-squared: 0.9251, Adjusted R-squared: 0.913
F-statistic: 76.54 on 10 and 62 DF, p-value: < 2.2e-16
>
> x1 <- x^.69517
>
> btreg.3 <- lm(y ~ family + x1 -1)
>
> summary(btreg.3)
Call:
lm(formula = y ~ family + x1 - 1)
Residuals:
Min 1Q Median 3Q Max
-3.69697 -0.55882 0.08404 0.39810 5.18303
Coefficients:
Estimate Std. Error t value Pr(>|t|)
family1 0.68157 0.45222 1.507 0.1368
family2 -0.23066 0.45222 -0.510 0.6118
family3 1.14046 0.45222 2.522 0.0142 *
family4 0.11490 0.45222 0.254 0.8003
family5 0.50379 0.45222 1.114 0.2695
family6 0.06046 0.45222 0.134 0.8941
family7 0.01268 0.45222 0.028 0.9777
family8 -0.14621 0.45222 -0.323 0.7475
x1 0.17512 0.01140 15.364 <2e-16 ***
---
Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
Residual standard error: 1.219 on 63 degrees of freedom
Multiple R-squared: 0.925, Adjusted R-squared: 0.9143
F-statistic: 86.39 on 9 and 63 DF, p-value: < 2.2e-16
>
> z1 <- 0.17512*x1*log(x1)
>
> btreg.4 <- lm(y ~ family + x1 +z1 -1)
>
> summary(btreg.4)
Call:
lm(formula = y ~ family + x1 + z1 - 1)
Residuals:
Min 1Q Median 3Q Max
-3.71333 -0.54549 0.09222 0.38316 5.16667
Coefficients:
Estimate Std. Error t value Pr(>|t|)
family1 0.71028 0.50434 1.408 0.1640
family2 -0.20194 0.50434 -0.400 0.6902
family3 1.16917 0.50434 2.318 0.0238 *
family4 0.14362 0.50434 0.285 0.7768
family5 0.53251 0.50434 1.056 0.2951
family6 0.08917 0.50434 0.177 0.8602
family7 0.04140 0.50434 0.082 0.9348
family8 -0.11749 0.50434 -0.233 0.8166
x1 0.15708 0.13613 1.154 0.2530
z1 0.02917 0.21934 0.133 0.8946
---
Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
Residual standard error: 1.229 on 62 degrees of freedom
Multiple R-squared: 0.9251, Adjusted R-squared: 0.913
F-statistic: 76.54 on 10 and 62 DF, p-value: < 2.2e-16
>
> x2 <- x1^1.02917
>
> btreg.5 <- lm(y ~ family + x2 -1)
>
> summary(btreg.5)
Call:
lm(formula = y ~ family + x2 - 1)
Residuals:
Min 1Q Median 3Q Max
-3.7133 -0.5455 0.0922 0.3832 5.1667
Coefficients:
Estimate Std. Error t value Pr(>|t|)
family1 0.70933 0.45136 1.572 0.1211
family2 -0.20289 0.45136 -0.450 0.6546
family3 1.16822 0.45136 2.588 0.0120 *
family4 0.14267 0.45136 0.316 0.7530
family5 0.53156 0.45136 1.178 0.2434
family6 0.08822 0.45136 0.195 0.8457
family7 0.04045 0.45136 0.090 0.9289
family8 -0.11844 0.45136 -0.262 0.7939
x2 0.15800 0.01028 15.367 <2e-16 ***
---
Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
Residual standard error: 1.219 on 63 degrees of freedom
Multiple R-squared: 0.9251, Adjusted R-squared: 0.9144
F-statistic: 86.41 on 9 and 63 DF, p-value: < 2.2e-16
>
> z2 <- 0.15800*x2*log(x2)
>
> btreg.6 <- lm(y ~ family + x2 +z2 -1)
>
> summary(btreg.6)
Call:
lm(formula = y ~ family + x2 + z2 - 1)
Residuals:
Min 1Q Median 3Q Max
-3.71333 -0.54549 0.09222 0.38316 5.16667
Coefficients:
Estimate Std. Error t value Pr(>|t|)
family1 7.094e-01 4.988e-01 1.422 0.1600
family2 -2.028e-01 4.988e-01 -0.407 0.6857
family3 1.168e+00 4.988e-01 2.342 0.0224 *
family4 1.427e-01 4.988e-01 0.286 0.7757
family5 5.316e-01 4.988e-01 1.066 0.2907
family6 8.828e-02 4.988e-01 0.177 0.8601
family7 4.050e-02 4.988e-01 0.081 0.9355
family8 -1.184e-01 4.988e-01 -0.237 0.8132
x2 1.580e-01 1.230e-01 1.284 0.2039
z2 6.179e-05 2.140e-01 0.000289 0.9998
---
Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
Residual standard error: 1.229 on 62 degrees of freedom
Multiple R-squared: 0.9251, Adjusted R-squared: 0.913
F-statistic: 76.54 on 10 and 62 DF, p-value: < 2.2e-16
>
> dev.off()
null device
1
>
>