473,320 Members | 1,766 Online
Bytes | Software Development & Data Engineering Community
Post Job

Home Posts Topics Members FAQ

Join Bytes and contribute your articles to a community of 473,320 developers and data experts.

Lesson 5 (part2) – r programming structures

nbiswas
149 100+
We have been using many built-in R functions throughout the lessons.
After you start using R more, you will reach the point when you will want to create reusable sections of statements.
Functions are simply a named group of R expressions that are considered an R object of type function().
Input values or function arguments are passed to functions by value. Each argument is either matched with the function definition by position or by name.
We have used built-in functions that are also considered generic, such as the plot() function. Generic functions are designed to accept a variety of different arguments. R functions can accept a variable length of input arguments and they can even accept references to other functions.
When functions are called the actual arguments are assigned to local variables within the function. The statements in the function body are evaluated based on the input data. Control will be returned to the invoking function when the return() statement is reached or when the final expression in the function body is reached.
If there is no explicit data returned to the calling function, the output generated by the last expression will be returned.
Let's take a look at an example of a simple R function.
We want find that we are often writing code to determine the number of students who passed a course so we decide to define a function to simply our code.
The function is defined using the function() expression. The name of this function has been specified as numPassed and the function takes a single argument.
The comments on line 3 through 6 help explain the usage of the function includes the input and output datatypes expected.
The body or content of the function is contained within the brace brackets ({}). We decided to use the vector enabled sum() function with a conditional expression to calculate the number of students with marks of 50 or higher. The pass object on line 8 is actually a vector of boolean values based on the criteria defined in line 7.
If you want to examine the source code of a function you can type the name of the function, without specifying any arguments, in the R Console.
This function will accepts a single dataframe or a list and it will perform a conditional expression tests for each element of the data structure. If the student has obtained a mark of 50 or more a logical value of TRUE is generated. These logical values are collected in a local vector called pass.
We then use the sum() function to count the number of TRUE values and return the data to the calling function as an integer.
Finally, the function is being called within the concatenate or cat() function.
One of the most used features of R is the *apply() family of functions.
Here we are using the lapply() or List Apply function. The function that we are invoking is the mean or average function within the base package.
The mean function will be executed across each column of the data.frame called m.
The lapply() function always returns a list as its ouput. In our scenario, the list with the mean or average by subject is provided.
In the box of the right side we see how the sapply() or simplified apply function can be used. sapply() basically performs the same task as lapply(), but it coerces the output into a vector or a matrix. Notice how sapply generated a vector of double values where each element has a name associated with it. The final expression demonstrates how to retrieve the element named "Math".
Here we are using the same approach of using lapply() to invoke a function across a dataset.
In the previous example we used the built-in function mean and in this example we are using the function numPassedCourse, that we created, to find the number of students who passed the course.
These examples of using the apply() function are quite simple, but you can probably see how a single line of R code can actually be performing a complex operation including iteration.
Now that we have learned about flow of control and functions, let's examine data relationships across multiple variables.
In our case study we want to determine if there is a relationship between the height of a person and their shoe size.
We will use R to simulate data sets for our analysis and we will use quantitative measures and visualizations to perform our data analysis.
There are a few different quantitative measurements that can be used to measure the degree of varying between multiple variables.
A positive covariance would indicate that there is a positive linear relationship between the two variables, and a negative covariance would indicate the opposite.
Correlation is a similar statistical measurement of the degree of linear relationship between variables, but it has the added benefit that it has a well defined range of values. The coefficient range is from -1 to 1.
For example, a correlation coefficient of 0.8 would indicate that there is a strong linear relationship between the two variables. Therefore, if variable 1 increases there is a similar observed increase in variable 2. If the correlation coefficient is close to zero (0) there is minimal relationship between the changes in values of the variables. Negative coefficients would indicate that the direction of change in value of the 2 variables are opposite one another.
Here we are creating 2 uniform distributions of heights and shoe size for our analysis.
We use the covariance or cov() function to compute the covariance measurement and we use the correlation or cor() function to compute the correlation.
These values are stored in variables as they will be used in our final visualization.
We are using the base graphics package plot() function to create a scatterplot of the generated data. Visually it seems like there is no relationship between these two measurements. The values of correlation and covariance support our visual understanding as their values are close to (0) zero.
On line 22 a linear regression model is created and stored in a variable called lm1.
The linear model function or lm() can be used to perform an analysis of multiple data sets and determine an linear approximation of the relationship between the variables. Notice that the lm() function is provided a formula. In R a formula consists of a tilde (~) character where in this case the variable on the left is considered the response and the variable on the right of the tilde is the term under consideration.
The object lm1 contains quite a large amount of detailed information. In this scenario we use the object to add a line to our graph to visualize the relationship.
We use a math expression on line 26 to adjust the shoe size for our next analysis.
As an aside, the expression uses integer division to introduce a relationship between the variables.
Again we determine the covariance and correlation coefficients of the revised variables.
In our initial simulation of randomly generated uniform shoe sizes and heights there was very little correlation.
Here we see visually that there is a strong linear correlation and the value of the correlation coefficient agrees with our visual analysis.
We have also performed a revised linear regression model and plotted the result on the graph.
Here we see a summary of the before and after of our measurements for heights and shoe sizes.
The strong correlation is quite evident in the revised data for shoe sizes shown on the right.
There are many lab exercises that you can work with following this lesson including: dice simulations, correlations, and recoding of datasets using the apply family of functions in R.
Sep 4 '14 #1
2 5373
zmbd
5,501 Expert Mod 4TB
nbiswas:
Please provide proper citations for these articles.
-z
Sep 29 '14 #2
Thanks for the information
May 12 '15 #3

Sign in to post your reply or Sign up for a free account.

Similar topics

15
by: christopher diggins | last post by:
I have written an article on how to do Aspect Oriented Programming in vanilla C++ (i.e. without language extensions or other tools such as AspectC++). The article is available at...
65
by: Roger Smythe | last post by:
A means for the progressive decomposition a problem space into increasingly simpler component parts such that these component parts represent higher levels of conceptual abstraction, and are...
25
by: Justin Robbs | last post by:
I am trying to write the communcations part of a Point of Sale program for the Convenience Store industry. The setup in each store will have varying numbers of registers. There could be as few...
6
by: Ken Allen | last post by:
OK, I admit that I have been programming since before C++ was invented, and I have developed more than my share of assembly language systems, and even contributed to operating system and compiler...
15
by: UrsusMaximus | last post by:
While preparing a Python411 podcast about classes and OOP, my mind wondered far afield. I found myself constructing an extended metaphor or analogy between the way programs are organized and...
3
by: Juan R. | last post by:
In http://canonicalscience.blogspot.com/2006/04/scientific-language-canonml-is.html] I presented some generic requirements for a markup language for science and mathematics. Basic features of...
151
by: istillshine | last post by:
There are many languages around: C++, JAVA, PASCAL, and so on. I tried to learn C++ and JAVA, but ended up criticizing them. Is it because C was my first programming language? I like C...
1
nbiswas
by: nbiswas | last post by:
In this lesson we will initially learn about the features and uses of R. R is a software environment that is excellent for data analysis and graphics. It was initially created in 1993 by Ross...
1
nbiswas
by: nbiswas | last post by:
Welcome to the lesson on R data structures. To perform any meaningful data analysis we need to collect our data into R data structures. In this lesson we will explore the most frequently used...
1
nbiswas
by: nbiswas | last post by:
The programming structures that we will examine include: control flow statements and user defined functions. In previous lessons we were either performing summary or descriptive statistics across...
0
isladogs
by: isladogs | last post by:
The next Access Europe meeting will be on Wednesday 6 Mar 2024 starting at 18:00 UK time (6PM UTC) and finishing at about 19:15 (7.15PM). In this month's session, we are pleased to welcome back...
1
isladogs
by: isladogs | last post by:
The next Access Europe meeting will be on Wednesday 6 Mar 2024 starting at 18:00 UK time (6PM UTC) and finishing at about 19:15 (7.15PM). In this month's session, we are pleased to welcome back...
0
by: Vimpel783 | last post by:
Hello! Guys, I found this code on the Internet, but I need to modify it a little. It works well, the problem is this: Data is sent from only one cell, in this case B5, but it is necessary that data...
0
by: jfyes | last post by:
As a hardware engineer, after seeing that CEIWEI recently released a new tool for Modbus RTU Over TCP/UDP filtering and monitoring, I actively went to its official website to take a look. It turned...
0
by: ArrayDB | last post by:
The error message I've encountered is; ERROR:root:Error generating model response: exception: access violation writing 0x0000000000005140, which seems to be indicative of an access violation...
1
by: PapaRatzi | last post by:
Hello, I am teaching myself MS Access forms design and Visual Basic. I've created a table to capture a list of Top 30 singles and forms to capture new entries. The final step is a form (unbound)...
1
by: CloudSolutions | last post by:
Introduction: For many beginners and individual users, requiring a credit card and email registration may pose a barrier when starting to use cloud servers. However, some cloud server providers now...
0
by: Faith0G | last post by:
I am starting a new it consulting business and it's been a while since I setup a new website. Is wordpress still the best web based software for hosting a 5 page website? The webpages will be...
0
isladogs
by: isladogs | last post by:
The next Access Europe User Group meeting will be on Wednesday 3 Apr 2024 starting at 18:00 UK time (6PM UTC+1) and finishing by 19:30 (7.30PM). In this session, we are pleased to welcome former...

By using Bytes.com and it's services, you agree to our Privacy Policy and Terms of Use.

To disable or enable advertisements and analytics tracking please visit the manage ads & tracking page.