R Exercise Practical 2. Follow the instructions described in this sheet. Type or paste any answers you are asked for into a Microsoft Word document.

Size: px

Start display at page:

Download "R Exercise Practical 2. Follow the instructions described in this sheet. Type or paste any answers you are asked for into a Microsoft Word document."

Earl Stewart
5 years ago
Views:

1 R Exercise Practical 2 Follow the instructions described in this sheet. Type or paste any answers you are asked for into a Microsoft Word document. Note: You should have this sheet on your computer for reference and to copy and paste some text. 1. Obtain and attach the data set cats by typing the following: >library(mass) >data(cats) >attach(cats) Have a look at the data by typing cats and also use the?cats command. Plot Hwt (Heart Weight) against Bwt (Body Weight). Choose Bwt to be on the horizontal axis. >plot(hwt~bwt) Obtain the best straight line through the points using >abline(lm(hwt~bwt)) So far, the analysis has not looked at all at the gender of the cat. Type in the R code >boxplot(hwt~sex, xlab='gender', ylab='heart weight') which clearly shows that heart weight depends on gender. (use the subcommand horizontal =T if you want a different orientation). Similarly, the R code > plot(hwt~bwt, xlab='body weight (kg)', ylab='heart weight (g)', type='n') > points(hwt[sex=='m']~bwt[sex=='m']) > points(hwt[sex=='f']~bwt[sex=='f'], pch=3) produces a plot of heart weight against body weight for male cats (denoted by o) and female cats (denoted by +). 2. You can use R to plot some well known mathematical functions. Try the following: (a) > curve(sin(x),from =0, to =2*pi) (b) > curve(2*x^3+x^2-2*x+10,from =-2, to =2) (c) > curve(log(x),from =0, to =10,col="magenta") (d) >curve(5-x,from =0, to =10,col="red", add=t)

2 3. Download a new library to your R package using the command: >install.packages( UsingR ) Choose any mirror site you like (St. Andrews worked for me). Obtain and attach the data set cfb by typing the following: >library(usingr) >data(cfb) >attach(cfb) One of the vectors in this data frame is called INCOME. It represents the yearly income figures by households in a survey of consumer finances is the US in (a) Obtain a histogram of the INCOME data. What does it tell you? (b) Calculate the mean income using >mean(income) How good is this statistic in describing the data? (c) A subcommand of the mean( ) command is trim. Use?mean to investigate what this subcommand does, and experiment with various values of trim. (d) Calculate the median of the INCOME data using >median(income) Does this give a better feel for the data set than the mean? (e) Calculate the coefficients of symmetry and kurtosis of the INCOME data by defining the functions discussed in the lectures (See lecture 3 on Vision). 4. Consider the following three random samples from a population. Enter each dataset into R using the scan command). (i) The data below are the number of hours between charging on a particular type of mobile phone (ii) The data below are the number of hours spent watching television per week for a sample of 34 households

3 (iii) The data below are the number of miles travelled to and from work each day by a sample of 12 company employees (a) For each data set find the mean, median, standard deviation, quartiles and interquartile range. Use the commands >summary(hours) >sd(hours) (b) (c) Obtain a boxplot for each dataset. Obtain a histogram for each dataset. (d) Obtain a smooth curve through the histogram using the commands; > hist(hours,freq=f) > lines(density(hours)) (e) Recall from your statistics lectures that when samples are taken the statistic measured varies from sample to sample. Use the command >t.test(hours) to find a 95% confidence interval for the population mean for each sample. 5. The built in data set islands contains the size of the world s land masses that exceed square miles. Make a stem-and-leaf plot, then compare the mean, median and 25% trimmed mean. Are they similar? Hints: >data(islands) >mean(islands, 0.25) 6. The median absolute deviation is defined as mad(x)=1.4826median( x i median(x) ) This is a resistant measure of spread and is implemented using the mad( ) function. Explain in words what it measures. Compare the values of the sample standard deviation, interquartile range and median absolute deviation for the exec.pay data set in the Using R library. Hints: >library(usingr) >data(exec.pay)

4 7. Consider the results below that give data on highway traffic flow. It shows the relationship between traffic flow, x (in thousands of vehicles per day) and lead content, y (micrograms per gram) of the bark of nearby trees. (a) Enter these data into two columns in R, using for example, x=c(8.3,8.3,12.1, ) (b) Plot the data using plot(y~x), labelling the diagram appropriately. (c) You now wish to find the regression line of y on x. This has the equation y= a + bx The values a and b can be obtained in R using the lm command >lm(y~x) This gives the equation for y in terms of x (or at least the gradient b and point that it cuts the y axis, a). Assign this to a vector called regress. >regress=lm(y~x) Write down the equation of the best straight line and copy and paste the output from the >summary(regress) (e) The calculated straight-line equation can be drawn on the graph using >plot(y~x) >abline(lm(y~x)) Copy and paste this graph into your Word document (d) Find the residuals using the command >residuals(regress) It is often useful to look at a plot of the residuals against x. Do this using >plot(residuals(regress)~x) Do they follow any pattern? (e) Obtain the Product-Moment Correlation Coefficient, r, and comment on what its value means. >cor(x,y)

5 8. A series of numbers x 1, x 2, x 3, can be created from constants a, b and c using the formula x n+1 = (ax n + c) mod b This is called the Linear Congruential Method of random number generation. Note that x mod y is the remainder obtained when x is divided by y. So 20mod3 =2 Use a for loop to generate the first 100 numbers of such a series for the values (a) x 1 =1, a=5, b=94 and c=0 (b) x 1 =2, a=4, b=99 and c=3 Hint x(mody) is obtained in R by x%%y Which is the more efficient at producing random numbers?

Further Maths Notes. Common Mistakes. Read the bold words in the exam! Always check data entry. Write equations in terms of variables

Further Maths Notes. Common Mistakes. Read the bold words in the exam! Always check data entry. Write equations in terms of variables Further Maths Notes Common Mistakes Read the bold words in the exam! Always check data entry Remember to interpret data with the multipliers specified (e.g. in thousands) Write equations in terms of variables