STAT UN2102 Homework
STAT UN2102 Homework 4 [100 pts]
Due 11:59pm Monday, May 6th on Canvas
Your homework should be submitted on Canvas as an R Markdown file. Please submit
the knitted .pdf or .html file along with the .Rmd file. We will not (and cannot) accept
any other formats. Please clearly label the questions in your responses and support your
answers by textual explanations and the code you use to produce the result. We may print
out your homeworks. Please do not waste paper by printing the dataset or any vector over,
say, length 20.
Goals: Simulating probability distributions using the accept-reject method, simulating a
sampling distribution related to the linear regression model.
1 Reject-Accept Method
Let random variable X denote the temperature at which a certain chemical reaction takes
place. Suppose that X has probability density function
Perform the following tasks:
1. Determine the maximum of f(x). Find an envelope function e(x) by using a uniform
distribution for g(x) and setting e(x) = maxx{f(x)}.
2. Using the Accept-Reject Algorithm, write a program that simulates 1000 draws
from the probability density function f(x) from Equation 1.
3. Plot a histogram of your simulated data with the density function f overlayed in the
graph. Label your plot appropriately.
2 Regression and Empirical Size
2.1 Regression
We work with the grocery retailer dataset from Canvas. The description follows:
1
A large national grocery retailer tracks productivity and costs of its facilities closely. Consider
a data set obtained from a single distribution center for a one-year period. Each data
point for each variable represents one week of activity. The variables included are number
of cases shipped in thousands (X1), the indirect costs of labor as a percentage of total
costs (X2), a qualitative predictor called holiday that is coded 1 if the week has a holiday
and 0 otherwise (X3), and total labor hours (Y ). Consider the multiple linear regression
model
(2) Yi = β0 + β1 Xi1 + β2 Xi2 + β3 Xi3 + i, i = 1, 2, . . . , 52,
and iid~ N(0, σ2).
Perform the following tasks:
4. Read in the grocery retailer dataset. Name the dataset grocery.
5. Use the least squares equation = (XTX)
1XTY to estimate regression model (2).
To estimate the model, use the linear model function in R, i.e., use lm().
6. Use R to estimate σ2, i.e., compute MSE =1
. To perform this task,
use the residuals function.
2.2 Test for Slope
STAT UN2102作业代做、代做R Markdown file作业、代写R课程作业
Now consider investigating if the number of cases shipped (X1) is statistically related to
total labor hours (Y ). To investigate the research question, we run a t-test on the coefficient
corresponding to X1, i.e., we test the null alternative pair
(3) H0 : β1 = 0 versus HA : β1 6= 0.
To run the hypothesis testing procedure, we use the t-statistic
1 is the second element of the least squares estimator β= (XTX)
1XTY and
SE(β1) is the standard error of β?
1. The least squares estimates, estimated standard errors,
t-statistics and p-values for all coefficients β0, β1, β2, β3 are nicely organized in the standard
linear regression output displayed in Table 1. To get this output in R, use the summary()
function on your model.
Test the manager’s claim in (3) using the R functions lm() and summary().
2
Table 1: Standard Multiple Linear Regression Output
Estimate Std. Error t value Pr(> |t|) or Sig
(Intercept) β
2.3 Sampling Distribution
Under model (2) and under the null hypothesis H0 : β1 = 0, the test statistic (4) has a
student’s t-distribution with n 4 degrees of freedom, i.e.,
The goal of this section is to simulate the sampling distribution of the t-statistic.
Perform the following tasks:
5. Write a loop that simulates the sampling distribution of the t-statistic under null
hypothesis (3) with the multiple linear regression model (2). To accomplish this task:
i. Assume the true model relating Y with X1, X2, X3 is
(5) Yi = 4200 + β1Xi1 ? 15X2 + 620X3 + i, i = 1, 2, . . . , 52,i
iid~ N(0, 20500).
ii. Assuming H0 : β1 = 0 is true, simulate 10,000 draws from model (5) using the
fixed covariates X2, X3.
iii. For each iteration of the loop, fit the full model
using the simulated Y and fixed covariates X1, X2, X3.
iv. For each iteration of the loop, also compute the t-statistic from equation (4).
Store these values in a vector t.stat. Hint: Use the summary function in R and
extract the actual summary table using the code summary(model)[[4]]. Then
extract the relevant t-statistic from the table.
v. Display the first six elements of your simulated t-values.
3
7. Plot a histogram of the simulated sampling distribution. Overlay the correct t-density
on this histogram, i.e., overlay the density t(df = 52 ? 4). Plot the density in green
and set breaks=40 in the histogram. Make sure to label the plot appropriately. You
can use base R or ggplot.
8. Recall that the significance level of a testing procedure is defined as
P(Type I error) = P(Rejecting H0 when H0 is true) = α.
The significance level is often called the size of the testing procedure. Based on
significance levels α = 0.10, 0.05, 0.01, compute the sample proportion of simulated
t-values that fell in the rejection region. The proportion of simulated rejected t-values
under the null is called the empirical size of a test. The three values should be close
to the actual α levels.
因为专业,所以值得信赖。如有需要,请加QQ:99515681 或邮箱:99515681@qq.com
微信:codinghelp
STAT UN2102 Homework的更多相关文章
- bzoj 4320: ShangHai2006 Homework
4320: ShangHai2006 Homework Time Limit: 10 Sec Memory Limit: 128 MB Description 1:在人物集合 S 中加入一个新的程序员 ...
- HDU 1789 Doing Homework again(贪心)
Doing Homework again 这只是一道简单的贪心,但想不到的话,真的好难,我就想不到,最后还是看的题解 [题目链接]Doing Homework again [题目类型]贪心 & ...
- hdu-1789-Doing Homework again
/* Doing Homework again Time Limit: 1000/1000 MS (Java/Others) Memory Limit: 32768/32768 K (Java/Oth ...
- HDU 1789 Doing Homework again (贪心)
Doing Homework again http://acm.hdu.edu.cn/showproblem.php?pid=1789 Problem Description Ignatius has ...
- Doing Homework 状态压缩DP
Doing Homework 题目抽象:给出n个task的name,deadline,need. 每个任务的罚时penalty=finish-deadline; task不可以同时做.问按怎样的 ...
- 机器学习 —— 概率图模型(Homework: Exact Inference)
在前三周的作业中,我构造了概率图模型并调用第三方的求解器对器进行了求解,最终获得了每个随机变量的分布(有向图),最大后验分布(双向图).本周作业的主要内容就是自行编写概率图模型的求解器.实际上,从根本 ...
- hdoj 1789 Doing Homework again
Doing Homework again Time Limit: 1000/1000 MS (Java/Others) Memory Limit: 32768/32768 K (Java/Oth ...
- homework做了些什么?
第一步:get_new_guid_uid_pairs_{$ymd} 参数是时间和100上的文件. 那么100上的文件是从哪里来的呢? 我们进入到100机器上,打开root权限下的cron,看到如下内容 ...
- HDU 1074 Doing Homework (dp+状态压缩)
题目链接:http://acm.hdu.edu.cn/showproblem.php?pid=1074 题目大意:学生要完成各科作业, 给出各科老师给出交作业的期限和学生完成该科所需时间, 如果逾期一 ...
随机推荐
- JavaSE 可变参数的方法重载
/** * 可变参数的方法重载 */ class A { public void test(int a, int b) { System.out.println(a+", "+b) ...
- 2017-2018 ACM-ICPC Latin American Regional Programming Contest GYM101889
挺有意思的一套题,题也没有啥毒瘤了,本来是队切的结果种种原因大家全挂机了. 只补了百人题,一共7个,其他的暂时先不补了,,也不会嘛qwq H:签到 #include <bits/stdc++.h ...
- mac homebrew
brew 全称Homebrew 是Mac OSX上的软件包管理工具 Homebrew 安装和卸载工具 只用一行命令就能完成 官方地址: http://brew.sh/index.html ...
- java个人博客源码
初入博客园,请各位多关照,来而不往非礼也. 如需要源码以及学习内容,qq:1397617269,我看到就回你们资源. 直接给链接: 链接:https://pan.baidu.com/s/1S_awtg ...
- SQL Server 2008作业失败:无法确定所有者是否有服务器访问权限
转自:http://www.mzwu.com/article.asp?id=3592 下午准备将服务器数据导一份到本地,连上服务器找数据库备份文件发现最后一次备份是7月1日,竟然十几天没生成新备份,查 ...
- Codeforces 1089K - King Kog's Reception - [线段树][2018-2019 ICPC, NEERC, Northern Eurasia Finals Problem K]
题目链接:https://codeforces.com/contest/1089/problem/K time limit per test: 2 seconds memory limit per t ...
- 程序员调 Bug 的样子,非常真实
程序员调 Bug 的样子,非常真实
- LeetCode 1013 Partition Array Into Three Parts With Equal Sum 解题报告
题目要求 Given an array A of integers, return true if and only if we can partition the array into three ...
- POJ3461 Oulipo 字符串
正解:kmp/哈希 解题报告: 传送门! 这题其实就kmp板子,,,用来复习下kmp的太久没打了QAQ 所以kmp做法就不港了放个代码就是了QAQ #include<algorithm> ...
- Jmeter学习之-从数据库取出数据并且校验数据是否准确
https://www.cnblogs.com/wuyonghuan/p/7479582.html 应用场景:调用某个接口像数据库中插入数据,需要在接口调用完成后查看数据更新或插入的数据是否正确的时候 ...