adaptcw
各位大虾:
这是我用matlab计算出的结果:
Estimate t value Pr(>|t|)
-0.842533 -0.214888 0.385611
2.77996 2.61703 0.0169656
0.501364 1.87157 0.0718219
0.434615 2.29138 0.0331347
-0.542054 -3.02443 0.00685769
这是r语言计算出的结果:
Estimate Std. Error t value Pr(>|t|)
(Intercept) -0.8423 3.9208 -0.215 0.83164
thirdrate 2.7799 1.0623 2.617 0.01484 *
rain 0.5014 0.2679 1.872 0.07301 .
tem 0.4346 0.1897 2.291 0.03064 *
popdis -0.5421 0.1792 -3.024 0.00569 **
常数项的Pr(>|t|) 相差很多,其他的都差不多,是什么原因呢?
是因为常数项的自由度和别的不一样吗?
bjt
标题太简单了,改改吧。
进来一看原来是 关于 matlab 的
adaptcw
该改成什么呢,其实我的目的就是问问常数项的自由度
yihui
各个系数t统计量的自由度都一样。
这两个结果差异也太离谱了,能否贴出原始数据?
xifan
这个还蛮逗的,俺上次发现excel算出的数两台电脑上不一样,我的电脑上的是错的
ddtd
原始数据贴出来看看吧…………
yihui
Excel有时候确实会抽风……
adaptcw
以下是原始数据,各位帮帮忙,谢谢:
index thirdrate rain tem popdis
5.19645 3.31874 -0.9143780 -0.3867570 6.77567
4.88410 3.63093 -1.2435400 -0.5207740 7.38728
5.33514 3.47513 0.4424980 0.4748460 7.03123
2.94296 3.59006 -1.2623000 -0.6625880 7.55810
3.03704 3.54279 -1.1425400 -0.6169150 7.58194
3.75322 3.37406 -0.9813460 -0.3615100 7.44702
5.77013 3.74744 -1.3330400 -0.7842770 7.77846
5.96275 3.46841 -0.0270493 0.0622238 5.08090
4.78057 3.33679 -0.6238570 -0.3144510 6.55534
5.57233 3.38034 -0.0987228 1.3464900 7.00288
4.70231 3.46197 -0.8795290 -0.3680310 6.97937
3.72059 3.51136 -1.0538600 -0.4487600 7.27216
5.69151 3.49742 -0.2772000 0.6698870 6.76713
4.51710 3.51925 -0.8669980 -0.3525630 6.74034
4.65238 3.53201 -1.1297500 -0.5027710 7.22100
5.01177 3.63019 -0.3082600 0.2053370 6.27553
6.26186 3.42406 0.4782370 0.8186960 5.87224
7.03251 3.60559 0.4119920 0.2326550 6.82934
6.94610 3.68977 0.5267760 3.8208900 7.23124
4.42417 3.52447 -0.4156290 -0.2286050 6.09310
5.64821 3.55205 -0.0668519 0.1698070 6.14830
5.13792 3.64310 -0.3572780 -0.0293571 6.79693
5.06867 3.81181 -1.0148700 -0.4142390 6.98886
4.08945 3.43173 -0.6073550 0.1405410 7.30870
7.04510 3.75087 -0.1215270 -0.1659990 4.71146
6.85623 3.70497 0.2733040 2.2977000 8.02337
5.30982 3.56591 1.4652500 0.4654830 8.00711
3.02700 3.44311 -0.8934410 -0.3665120 7.65683
4.92527 3.48800 -1.1235000 -0.4204090 7.13565
4.70903 3.57289 -0.8282970 -0.3865540 7.23980
yihui
毫无疑问,自由度是30-4-1=25,上面MatLab的结果不对,但我怀疑并不是MatLab本身的问题(它不可能连回归都算不对),而是你自己的问题,比如数据的问题等。
adaptcw
结果相差的是不多的,括号中的是我的计算结果,前面系数的估计和t value值都相差不大
p值稍微有点偏差,可能是小数位数的关系,但但第一项差太多了(加粗显示的部分):
Estimate t value Pr(>|t|)
(Intercept) -0.8423(-0.842533) -0.215(-0.214888) 0.83164 (0.385611)
thirdrate 2.7799(2.77996) 2.617(2.61703) 0.01484 * (0.0169656)
rain 0.5014(0.501364) 1.872(1.87157) 0.07301 . (0.0718219)
tem 0.4346(0.434615) 2.291(2.29138) 0.03064 * (0.0331347)
popdis -0.5421(-0.542054) -3.024(-3.02443) 0.00569 **(0.00685769)
adaptcw
自己解答一下吧,是我函数用错了,matlab中求p值得函数为tcdf,我用成tpdf了
yihui
我就说MatLab不可能连回归都算不对嘛
xifan
貌似大家毕业论文都写完了