楼主: oceanlover
5287 7

[问答] 请大家帮忙看看这段代码为啥总是出问题 [推广有奖]

  • 8关注
  • 4粉丝

讲师

86%

还不是VIP/贵宾

-

0%
威望
0 级
论坛币
1272 个
通用积分
10.2974
学术水平
19 点
热心指数
30 点
信用等级
17 点
经验
354 点
帖子
312
精华
0
在线时间
898 小时
注册时间
2006-11-1
最后登录
2024-3-22
楼主
oceanlover 发表于 2014-8-22 16:14:33 |只看作者 |坛友微信交流群|倒序 |AI写论文

+2 论坛币
k人 参与回答

经管之家送您一份

应届毕业生专属福利!

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

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

我在运行Christoph写的关于做Toda-Yamamoto 协整的R代码。具体代码见：

http://www.christophpfeiffer.org/2012/11/07/toda-yamamoto-implementation-in-r/

运行到vcov(lm1)的时候，出现错误提示：

错误于if (is.finite(resvar) && resvar < (mean(f)^2 + var(f)) * 1e-30) warning("essentially perfect fit: summary may be unreliable") :

  需要TRUE/FALSE值的地方不可以用缺少值


应该就是做了回归之后，无法使用summary(), vcov()等函数。但是作者自己解释说，这个程序是没有问题的，能做出结果。

我已经按照作者说的：安装了所有需要的package; 同时重新安装了R程序；但是，这个错误依然存在。


希望大家帮忙看一下，问题在哪里，我应该如何修正？

非常感谢！

二维码

扫码加我 拉你入群

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

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

分享0
关键词：Implementa implement Essential Yamamoto Reliable perfect warning 程序

相关帖子

幸福就是找个温暖的人过一辈子

使用道具

沙发
ywh19860616 发表于 2014-8-22 17:26:35 |只看作者 |坛友微信交流群
  1. > vcov(lm2)
  2.                              (Intercept) arab.l[, 2:8]Arabica.lag-1
  3. (Intercept)                19.5748580946               1.815201e-02
  4. arab.l[, 2:8]Arabica.lag-1  0.0181520084               3.713056e-03
  5. arab.l[, 2:8]Arabica.lag-2 -0.0156665475              -4.506019e-03
  6. arab.l[, 2:8]Arabica.lag-3  0.0013505734               9.849620e-04
  7. arab.l[, 2:8]Arabica.lag-4  0.0003868079              -9.099862e-05
  8. arab.l[, 2:8]Arabica.lag-5  0.0035104653               2.090135e-05
  9. arab.l[, 2:8]Arabica.lag-6 -0.0042234735              -2.314729e-04
  10. arab.l[, 2:8]Arabica.lag-7 -0.0051147654               2.821734e-04
  11. robu.l[, 2:8]Robusta.lag-1 -0.0481667612              -4.153139e-03
  12. robu.l[, 2:8]Robusta.lag-2  0.0275693685               5.452330e-03
  13. robu.l[, 2:8]Robusta.lag-3 -0.0078572398              -1.997996e-03
  14. robu.l[, 2:8]Robusta.lag-4 -0.0058464032               2.594802e-04
  15. robu.l[, 2:8]Robusta.lag-5 -0.0021254251               4.325110e-04
  16. robu.l[, 2:8]Robusta.lag-6  0.0174339934               6.108248e-04
  17. robu.l[, 2:8]Robusta.lag-7 -0.0292793644              -8.334505e-04
  18. index(arab)                -0.0398738302              -1.050853e-04
  19.                            arab.l[, 2:8]Arabica.lag-2
  20. (Intercept)                             -1.566655e-02
  21. arab.l[, 2:8]Arabica.lag-1              -4.506019e-03
  22. arab.l[, 2:8]Arabica.lag-2               9.112665e-03
  23. arab.l[, 2:8]Arabica.lag-3              -5.700914e-03
  24. arab.l[, 2:8]Arabica.lag-4               1.261055e-03
  25. arab.l[, 2:8]Arabica.lag-5              -6.507249e-05
  26. arab.l[, 2:8]Arabica.lag-6              -1.340295e-05
  27. arab.l[, 2:8]Arabica.lag-7              -1.168445e-04
  28. robu.l[, 2:8]Robusta.lag-1               5.151479e-03
  29. robu.l[, 2:8]Robusta.lag-2              -1.082615e-02
  30. robu.l[, 2:8]Robusta.lag-3               7.907287e-03
  31. robu.l[, 2:8]Robusta.lag-4              -2.496323e-03
  32. robu.l[, 2:8]Robusta.lag-5              -3.202493e-04
  33. robu.l[, 2:8]Robusta.lag-6               2.512433e-04
  34. robu.l[, 2:8]Robusta.lag-7               4.074080e-04
  35. index(arab)                              4.320115e-05
  36.                            arab.l[, 2:8]Arabica.lag-3
  37. (Intercept)                              1.350573e-03
  38. arab.l[, 2:8]Arabica.lag-1               9.849620e-04
  39. arab.l[, 2:8]Arabica.lag-2              -5.700914e-03
  40. arab.l[, 2:8]Arabica.lag-3               9.474094e-03
  41. arab.l[, 2:8]Arabica.lag-4              -5.923738e-03
  42. arab.l[, 2:8]Arabica.lag-5               1.271174e-03
  43. arab.l[, 2:8]Arabica.lag-6               2.741945e-05
  44. arab.l[, 2:8]Arabica.lag-7              -1.231262e-04
  45. robu.l[, 2:8]Robusta.lag-1              -1.600254e-03
  46. robu.l[, 2:8]Robusta.lag-2               7.264122e-03
  47. robu.l[, 2:8]Robusta.lag-3              -1.167490e-02
  48. robu.l[, 2:8]Robusta.lag-4               8.218547e-03
  49. robu.l[, 2:8]Robusta.lag-5              -2.403915e-03
  50. robu.l[, 2:8]Robusta.lag-6              -4.343963e-04
  51. robu.l[, 2:8]Robusta.lag-7               6.154511e-04
  52. index(arab)                             -5.174952e-06
  53.                            arab.l[, 2:8]Arabica.lag-4
  54. (Intercept)                              3.868079e-04
  55. arab.l[, 2:8]Arabica.lag-1              -9.099862e-05
  56. arab.l[, 2:8]Arabica.lag-2               1.261055e-03
  57. arab.l[, 2:8]Arabica.lag-3              -5.923738e-03
  58. arab.l[, 2:8]Arabica.lag-4               9.525807e-03
  59. arab.l[, 2:8]Arabica.lag-5              -6.001868e-03
  60. arab.l[, 2:8]Arabica.lag-6               1.460259e-03
  61. arab.l[, 2:8]Arabica.lag-7              -2.089740e-04
  62. robu.l[, 2:8]Robusta.lag-1               6.868756e-04
  63. robu.l[, 2:8]Robusta.lag-2              -2.704803e-03
  64. robu.l[, 2:8]Robusta.lag-3               7.797930e-03
  65. robu.l[, 2:8]Robusta.lag-4              -1.176794e-02
  66. robu.l[, 2:8]Robusta.lag-5               8.283075e-03
  67. robu.l[, 2:8]Robusta.lag-6              -2.792592e-03
  68. robu.l[, 2:8]Robusta.lag-7               4.744160e-04
  69. index(arab)                             -7.042079e-06
  70.                            arab.l[, 2:8]Arabica.lag-5
  71. (Intercept)                              3.510465e-03
  72. arab.l[, 2:8]Arabica.lag-1               2.090135e-05
  73. arab.l[, 2:8]Arabica.lag-2              -6.507249e-05
  74. arab.l[, 2:8]Arabica.lag-3               1.271174e-03
  75. arab.l[, 2:8]Arabica.lag-4              -6.001868e-03
  76. arab.l[, 2:8]Arabica.lag-5               9.485777e-03
  77. arab.l[, 2:8]Arabica.lag-6              -5.862175e-03
  78. arab.l[, 2:8]Arabica.lag-7               1.195262e-03
  79. robu.l[, 2:8]Robusta.lag-1              -1.898241e-04
  80. robu.l[, 2:8]Robusta.lag-2               8.748873e-04
  81. robu.l[, 2:8]Robusta.lag-3              -2.775130e-03
  82. robu.l[, 2:8]Robusta.lag-4               7.845980e-03
  83. robu.l[, 2:8]Robusta.lag-5              -1.165310e-02
  84. robu.l[, 2:8]Robusta.lag-6               7.987624e-03
  85. robu.l[, 2:8]Robusta.lag-7              -2.150191e-03
  86. index(arab)                             -1.678884e-05
  87.                            arab.l[, 2:8]Arabica.lag-6
  88. (Intercept)                             -4.223473e-03
  89. arab.l[, 2:8]Arabica.lag-1              -2.314729e-04
  90. arab.l[, 2:8]Arabica.lag-2              -1.340295e-05
  91. arab.l[, 2:8]Arabica.lag-3               2.741945e-05
  92. arab.l[, 2:8]Arabica.lag-4               1.460259e-03
  93. arab.l[, 2:8]Arabica.lag-5              -5.862175e-03
  94. arab.l[, 2:8]Arabica.lag-6               8.928939e-03
  95. arab.l[, 2:8]Arabica.lag-7              -4.425125e-03
  96. robu.l[, 2:8]Robusta.lag-1               3.708556e-04
  97. robu.l[, 2:8]Robusta.lag-2              -2.917550e-04
  98. robu.l[, 2:8]Robusta.lag-3               7.957030e-04
  99. robu.l[, 2:8]Robusta.lag-4              -2.873956e-03
  100. robu.l[, 2:8]Robusta.lag-5               7.474386e-03
  101. robu.l[, 2:8]Robusta.lag-6              -1.057110e-02
  102. robu.l[, 2:8]Robusta.lag-7               5.234833e-03
  103. index(arab)                              3.756914e-05
  104.                            arab.l[, 2:8]Arabica.lag-7
  105. (Intercept)                             -5.114765e-03
  106. arab.l[, 2:8]Arabica.lag-1               2.821734e-04
  107. arab.l[, 2:8]Arabica.lag-2              -1.168445e-04
  108. arab.l[, 2:8]Arabica.lag-3              -1.231262e-04
  109. arab.l[, 2:8]Arabica.lag-4              -2.089740e-04
  110. arab.l[, 2:8]Arabica.lag-5               1.195262e-03
  111. arab.l[, 2:8]Arabica.lag-6              -4.425125e-03
  112. arab.l[, 2:8]Arabica.lag-7               3.737247e-03
  113. robu.l[, 2:8]Robusta.lag-1              -5.105251e-04
  114. robu.l[, 2:8]Robusta.lag-2               3.265068e-04
  115. robu.l[, 2:8]Robusta.lag-3              -7.943013e-05
  116. robu.l[, 2:8]Robusta.lag-4               7.900997e-04
  117. robu.l[, 2:8]Robusta.lag-5              -1.881306e-03
  118. robu.l[, 2:8]Robusta.lag-6               5.115619e-03
  119. robu.l[, 2:8]Robusta.lag-7              -4.114939e-03
  120. index(arab)                             -7.906184e-05
  121.                            robu.l[, 2:8]Robusta.lag-1
  122. (Intercept)                             -0.0481667612
  123. arab.l[, 2:8]Arabica.lag-1              -0.0041531385
  124. arab.l[, 2:8]Arabica.lag-2               0.0051514790
  125. arab.l[, 2:8]Arabica.lag-3              -0.0016002538
  126. arab.l[, 2:8]Arabica.lag-4               0.0006868756
  127. arab.l[, 2:8]Arabica.lag-5              -0.0001898241
  128. arab.l[, 2:8]Arabica.lag-6               0.0003708556
  129. arab.l[, 2:8]Arabica.lag-7              -0.0005105251
  130. robu.l[, 2:8]Robusta.lag-1               0.0071274451
  131. robu.l[, 2:8]Robusta.lag-2              -0.0095975455
  132. robu.l[, 2:8]Robusta.lag-3               0.0038516554
  133. robu.l[, 2:8]Robusta.lag-4              -0.0010192720
  134. robu.l[, 2:8]Robusta.lag-5              -0.0002690005
  135. robu.l[, 2:8]Robusta.lag-6              -0.0008009697
  136. robu.l[, 2:8]Robusta.lag-7               0.0011112874
  137. index(arab)                              0.0001761388
  138.                            robu.l[, 2:8]Robusta.lag-2
  139. (Intercept)                              2.756937e-02
  140. arab.l[, 2:8]Arabica.lag-1               5.452330e-03
  141. arab.l[, 2:8]Arabica.lag-2              -1.082615e-02
  142. arab.l[, 2:8]Arabica.lag-3               7.264122e-03
  143. arab.l[, 2:8]Arabica.lag-4              -2.704803e-03
  144. arab.l[, 2:8]Arabica.lag-5               8.748873e-04
  145. arab.l[, 2:8]Arabica.lag-6              -2.917550e-04
  146. arab.l[, 2:8]Arabica.lag-7               3.265068e-04
  147. robu.l[, 2:8]Robusta.lag-1              -9.597546e-03
  148. robu.l[, 2:8]Robusta.lag-2               1.992425e-02
  149. robu.l[, 2:8]Robusta.lag-3              -1.473979e-02
  150. robu.l[, 2:8]Robusta.lag-4               5.397486e-03
  151. robu.l[, 2:8]Robusta.lag-5              -5.579183e-04
  152. robu.l[, 2:8]Robusta.lag-6               9.589586e-05
  153. robu.l[, 2:8]Robusta.lag-7              -7.068956e-04
  154. index(arab)                             -8.355191e-05
  155.                            robu.l[, 2:8]Robusta.lag-3
  156. (Intercept)                             -7.857240e-03
  157. arab.l[, 2:8]Arabica.lag-1              -1.997996e-03
  158. arab.l[, 2:8]Arabica.lag-2               7.907287e-03
  159. arab.l[, 2:8]Arabica.lag-3              -1.167490e-02
  160. arab.l[, 2:8]Arabica.lag-4               7.797930e-03
  161. arab.l[, 2:8]Arabica.lag-5              -2.775130e-03
  162. arab.l[, 2:8]Arabica.lag-6               7.957030e-04
  163. arab.l[, 2:8]Arabica.lag-7              -7.943013e-05
  164. robu.l[, 2:8]Robusta.lag-1               3.851655e-03
  165. robu.l[, 2:8]Robusta.lag-2              -1.473979e-02
  166. robu.l[, 2:8]Robusta.lag-3               2.202823e-02
  167. robu.l[, 2:8]Robusta.lag-4              -1.546093e-02
  168. robu.l[, 2:8]Robusta.lag-5               5.227887e-03
  169. robu.l[, 2:8]Robusta.lag-6              -4.946695e-04
  170. robu.l[, 2:8]Robusta.lag-7              -3.609396e-04
  171. index(arab)                              2.497706e-05
  172.                            robu.l[, 2:8]Robusta.lag-4
  173. (Intercept)                             -5.846403e-03
  174. arab.l[, 2:8]Arabica.lag-1               2.594802e-04
  175. arab.l[, 2:8]Arabica.lag-2              -2.496323e-03
  176. arab.l[, 2:8]Arabica.lag-3               8.218547e-03
  177. arab.l[, 2:8]Arabica.lag-4              -1.176794e-02
  178. arab.l[, 2:8]Arabica.lag-5               7.845980e-03
  179. arab.l[, 2:8]Arabica.lag-6              -2.873956e-03
  180. arab.l[, 2:8]Arabica.lag-7               7.900997e-04
  181. robu.l[, 2:8]Robusta.lag-1              -1.019272e-03
  182. robu.l[, 2:8]Robusta.lag-2               5.397486e-03
  183. robu.l[, 2:8]Robusta.lag-3              -1.546093e-02
  184. robu.l[, 2:8]Robusta.lag-4               2.213537e-02
  185. robu.l[, 2:8]Robusta.lag-5              -1.538152e-02
  186. robu.l[, 2:8]Robusta.lag-6               5.455078e-03
  187. robu.l[, 2:8]Robusta.lag-7              -1.088530e-03
  188. index(arab)                              2.188369e-05
  189.                            robu.l[, 2:8]Robusta.lag-5
  190. (Intercept)                             -2.125425e-03
  191. arab.l[, 2:8]Arabica.lag-1               4.325110e-04
  192. arab.l[, 2:8]Arabica.lag-2              -3.202493e-04
  193. arab.l[, 2:8]Arabica.lag-3              -2.403915e-03
  194. arab.l[, 2:8]Arabica.lag-4               8.283075e-03
  195. arab.l[, 2:8]Arabica.lag-5              -1.165310e-02
  196. arab.l[, 2:8]Arabica.lag-6               7.474386e-03
  197. arab.l[, 2:8]Arabica.lag-7              -1.881306e-03
  198. robu.l[, 2:8]Robusta.lag-1              -2.690005e-04
  199. robu.l[, 2:8]Robusta.lag-2              -5.579183e-04
  200. robu.l[, 2:8]Robusta.lag-3               5.227887e-03
  201. robu.l[, 2:8]Robusta.lag-4              -1.538152e-02
  202. robu.l[, 2:8]Robusta.lag-5               2.179263e-02
  203. robu.l[, 2:8]Robusta.lag-6              -1.464328e-02
  204. robu.l[, 2:8]Robusta.lag-7               3.916221e-03
  205. index(arab)                              1.817761e-05
  206.                            robu.l[, 2:8]Robusta.lag-6
  207. (Intercept)                              1.743399e-02
  208. arab.l[, 2:8]Arabica.lag-1               6.108248e-04
  209. arab.l[, 2:8]Arabica.lag-2               2.512433e-04
  210. arab.l[, 2:8]Arabica.lag-3              -4.343963e-04
  211. arab.l[, 2:8]Arabica.lag-4              -2.792592e-03
  212. arab.l[, 2:8]Arabica.lag-5               7.987624e-03
  213. arab.l[, 2:8]Arabica.lag-6              -1.057110e-02
  214. arab.l[, 2:8]Arabica.lag-7               5.115619e-03
  215. robu.l[, 2:8]Robusta.lag-1              -8.009697e-04
  216. robu.l[, 2:8]Robusta.lag-2               9.589586e-05
  217. robu.l[, 2:8]Robusta.lag-3              -4.946695e-04
  218. robu.l[, 2:8]Robusta.lag-4               5.455078e-03
  219. robu.l[, 2:8]Robusta.lag-5              -1.464328e-02
  220. robu.l[, 2:8]Robusta.lag-6               1.966827e-02
  221. robu.l[, 2:8]Robusta.lag-7              -9.515665e-03
  222. index(arab)                             -7.747491e-05
  223.                            robu.l[, 2:8]Robusta.lag-7   index(arab)
  224. (Intercept)                             -0.0292793644 -3.987383e-02
  225. arab.l[, 2:8]Arabica.lag-1              -0.0008334505 -1.050853e-04
  226. arab.l[, 2:8]Arabica.lag-2               0.0004074080  4.320115e-05
  227. arab.l[, 2:8]Arabica.lag-3               0.0006154511 -5.174952e-06
  228. arab.l[, 2:8]Arabica.lag-4               0.0004744160 -7.042079e-06
  229. arab.l[, 2:8]Arabica.lag-5              -0.0021501906 -1.678884e-05
  230. arab.l[, 2:8]Arabica.lag-6               0.0052348327  3.756914e-05
  231. arab.l[, 2:8]Arabica.lag-7              -0.0041149390 -7.906184e-05
  232. robu.l[, 2:8]Robusta.lag-1               0.0011112874  1.761388e-04
  233. robu.l[, 2:8]Robusta.lag-2              -0.0007068956 -8.355191e-05
  234. robu.l[, 2:8]Robusta.lag-3              -0.0003609396  2.497706e-05
  235. robu.l[, 2:8]Robusta.lag-4              -0.0010885304  2.188369e-05
  236. robu.l[, 2:8]Robusta.lag-5               0.0039162207  1.817761e-05
  237. robu.l[, 2:8]Robusta.lag-6              -0.0095156655 -7.747491e-05
  238. robu.l[, 2:8]Robusta.lag-7               0.0071334375  1.465306e-04
  239. index(arab)                              0.0001465306  1.452517e-04
  240. > wald.test(b=coef(lm2), Sigma=vcov(lm2), Terms= c(2:7),df=6)
  241. Wald test:
  242. ----------

  244. Chi-squared test:
  245. X2 = 12.3, df = 6, P(> X2) = 0.056

  247. F test:
  248. W = 2.0, df1 = 6, df2 = 6, P(> W) = 0.2
复制代码


利用R 3.0.2运行作者提供的例子，可以正常运行。

使用道具

藤椅
Aharach 发表于 2014-8-24 06:39:28 |只看作者 |坛友微信交流群
在判断条件里面加个 !is.null(resvar)

使用道具

板凳
oceanlover 发表于 2014-8-24 19:30:04 来自手机 |只看作者 |坛友微信交流群
ywh19860616 发表于 2014-8-22 17:26
利用R 3.0.2运行作者提供的例子，可以正常运行。
谢谢！我现在用的是cran网站是的当前版本，也按照提示安装了。为什么运行时就出错？

使用道具

报纸
oceanlover 发表于 2014-8-24 19:31:42 来自手机 |只看作者 |坛友微信交流群
Aharach 发表于 2014-8-24 06:39
在判断条件里面加个 !is.null(resvar)
谢谢，但不知要加到哪里？是vcov()函数里面么？

使用道具

地板
Aharach 发表于 2014-8-24 21:35:24 |只看作者 |坛友微信交流群
oceanlover 发表于 2014-8-24 19:31
谢谢，但不知要加到哪里？是vcov()函数里面么？
  1. if (!is.null(resvar) && is.finite(resvar) && resvar < (mean(f)^2 + var(f)) * 1e-30)
复制代码

使用道具

7楼
yoghurtim 发表于 2014-9-16 14:14:43 |只看作者 |坛友微信交流群
ywh19860616 发表于 2014-8-22 17:26
利用R 3.0.2运行作者提供的例子，可以正常运行。
我也遇到这个问题了，请问楼主怎么解决的？

使用道具

8楼
chongchong719 发表于 2018-3-15 18:14:47 |只看作者 |坛友微信交流群
我也遇到同样问题, 后来发现, 是自己的数据中有NA

使用道具

返回列表
高级模式
B Color Image Link Quote Code Smilies
您需要登录后才可以回帖 登录 | 我要注册

本版微信群
加好友,备注cda
拉您进交流群
  • 手机版  |
  • 意见反馈  |
  • 帮助  |
  • 新手入门  |
  • 用户手册  |
  • 友情链接  |

    • 如有投资本站、合作意向或投放广告，请联系：13661292478（刘老师）

    联系客服

    邮箱：service@pinggu.org 投诉或不良信息处理：（010-68466864）

    010-68466864

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

    GMT+8, 2024-5-10 00:31