In this part of the project, we will be using R to run simulations of the...
Fantastic news! We've Found the answer you've been seeking!
Question:
Transcribed Image Text:
In this part of the project, we will be using R to run simulations of the model that is, the model proposed by J. Chang in his paper "Recent common ancestors of all present-day individuals". Before proceeding, make sure you fully understand that model. 1.1 Function to simulate the TMRCA Write an R. function that receives as input the size n of the population, runs a simulation of that model having that population size, and returns the number of generations it took to reach a MRCA. (Make sure you test this function is free of bugs by testing it for "small" values of n.) 1.2 Simulating the TMRCA Now we are going to "collect our data". That is, we will choose a "large" value of n and call the function defined above many times using that value of n. Make sure you collect the data (TMRCA) for each run of that function. Repeat this procedure for different values of n. We would like n to be very large, as close as possible to the size of the actual human population today. As you will see, if n is that large, the R code will run for too long. So you will have to find a suitable value of n that is not too small but at the same time allows you to actually run that function in a reasonable time. How many times should you call the function? Good question! If you read Chang's paper, Table 1 shows he called the function 25 times. You should also have at least 25 simulations (for each value of n). How small could n be? Another good question. Again, if you read Chang's paper, you will see in the same table he chose the largest n to be 4000. Your largest n should be at least 4000. 1.3 Presenting the results of the simulation Present the results of your simulations in a "nice" way. The bare minimum is a table like the one in Chang's paper. For full marks, though, your presentation of the simulations will have to be "nicer" both in terms of its form and of its content. This is the time to use your knowledge of descriptive statistics and data organization skills. Include also how much time it took to run the simulations. (In R. you can get the current time by calling Sys.time ().) 1.4 (Bonus) Further exploration If you want to play more with these simulations, this is your chance! Let your imagination ask interesting questions and then seek the answers! If you can come up with your own ideas, it is better but, if not, here are some ideas: You could find the time to "a generation in which each individual is either a CA of all present-day individuals or an ancestor of no present-day individual" (this is what Chang calls Un). • Can you make the model slightly more realist and then write an R code to run simulations of those modified versions of the model? • At the generation of the MRCA, are there usually more than one MRCA? If so, how many? You could do some meaningful statistics (confidence interval, hypothesis testing, ...) with the data collected. In this part of the project, we will be using R to run simulations of the model that is, the model proposed by J. Chang in his paper "Recent common ancestors of all present-day individuals". Before proceeding, make sure you fully understand that model. 1.1 Function to simulate the TMRCA Write an R. function that receives as input the size n of the population, runs a simulation of that model having that population size, and returns the number of generations it took to reach a MRCA. (Make sure you test this function is free of bugs by testing it for "small" values of n.) 1.2 Simulating the TMRCA Now we are going to "collect our data". That is, we will choose a "large" value of n and call the function defined above many times using that value of n. Make sure you collect the data (TMRCA) for each run of that function. Repeat this procedure for different values of n. We would like n to be very large, as close as possible to the size of the actual human population today. As you will see, if n is that large, the R code will run for too long. So you will have to find a suitable value of n that is not too small but at the same time allows you to actually run that function in a reasonable time. How many times should you call the function? Good question! If you read Chang's paper, Table 1 shows he called the function 25 times. You should also have at least 25 simulations (for each value of n). How small could n be? Another good question. Again, if you read Chang's paper, you will see in the same table he chose the largest n to be 4000. Your largest n should be at least 4000. 1.3 Presenting the results of the simulation Present the results of your simulations in a "nice" way. The bare minimum is a table like the one in Chang's paper. For full marks, though, your presentation of the simulations will have to be "nicer" both in terms of its form and of its content. This is the time to use your knowledge of descriptive statistics and data organization skills. Include also how much time it took to run the simulations. (In R. you can get the current time by calling Sys.time ().) 1.4 (Bonus) Further exploration If you want to play more with these simulations, this is your chance! Let your imagination ask interesting questions and then seek the answers! If you can come up with your own ideas, it is better but, if not, here are some ideas: You could find the time to "a generation in which each individual is either a CA of all present-day individuals or an ancestor of no present-day individual" (this is what Chang calls Un). • Can you make the model slightly more realist and then write an R code to run simulations of those modified versions of the model? • At the generation of the MRCA, are there usually more than one MRCA? If so, how many? You could do some meaningful statistics (confidence interval, hypothesis testing, ...) with the data collected.
Expert Answer:
Answer rating: 100% (QA)
To fulfill the requirements of the project and conduct the simulations as described it is advisable to have a programming environment with R installed ... View the full answer
Related Book For
Income Tax Fundamentals 2013
ISBN: 9781285586618
31st Edition
Authors: Gerald E. Whittenburg, Martha Altus Buller, Steven L Gill
Posted Date:
Students also viewed these accounting questions
-
Read the case study "Southwest Airlines," found in Part 2 of your textbook. Review the "Guide to Case Analysis" found on pp. CA1 - CA11 of your textbook. (This guide follows the last case in the...
-
The following additional information is available for the Dr. Ivan and Irene Incisor family from Chapters 1-5. Ivan's grandfather died and left a portfolio of municipal bonds. In 2012, they pay Ivan...
-
List three specific parts of the Case Guide, Objectives and Strategy Section (See below) that you had the most difficulty understanding. Describe your current understanding of these parts. Provide...
-
2. Evaluate the following definite integrals (a) (2x - 6x) dx /6 (b) 16 cos 3t + 2 sin 3t dt x+1 2 da
-
A vertical lock gate is 4 m wide and separates 20C water levels of 2 m and 3 m, respectively. Find the moment about the bottom required to keep the gate stationary.
-
At December 31, 2020, Higley Corporation has one temporary difference which will reverse and cause taxable amounts in 2021. In 2020, a new tax act set taxes equal to 35% for 2020, 30% for 2021, and...
-
Migami Computer Parts Inc. is in the process of setting a selling price on a new component it has just designed and developed. The following cost estimates for this new component have been provided...
-
Chloe and Emma start a new business, Cement Sidewalks and Accessories (CSA), during the current year. CSA is organized as a partnership. Chloe owns 40% of CSA; Emma owns the remaining 60%. Chloe and...
-
Question 1 Points 2 Find the surface area of a sphere with a radius of 22 cm. 6,084.57 cm 1,024.57 cm 2,585.5 cm O 245.67 cm
-
Problems 8 and 15 in Chapter 13 presented a scenario for Bangs Leisure Chairs. The Solver Sensitivity Report for the optimal solution (Problem 15a only) is shown in Figure 15.38. Clearly explain the...
-
Accounting Beachfound usen three processes to marwitacture lits for personal walercraft forming The cost transferred into Finished Goods inventary is the cost of the ifis. a Iftin parts from...
-
Nancy has an individual RESP for her son, Matthew which was opened 8 years ago. Matthew, age 20 passed away in a motorcycle accident this year. Current RESP breakdown: Contribution: $20,000 Canada...
-
Virat landed his first job at Telus upon graduating from Concordia in December 2023. After his 3 months of probation, Virat was offered to participate in Telus' Group Registered Retirement Savings...
-
ABC Company manufactures rubber shoes. It gets its raw materials from a supplier that grants it 30 days credit. ABC purchases P2.4 million in raw materials per year. Answer the questions that follow...
-
Malfoy Industries, who uses the FIFO method of process costing, manufactures products in three sequential processes: Molding, Plating, and Finishing. In the Plating process, direct materials are...
-
I am working on a project to use a renewable fuel for powering a village in india . I am using soybean oil for my diesel generators. Q . If plant A has a capital cost of $ 2 0 million for production...
-
Activity 2B Refferto Case Study A: Sustainability Policy fforQuestion 1 1. Identiffy three resources that may be required to implement the sustainability policy detailed in the case study. 2. Using...
-
Match each of the key terms with the definition that best fits it. _______________ A record of the sequence of data entries and the date of those entries. Here are the key terms from the chapter. The...
-
The following additional information is available for the Dr. Ivan and Irene Incisor family from Chapters 1 and 2. On September 1, Irene opened a retail store that specializes in sports car...
-
Tom has a successful business with $100,000 of income in 2012. He purchases one new asset in 2012, a new machine which is 7-year MACRS property and costs $25,000. If you are Tom's tax advisor, how...
-
Jan has two jobs during 2012. One employer withheld and paid FICA taxes on $66,600 of Jan's salary, and the other employer withheld and paid FICA taxes on $44,400 in salary paid to Jan. Calculate the...
-
The Mendoza Company is installing an absorption standard-costing system and a flexible-overhead budget. Standard costs have recently been developed for its only product and are as follows: Expected...
-
The expected activity of the paper-making plant of Leventhal Paper Company was 45,000 machine hours per month. Practical capacity was 60,000 machine hours per month. The standard machine hours...
-
Precision AutoParts Company has a highly automated manufacturing process for producing a variety of auto parts. Through the use of computer-aided manufacturing and robotics, the company has reduced...
Study smarter with the SolutionInn App