#### 联系方式

• QQ：99515681
• 邮箱：99515681@qq.com
• 工作时间：8:00-23:00
• 微信：codinghelp

#### 您当前位置：首页 >> Matlab编程Matlab编程

###### 日期：2021-03-12 10:16

Statistical Methods for Data Science

DATA7202

Semester 1, 2021

Assignment 1 (Weight: 25%)

Assignment 1 is due on 22 Mar 2021 17:00).

Please answer the questions below. For theoretical questions, you should present rigorous proofs

and appropriate explanations. Your report should be visually appealing and all questions should

be answered in the order of their appearance. For programming questions, you should present your

analysis of data using Python, Matlab, or R, as a short report, clearly answering the objectives

and justifying the modeling (and hence statistical analysis) choices you make, as well as discussing

your conclusions. Do not include excessive amounts of output in your reports. All the code should

be copied into the appendix and the sources should be packaged separately and submitted on the

blackboard in a zipped folder with the name:

"student_last_name.student_first_name.student_id.zip".

For example, suppose that the student name is John Smith and the student ID is 123456789.

Then, the zipped file name will be John.Smith.123456789.zip.

(a) The data is generated via the following data generation mechanism. Xi ～ U(0, 1) for

i ∈ {1, 2, 3}; here U(0, 1) stands for the continuous uniform distribution over the [0, 1] set.

However, we require that X1 + X2 + X3 = 1, that is, the explanatory variables stand for

a percentage of the budget.

(b) In addition, the model for y is as follow:

Y = 0.5X1 + 3X2 + 5X3 + 5X2X3 + 2X1X2X3 + W, (1)

where W ～ U(0, 1).

Similar to the original example, generate train and test sets of size N = 1000. Fit the linear regression

and the random forest models to the data. For the linear regression, make an inference

types to sales. Use the linear model and the RF (with 500 trees), to make a prediction (using

the test set), and report the corresponding mean squared errors.

When constructing datasets, please use “1” and “2” seeds for the train and the test sets,

respectively.

2. [10 Marks] Consider the following variant of the cross-validation procedure.

(i) Using the available data, find a subset of “good” predictors that show correlation with

the response variable.

(ii) Using these predictors, construct a model (for regression or classification).

(iii) Use cross-validation to estimate the model prediction error.

1

Is this a good method? Do you expect to obtain the true prediction error? Explain your

3. [5 Marks] Suppose that we observe X1, . . . , Xn ～ F. We model F as a normal distribution

with mean μ and standard deviation of σ. For this problem, determine the hypothesis class

H = {f(x, θ); θ ∈ Θ}.

and state explicitly what is θ and Θ.

4. [15 Marks] Let H be a class of binary classifiers over a set Z. Let D be an unknown distribution

over X , and let g be a target hypothesis in H. F Show that the expected value of LossT (g)

over the choice of T equals LossD(g), namely,

ET LossD(g) = LossD(g).

5. [15 Marks (see details below)] Consider the following dataset.

Now, suppose that we would like to consider two models.

Model1 : y = β1x1 + ε,

and

Model2 : y = β0 + β1x1 + ε,

where ε ～ N(0, 1). That is, we consider two linear models with and without the intercept.

(a) [5 Marks)] Fit these models tot the data and write the corresponding coefficients. Namely,

fill the following table:

Model β0 β1

Model1 0

Model2

(b) [5 Marks)] Consider the squared error loss, the absolute error loss, and the L1.5 loss. Find

the average loss for each model. Namely, fill the following table:

Model squared error loss absolute error loss L1.5 loss

Model1

Model2

(c) [5 Marks)] Draw a conclusion from the obtained results.

6. [30 Marks (see details below)] Consider the Hitters data-set (given in Hitters.csv). Our

objective is to predict a hitter’s salary via linear models.

(a) [5 Marks)] Load the data-set and replace all categorical values with numbers. (You can

use the LabelEncoder object in Python).

2

(b) [5 Marks)] Generally, it is better to use OneHotEncoder when dealing with categorical

variables. Justify the usage of LabelEncoder in (a).

(c) [20 Marks)] Fit linear regression and report 10-Fold Cross-Validation mean squared error.

7. [10 Marks)] Consider a function(2)

Suppose that a = 1, b = 2, and c = 3, and write a Crude Monte Carlo algorithm for the

estimation of ` using N = 10000 sample size. Deliver the 95% confidence interval. Compare

the obtained estimation with the true value ` as given in (2).

3