1 2 Windows 7 3 Windows 4 R R Console R R Console ˆ R GUI R R R *5 R 2 R R R 6.1 ˆ 2 ˆ 2 ˆ Graphics Device 1 Rcmdr R Console R Rconsole R --sdi R M

Size: px

Start display at page:

Download "1 2 Windows 7 *3 Windows * 4 R R Console R R Console ˆ R GUI R R R *5 R 2 R R R 6.1 ˆ 2 ˆ 2 ˆ Graphics Device 1 Rcmdr R Console R Rconsole R --sdi R M"

あつみねてっちがわら
7 years ago
Views:

1 R John Fox and Milan Bouchet-Valat Version R Fox 2005 R R Core Team, 2013 GUI R R R R R R R R R the Comprehensive R Archive Network (CRAN) R CRAN 6.4 R Windows R Rcmdr Mac OS X Linux Unix R OS R R < <tinyurl.com/rcmdr> Windows R Rcmdr GUI R * 1 2 R R R Console library(rcmdr) Rcmdr R GUI Windows R R SDI *2 R Console R Fox(2005) [email protected] R Rcmdr [email protected] *1 ( ) *2 R Windows MDIR Console R 1 SDIR Console R SDI R etc 1

2 1 2 Windows 7 *3 Windows * 4 R R Console R R Console ˆ R GUI R R R *5 R 2 R R R 6.1 ˆ 2 ˆ 2 ˆ Graphics Device 1 Rcmdr R Console R Rconsole R --sdi R MDI R *3 ( ) Windows 8 *4 Rcmdr R Rcmdr Rcmdr Comprehensive R Archive Network (CRAN) < R CD-ROMWindows R GUI Rcmdr Rcmdr install.packages Rcmdr dependencies = TRUE Dirk Eddelbuettel Debian Linux $ apt-get install r-cran-rcmdr Rcmdr Linux Rcmdr Mac OS X Rcmdr tcltk X-Windows Tcl/Tk X-Windows R *5 R Console > R 2

3 2 R R Rcmdr Rcmdr Fox, 2007 Fox and Carvalho, R R R

4 * 6 Rcmdr R R R R R R R knitr knitr *6 R R Windows R dividers... 4

5 Rcmdr URL SPSS SAS Minitab STATA Excel Access dbase [32-bit Windows ] Excel [64-bit Windows ] () 5

6 Rcmdr F k- 6

7 Rcmdr QQ 3 3 PDF/Postscript/EPS RGL 7

8 Rcmdr AIC BIC RESET QQ 8

9 Rcmdr F F F F F 9

10 Rcmdr Rcmdr Rcmdr Rcmdr Commander R R Rcmdr R R R R 10

11 ˆ 2 R 2 * 7 2 R ˆ R GUI R Ctrl-r *8 Ctrl-Tab R Ctrl-a Ctrl-s ˆ R R R Console ˆ Rcmdr R Console R R R... R * 9 Rcmdr R 3 R * 10 1 R *7 David Firth relimp showdata 100R View R View 0 R *8 Ctrl Control r *9... GUI *10 11

12 R * 11 ˆ Mac OS X... ˆ ascii URL Minitab SPSS StataWindowsExcel Access dbase ˆ R 3.1 Nations.txt * 12 TFR contraception infant.mortality GDP region Afghanistan 6.90 NA Asia Albania 2.60 NA Europe Algeria Africa American-Samoa NA NA 11 NA Oceania Andorra NA NA NA NA Europe Angola 6.69 NA Africa Antigua NA Americas Argentina 2.62 NA Americas Armenia Europe Australia Oceania... ˆ 1 TFR 1 contraception ( ) infant.mortality 1000 GDP 1 US region ˆ R read.table ˆ R NA not available ˆ TFR contraception infant.mortality GDP region R region R Nations.txt R R URL... 3 *11 *12 Rcmdr etc 12

13 URL Nations 3 R. 0 9 R nations Nations NATIONS OK 4 Nations.txt R 5 R 5 Nations read.table showdatar R relimp library read.table R R R 3.2 R * 13 (2000) Problem 2.44 Moore ˆ R... Problem2.44 OK R *13 R Mac OS X URL... 13

14 4 5 ˆ 2 Enter 6 ˆ 1 var1 7 ˆ age Enter 2 height 14

15 8 ˆ R R * 14 9 * 15 R *14 R *15 R 15

16 8 2 9 car Prestige 4 R GUI R Nations Moore (2000) 5 car Prestige R 10 TFR contraception infant.mortality GDP 1 3 region 10 R R 11 region infant.mortality * 16 2 *

17 12 OK > numsummary(nations[,"infant.mortality"], statistics=c("mean", "sd", "IQR", + "quantiles"), quantiles=c(0,.25,.5,.75,1)) mean sd IQR 0% 25% 50% 75% 100% n NA sdiqr n NA R 11 OK * 17 R R OK R OK Nations 1 region OK... region 14 2 GDP infant.mortality OK > numsummary(nations[,c("gdp", "infant.mortality")], groups=nations$region, + statistics=c("mean", "sd", "IQR", "quantiles"), quantiles=c(0,.25,.5,.75,1)) Variable: GDP mean sd IQR 0% 25% 50% 75% 100% n NA Africa Americas Asia Europe Oceania Variable: infant.mortality mean sd IQR 0% 25% 50% 75% 100% n NA Africa Americas Windows Shift Ctrl *17 17

Asia 45.65854 32.980001 50.0 5 22.00 37.0 72.00 154 41 0 Europe 11.85366 7.122363 10.0 5 6.00 8.0 16.

18 Asia Europe Oceania R R R infant.mortality OK 16 18

19 region 2 1 * 18 *18 Windows Windows R Page Up Page Down Windows RGL Fox, 2003 Fox and Hong

20 15 16 Nations infant.mortality 5 R Venables and Ripley (2002) 2 nnet MASS 17 * 19 *19 R R 20

21 ... Prestige Prestige Prestige 3.3 car 17 ˆ ˆ ˆ ˆ log(income) ˆ LinearModel.1 ˆ R lm subset 1 TRUE FALSE type!= "prof" Prestige prof OK LinearModel.1 > LinearModel.1 <- lm(prestige ~ (education + log(income ))*type, data=prestige) > summary(linearmodel.1) Introduction to R R Console PDF 21

22 Call: lm(formula = prestige ~ (education + log(income)) * type, data = Prestige) Residuals: Min 1Q Median 3Q Max Coefficients: Estimate Std. Error t value Pr(> t ) (Intercept) e-08 *** education * log(income) e-08 *** type[t.prof] ** type[t.wc] education:type[t.prof] education:type[t.wc] * log(income):type[t.prof] * log(income):type[t.wc] Signif. codes: 0 *** ** 0.01 * Residual standard error: on 89 degrees of freedom (4 observations deleted due to missingness) Multiple R-squared: 0.871,Adjusted R-squared: F-statistic: on 8 and 89 DF, p-value: < 2.2e-16 Type II > Anova(LinearModel.1, type="ii") Anova Table (Type II tests) Response: prestige Sum Sq Df F value Pr(>F) education e-07 *** log(income) e-09 *** type ** education:type log(income):type * Residuals Signif. codes: 0 *** ** 0.01 *

23 6 6.1 R R R 18 R R R HTML R R * 20 R R R * 21 {r}r R R R {r echo=false}) R knitr Xie R  Replace with Main Title ======================= ### Your Name ### r as.character(sys.date()) {r echo=false} # include this code chunk as-is to set options opts_chunk$set(comment=na, prompt=true, out.width=750, fig.height=8, fig.width=8) library(rcmdr) {r} Nations <- read.table("c:/r/r-3.0.1patched/library/rcmdr/etc/nations.txt", header=true, sep="", na.strings="na", dec=".", strip.white=true) *20 *21 R 23

24 ... {r} data(prestige, package="car")... Let us regress occupational prestige on the education and income levels of the occupations, transforming income to linearize its relationship to prestige: {r} LinearModel.1 <- lm(prestige ~ (education + log(income))*type, data=prestige) summary(linearmodel.1) Anova(LinearModel.1, type="ii") Your Name Replace with Main Title R {r} R * *this is important* Let us regress occupational prestige... prestige... R R HTML HTML R R ( 18 R R R R R Ctrl-E R HTML OK R 6.2 R R R R 24

25 18 R R Word OpenOffice Writer Windows WordPad Ctrl-c Ctrl-v 1 R Courier New R Ctrl-v Ctrl-w R * 22 R R R R *22 Windows 25

26 6.3 R R R R R Windows Mac OS X RStudio < * R R R R CRAN R R R R Rcmdr R 6.5 R R R R R R Console R R [1] Firth, D. (2011). relimp: Relative Contribution of Effects in a Regression Model. R package version [2] Fox, J. (2003). Effect displays in R for generalised linear models. Journal of Statistical Software, 8(15):1-27. [3] Fox, J. (2005). The R Commander: A basic-statistics graphical user interface to R. Journal of Statistical Software, 19(9):1-42. [4] Fox, J. (2007). Extending the Rcmdr by plug-in Packages. R News, 7(3): [5] Fox, J. and Sá Carvalho, M. (2012). The RcmdrPlugin.survival package: Extending the R Commander to survival analysis. Journal of Statistical Software, 49(7):1-32. [6] Moore, D. S. (2000). The Basic Practice of Statistics. Freeman, New York, second edition. [7] Venables, W. N. and Ripley, B. D. (2002). Modern Applied Statistics with S. Springer, New York, fourth edition. *23 R RStudio R RStudio R RStudio 26

27 [8] Xie, Y. (2013). knitr: A general-purpose package for dynamic report generation in R. R package version

R John Fox R R R Console library(rcmdr) Rcmdr R GUI Windows R R SDI 1 R Console R 1 2 Windows XP Windows 2 R R Console R ˆ R

R John Fox R R R Console library(rcmdr) Rcmdr R GUI Windows R R SDI *1 R Console R 1 2 Windows XP Windows * 2 R R Console R ˆ R R John Fox 2006 8 26 2008 8 28 1 R R R Console library(rcmdr) Rcmdr R GUI Windows R R SDI *1 R Console R 1 2 Windows XP Windows * 2 R R Console R ˆ R GUI R R R Console > ˆ 2 ˆ Fox(2005) [email protected]