Getting Started with R

Some Commands

Misc

data     <- read.table(filename,header=TRUE)           # Read tab or space delimited file with a header
csv_data <- read.table(filename,header=TRUE,sep=',')   # Read csv files with a header

v1       <- c(1,2,3,4,5,6,7)                           # Create a vector from these values
v2       <- c(1:10)                                    # Create a vector from 1 to 10
v3       <- c(v1,v2)                                   # Concatenate the vectors

matrix1  <- cbind(v1,v2)                               # Column bind into a 2 x n matrix

ls                                                     # List objects

Mean, SD, etc.

> sapply(h, length)
WList_Len 
     4599 
> sapply(h, mean)
WList_Len 
  15.6845 
> sapply(h, sd)
WList_Len 
 15.74873 
> 

Installing Modules

> install.packages('vars',,'http://cran.ms.unimelb.edu.au')
also installing the dependencies ‘zoo’, ‘strucchange’, ‘urca’, ‘lmtest’, ‘sandwich’

trying URL 'http://cran.ms.unimelb.edu.au/bin/windows/contrib/3.3/zoo_1.7-13.zip'
Content type 'application/zip' length 900080 bytes (878 KB)
downloaded 878 KB

trying URL 'http://cran.ms.unimelb.edu.au/bin/windows/contrib/3.3/strucchange_1.5-1.zip'
Content type 'application/zip' length 789069 bytes (770 KB)
downloaded 770 KB

trying URL 'http://cran.ms.unimelb.edu.au/bin/windows/contrib/3.3/urca_1.2-9.zip'
Content type 'application/zip' length 935642 bytes (913 KB)
downloaded 913 KB

trying URL 'http://cran.ms.unimelb.edu.au/bin/windows/contrib/3.3/lmtest_0.9-34.zip'
Content type 'application/zip' length 284171 bytes (277 KB)
downloaded 277 KB

trying URL 'http://cran.ms.unimelb.edu.au/bin/windows/contrib/3.3/sandwich_2.3-4.zip'
Content type 'application/zip' length 504789 bytes (492 KB)
downloaded 492 KB

trying URL 'http://cran.ms.unimelb.edu.au/bin/windows/contrib/3.3/vars_1.5-2.zip'
Content type 'application/zip' length 559244 bytes (546 KB)
downloaded 546 KB

package ‘zoo’ successfully unpacked and MD5 sums checked
package ‘strucchange’ successfully unpacked and MD5 sums checked
package ‘urca’ successfully unpacked and MD5 sums checked
package ‘lmtest’ successfully unpacked and MD5 sums checked
package ‘sandwich’ successfully unpacked and MD5 sums checked
package ‘vars’ successfully unpacked and MD5 sums checked

Some Links

• http://www.statmethods.net/interface/workspace.html
• http://www.rstudio.com/ide/docs/using/history

• https://www.zoology.ubc.ca/~schluter/R/calculate/
• http://www.r-tutor.com/elementary-statistics/numerical-measures/mean
• http://stat.ethz.ch/R-manual/R-patched/library/stats/html/median.html
• http://www.cyclismo.org/tutorial/R/index.html
• http://answers.oreilly.com/topic/2597-how-to-compute-basic-statistics-with-r/
• http://grokbase.com/g/r/r-help

Personality Project

• http://www.personality-project.org/r/r.commands.html
• http://www.personality-project.org/r/r.guide.html
• http://personality-project.org/r/r.short.html
• http://www.personality-project.org/r/psych/HowTo/

Neil Gunther's 911 Help

See File:911r.pdf for the PDF document.

# Emergency help for R newbies
# Created by NJG on Fri Aug 15, 2011
# Updated by NJG on Mon Jan  9  2012

# 1. Set YOUR path to this 911.r script in the following variable:

# LUNUX:   Command needs to gnome-open
# WINDOWS: C:\\You path\\ needs "\\" to esc '\' and '9' in 991.r
# MACOSX:  Put in explicit full path name
path911 <- "/Users/njg/Desktop//911/"

# 2. Now run this script in the R Console.
# It is assumed that Acrobat Reader is installed.

help.911 <- function() {
	ost <- "unknown"
	switch(Sys.info()[['sysname']],
		Darwin  = {
			ost <- "Mac OS X"
			cmd <- "open"
		},
		Linux = {
			ost <- "Linux"
			cmd <- "gnome-open"
		},
		Windows = {
			ost <- "Windows"
			cmd <- "open"
		}
	)
	if(ost == "unknown") { stop("Unknown operating system") }
	setwd(path911)
	t <- try(system(paste(cmd, "911r.pdf")))
	if(t != 0) { stop(sprintf("system() call problem running on %s system", ost)) }
}

# 3. Usage
cat("Type 'help.911()' at anytime at the R Console prompt\n")
cat("to open the 911r.pdf file in Acrobat Reader.\n")
cat("Like this:\n")
cat("> help.911()\n")

help.911()

Command History

files <- list.files(pattern="*.csv")
files

data <- lapply(files, read.table, sep=",", header=T)
data

stuff <- as.data.frame(do.call(rbind, data))
stuff

dim(stuff)


print(dframe)
print(stuff)

mean(stuff)
sd(stuff)

sapply(stuff, length)
sapply(stuff, mean)
sapply(stuff, sd)

getwd()
ls()
ls
options()
history()
history(max.show=Inf)