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 > >