楼主: 子衿1219
1033 2

[问答] 线性模型的问题。 [推广有奖]

  • 1关注
  • 4粉丝

讲师

70%

还不是VIP/贵宾

-

威望
0
论坛币
60 个
通用积分
1.0023
学术水平
2 点
热心指数
0 点
信用等级
1 点
经验
13521 点
帖子
326
精华
0
在线时间
431 小时
注册时间
2012-3-19
最后登录
2018-7-8

+2 论坛币
k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群
赵安豆老师微信:zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

感谢您参与论坛问题回答

经管之家送您两个论坛币!

+2 论坛币
程序如下:
library(mvtnorm)
nm<-100
p<-10
rmat<-diag(p)

x<-rmvnorm(nm,mean=rep(0,p),sigma=rmat)
eps<-rnorm(nm,mean=0,1)
y<-x[,1]+x[,2]+x[,3]+x[,4]+x[,5]+eps


a<-lm(y~0+x)
结果如下:
Call:
lm(formula = y ~ 0 + x)

Coefficients:
       x1         x2         x3         x4         x5         x6         x7  
1.164506   0.938315   0.766277   1.091936   0.994549   0.184499  -0.102038  
       x8         x9        x10  
0.194357   0.017017  -0.006841

问题如下:
真实值只与x1-x5有关,为何差别那么大?
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

关键词:线性模型 Library BRARY Norm bra 模型

沙发
水天一色DIY 在职认证  发表于 2013-4-24 16:44:55 |只看作者 |坛友微信交流群
  1. Call:
  2. lm(formula = y ~ 0 + x)

  3. Coefficients:
  4.       x1        x2        x3        x4        x5        x6        x7        x8        x9       x10  
  5. 0.80940   0.95171   0.90152   0.91000   1.06887   0.08463   0.01374   0.07429   0.19484  -0.14911  
复制代码
Coefficients:
    Estimate Std. Error t value Pr(>|t|)   
x1   0.80940    0.11192   7.232 1.51e-10 ***
x2   0.95171    0.10420   9.133 1.82e-14 ***
x3   0.90152    0.10386   8.680 1.59e-13 ***
x4   0.91000    0.09277   9.809 7.12e-16 ***
x5   1.06887    0.11291   9.467 3.67e-15 ***
x6   0.08463    0.11433   0.740   0.4611   
x7   0.01374    0.11019   0.125   0.9010   
x8   0.07429    0.09914   0.749   0.4556   
x9   0.19484    0.09991   1.950   0.0543 .  
x10 -0.14911    0.11487  -1.298   0.1976
结果并不是相差很大,并且其它变量t值没有通过。

使用道具

藤椅
子衿1219 发表于 2013-4-25 09:20:12 |只看作者 |坛友微信交流群
水天一色DIY 发表于 2013-4-24 16:44
Coefficients:
    Estimate Std. Error t value Pr(>|t|)   
x1   0.80940    0.11192   7.232 1.51e-1 ...
奥,在这里啊,谢谢了。

使用道具

您需要登录后才可以回帖 登录 | 我要注册

本版微信群
加好友,备注cda
拉您进交流群

京ICP备16021002-2号 京B2-20170662号 京公网安备 11010802022788号 论坛法律顾问:王进律师 知识产权保护声明   免责及隐私声明

GMT+8, 2024-5-12 21:59