- “2016年9月19日晚”, 穿越吗?
第一类错误是弃真,即$H_0$真,却拒绝$H_0$。
故答案为B.
当然,接受$H_1$的说法也是不太妥当。你的R的工作目录。
用广义最小二乘法,但是关键还是相关系数<bblatex>\rho</bblatex>应该为多少,这不是R的事情,而是计量的事情。
广义最小二乘使用nlme中的gls命令。可惜里面的有个参数我还是没有搞明白。顶起来!希望高手给予关注,当然也是为了R在计量中的发展努力!
DW检验使用
library(lmtest)
dwtest()
即可,无需编程。但只能检验一阶序列相关。
高阶的可用lmtest中
bgtest()即可!李子奈教授计量经济学第三版第四章例题4.2.1中
前面部分关于回归一元、二元很容易实现,后面的对4.2.24式进行序列相关性的lagrange乘数检验即
<bblatex>e_t=\alpha+\beta X+\beta_1 T^2+\beta_2 e_{t-1}</bblatex>的回归如何实现?
以及广义差分的估计如何实现?
我采用dynlm包中的dynlm命令得出的结论和教材完全不一样,而且感觉也不对。利用arima模型好像结果也和教材有较大差异。现奉上我前面的程序及所需要的数据,望高手给出R程序的实现。(还是李的例题有问题?)
library(bstats)
library(dynlm)
bbc<-read.csv("2.6.2.csv",header=TRUE)
attach(bbc)
lm.1<-lm(Y~X)
summary(lm.1)
dw.test(lm.1)
rs<-(resid(lm.1))
plot(bbc$year,rs,type="l",col=2)
abline(h=0,col=4)
rs1<-rs[-1]
rs2<-rs[-length(rs)]
plot(rs2,rs1,col=4,xlab="e(t-1)",ylab="e(t)")
abline(h=0,col=2)
abline(v=0,col=2)
T<-bbc$year-1977
lm.2<-lm(Y~X+I(T^2),data=bbc)
summary(lm.2)
dw.test(lm.2)
res<-resid(lm.2)
lm.3<-dynlm(res~X+I(T^2)+L(res),data=bbc)
summary(lm.3)
结果lm.3的结果与教材完全不同。虽然有Rstudio这么好的编辑器,但是对使用latex+sweave的,用rwindet也不错。
我参考http://cos.name/cn/topic/102932 成功。看了上面各位的讨论,受益不少!谢谢你们!
确实用R2.13.0编译中文没有问题,但是仍然希望哪位高人能给几个比较理想的详细的模板(如windows下中文xetex+加sweave+article;中文xetex+加sweave+beamer等)!
想学ubuntu的,担心太费时间。我在利用Sweave处理包含中文的.Rnw文件(谢益辉提供的)时总是显示
"has unknown encoding: assuming Latin-1"
不通过.后来查到要在Rprofile.site中增加 options(encoding='UTF-8')
可是我加了以后还是问题依旧。
我的是xp sp3, r 2.13.1.而我的另一台同样的系统的本本却可以通过Sweave[s:12]
我卸掉R后重装问题依旧,就差重做系统了,太费时间,不敢!
问题何在?各位给个建议,如何处理!另强烈建议高手(如谢益辉同学等[s:13])抽时间给大家做个windows下ctex版本的带中文的sweave的beamer模板。毕竟windows用的更多!
试过lyx,但是最后还是选择winedt。谢谢bioshaw。
真的很好用!我有一较大的数据集,几万行,有些行数据缺失严重,想删掉该行,如何在r中删除?手工的话工作量太大。
我编的条件循环,如下(其中cbb已经是一数据框)
for(i in 1:24000){
for(j in 1:10){
if(cbb[i,j]==cbb[i,j+1]==cbb[i,j+2]== NA) na.omit(cbb[i,])}}
总是出问题错误: 意外的'==' 于
" for(j in 1:10){
if(cbb[i,j]==cbb[i,j+1]=="改成=也不行。
各位帮忙看看如何实现。谢谢!江苏大学的统计系也可以考虑。
- 大家,值得认真研读。
- 我怎么看不到?威望不够?