苹果/安卓/wp
积分 70, 距离下一级还需 15 积分
道具: 彩虹炫, 涂鸦板, 雷达卡, 热点灯, 金钱卡, 显身卡下一级可获得
权限: 自定义头衔
购买后可立即获得
权限: 隐身
道具: 金钱卡, 彩虹炫, 雷达卡, 热点灯, 涂鸦板
群里有大牛，了解熟悉Bai和Perron(2002)结构性断点检验的方法么？国外学者已做出相应部分的R程序（sturcchange)，但在具体结合自己数据运用上仍存在问题。可否透露下代码怎么写么？或数据录入通过哪种形势能达到和包自带数据库效果一样。。。求助
支持楼主：、
购买后，论坛将把您花费的资金全部奖励给楼主，以表示您对TA发好贴的支持
载入中......
无限扩大经管职场人脉圈！每天抽选10位免费名额，现在就扫& 论坛VIP& 贵宾会员& 可免费加入
加入我们,立即就学扫码下载「就学」app& Join us!& JoinLearn&
&nbsp&nbsp|
&nbsp&nbsp|
&nbsp&nbsp|
&nbsp&nbsp|
&nbsp&nbsp|
&nbsp&nbsp|
如有投资本站或合作意向，请联系（010-）；
邮箱：service@pinggu.org
投诉或不良信息处理：（010-）
京ICP证090565号
京公网安备号
论坛法律顾问：王进律师苹果/安卓/wp
积分 37, 距离下一级还需 8 积分
道具: 彩虹炫, 涂鸦板, 雷达卡, 热点灯, 金钱卡下一级可获得
道具: 显身卡
购买后可立即获得
权限: 隐身
道具: 金钱卡, 彩虹炫, 雷达卡, 热点灯, 涂鸦板
开心签到天数: 130 天连续签到: 1 天[LV.7]常住居民III
不知道用R应该怎么做断点检验？
载入中......
无限扩大经管职场人脉圈！每天抽选10位免费名额，现在就扫& 论坛VIP& 贵宾会员& 可免费加入
加入我们,立即就学扫码下载「就学」app& Join us!& JoinLearn&
&nbsp&nbsp|
&nbsp&nbsp|
&nbsp&nbsp|
&nbsp&nbsp|
&nbsp&nbsp|
&nbsp&nbsp|
如有投资本站或合作意向，请联系（010-）；
邮箱：service@pinggu.org
投诉或不良信息处理：（010-）
京ICP证090565号
京公网安备号
论坛法律顾问：王进律师苹果/安卓/wp
积分 29, 距离下一级还需 16 积分
道具: 彩虹炫, 涂鸦板, 雷达卡, 热点灯, 金钱卡下一级可获得
道具: 显身卡
购买后可立即获得
权限: 隐身
道具: 金钱卡, 彩虹炫, 雷达卡, 热点灯, 涂鸦板
尴尬签到天数: 1 天连续签到: 1 天[LV.1]初来乍到
我用了adf.test和pp.test函数分别对同一组数据做了单位根检验，但是两个结果完全相反。很疑惑，不知道该用哪种，求大神解答~
& pp.test(pm25)
& & & & Phillips-Perron Unit Root Test
data:&&pm25
Dickey-Fuller Z(alpha) = -69.7489, Truncation lag parameter = 4, p-value = 0.01
alternative hypothesis: stationary
警告信息：
In pp.test(pm25) : p-value smaller than printed p-value
& adf.test(pm25)
& & & & Augmented Dickey-Fuller Test
data:&&pm25
Dickey-Fuller = -1.8316, Lag order = 4, p-value = 0.6465
alternative hypothesis: stationary
载入中......
ADF的P值大，无法拒接原假设，应该使用PP检验的结果
shaoqinglong11 发表于
ADF的P值大，无法拒接原假设，应该使用PP检验的结果我也发现adf.test的p值偏大，但是两个方法相差太多了，接受pp.test的结果，潜意识不是接受了这个序列平稳的假设吗？
noisyman 发表于
我也发现adf.test的p值偏大，但是两个方法相差太多了，接受pp.test的结果，潜意识不是接受了这个序列平稳 ...那你的本意是不是觉得这个序列不可能是平稳的？要不再用DF检验试试吧，我得到的结果一般这两种方法是相同的
我估计你的数据可能有离群值，导致两个检验结果不同，你用dfao命令试试
先排除确定性趋势的可能，之后按照trend-intercept-none顺序做ADF检验。最后别忘了做一下perron检验。
shaoqinglong11 发表于
我估计你的数据可能有离群值，导致两个检验结果不同，你用dfao命令试试被大神说中了，确实有几个离群值。。。我其实在做上海市pm2.5的分析预测，去年年底那几天爆表的时候确实是异常的。但我不知道该怎么处理，是去除了好还是留下好……dfao命令没有搜到，是什么意思呢？
本帖最后由 noisyman 于
21:28 编辑 nuomin 发表于
先排除确定性趋势的可能，之后按照trend-intercept-none顺序做ADF检验。最后别忘了做一下perron检验。求问如何排除确定性趋势呢，是decompose()之后取什么分量吗？
trend-intercept-none顺序是做3组adf检验吗？求大神解答，不甚感激^_^
noisyman 发表于
被大神说中了，确实有几个离群值。。。我其实在做上海市pm2.5的分析预测，去年年底那几天爆表的时候确实是 ...不用去除离群值，设置虚拟变量就可以了。dfao命令在STATA里面有的，我一般用STATA做计量，不晓得你用什么软件。
shaoqinglong11 发表于
不用去除离群值，设置虚拟变量就可以了。dfao命令在STATA里面有的，我一般用STATA做计量，不晓得你用什么 ...我用的是R……可以再多问一句，怎么设置虚拟变量呢？
无限扩大经管职场人脉圈！每天抽选10位免费名额，现在就扫& 论坛VIP& 贵宾会员& 可免费加入
加入我们,立即就学扫码下载「就学」app& Join us!& JoinLearn&
&nbsp&nbsp|
&nbsp&nbsp|
&nbsp&nbsp|
&nbsp&nbsp|
&nbsp&nbsp|
&nbsp&nbsp|
如有投资本站或合作意向，请联系（010-）；
邮箱：service@pinggu.org
投诉或不良信息处理：（010-）
京ICP证090565号
京公网安备号
论坛法律顾问：王进律师苹果/安卓/wp
积分 37, 距离下一级还需 8 积分
道具: 彩虹炫, 涂鸦板, 雷达卡, 热点灯, 金钱卡下一级可获得
道具: 显身卡
购买后可立即获得
权限: 隐身
道具: 金钱卡, 彩虹炫, 雷达卡, 热点灯, 涂鸦板
开心签到天数: 130 天连续签到: 1 天[LV.7]常住居民III
请问R要怎么做断点检验？
支持楼主：、
购买后，论坛将把您花费的资金全部奖励给楼主，以表示您对TA发好贴的支持
载入中......
在gap包里有chow.test复制代码
jiangbeilu 发表于
在gap包里有(=@__@=)请问语句应该怎么写？例如f1为无约束方程，f2、f3为受约束方程
初级热心勋章
初级热心勋章
中级热心勋章
中级热心勋章
初级信用勋章
初级信用勋章
初级学术勋章
初级学术勋章
中级信用勋章
中级信用勋章
高级热心勋章
高级热心勋章
无限扩大经管职场人脉圈！每天抽选10位免费名额，现在就扫& 论坛VIP& 贵宾会员& 可免费加入
加入我们,立即就学扫码下载「就学」app& Join us!& JoinLearn&
&nbsp&nbsp|
&nbsp&nbsp|
&nbsp&nbsp|
&nbsp&nbsp|
&nbsp&nbsp|
&nbsp&nbsp|
如有投资本站或合作意向，请联系（010-）；
邮箱：service@pinggu.org
投诉或不良信息处理：（010-）
京ICP证090565号
京公网安备号
论坛法律顾问：王进律师苹果/安卓/wp
积分 11681, 距离下一级还需 6619 积分
权限: 自定义头衔, 签名中使用图片, 隐身, 设置帖子权限, 设置回复可见, 签名中使用代码
道具: 彩虹炫, 涂鸦板, 雷达卡, 热点灯, 金钱卡, 显身卡, 匿名卡, 抢沙发, 提升卡, 沉默卡, 千斤顶, 变色卡下一级可获得
道具: 置顶卡
购买后可立即获得
权限: 隐身
道具: 金钱卡, 彩虹炫, 雷达卡, 热点灯, 涂鸦板
开心签到天数: 4 天连续签到: 1 天[LV.2]偶尔看看I
本帖最后由 耕耘使者 于
11:49 编辑
这是非参数统计学中的，KS检验是单样本的科尔莫哥罗夫-斯米尔诺夫检验，都是用来检验分布类型的。
具体一些：
& & 1、关于单样本卡方分布，我想问的关键点是连续分布的检验
& & 如，给定一个连续型变量样本值100个，比如零件长度，80，84.2，85.7……
　若用卡方检验，检验它是否服从某一类型分布（如正态分布、t分布等等），该如何用R做？麻烦的办法我会：自己划分一个个小区间，自己统计观察频数，再自己计算期望频数，之后再把这两个变量录入R，进行操作。如果只能这样，我宁可放弃单样本卡方检验，和完全手算差不了多少。
&&2、关于单样本KS分布
& & 查阅了帮助，函数是ks.test，但参数中的分布类型没有列出，案例中，只给了一下&pgamma&, 即，若要检验是否服从gamma分布，需写上这个参数。但其它分布写什么？如正态分布，写&norm&？t分布，写&t&？帮助上没写明，只能依靠猜吗？
& & 请大家解惑，谢谢！
wilcox.test
wilcox检验
binom.test
贝努力试验检验
chisq.test
fisher精确检验
科尔莫哥罗夫-斯米尔诺夫检验
shapiro.test
shapio-wilk正态分布检验
phillips-perron检验
quada.test
friedman.test
friedman秩和检验
更多搜索 help.search("test")
载入中......
t.test&&t检验
wilcox.test&&wilcox检验
prop.test&&
binom.test&&贝努力试验检验
chisq.test&&卡方检验
fish.test&&fisher精确检验
ks.test& & 科尔莫哥罗夫-斯米尔诺夫检验
shapiro.test& & shapio-wilk正态分布检验
pp.test.& &phillips-perron检验
quada.test& &quade检验
friedman.test& &friedman秩和检验
更多搜索 help.search(&test&)
热心帮助其他会员
总评分:&学术水平 + 6&
热心指数 + 6&
信用等级 + 5&
是什么给了你自信
chisq.test
热心帮助其他会员
总评分:&学术水平 + 1&
热心指数 + 1&
ls这些差不多够用了，我再补充下吧，如果是针对正态检验的话还可以考虑nortest包，简单说明一下lillie.test()更精确的Kolmogorov-Smirnov检验。
ad.test()进行Anderson-Darling正态性检验。
cvm.test()进行Cramer-von Mises正态性检验。
pearson.test()进行Pearson卡方正态性检验。
sf.test()进行Shapiro-Francia正态性检验复制代码
热心帮助其他会员
热心帮助其他会员
总评分:&论坛币 + 100&
学术水平 + 6&
热心指数 + 6&
信用等级 + 5&
个人主页，欢迎订阅/
wtxhpx1991 发表于
chisq.test非常感谢，请看主贴。
yywan0913 发表于
t.test&&t检验
wilcox.test&&wilcox检验
prop.test非常感谢，请看主贴。
求证1加1 发表于
ls这些差不多够用了，我再补充下吧，如果是针对正态检验的话还可以考虑nortest包，简单说明一下辛苦版主了，请看主贴。
初级热心勋章
初级热心勋章
初级信用勋章
初级信用勋章
中级热心勋章
中级热心勋章
无限扩大经管职场人脉圈！每天抽选10位免费名额，现在就扫& 论坛VIP& 贵宾会员& 可免费加入
加入我们,立即就学扫码下载「就学」app& Join us!& JoinLearn&
&nbsp&nbsp|
&nbsp&nbsp|
&nbsp&nbsp|
&nbsp&nbsp|
&nbsp&nbsp|
&nbsp&nbsp|
如有投资本站或合作意向，请联系（010-）；
邮箱：service@pinggu.org
投诉或不良信息处理：（010-）
京ICP证090565号
京公网安备号
论坛法律顾问：王进律师}