MAT2001-STATISTICS FOR ENGINEERS-

EMBEDDED LAB
R-BASICS-LAB-I & LAB-2

“There are three kinds of lies: lies, damned lies, and

statistics.” (B.Disraeli)
What is R?
 Introduction to R
 R is a language and environment for statistical computing and
graphics. It is a GNU project which is similar to the S
language and environment which was developed at Bell
Laboratories (formerly AT&T, now Lucent Technologies) by
John Chambers and colleagues. R can be considered as a
different implementation of S. There are some important
differences, but much code written for S runs unaltered
under R.
 R provides a wide variety of statistical (linear and nonlinear
modelling, classical statistical tests, time-series analysis,
classification, clustering, …) and graphical techniques, and is
highly extensible. The S language is often the vehicle of
choice for research in statistical methodology, and R provides
an Open Source route to participation in that activity
 One of R’s strengths is the ease with which well-designed
publication-quality plots can be produced, including
mathematical symbols and formulae where needed. Great
care has been taken over the defaults for the minor design
choices in graphics, but the user retains full control.
 R is available as Free Software under the terms of the Free
Software Foundation’s GNU General Public License in
source code form
The R environment
 R is an integrated suite of software facilities for data manipulation,
calculation and graphical display. It includes
 an effective data handling and storage facility,
 a suite of operators for calculations on arrays, in particular
matrices,
 a large, coherent, integrated collection of intermediate tools for
data analysis,
 graphical facilities for data analysis and display either on-screen or
on hardcopy, and
 a well-developed, simple and effective programming language
which includes conditionals, loops, user-defined recursive
functions and input and output facilities.
Why should you adopt R
Programming Language?
 R, SAS, and SPSS are three statistical languages. Of these
three statistical languages, R is the only an open source. SAS
is the most important private software business in the world.
SPSS is now overseen by IBM. R Programming is extensible
and hence, R groups are noted for its energetic
contributions. Lots of Rs typical features can be written in R
itself and hence, R has gotten faster over time and serves as a
glue language.
 Installation
 R can be downloaded from one of the mirror sites
in http://cran.r-project.org/mirrors.html. You should pick
your nearest location.
 Using External Data
 R offers plenty of options for loading external data, including
Excel, Minitab and SPSS files. We have included a tutorial
titled Data Import on the subject for the purpose.
Command Line vs Scripts
What is command line?
 Execution of commands in R is not menu driven. (Not like
Clicking over buttons to get outcome)

 We need to type the commands.

 Single line and multi line commands are possible to write.

 When writing multi-line programs, it is useful to use a text

editor rather than execute everything directly at the
command line.
BASICS
 > is the prompt sign in R. you can enter numbers
and perform calculations.
 The assignment operators are the left arrow with
dash <- ,->and equal sign =.
 > x = 30 #assigns the value 30 to x.
 >y=x*2 #assigns the value 2*x to y.
 >z=x+y #assigns the value x + y to z.

Output:-
 # : The character # marks the beginning of a
comment.All characters until the end of the line are
ignored.
> # mu is the mean
> # x <- 20 is treated as comment only
 Capital and small letters are different.
 The command c(1,2,3,4,5) combines the numbers 1,2,3,4
and 5 to a vector.
R as a calculator
 Integer Division: Division in which the fractional part
(remainder) is discarded

 x mod y : modulo operation finds the remainder

after division of one number by another
Maximum: max & Minimum: min
Overview Over Further Functions
 log() function computes natural logarithms (Ln) for
a number or vector. log10 computes common
logarithms (Lg).log2 computes binary logarithms
(Log2). log(x,b) computes logarithms with base b.
 Complex functions
Matrix
 Matrices are important objects in any calculation. A
matrix is a rectangular array with p rows and n
columns.
 Note:-
 The parameter nrow defines the row number of a matrix.
 The parameter ncol defines the column number of a
matrix.
 The parameter data assigns specified values to the matrix
elements.
 In R, a 4 × 2-matrix X can be created with a following
command:

 One can access a single element of a matrix with x[i,j]:

 In case, the data has to be entered row wise, then a 4 × 2-
matrix X can be created with
 We can get specific properties of a matrix:
 Assigning a specified number to all matrix elements:
 Construction of a diagonal matrix, here the identity
matrix of a dimension 2:
 Transpose of a matrix X: X’
 Multiplication of a matrix with a constant

 Matrix multiplication: operator %*%

 Cross product of a matrix X, X’X, with a function
crossprod

 Note: Command crossprod() executes the multiplication

faster than the conventional method with t(x)%*%x
 Addition and subtraction of matrices (of same
dimensions!) can be executed with the usual operators +
and –
 Access to rows, columns or submatrices:
 Inverse of a matrix:
solve( ) finds the inverse of a positive definite matrix
 5x = 10, what's x?

 Let's see two variables examples:

 3x + 2y = 8 , x + y =2; What's x and y?
 Eigen Values and Eigen Vectors:
eigen() finds the eigen values and eigen vectors of a
positive definite matrix
 Changing the matrix's elements
 Sequences
 A sequence is a set of related numbers, events, movements, or
items that follow each other in a particular order. The regular
sequences can be generated in R.
 Sequence with constant increment and decrement:
 Generating sequences of dates
 Generating sequences of letters-lower case alphabets
 sequence of uppercase alphabets
 Repeats-the command rep
 Repetition of characters
 Sorting
sort function sorts the values of a vector in ascending
order (by default) or descending order.
 Length
length() function gets or sets the length of a vector (list) or other
objects
 Data Structures in R
 Extracting vector elements by the element index(es)
 Extracting vector elements by a logical expression
 set operations (union, intersection, asymmetric difference,
equality and membership) on two vectors.
 Union() is not the same as concatenation c() because c() will
duplicate values that are common to both vectors
 Missing data
R represents missing observations through the data value NA
We can detect missing values using is.na

 Now try a vector to know if any value is missing?

 Logical Operators and Comparisons

 TRUE and FALSE are reserved words denoting logical

constants.
 Conditional execution
ifelse(test, yes, no)

 Interpretation
If x < 6 (TRUE), then x = x2 (YES) .
If x ≥ 6 (FALSE), then x = x + 1 (NO).
 The for loop
 If the number of repetitions is known in advance a for()
loop can be used.
 The while() loop
If the number of loops is not known in before, e.g. when
an iterative algorithm to maximize a likelihood function
is used, one can use a while() loop.
 The repeat loop
The repeat loop doesn’t test any condition — in contrast
to the while() loop — before entering the loop and also
not during the execution of the loop.
 Data Frames
 In a data frame, we can combine variables of equal
length,with each row in the data frame containing
observations on the same unit.
 Data frames contain complete data sets that are mostly
created with other programs (spreadsheet-files, software
SPSS-files, Excel-files etc.).
 Variables in a data frame may be numeric (numbers) or
categorical (characters or factors).
 Package ―MASS‖ describes functions and datasets to support
Venables and Ripley, ``Modern Applied Statistics with S'' (4th
edition 2002)
 Structure of the data:
The result of str gives the dimension as well as the name
and type of each variable.
 Extract a variable from data frame using $
Suppose we want to extract information on variable School from the
data set painters.
 Subsets of a data frame can be obtained with subset() or
with the second equivalent command:
 The command split partitions the data set by values of a
specific variable. This should preferably be a factor
variable.
 Practice
1. x=23, y=x^2, z=y^3+x^2

2. 2+3-4*6/3-2^3/4+7-6*2**3/2+2-4+5 ?

3. -478/365-12^3/47-6**2**3-27^-(1/3)-7^6*3

4. c(1,2,3,4)*c(1,2,2,4)^c(2,1,3,2)**c(1,2,3,3)-
c(2,3,402653185,262145) ?

5. c(2,3,4,5)^c(2,3)+c(12,23,14,25)^c(3,2)-c(5,6,7,8)*c(2,3)?
 max(c(62,83,44,75)^-c(9,-3)) / min(c(52,62,71,85) ^
c(2,3) ) -prod(c(1,2,1,2) ^c(1,2)) +max(c(12,13,14,15) ^
c(2,3)) ?

 X1 <- c(123,258,318,624), X2 <- sqrt(X1^3)+X1/X1^2-

X1**(1/2)

 X<-matrix(nrow=3,ncol=3,data=
c(10,20,30,40,50,60,70,80,90),byrow=F) then X[ ,2]?

 X<-matrix(nrow=3, ncol=3, data =

c(10,20,30,40,50,60,70,80,90) , byrow=F) then
X[2:3,2:3] ?
 x = 33:53 then specifies the outcomes of the following
statements:
x[(x > 28) & (x < 59)] and x[(x > 28) || (x < 59)]?

 sqrt(abs(seq(-6,6, by = 3)))?

 letter[1:3]?

 LETTERS[18:12]

 rep(2:5,4)
 rep(70:65, times=3)?

 rep(20:25, each=3)?

 x <- matrix(nrow=2, ncol=2, data=1:4, byrow=T)

rep(x, each=3)?

 sort(c(20,50, 10, 30, 90,70, 80), decreasing = FALSE)?

 sort (c(20,50, 10, 30, 90,70, 80), increasing = TRUE)?

 order(c(20,50, 10, 30, 90,70, 80), decreasing = FALSE)?

 mode(c(1, 2, "3", 8+9, "7+9", 6.7, 110*45))?

 x <- c(10, 75, 20, 35, 30, 40, 180, 50, 60, 27, 70, 67, 80, 50,
39, 120)
>x[(x>50)]?
>x[(x - 20 > 40)]
>x[(x^2 + 10 > 50)]

• factor( c(rep("male",2), rep("female", 3)))?