[R] Coefficient of Partial Determination
Nordlund, Dan (DSHS/RDA)
NordlDJ at dshs.wa.gov
Wed Mar 1 21:31:20 CET 2017
The rms and Hmisc packages can be found on CRAN. Just run install packages from the console menu.
Dan
Daniel Nordlund, PhD
Research and Data Analysis Division
Services & Enterprise Support Administration
Washington State Department of Social and Health Services
> -----Original Message-----
> From: R-help [mailto:r-help-bounces at r-project.org] On Behalf Of AbouEl-
> Makarim Aboueissa
> Sent: Wednesday, March 01, 2017 9:57 AM
> To: R-help at r-project.org
> Subject: [R] Coefficient of Partial Determination
>
> Dear All:
>
> Can the *Coefficient of Partial Determination* in multiple linear
> regression be computed in R? If so, could you please let me know how?
>
>
> fullmodel <- lm(Price ~ Size + Lotsize + Bedrooms + Bathrooms)
>
>
> *I found this in the internet, but I could find the package "rms"*
>
> *library(rms) # will also load Hmisc*
>
> *fit <- ols(y ~ x1 + x2, data=bf.dat) *
>
> *plt <- plot(anova(fit), what='partial R2')*
>
> *plt*
>
>
>
> Here is part of the data as an example:
> ---------------------------------------------------
>
> Taxes Bedrooms Bathrooms Price Size Lotsize
> 296 3 3 795000 2371 5850
> 242 4 3 399000 2818 4000
> 242 4 3 545000 3032 3060
> 222 4 4 909000 3540 6650
> 222 3 1 109900 1249 6360
> 222 3 3 324900 1800 4160
> 311 4 2 192900 1603 3880
> 311 3 2 215000 1450 4160
> 311 4 3 999000 3360 4800
> 311 3 2 319000 1323 5500
> 311 3 2 350000 1750 7200
> 311 3 2 249000 1400 3000
> 311 2 2 299000 1257 1700
> 307 3 2 235900 1400 2880
> 307 3 2 348000 1600 3600
> 307 4 3 314000 1794 3185
> 307 4 2 399000 1850 3300
> 307 3 3 599000 2950 5200
> 307 3 2 299000 1719 3450
> 307 3 3 425000 1472 3986
> 307 4 3 1100000 4168 4785
> 307 3 3 1500000 3880 4510
> 307 2 1 110000 1000 4000
> 307 3 2 200000 1139 3934
> 307 3 1 134900 1080 4960
> 307 4 3 250000 2000 3000
> 307 3 4 950000 1920 3800
> 307 4 2 239950 1348 4960
> 307 3 2 170000 1280 3000
> 307 3 2 285000 2400 4500
> 307 3 3 279000 1700 3500
> 307 3 2 219000 1600 3500
> 307 3 2 155000 1050 4000
> 307 3 2 389000 1415 4500
> 307 3 1 340000 1110 6360
> 279 2 1 95000 797 4500
> 279 2 2 140000 1100 4032
> 279 3 3 1100000 2602 5170
> 279 4 3 360000 2351 5400
> 252 3 3 415000 1350 3150
> 252 4 2 250000 1206 3745
> 233 3 3 559000 2628 4520
> 233 3 3 525000 2365 4640
> 233 3 3 779000 2990 8580
> 233 3 2 595000 1750 2000
> 233 4 5 1150000 5500 2160
> 233 3 2 550000 1852 3040
> 233 3 2 500000 2100 3090
> 233 4 3 279000 2580 4960
> 233 4 2 375000 1963 3350
> 243 3 3 330000 1900 5300
> 243 3 3 199000 1450 4100
> 243 2 2 165000 1000 9166
> 243 4 3 1399000 6500 4040
> 469 3 2 255000 1218 3630
> 226 2 2 325000 893 3620
>
>
>
> Thank you very much for your help and support
>
> abou
> ______________________
> AbouEl-Makarim Aboueissa, PhD
> University of Southern Maine
> Department of Mathematics and Statistics
>
> [[alternative HTML version deleted]]
>
> ______________________________________________
> R-help at r-project.org mailing list -- To UNSUBSCRIBE and more, see
> https://stat.ethz.ch/mailman/listinfo/r-help
> PLEASE do read the posting guide http://www.R-project.org/posting-
> guide.html
> and provide commented, minimal, self-contained, reproducible code.
More information about the R-help
mailing list