質問編集履歴
4
maxitの数値を変更した結果、別なエラーが発生した。
test
CHANGED
|
File without changes
|
test
CHANGED
|
@@ -1,6 +1,14 @@
|
|
|
1
|
+
長くなっておりますが、今後同様のエラーが発生した人のためにこの質問の中で試行錯誤しておきます。
|
|
2
|
+
|
|
3
|
+
追記内容については、時間が経ちましたら質問内容下部に移動します。
|
|
4
|
+
|
|
5
|
+
|
|
6
|
+
|
|
7
|
+
|
|
8
|
+
|
|
1
|
-
＝＝追記＝＝
|
|
9
|
+
＝＝11/27 追記＝＝
|
|
2
|
-
|
|
10
|
+
|
|
3
|
-
|
|
11
|
+
### 解決の見込み
|
|
4
12
|
|
|
5
13
|
Rに詳しい知人にエラーについて簡単に調べてもらったところ、glm関数のオプションmaxitの値をいじくることでエラー発生を抑えることができるかもしれないとのことです。
|
|
6
14
|
|
|
@@ -14,6 +22,44 @@
|
|
|
14
22
|
|
|
15
23
|
|
|
16
24
|
|
|
25
|
+
### 結果
|
|
26
|
+
|
|
27
|
+
##### maxit<=19の場合
|
|
28
|
+
|
|
29
|
+
```R
|
|
30
|
+
|
|
31
|
+
#説明変数:距離対数、輸出年度、FTA/WTO、ブランチ、輸出相手国
|
|
32
|
+
|
|
33
|
+
res4 <-glm(Value_USD_r~ log(distw)+year+fta_wto+branch+ImporterISO3,
|
|
34
|
+
|
|
35
|
+
data = dat_est_0201, family = poisson, maxit=19)
|
|
36
|
+
|
|
37
|
+
```
|
|
38
|
+
|
|
39
|
+
##### 別なエラー発生
|
|
40
|
+
|
|
41
|
+
```R
|
|
42
|
+
|
|
43
|
+
エラー: inner loop 1; cannot correct step size
|
|
44
|
+
|
|
45
|
+
追加情報: 警告メッセージ:
|
|
46
|
+
|
|
47
|
+
step size truncated due to divergence
|
|
48
|
+
|
|
49
|
+
```
|
|
50
|
+
|
|
51
|
+
|
|
52
|
+
|
|
53
|
+
このエラーについて調べてみてもよく分かりませんでした。。。
|
|
54
|
+
|
|
55
|
+
|
|
56
|
+
|
|
57
|
+
### 今後試すこと
|
|
58
|
+
|
|
59
|
+
ImpoterISO3の変数の数が多いということがエラーの原因であると考えられるため、変数の数を少なくしてみてプログラムが正常に走るかを検証してみたいと思います。
|
|
60
|
+
|
|
61
|
+
|
|
62
|
+
|
|
17
63
|
＝＝＝＝＝＝
|
|
18
64
|
|
|
19
65
|
### 前提・実現したいこと
|
|
@@ -80,7 +126,7 @@
|
|
|
80
126
|
|
|
81
127
|
#説明変数：輸出年、距離対数、貿易協定、ブランチ、輸出国
|
|
82
128
|
|
|
83
|
-
res4 <-glm(Value_USD_r~
|
|
129
|
+
res4 <-glm(Value_USD_r~ log(distw)+year+fta_wto+branch+ImporterISO3,
|
|
84
130
|
|
|
85
131
|
data = dat_est_0201, family = poisson)
|
|
86
132
|
|
|
@@ -236,9 +282,7 @@
|
|
|
236
282
|
|
|
237
283
|
```
|
|
238
284
|
|
|
239
|
-
==追記==
|
|
240
|
-
|
|
241
|
-
11/24
|
|
285
|
+
==11/24 追記==
|
|
242
286
|
|
|
243
287
|
#### 各変数データの詳細
|
|
244
288
|
|
3
解決方法の見込みを共有
test
CHANGED
|
File without changes
|
test
CHANGED
|
@@ -1,3 +1,21 @@
|
|
|
1
|
+
＝＝追記＝＝
|
|
2
|
+
|
|
3
|
+
**解決の見込み**
|
|
4
|
+
|
|
5
|
+
Rに詳しい知人にエラーについて簡単に調べてもらったところ、glm関数のオプションmaxitの値をいじくることでエラー発生を抑えることができるかもしれないとのことです。
|
|
6
|
+
|
|
7
|
+
https://stackoverflow.com/questions/20814590/r-na-nan-inf-in-x-error
|
|
8
|
+
|
|
9
|
+
|
|
10
|
+
|
|
11
|
+
今日の午後に走らせてみて、結果を報告します。
|
|
12
|
+
|
|
13
|
+
またmaxitをいじくることによる統計的な影響などについて詳しい方がいましたら、ご教授いただきたいです。
|
|
14
|
+
|
|
15
|
+
|
|
16
|
+
|
|
17
|
+
＝＝＝＝＝＝
|
|
18
|
+
|
|
1
19
|
### 前提・実現したいこと
|
|
2
20
|
|
|
3
21
|
|
|
@@ -220,7 +238,7 @@
|
|
|
220
238
|
|
|
221
239
|
==追記==
|
|
222
240
|
|
|
223
|
-
|
|
241
|
+
11/24(Sunday)
|
|
224
242
|
|
|
225
243
|
#### 各変数データの詳細
|
|
226
244
|
|
2
ImpoterISO3を除去した場合に正常に走ったプログラム結果について追記しました。
test
CHANGED
|
File without changes
|
test
CHANGED
|
@@ -176,6 +176,8 @@
|
|
|
176
176
|
|
|
177
177
|
**→ImporterISO3を説明変数とすることに問題あり**
|
|
178
178
|
|
|
179
|
+
|
|
180
|
+
|
|
179
181
|
|
|
180
182
|
|
|
181
183
|
#### 説明変数の数がデータ数を上回っていないことを確認
|
1
各変数データの詳細について追記しました。
test
CHANGED
|
File without changes
|
test
CHANGED
|
@@ -215,3 +215,63 @@
|
|
|
215
215
|
|
|
216
216
|
|
|
217
217
|
```
|
|
218
|
+
|
|
219
|
+
==追記==
|
|
220
|
+
|
|
221
|
+
|
|
222
|
+
|
|
223
|
+
#### 各変数データの詳細
|
|
224
|
+
|
|
225
|
+
distwのデータの中身についての質問がありましたので、各データ内容の詳細について追記しておきます。
|
|
226
|
+
|
|
227
|
+
（データの分かりやすい見せ方が思いつかなかったので、Excelにて失礼します。）
|
|
228
|
+
|
|
229
|
+
|
|
230
|
+
|
|
231
|
+
##### distw（貿易国間の距離）
|
|
232
|
+
|
|
233
|
+
0データは含まれていません。
|
|
234
|
+
|
|
235
|
+
2267データのそれぞれに、貿易国間のnumデータがあります。
|
|
236
|
+
|
|
237
|
+
|
|
238
|
+
|
|
239
|
+
|
|
240
|
+
|
|
241
|
+
##### Value_USD_r（輸出額）
|
|
242
|
+
|
|
243
|
+
貿易していない国同士も含まれていますので、0データがあります。
|
|
244
|
+
|
|
245
|
+
|
|
246
|
+
|
|
247
|
+
|
|
248
|
+
|
|
249
|
+
##### ImpoterISO3（輸出相手国）
|
|
250
|
+
|
|
251
|
+
198種類の国名コードが含まれています。
|
|
252
|
+
|
|
253
|
+
|
|
254
|
+
|
|
255
|
+
|
|
256
|
+
|
|
257
|
+
##### year（貿易年）
|
|
258
|
+
|
|
259
|
+
2004年から2015年までのデータが含まれています。
|
|
260
|
+
|
|
261
|
+
|
|
262
|
+
|
|
263
|
+
|
|
264
|
+
|
|
265
|
+
##### fta_wto
|
|
266
|
+
|
|
267
|
+
貿易協定の有無ですので、二値データとなっています。
|
|
268
|
+
|
|
269
|
+
|
|
270
|
+
|
|
271
|
+
|
|
272
|
+
|
|
273
|
+
##### branch
|
|
274
|
+
|
|
275
|
+
ブランチの有無ですので、二値データとなっています。
|
|
276
|
+
|
|
277
|
+
|