应用统计学第11周作业(多元线性回归、正则方程组)

1.

在这里插入图片描述
相关系数 = 0.95,R Square = 0.897,F检验值为46.7,p-value = 3.88e-08<<0.05,F检验通过,整体回归是显著的。
三个自变量房产估价的p-value = 0.0013,t检验通过,说明房产估价对销售价格的解释性很强。
自变量地产估价和使用面积的p-value>0.05,对销售价格的解释性不是很强。

2.

在这里插入图片描述
使用逐步回归的方法确定最后留下的变量,最后留下的是各项贷款余额和本年固定资产投资额。
在这里插入图片描述
复相关系数R=0.872,表明模型中的两个自变量与因变量线性相关密切;决定系数 R 2 = 0.761 R^2=0.761 R2=0.761,说明因变量固体垃圾排放量变异的76.1%可由预测变量各项贷款余额和本年固定资产投资额来解释,表明模型对数据的拟合程度较好;调整的决定系数 A d j R 2 = 0.739 AdjR^2=0.739 AdjR2=0.739,也表明模型的拟合效果较好.
在这里插入图片描述
显著性均为0,说明在SPSS默认的引入、剔除(0.05、0.1)自变量的统计学标准下,所拟合的模型具有统计学意义。
在这里插入图片描述
偏回归系数解释自变量对因变量的影响,而标准化回归系数可以比较各自变量对因变量影响的大小,标准化回归系数的绝对值越大,自变量对的影响越大,由此可知,x1的影响大于x2对 y ^ \hat{y} y^的影响。

在这里插入图片描述
逐步回归过程中每一步被剔除的变量情况及其检验结果。以本研究的最优回归模型2为例,本年累计应收贷款和贷款项目个数是未进入模型的两个自变量,看其显著性值可知,本年累计应收贷款为0.066(p>0.05),贷款项目个数为0.788(p>0.05),均不满足进入标准(p≦0.05),因此被排除。

3.

正 则 方 程 组 : ∂ ( ∑ e ^ i 2 ) ∂ b 0 = − 2 ∑ i = 1 n ( y i − b 0 − b 1 x i 1 − ⋯ − b k x i k ) = 0 ∂ ( ∑ e ^ i 2 ) ∂ b 1 = − 2 ∑ i = 1 n x i 1 ( y i − b 0 − b 1 x i 1 − ⋯ − b k x i k ) = 0 … … … . . ∂ ( ∑ e ^ i 2 ) ∂ b k = − 2 ∑ i = 1 n x i k ( y i − b 0 − b 1 x i 1 − ⋯ − b k x i k ) = 0  解:  b 0 , b 1 , ⋯   , b k 正则方程组:\\\begin{array}{l} \frac{\partial\left(\sum \hat{e}_{i}^{2}\right)}{\partial b_{0}}=-2 \sum_{i=1}^{n}\left(y_{i}-b_{0}-b_{1} x_{i 1}-\cdots-b_{k} x_{i k}\right)=0 \\ \frac{\partial\left(\sum \hat{e}_{i}^{2}\right)}{\partial b_{1}}=-2 \sum_{i=1}^{n} x_{i 1}\left(y_{i}-b_{0}-b_{1} x_{i 1}-\cdots-b_{k} x_{i k}\right)=0 \\ \ldots \ldots \ldots . . \\ \frac{\partial\left(\sum \hat{e}_{i}^{2}\right)}{\partial b_{k}}=-2 \sum_{i=1}^{n} x_{i k}\left(y_{i}-b_{0}-b_{1} x_{i 1}-\cdots-b_{k} x_{i k}\right)=0 \\ \text { 解: } \quad b_{0}, b_{1}, \cdots, b_{k} \end{array} b0(e^i2)=2i=1n(yib0b1xi1bkxik)=0b1(e^i2)=2i=1nxi1(yib0b1xi1bkxik)=0..bk(e^i2)=2i=1nxik(yib0b1xi1bkxik)=0 b0,b1,,bk
证 明 : 由 k 个 正 则 方 程 式 可 得 ∑ i = 1 n y i = b 0 + b 1 x i 1 + ⋯ + b k x i k ∑ i = 1 n x i 1 y i = b 0 ∑ i = 1 n 1 ⋅ x i 1 + b 1 ∑ i = 1 n x i 1 2 + ⋯ + b k ∑ i = 1 n x i 1 x i k … … … . . ∑ i = 1 n x i k y i = b 0 ∑ i = 1 n x i k + b 1 ∑ i = 1 n x i k x i 1 + ⋯ + b k ∑ i = 1 n x i k 2 证明:由k个正则方程式可得\\ \begin{array}{l} \sum_{i=1}^{n} y_{i} = b_{0}+b_{1} x_{i 1}+\cdots+b_{k} x_{i k} \\ \sum_{i=1}^{n} x_{i 1} y_{i} = b_{0}\sum_{i=1}^n1·x_{i1}+b_{1} \sum_{i=1}^nx_{i 1}^2+\cdots+b_{k} \sum_{i=1}^nx_{i1}x_{i k} \\ \ldots \ldots \ldots . . \\ \sum_{i=1}^{n} x_{i k} y_{i}=b_{0}\sum_{i=1}^nx_{ik}+b_{1} \sum_{i=1}^nx_{ik}x_{i 1}+\cdots+b_{k}\sum_{i=1}^n x_{i k}^2\\ \end{array} ki=1nyi=b0+b1xi1++bkxiki=1nxi1yi=b0i=1n1xi1+b1i=1nxi12++bki=1nxi1xik..i=1nxikyi=b0i=1nxik+b1i=1nxikxi1++bki=1nxik2
所 以 X ′ Y ( k + 1 ) × 1 = ( X ′ X ) ( k + 1 ) × ( k + 1 ) B o l s ( k + 1 ) × 1 即 两 边 同 乘 以 ( X ′ X ) − 1 , 最 终 可 得 B o l s = ( X ′ X ) − 1 X ′ Y 所以 X'Y_{(k+1)\times1} = (X'X)_{(k+1)\times (k+1)}B{ols} _{(k+1)\times 1}\\ 即两边同乘以 (X'X)^{-1},最终可得\\ B_{ols} = (X'X)^{-1}X'Y XY(k+1)×1=(XX)(k+1)×(k+1)Bols(k+1)×1(XX)1Bols=(XX)1XY

4.

帽子矩阵 H = X ( X ′ X ) − 1 X ′ H = X(X'X)^{-1}X' H=X(XX)1X
H 2 = H ⋅ H = X ( X ′ X ) − 1 X ′ ⋅ X ( X ′ X ) − 1 X ′ = X ( X ′ X ) − 1 ( X ′ X ) ( X ′ X ) − 1 X ′ = X ( X ′ X ) − 1 X ′ = H H^{2} = H·H = X(X'X)^{-1}X'·X(X'X)^{-1}X'\\ = X(X'X)^{-1}(X'X)(X'X)^{-1}X'\\ = X(X'X)^{-1}X' = H H2=HH=X(XX)1XX(XX)1X=X(XX)1(XX)(XX)1X=X(XX)1X=H

  • 2
    点赞
  • 1
    收藏
    觉得还不错? 一键收藏
  • 打赏
    打赏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

Cachel wood

你的鼓励将是我创作的最大动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付:¥1
获取中
扫码支付

您的余额不足,请更换扫码支付或充值

打赏作者

实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值