讲解:Econ 325、Limit Theorem、R、RJava|Prolog

Econ 325 (004)Winter Session, Term 1, 2019M. VaneyLab 2 - Demonstration of the Central Limit TheoremDue: Monday November 25. Submit your work online.PurposeIn this lab R is used to demonstrate the Central Limit Theorem, a theorem that provides atheoretical basis for estimation and inference even for underlying populations that are notnormally distributed. The lab reinforces the use of .do Öles as an e¢ cient way to execute aseries of commands and the use of loops to automate repetitive tasks. The lab also introducesa few additional R commands.Central Limit TheoremGiven a random sample of size n from underlying distribution f(x) with 1(Önite mean) and 0 approximately normal with X. This can also be expressed as limn!1 X.One implication of this for estimation is that even if the underlying distribution is not normallydistributed, by appealing to the Central Limit Theorem we may treat the sample mean,X�n; as an approximately normally distributed random variable. The following Ögure showsthe underlying distribution of a random variable X as a solid line. Clearly X is not normallydistributed. The random variable X has realizations only over the interval [0; 3] rather than(1;1);X is not symmetric, X is not uni-modal. However, taking random samples ofsize n and computing the sample mean for each di§erent random sample we see that thedistribution of the sample mean (red dashed line) has many of the features characteristic ofa normally distributed random variable (uni-modal, symmetric, bell-shaped).How closely the sample mean conforms to a normal distribution will depend on features ofthe underlying distribution and the sample size. The larger the sample size the more closelythe distribution of the sample mean will resemble a normally distributed random variable.Data and MethodologyA number of ëpopulationsíare provided. In order to demonstrate the CLT it will be necessaryto describe the distribution of the sample mean for each of the populations.DataThe Öle lab2-variables.csv contains N = 700 observations for each of 5 random variables(called x1; : : : ; x5). Each of these can be thought of as a di§erent Population with a givenunderlying distribution f(x1); g(x2); : : : ; k(x5).MethodsUse R to carry out the following tasks:1. (a) Generate summary statistics and create histograms for each of the 5 variables.(b) Draw 1000 random samples of size n = 4; 25 and 144 for each of the randomvariables (without replacement). Compute the sample mean for each randomsample and construct a histogram of the sample means..R commandsThis lab will make use of some commands that are found two additional packages availablein R: dplyr and ggplot2. Both of these packages must be loaded in R. You can check tosee which packages are loaded by selecting the packages tab in the lower right corner of thescreen. If a package has not been installed in the console the following command can beentered:install.packages(ggplot2)the ggplot2 package will be installed (it may take a minute or two)In order to make use of the additional commands available in a package your script Ölemust refer to the packages through a library commnad. It is best to start the script withspeciÖcation of the required packages:library(ggplot2)library(dplyr)The dplyr package has a number of commands that are useful for re-organizing data. Thecommand that we will use in this lab is sample_n(data, sample size)The ggplot2 package is used for making various graphs and Ögures. A very useful resourcefor creating histograms in ggplot2 can be found at the link provided in the Lab folder onCanvas.The sample_n() command will draw a single random sample (of rows of a dataset) of aspeciÖc size, n. To generate 1000 random samples, sample_n() command along with a2command to take the mean can be embedded in the command replicate() which will repeatthese commands a speciÖed number of times.Results and DiscussionPresent and provide some discussion of the following:Submit your .do Öle for this lab. Do not submit raw data.1. (a) Consider the summary statistics and graphics for the underlying populations. Dothe underlying distributions appear to be Normally distributed? Comment onthe apparent distributions of each of the variables (symmetric, skewed, numberof modes,di§erence between mean and median, etc.).(b) Discuss how changing the size of the sample alters the distribution of the samplemean for each of the di§erent variables. Do the results conform with the predictionof the Central Limit Theorem?3转自:http://www.daixie0.com/contents/18/4349.html

©著作权归作者所有,转载或内容合作请联系作者
  • 序言:七十年代末,一起剥皮案震惊了整个滨河市,随后出现的几起案子,更是在滨河造成了极大的恐慌,老刑警刘岩,带你破解...
    沈念sama阅读 216,651评论 6 501
  • 序言:滨河连续发生了三起死亡事件,死亡现场离奇诡异,居然都是意外死亡,警方通过查阅死者的电脑和手机,发现死者居然都...
    沈念sama阅读 92,468评论 3 392
  • 文/潘晓璐 我一进店门,熙熙楼的掌柜王于贵愁眉苦脸地迎上来,“玉大人,你说我怎么就摊上这事。” “怎么了?”我有些...
    开封第一讲书人阅读 162,931评论 0 353
  • 文/不坏的土叔 我叫张陵,是天一观的道长。 经常有香客问我,道长,这世上最难降的妖魔是什么? 我笑而不...
    开封第一讲书人阅读 58,218评论 1 292
  • 正文 为了忘掉前任,我火速办了婚礼,结果婚礼上,老公的妹妹穿的比我还像新娘。我一直安慰自己,他们只是感情好,可当我...
    茶点故事阅读 67,234评论 6 388
  • 文/花漫 我一把揭开白布。 她就那样静静地躺着,像睡着了一般。 火红的嫁衣衬着肌肤如雪。 梳的纹丝不乱的头发上,一...
    开封第一讲书人阅读 51,198评论 1 299
  • 那天,我揣着相机与录音,去河边找鬼。 笑死,一个胖子当着我的面吹牛,可吹牛的内容都是我干的。 我是一名探鬼主播,决...
    沈念sama阅读 40,084评论 3 418
  • 文/苍兰香墨 我猛地睁开眼,长吁一口气:“原来是场噩梦啊……” “哼!你这毒妇竟也来了?” 一声冷哼从身侧响起,我...
    开封第一讲书人阅读 38,926评论 0 274
  • 序言:老挝万荣一对情侣失踪,失踪者是张志新(化名)和其女友刘颖,没想到半个月后,有当地人在树林里发现了一具尸体,经...
    沈念sama阅读 45,341评论 1 311
  • 正文 独居荒郊野岭守林人离奇死亡,尸身上长有42处带血的脓包…… 初始之章·张勋 以下内容为张勋视角 年9月15日...
    茶点故事阅读 37,563评论 2 333
  • 正文 我和宋清朗相恋三年,在试婚纱的时候发现自己被绿了。 大学时的朋友给我发了我未婚夫和他白月光在一起吃饭的照片。...
    茶点故事阅读 39,731评论 1 348
  • 序言:一个原本活蹦乱跳的男人离奇死亡,死状恐怖,灵堂内的尸体忽然破棺而出,到底是诈尸还是另有隐情,我是刑警宁泽,带...
    沈念sama阅读 35,430评论 5 343
  • 正文 年R本政府宣布,位于F岛的核电站,受9级特大地震影响,放射性物质发生泄漏。R本人自食恶果不足惜,却给世界环境...
    茶点故事阅读 41,036评论 3 326
  • 文/蒙蒙 一、第九天 我趴在偏房一处隐蔽的房顶上张望。 院中可真热闹,春花似锦、人声如沸。这庄子的主人今日做“春日...
    开封第一讲书人阅读 31,676评论 0 22
  • 文/苍兰香墨 我抬头看了看天上的太阳。三九已至,却和暖如春,着一层夹袄步出监牢的瞬间,已是汗流浃背。 一阵脚步声响...
    开封第一讲书人阅读 32,829评论 1 269
  • 我被黑心中介骗来泰国打工, 没想到刚下飞机就差点儿被人妖公主榨干…… 1. 我叫王不留,地道东北人。 一个月前我还...
    沈念sama阅读 47,743评论 2 368
  • 正文 我出身青楼,却偏偏与公主长得像,于是被迫代替她去往敌国和亲。 传闻我的和亲对象是个残疾皇子,可洞房花烛夜当晚...
    茶点故事阅读 44,629评论 2 354

推荐阅读更多精彩内容