presentation 1 r
play

Presentation 1: R Murray Logan 10 Apr 2016 Schedule - 1 week - PowerPoint PPT Presentation

Presentation 1: R Murray Logan 10 Apr 2016 Schedule - 1 week course Day AM PM Monday Intro to R Intro to hypothesis testing Tuesday Simple Multiple regression regression Wednesday One-way ANOVA Factorial ANOVA Thursday Mixed


  1. [628] "xkcd.ttf" "pres.7.5a_notes.log" [412] "pres.7.4a.toc" "pres.7.4a.vrb" "pres.7.5a.aux" [415] "pres.7.5a.html" "pres.7.5a.log" "pres.7.5a.nav" [418] "pres.7.5a_notes.aux" "pres.7.5a_notes.out" "pres.7.4a.tex~" [421] "pres.7.5a_notes.pdf" "pres.7.5a_notes.tex" "pres.7.5a_notes.toc" [424] "pres.7.5a.out" "pres.7.5a.pdf" "pres.7.5a.R" [427] "#pres.7.5a.Rmd#" "pres.7.4a_tmp.md" [409] "pres.7.4a.tex" "pres.7.5a.Rmd~" "pres.7.4a_notes.pdf" [391] "pres.7.4a.aux" "pres.7.4a.html" "pres.7.4a.log" [394] "pres.7.4a.nav" "pres.7.4a_notes.aux" "pres.7.4a_notes.log" [397] "pres.7.4a_notes.out" "#pres.7.4a_notes.tex#" "pres.7.4a.snm" [400] "pres.7.4a_notes.tex" "pres.7.4a_notes.tex~" "pres.7.4a_notes.toc" [403] "pres.7.4a.out" "pres.7.4a.pdf" "pres.7.4a.R" [406] "pres.7.4a.Rmd" "pres.7.4a.Rmd~" "pres.7.5a.Rmd" [430] "pres.7.5a.snm" "pres.7.3a.toc" "pres.8.2a.nav" "pres.7.6a_tmp.md" "pres.7.6a.toc" [454] "pres.7.6a.vrb" "pres.8.2a.aux" "pres.8.2a.html" [457] "pres.8.2a.log" "pres.8.2a.md" [460] "pres.8.2a_notes.aux" "pres.7.6a.snm" "pres.8.2a_notes.log" "pres.8.2a_notes.out" [463] "pres.8.2a_notes.pdf" "pres.8.2a_notes.tex" "pres.8.2a_notes.tex~" [466] "pres.8.2a_notes.toc" "pres.8.2a.out" [451] "pres.7.6a.tex" "pres.7.6a.Rmd~" "pres.7.5a.tex" [439] "pres.7.6a_notes.aux" "pres.7.5a_tmp.md" [433] "pres.7.5a.toc" "pres.7.5a.vrb" "pres.7.6a.aux" [436] "pres.7.6a.html" "pres.7.6a.log" "pres.7.6a.nav" "pres.7.6a_notes.log" [448] "pres.7.6a.Rmd" "pres.7.6a_notes.out" [442] "pres.7.6a_notes.pdf" "pres.7.6a_notes.tex" "pres.7.6a_notes.toc" [445] "pres.7.6a.out" "pres.7.6a.pdf" "pres.7.6a.R" "pres.7.3a.vrb" [388] "pres.7.3a_tmp.md" [469] "pres.8.2a.R" "pres.5.2.R" [331] "pres.5.2_notes.log" "pres.5.2_notes.out" "pres.5.2_notes.pdf" [334] "pres.5.2_notes.tex" "pres.5.2_notes.toc" "pres.5.2.out" [337] "pres.5.2.pdf" "#pres.5.2.Rmd#" "pres.5.2.nav" [340] "pres.5.2.Rmd" "pres.5.2.Rmd~" "pres.5.2.snm" [343] "pres.5.2.tex" "pres.5.2.tex~" "pres.5.2_tmp.md" [346] "pres.5.2.toc" "pres.5.2_notes.aux" [328] "pres.5.2.md" "pres.7.2a.aux" "pres.4.pdf" "pres.3.2.R" "pres.3.2.Rmd" [313] "pres.3.2.Rmd~" "pres.4.html" "pres.4.html~" [316] "pres.4_notes.pdf" "pres.4_notes.tex~" [319] "pres.4pdf_knit_.Rmd~" "pres.5.2.log" "pres.4.R" "pres.4.Rmd" [322] "pres.4.Rmd~" "pres.4.tex~" "pres.4_tmp.md~" [325] "pres.5.2.aux" "pres.5.2.html" "pres.5.2.vrb" [349] "pres.7.2a.html" "pres.7.3a.tex" "pres.7.3a_notes.tex" "pres.7.3a.html" "pres.7.3a.log" [373] "pres.7.3a.nav" "pres.7.3a_notes.aux" "pres.7.3a_notes.log" [376] "pres.7.3a_notes.out" "pres.7.3a_notes.pdf" [379] "pres.7.3a_notes.toc" "pres.7.2a.vrb" "pres.7.3a.out" "pres.7.3a.pdf" [382] "pres.7.3a.R" "#pres.7.3a.Rmd#" "pres.7.3a.Rmd" [385] "pres.7.3a.Rmd~" "pres.7.3a.snm" [370] "pres.7.3a.aux" "pres.7.2a.toc" "pres.7.2a.log" [358] "pres.7.2a.out" "pres.7.2a.nav" [352] "pres.7.2a_notes.aux" "pres.7.2a_notes.log" "pres.7.2a_notes.out" [355] "pres.7.2a_notes.pdf" "pres.7.2a_notes.tex" "pres.7.2a_notes.toc" "pres.7.2a.pdf" [367] "pres.7.2a_tmp.md" "pres.7.2a.R" [361] "#pres.7.2a.Rmd#" "pres.7.2a.Rmd" "pres.7.2a.Rmd~" [364] "pres.7.2a.snm" "pres.7.2a.tex" "pres.7.2a.tex~" "pres.8.2a.pdf" "#pres.8.2a.Rmd#" "pres.3.2.html" "pres.9.4a.tex" "pres.9.4a_notes.toc" "pres.9.4a.out" [574] "pres.9.4a.pdf" "pres.9.4a.R" "#pres.9.4a.Rmd#" [577] "pres.9.4a.Rmd" "pres.9.4a.snm" [580] "pres.9.4a_tmp.md" "pres.9.4a_notes.pdf" "pres.9.4a.toc" "pres.9.4a.vrb" [583] "pres.b1.aux" "pres.b1.html" "pres.b1.log" [586] "pres.b1_notes.aux" "pres.b1_notes.log" [571] "pres.9.4a_notes.tex" "pres.9.4a_notes.out" [589] "pres.b1_notes.pdf" [559] "pres.9.3a.tex" "pres.9.3a.out" [553] "pres.9.3a.pdf" "pres.9.3a.R" "#pres.9.3a.Rmd#" [556] "pres.9.3a.Rmd" "pres.9.3a.Rmd~" "pres.9.3a.snm" "pres.9.3a_tmp.md" [568] "pres.9.4a_notes.log" "pres.9.3a.toc" [562] "pres.9.3a.vrb" "pres.9.4a.aux" "pres.9.4a.html" [565] "pres.9.4a.log" "pres.9.4a.nav" "pres.9.4a_notes.aux" "pres.b1_notes.out" "pres.b1_notes.tex" [550] "pres.9.3a_notes.tex" [619] "quiz.toc" "quiz.pdf" [613] "quiz.R" "quiz.Rmd" "quiz.Rmd~" [616] "quiz.snm" "quiz.tex" "quiz_tmp.md" "*shell*" [610] "quiz_notes.toc" "*shell*~" [622] "#*shell*#" "#test.bugs#" "test.bugs" [625] "test.png" "texput.log" "xkcd.png" "quiz.out" "quiz_notes.tex" "pres.b1_notes.toc" "pres.b1_tmp.md" [592] "pres.b1.out" "pres.b1.pdf" "pres.b1.R" [595] "#pres.b1.Rmd#" "pres.b1.Rmd" "pres.b1.Rmd~" [598] "pres.b1.tex" "pres.b1.vrb" "quiz_notes.pdf" [601] "quiz.aux" "quiz.html" "quiz.log" [604] "quiz.nav" "quiz_notes.aux" "quiz_notes.log" [607] "quiz_notes.out" "pres.9.3a_notes.toc" "pres.9.3a_notes.pdf" "pres.8.2a.Rmd" [499] "pres.8.3a.vrb" "#pres.8.3a.Rmd#" [493] "pres.8.3a.Rmd" "pres.8.3a.Rmd~" "pres.8.3a.snm" [496] "pres.8.3a.tex" "pres.8.3a_tmp.md" "pres.8.3a.toc" "pres.9.1a.Rmd~" [490] "pres.8.3a.pdf" "pres.9.1.aux" [502] "pres.9.1.html" "pres.9.1.log" "pres.9.1.md" [505] "pres.9.1.nav" "pres.9.1_notes.aux" "pres.9.1_notes.log" "pres.8.3a.R" "pres.8.3a.out" "pres.9.1_notes.pdf" "pres.8.3a.html" [472] "pres.8.2a.Rmd~" "pres.8.2a.snm" "pres.8.2a.tex" [475] "pres.8.2a_tmp.md" "pres.8.2a.toc" "pres.8.2a.vrb" [478] "pres.8.3a.aux" "pres.8.3a.log" "pres.8.3a_notes.toc" [481] "pres.8.3a.md" "pres.8.3a.nav" "pres.8.3a_notes.aux" [484] "pres.8.3a_notes.log" "pres.8.3a_notes.out" "pres.8.3a_notes.pdf" [487] "pres.8.3a_notes.tex" [508] "pres.9.1_notes.out" "pres.9.1_notes.tex" "pres.9.3a_notes.out" "pres.9.2a_tmp.md" [532] "pres.9.2a.out" "pres.9.2a.pdf" "pres.9.2a.R" [535] "pres.9.2a.Rmd" "pres.9.2a.Rmd~" "pres.9.2a.snm" [538] "pres.9.2a.tex" "pres.9.2a.toc" "pres.9.2a_notes.tex" [541] "pres.9.2a.vrb" "pres.9.3a.aux" "pres.9.3a.html" [544] "pres.9.3a.log" "pres.9.3a.nav" "pres.9.3a_notes.aux" [547] "pres.9.3a_notes.log" "pres.9.2a_notes.toc" [529] "pres.9.2a_notes.pdf" [511] "pres.9.1_notes.toc" "pres.9.1_tmp.md" "pres.9.1.out" "pres.9.1.pdf" [514] "pres.9.1.R" "pres.9.1.Rmd" "pres.9.1.Rmd~" [517] "pres.9.1.snm" "pres.9.1.tex" [520] "pres.9.1.toc" "pres.9.2a_notes.out" "pres.9.1.vrb" "pres.9.2a.aux" [523] "pres.9.2a.html" "pres.9.2a.log" "pres.9.2a.nav" [526] "pres.9.2a_notes.aux" "pres.9.2a_notes.log" [310] "pres.3.2.pdf" "pres.2.4.vrb" "pres.12.2.vrb" "pres.11.2a.md" "pres.10.6a_tmp.md" [100] "pres.10.6a.toc" "pres.10.6a.vrb" "pres.11.2a.aux" [103] "pres.11.2a.html" "pres.11.2a.log" [106] "pres.11.2a.nav" [97] "pres.10.6a.snm" "pres.11.2a_notes.aux" "pres.11.2a_notes.log" [109] "pres.11.2a_notes.out" "pres.11.2a_notes.pdf" "pres.11.2a_notes.tex" [112] "pres.11.2a_notes.toc" "pres.11.2a.out" "pres.11.2a.pdf" [115] "pres.11.2a.R" "pres.10.6a.tex" "pres.10.6a.Rmd~" "pres.11.2a.Rmd" [82] "pres.10.6a.html" [76] "pres.10.5a.snm" "pres.10.5a.tex" "pres.10.5a_tmp.md" [79] "pres.10.5a.toc" "pres.10.5a.vrb" "pres.10.6a.aux" "pres.10.6a.log" "pres.10.6a.Rmd" "pres.10.6a.nav" [85] "pres.10.6a_notes.aux" "pres.10.6a_notes.log" "pres.10.6a_notes.out" [88] "pres.10.6a_notes.pdf" "pres.10.6a_notes.tex" "pres.10.6a_notes.toc" [91] "pres.10.6a.out" "pres.10.6a.pdf" "pres.10.6a.R" [94] "#pres.10.6a.Rmd#" "#pres.11.2a.Rmd#" [118] "pres.11.2a.Rmd~" "pres.10.5a.Rmd" [145] "pres.13.aux" [139] "pres.12.2.Rmd" "pres.12.2.snm" "pres.12.2.tex" [142] "pres.12.2_tmp.md" "pres.12.2.toc" [307] "pres.2.4.toc" "pres.13.html" "pres.12.2.R" "pres.13.log" [148] "pres.13.nav" "pres.13_notes.aux" "pres.13_notes.log" [151] "pres.13_notes.out" "pres.13_notes.pdf" "#pres.12.2.Rmd#" [136] "pres.12.2.pdf" "pres.11.2a.snm" "pres.12.2.log" "pres.11.2a.tex" [121] "pres.11.2a_tmp.md" "pres.11.2a.toc" "pres.11.2a.vrb" [124] "pres.12.2.aux" "pres.12.2.html" [127] "pres.12.2.md" "pres.12.2.out" "pres.12.2.nav" "pres.12.2_notes.aux" [130] "pres.12.2_notes.log" "pres.12.2_notes.out" "pres.12.2_notes.pdf" [133] "pres.12.2_notes.tex" "pres.12.2_notes.toc" "pres.10.5a.Rmd~" [73] "#pres.10.5a.Rmd#" [154] "pres.13_notes.toc" [22] "courses_tmp.md" [16] "courses.pdf" "courses.R" "courses.Rmd" [19] "courses.Rmd~" "courses.snm" "courses.tex" "courses.toc" "courses_notes.toc" "end1.matter" [25] "figure" "getInits.R" "ggmapTemp.png" [28] "head.template" "images" "courses.out" [13] "courses_notes.tex" [31] "macnally.jpg" "beamerHeader.sty" [1] "/home/murray/tmp" #change the current working directory path > setwd ("/home/murray/tmp/") > list.files (path= getwd ()) #list all files (and directories) in the current working directory [1] "A_knit_.rmd~" "A.Rmd~" [4] "cache" "courses_notes.pdf" "courses.aux" "courses.html" [7] "courses.log" "courses.nav" "courses_notes.aux" [10] "courses_notes.log" "courses_notes.out" "junk2.dzslides" "Makefile" "pres.10.5a.R" "pres.10.4.vrb" "pres.10.4.Rmd" "pres.10.4.Rmd~" [55] "pres.10.4.snm" "pres.10.4.tex" "pres.10.4_tmp.md" [58] "pres.10.4.toc" "pres.10.5a.aux" "pres.10.4.pdf" [61] "pres.10.5a.html" "pres.10.5a.log" "pres.10.5a.nav" [64] "pres.10.5a_notes.aux" "pres.10.5a_notes.log" "pres.10.5a_notes.out" [67] "pres.10.5a_notes.pdf" "pres.10.5a_notes.tex" "pres.10.5a_notes.toc" [70] "pres.10.5a.out" "pres.10.5a.pdf" [52] "pres.10.4.R" "pres.10.4.out" "Makefile~" [40] "pres.10.4.aux" [34] "MeansParam.aux" "MeansParam.log" "MeansParam.pdf" [37] "#model.txt#" "model.txt" "NotebookPaper.jpg" "pres.10.4.html" [49] "pres.10.4_notes.toc" "pres.10.4.log" [43] "pres.10.4.nav" "pres.10.4_notes.aux" "pres.10.4_notes.log" [46] "pres.10.4_notes.out" "pres.10.4_notes.pdf" "pres.10.4_notes.tex" "pres.13_notes.tex" "pres.13.out" > getwd () "pres.2.1_notes.tex" "pres.2.1.log" "pres.2.1.nav" [253] "pres.2.1_notes.aux" "pres.2.1_notes.log" "pres.2.1_notes.out" [256] "pres.2.1_notes.pdf" "pres.2.1_notes.toc" "pres.2.1.aux" [259] "pres.2.1.out" "pres.2.1.pdf" "pres.2.1.R" [262] "pres.2.1.Rmd" "pres.2.1.Rmd~" "pres.2.1.snm" [250] "pres.2.1.html" "pres.1.vrb" "pres.2.1_tmp.md" [238] "pres.1old.vrb" [232] "pres.1old.Rmd" "pres.1old.Rmd~" "pres.1old.snm" [235] "pres.1old.tex" "pres.1old_tmp.md" "pres.1old.toc" "pres.1.out" [247] "pres.1.toc" "pres.1.pdf" [241] "pres.1.R" "pres.1.Rmd" "pres.1.Rmd~" [244] "pres.1.snm" "pres.1.tex" "pres.1_tmp.md" [265] "pres.2.1.tex" "pres.2.1.toc" "pres.1old.pdf" [298] "pres.2.4old.vrb" [292] "pres.2.4old.R" "pres.2.4old.Rmd" "pres.2.4old.snm" [295] "pres.2.4old.tex" "pres.2.4old_tmp.md" "pres.2.4old.toc" "pres.2.4.out" [289] "pres.2.4old_notes.toc" "pres.2.4old.out" "pres.2.4.pdf" [301] "pres.2.4.R" "pres.2.4.Rmd" "pres.2.4.Rmd~" [304] "pres.2.4.snm" "pres.2.4.tex" "pres.2.4_tmp.md" "pres.2.4old.pdf" [286] "pres.2.4old_notes.out" "pres.2.4old_notes.pdf" "pres.2.4old_notes.tex" [268] "pres.2.1.vrb" "pres.2.4_notes.log" "pres.2.4.aux" "pres.2.4.html" [271] "pres.2.4.log" "pres.2.4.md" "pres.2.4.nav" [274] "pres.2.4_notes.aux" "pres.2.4_notes.out" "pres.2.4old_notes.aux" "pres.2.4old_notes.log" [277] "pres.2.4_notes.pdf" "pres.2.4_notes.tex" "pres.2.4_notes.toc" [280] "pres.2.4old.aux" "pres.2.4old.html" "pres.2.4old.log" [283] "pres.2.4old.nav" "pres.1old.R" [229] "pres.1old.out" "pres.13.pdf" "pres.14.tex" "pres.14.out" "pres.14.pdf" [178] "pres.14.R" "pres.14.Rmd" "pres.14.Rmd~" [181] "pres.14.snm" "pres.14_tmp.md" "pres.14_notes.tex" [184] "pres.14.toc" "pres.14.vrb" "pres.15.aux" [187] "pres.15.html" "pres.15.log" "pres.15.nav" [175] "pres.14_notes.toc" "pres.14_notes.pdf" "pres.15_notes.log" [163] "pres.13.toc" [157] "pres.13.R" "pres.13.Rmd" "pres.13.Rmd~" [160] "pres.13.snm" "pres.13.tex" "pres.13_tmp.md" "pres.13.vrb" [172] "pres.14_notes.out" "pres.14.aux" [166] "pres.14.html" "pres.14.log" "pres.14.md" [169] "pres.14.nav" "pres.14_notes.aux" "pres.14_notes.log" [190] "pres.15_notes.aux" "pres.15_notes.out" "pres.1old_notes.toc" "pres.1old.html" "pres.1_notes.log" [214] "pres.1_notes.out" "pres.1_notes.pdf" "pres.1_notes.tex" [217] "pres.1_notes.toc" "pres.1old.aux" [220] "pres.1old_knit_.Rmd" [211] "pres.1.nav" "pres.1old.log" "pres.1old.nav" [223] "pres.1old_notes.aux" "pres.1old_notes.log" "pres.1old_notes.out" [226] "pres.1old_notes.pdf" "pres.1old_notes.tex" "pres.1_notes.aux" "pres.1.md" [193] "pres.15_notes.pdf" "pres.15.Rmd" "pres.15_notes.tex" "pres.15_notes.toc" [196] "pres.15.out" "pres.15.pdf" "pres.15.R" [199] "#pres.15.Rmd#" "pres.15.Rmd~" "pres.1.log" [202] "pres.15.snm" "pres.15.tex" "pres.15_tmp.md" [205] "pres.15.toc" "pres.15.vrb" "pres.1.aux" [208] "pres.1.html" #review the current working directory Basic usage y t o r r e c d i i n g o r k w T h e Note that R uses Unix directory slashes ( / ) NOT Windows directory slashes ( \ )

  2. > q () Basic usage n g t t i Q u i NOTE - do not put such a command in a script!

  3. Basic syntax o n s c t i F u n • collections of commands that expand the syntax of R • perform a single action • parameters that alter behaviour

  4. 9 10 > seq (from=2,to=10) 8 7 6 5 4 3 2 [1] Basic syntax s i o n n c t F u seq() - a function that generates sequences

  5. [1] 2 ...) function (from = 1, to = 1, by = ((to - from)/(length.out - 1)), length.out = NULL, along.with = NULL, 9 10 8 7 6 5 4 3 > seq (from=2,to=10) Basic syntax s i o n n c t F u seq() - a function that generates sequences • type seq( and then hit the TAB key

  6. 8 10 3 6 4 2 [1] > seq (from=2,to=10,by=2) 9 10 8 7 6 5 4 2 [1] > seq (from=2,to=10) Basic syntax s i o n n c t F u seq()

  7. 6 10 5 2 [1] > seq (from=2,to=10,length.out=3) 8 10 6 4 2 [1] > seq (from=2,to=10,by=2) 9 10 8 7 6 4 3 2 [1] > seq (from=2,to=10) Basic syntax s i o n n c t F u seq()

  8. 6 10 8 2 [1] > seq (from=2,to=10,l=3) 6 10 2 [1] > seq (from=2,to=10,length.out=3) 8 10 6 4 2 [1] > seq (from=2,to=10,by=2) 9 10 7 6 > seq (from=2,to=10) 5 4 3 2 [1] Basic syntax s i o n n c t F u seq() Unique parameter identifiers

  9. Your turn • generate a sequence of 10 numbers that increments by 2 and starting at 8

  10. > seq (from=8, len=10, by=2) [1] 8 10 12 14 16 18 20 22 24 26 Your turn • generate a sequence of 10 numbers that increments by 2 and starting at 8

  11. function (from, to, by, length.out = NULL, along.with = NULL, ...) > apropos ("^seq\\.") > str (seq.Date) ...) function (from = 1, to = 1, by = ((to - from)/(length.out - 1)), length.out = NULL, along.with = NULL, > str (seq.default) "seq.POSIXt" "seq.default" "seq.int" [1] "seq.Date" Basic syntax s i o n n c t F u d e d l o a v e r O seq()

  12. [9] "2000-04-25" > # create a sequence of dates spaced 7 days apart between 29th Feb 2000 and 30th Apr 2000 [7] "2000-04-11" "2000-04-18" [5] "2000-03-28" "2000-04-04" [3] "2000-03-14" "2000-03-21" [1] "2000-02-29" "2000-03-07" > sampleDates to= as.Date ("2000-04-30"), by="7 days") + > sampleDates <- seq (from= as.Date ("2000-02-29"), Basic syntax n s t i o u n c d F a d e r l o O v e seq() > # print (view) these dates

  13. [1] "2000-03-28" > mean ( c (1,2,3,4)) > mean (sampleDates) [1] 2.5 Basic syntax n s t i o u n c d F d e l o a v e r O

  14. > help (mean) Getting help

  15. > ?mean > help (mean) Getting help

  16. > ?mean > help (mean) > example (mean) Getting help

  17. #list all demos available on your system > demo (graphics) > demo () #run the graphics demo Getting help o n s a t i s t r m o n D e

  18. > apropos ('mea') Getting help n g c h i e a r S

  19. > apropos ('mea') > help.search ('mean') #search the local R manuals > help.start () #search the local HTML R manuals Getting help i n g r c h S e a

  20. > apropos ('mea') > help.search ('mean') #search the local R manuals > help.start () #search the local HTML R manuals Getting help i n g r c h S e a

  21. > args (mean) function (x, trim = 0, na.rm = FALSE, ...) #the arguments that apply to the mean function function (x, ...) NULL > args (mean.default) NULL > args (list.files) #the arguments that apply to the list.files function function (path = ".", pattern = NULL, all.files = FALSE, full.names = FALSE, recursive = FALSE, ignore.case = FALSE, include.dirs = FALSE, no.. = FALSE) NULL Getting help s e n t g u m a r i o n n c t F u and don฀t forget the auto-complete

  22. [10] "package:base" [1] ".GlobalEnv" "Autoloads" "package:methods" [7] "package:datasets" "package:grDevices" "package:utils" [4] "package:graphics" "package:stats" "package:knitr" > search () Package management g e s c k a P a Loaded packages

  23. "package:base" [1] ".GlobalEnv" [10] "Autoloads" "package:methods" "package:datasets" [7] "package:utils" "package:grDevices" "package:graphics" [4] "package:stats" "package:knitr" > library (MASS) "package:MASS" Package management g e s c k a P a Loading packages > search ()

  24. [10] "package:base" > detach ("package:MASS") "Autoloads" "package:methods" [7] "package:datasets" "package:grDevices" "package:utils" [4] "package:graphics" "package:stats" "package:knitr" [1] ".GlobalEnv" > library (MASS) Package management g e s c k a P a Loading packages > search () Unloading packages > search ()

  25. c ("Package", "LibPath", "Version", "Depends", "Built")) > installed.packages () + > installed.packages (fields = Package management s a g e a c k d p l l e s t a i n i n g s t L i

  26. > install.packages () > available.packages () Package management g e s c k a p a l e l a b v a i A g e s c k a p a i n g a l l n s t I

  27. > install.packages ('car') > available.packages () > install.packages () Package management e s k a g p a c l e l a b v a i A g e s c k a p a i n g a l l n s t I - For example, to install the car package

  28. Package management (Rstudio) s a g e a c k e p a b l a i l A v

  29. Package management (Rstudio) s a g e a c k g p l i n t a l I n s

  30. Data types (atomic modes) e s t y p t a d a r y i m a P r Type Description integer whole numbers numeric numbers (floating points) character words logical TRUE/FALSE (0/1)

  31. Data types e s t y p t a d a e d r i v D e Type Description factor categorical variable date (Date) number of days since 1970-01-01 POSIX (Time and Date) number of seconds since 1900-01-01 00:00:00

  32. > Date <- c ('2000-02-29','2002-08-20','2004-02-21') > (Dates<- as.Date (Date)) [1] "2000-02-29" "2002-08-20" "2004-02-21" > mean (Dates) [1] "2002-04-24" Dates • But what about other formats!

  33. > as.Date ('29/02/2000', format='%d/%m/%Y') [1] "2000-02-29" > as.Date ('29th Feb 2000', format='%dth %b %Y') [1] "2000-02-29" Dates check out the help for strptime

  34. [1] "2000-02-29 07:22:30 AEST" > as.POSIXct ('29/02/2000 07:22:30', format='%d/%m/%Y %H:%M:%S') Dates and times

  35. > Dates [1] "2000-02-29" "2002-08-20" "2004-02-21" > format (Dates,format='%d/%m/%Y') [1] "29/02/2000" "20/08/2002" "21/02/2004" > format (Dates,format='%b %Y') [1] "Feb 2000" "Aug 2002" "Feb 2004" > format (Dates,format='%Y') [1] "2000" "2002" "2004" Formatting dates and times

  36. > Date [1] "2000-02-29" "2002-08-20" "2004-02-21" > library (lubridate) [1] "2000-02-29 UTC" "2002-08-20 UTC" "2004-02-21 UTC" > dmy ('29/02/2000') [1] "2000-02-29 UTC" Dates and times packages: lubridate > ymd (Date)

  37. > Dates [1] "2000-02-29" "2002-08-20" "2004-02-21" > library (lubridate) [1] 2000 2002 2004 > class (Date) [1] "character" > decimal_date (Dates) [1] 2000.161 2002.633 2004.139 > week (Dates) [1] 9 34 8 Dates and times packages: lubridate > year (Dates)

  38. > ddays (100) [1] "8640000s (~100 days)" > Dates+100 [1] "2000-06-08" "2002-11-28" "2004-05-31" > Dates+ ddays (100) [1] "2000-06-08" "2002-11-28" "2004-05-31" > Dates+ dweeks (5) [1] "2000-04-04" "2002-09-24" "2004-03-27" Dates and times packages: lubridate

  39. > dmy_hms ('29/02/2000 07:22:30') [1] "2000-02-29 07:22:30 UTC" > dmy_hms ('29/02/2000 07:22:30')+100 [1] "2000-02-29 07:24:10 UTC" > dmy_hms ('29/02/2000 07:22:30')+ ddays (100) [1] "2000-06-08 07:22:30 UTC" > quarter (Dates) [1] 1 3 1 Dates and times packages: lubridate

  40. Storage types Type Description vector 1-d array (same type) matrix 2-d array (same type and length) list collection of vectors data frame 2-d array (any type, same length)

  41. > # a numeric vector > TEMPERATURE <- c (36.1, 30.6, 31, 36.3, 39.9, 6.5, 11.2, 12.8, 9.7, 15.9) > TEMPERATURE [1] 36.1 30.6 31.0 36.3 39.9 6.5 11.2 12.8 9.7 [10] 15.9 Data types o r s e c t V

  42. [1] "Fish" "Rock" "Tree" "Git" > TEMPERATURE > WORDS > WORDS<- c ('Fish', 'Rock', 'Tree', "Git") > # a character vector [10] 15.9 9.7 6.5 11.2 12.8 [1] 36.1 30.6 31.0 36.3 39.9 > TEMPERATURE <- c (36.1, 30.6, 31, 36.3, 39.9, 6.5, 11.2, 12.8, 9.7, 15.9) > # a numeric vector Data types r s t o V e c

  43. TRUE [10] 15.9 TRUE FALSE TRUE [1] > BOOL > BOOL<- c (TRUE, TRUE, FALSE, TRUE) > # a boolean vector [1] "Fish" "Rock" "Tree" "Git" > WORDS > WORDS<- c ('Fish', 'Rock', 'Tree', "Git") > # a character vector 9.7 6.5 11.2 12.8 [1] 36.1 30.6 31.0 36.3 39.9 > TEMPERATURE > TEMPERATURE <- c (36.1, 30.6, 31, 36.3, 39.9, 6.5, 11.2, 12.8, 9.7, 15.9) > # a numeric vector Data types s t o r V e c

  44. [1] 4 4 4 4 4 > rep (4,5) > rep ('Fish',5) [1] "Fish" "Fish" "Fish" "Fish" "Fish" Data types o r s e c t V Regular sequences

  45. "Q10" [1] "Q1" "Q9" [8] "Q8" "Q7" "Q6" "Q5" "Q4" "Q3" "Q2" > QUADRATS > QUADRATS <- c ("Q1","Q2","Q3","Q4","Q5","Q6","Q7","Q8","Q9","Q10") Data types s t o r V e c paste()

  46. "Q10" "Q7" "Q9" [8] "Q8" "Q7" "Q6" "Q5" "Q4" "Q3" "Q2" [1] "Q1" > QUADRATS > QUADRATS <- paste ("Q",1:10,sep="") "Q10" "Q9" [8] "Q8" "Q6" "Q5" "Q4" "Q3" "Q2" [1] "Q1" > QUADRATS > QUADRATS <- c ("Q1","Q2","Q3","Q4","Q5","Q6","Q7","Q8","Q9","Q10") Data types r s c t o V e paste()

  47. [1] 36.1 30.6 31.0 36.3 39.9 > TEMPERATURE 6.5 11.2 12.8 9.7 [10] 15.9 > names (TEMPERATURE) NULL Data types r s t o V e c properties

  48. 9.7 15.9 > names (TEMPERATURE) <- QUADRATS 6.5 11.2 12.8 36.1 30.6 31.0 36.3 39.9 Q10 Q9 Q8 Q7 Q6 Q5 Q4 Q3 Q2 Q1 > TEMPERATURE NULL > names (TEMPERATURE) [10] 15.9 9.7 6.5 11.2 12.8 [1] 36.1 30.6 31.0 36.3 39.9 > TEMPERATURE Data types s t o r V e c properties

  49. > SITE <- paste ( rep (LETTERS[1:5], each = 2), 1:2, sep = "") > SITE [1] "A1" "A2" "B1" "B2" "C1" "C2" "D1" "D2" "E1" [10] "E2" Data types r s c t o V e

  50. Your turn Calculate the area of circles with the following radii: • circle A=10cm • circle B=12.3cm • circle C=25.6cm and store the results in a vector with item names reflecting the circle names

  51. > RADIUS <- c (10,12.3,25.6) > names (RADIUS) <- paste ('circle',LETTERS[1:3]) > names (RADIUS) <- paste ('circle',LETTERS[1: length (RADIUS)]) > AREA <- pi*RADIUS^2 > AREA circle A circle B circle C 314.1593 475.2916 2058.8742 Your turn > names (RADIUS) <- c ('circle A', 'circle B', 'circle C') > #OR > #OR even better

  52. [7] "full" "full" "full" "full" > SHADE "full" "no" "no" "no" "no" [1] "no" > SHADE <- rep ( c ("no","full"),each=5) Data types o r s c t V e Factors

  53. Levels: full no "no" [10] full full full full full no no no no [1] no > SHADE > SHADE <- factor (SHADE) [7] "full" "full" "full" "full" "full" "no" "no" "no" [1] "no" > SHADE > SHADE <- rep ( c ("no","full"),each=5) Data types o r s e c t V Factors

  54. Levels: no full > SHADE <- factor (SHADE, levels= c ("no","full")) no no no full full full full [10] full Levels: full no > SHADE [1] no [1] no no no no no full full full full [10] full no > SHADE > SHADE > SHADE <- rep ( c ("no","full"),each=5) [7] "full" "full" "full" "full" "full" "no" "no" "no" "no" [1] "no" > SHADE <- factor (SHADE) Data types s t o r V e c Factors

  55. Levels: no full > SHADE [10] full full full full full no no no no [1] no > SHADE <- gl (2,5,10, c ("no","full")) Data types s t o r V e c Factors

  56. Your turn Create a categorical vector with: - three levels (A, B and C) - four replicates of each level

  57. > gl (3,4,12,lab=LETTERS[1:3]) [1] A A A A B B B B C C C C Levels: A B C Your turn But what if you needed to arrange such that there was only two replicates in a row??

  58. [1] A A B B C C A A B B C C > gl (3,2,12,lab=LETTERS[1:3]) Levels: A B C Your turn

  59. [,1] [,2] > matrix (SHADE,nrow=5) [1,] "no" "full" [2,] "no" "full" [3,] "no" "full" [4,] "no" "full" [5,] "no" "full" Data types s c e t r i M a

  60. [5,] 11.77 13.12 > XY [4,] 15.49 12.20 7.61 16.65 [3,] [2,] 24.03 12.93 8.37 [1,] 16.92 Y > X <- c (16.92,24.03,7.61,15.49,11.77) X Data types s i c e a t r M > Y<- c (8.37,12.93,16.65,12.2,13.12) > XY <- cbind (X,Y)

  61. 8.37 12.93 16.65 12.20 13.12 [2,] 24.03 12.93 Y 7.61 15.49 11.77 X 16.92 24.03 [,5] [,4] [,3] [,2] [,1] > rbind (X,Y) [5,] 11.77 13.12 [4,] 15.49 12.20 7.61 16.65 [3,] 8.37 [1,] 16.92 Y X > XY > X <- c (16.92,24.03,7.61,15.49,11.77) Data types s i c e a t r M > Y<- c (8.37,12.93,16.65,12.2,13.12) > XY <- cbind (X,Y)

  62. "E2" Q7 "B2" Q5 "39.9" "C1" Q6 "6.5" "C2" "11.2" Q4 "D1" Q8 "12.8" "D2" Q9 "9.7" "E1" Q10 "15.9" "36.3" "B1" > cbind (TEMPERATURE,SITE) "31" Q3 "A2" "30.6" Q2 "A1" "36.1" Q1 TEMPERATURE SITE Data types s i c e a t r M Beware of mixing data types

  63. 2 11.2 Q5 39.9 1 Q6 6.5 2 Q7 2 36.3 Q8 12.8 2 Q9 9.7 2 Q10 15.9 1 Q4 TEMPERATURE SHADE > cbind (TEMPERATURE,SHADE) 31.0 Q3 1 30.6 Q2 1 36.1 Q1 1 Data types s i c e a t r M Beware of mixing data types

  64. 9.7 15.9 Levels: no full [3,] 7.61 16.65 [4,] 15.49 12.20 [5,] 11.77 13.12 $SHADE [1] no no no no no full full full full [10] full $TEMPERATURE 8.37 Q1 Q2 Q3 Q4 Q5 Q6 Q7 Q8 Q9 Q10 36.1 30.6 31.0 36.3 39.9 6.5 11.2 12.8 [2,] 24.03 12.93 [1,] 16.92 $QUADRATS "Q2" COORDINATES = XY, SHADE = SHADE, + TEMPERATURE = TEMPERATURE) > EXPERIMENT $SITE [1] "A1" "A2" "B1" "B2" "C1" "C2" "D1" "D2" "E1" [10] "E2" Y [1] "Q1" "Q3" > EXPERIMENT <- list (SITE=SITE,QUADRATS = QUADRATS, "Q4" "Q5" "Q6" "Q7" [8] "Q8" "Q9" "Q10" $COORDINATES X + Data types s i s t L

Download Presentation
Download Policy: The content available on the website is offered to you 'AS IS' for your personal information and use only. It cannot be commercialized, licensed, or distributed on other websites without prior consent from the author. To download a presentation, simply click this link. If you encounter any difficulties during the download process, it's possible that the publisher has removed the file from their server.

Recommend


More recommend