作者boshings (bo)
看板Statistics
标题[程式] SAS与R R^2不同~
时间Mon May 25 23:49:15 2015
[软体程式类别]:
请填入软体程式类别，例如：SAS、SPSS、R、EVIEWS...等
SAS、R
[程式问题]:
资料处理、回归、叙述统计、logistic、probit...等
回归
[软体熟悉度]:
请把以下不需要的部份删除
中(3个月到1年)
[问题叙述]:
请问两个软体，同样的估计方法，为什麽R^2差异这麽大~~
[程式范例]:
虽然张贴程式很可怕，但基本上有些程式还是要张贴才能解决
[SAS CODE]
-----------------
DATA A;
INPUT firm year p c @@;
DATALINES;
1 1955 5.36598 1.14867 1 1960 6.03787 1.45185
1 1965 6.37673 1.52257 1 1970 6.93245 1.76627
2 1955 6.54535 1.35041 2 1960 6.69827 1.71109
2 1965 7.40245 2.09519 2 1970 7.82644 2.39480
;
RUN;
PROC SORT DATA=A;
BY firm year;
RUN;
PROC PANEL DATA=A PRINTFIXED;
MODEL p = c/ FIXTWO ;
ID firm year;
RUN;
[OUTPUT]
------------
SSE = 0.0724 DFE = 2
MSE = 0.0362 根 MSE = 0.1903
R 平方 = 0.9825
[R CODE]
-----------
library(plm)
data = read.csv('C:/data.csv',header=TRUE)
pdata = plm.data(data,index=c('firm','year'))
pm1 = plm(p~c, model='within', data=pdata, effect='twoways')
summary(pm1)
[OUTPUT]
-----------
Total Sum of Squares: 0.072469
Residual Sum of Squares: 0.072418
R-Squared : 0.00071469
Adj. R-Squared : 0.00017867
F-statistic: 0.00143041 on 1 and 2 DF, p-value: 0.97327
-----------------------------------------------------------------------------
--
※ 发信站: 批踢踢实业坊(ptt.cc), 来自: 59.115.198.160
※ 文章网址: https://webptt.com/cn.aspx?n=bbs/Statistics/M.1432568958.A.F89.html
※ 编辑: boshings (59.115.198.160), 05/25/2015 23:50:13
1F：→ luenchang: numbers are stored as 8 digit floating in SAS. 05/26 00:10
2F：→ luenchang: not sure how R stores numeric values 05/26 00:10
不是很了解你的意思，可能储存的位数有差，但是应该不会严重影响估计结果
※ 编辑: boshings (59.115.198.160), 05/26/2015 00:28:42
3F：推 tew: 不懂R 但你应该再看看系数值是不是一样 05/26 16:22
4F：→ tew: 说不定你根本是用了两个不同的模型 05/26 16:23
两者估计系数相同，SSE也相同 应该是同一个模型
SAS:
变数 自由度 估计值 标准误差 t 值 Pr > |t|
c 1 -0.02718 0.7188 -0.04 0.9733
R:
Estimate Std. Error t-value Pr(>|t|)
c -0.027184 0.718755 -0.0378 0.9733
※ 编辑: boshings (59.115.156.102), 05/26/2015 18:51:09
5F：推 asdiy: 是不是你的two way 变成因子实验 05/27 15:25
6F：→ asdiy: 觉得是effect的问题 05/27 15:37