Salmon cycle
安井翔太	
 /	
  Shota	
  Yasui	
  
Twitter:	
  @housecat442	
  
	
  
	
  
	
  
<けーれき>	
  

¡ 

日本のド文系経済学部	
  
アメリカで計量経済学1年	
  
ノルウェーで資源・環境経済学修士取得	
  
SNF研究所でデータ分析(環境税作成)	
  

¡ 

ネット広告代理店でデータ分析	
  

¡ 
¡ 
¡ 

§ 

アトリビューション分析	
  

§ 

マス広告評価	
  
動画広告評価	
  
Etc…	
  

§ 
§ 

	
  
Salmon cycle
!"2)-(,:)"&)2&
&%8/&(,-9/%

;/<=>/&
"/?%&="#$%&
,()$"%

!"#$%&>/<='=)"

!"2)-(,:)"&)2&
&%8/&(,-9/%

!"#$%&>/<='=)"

7$##*6&%)&%8/&(,-9/%

!"#$%&'()*%

!"#$%&'()*%

!"#$%&>/<='=)"

!"#$%&'()*%

;/<=>/&
"/?%&="#$%&
,()$"%

!"#$%&>/<='=)"

7$##*6&%)&%8/&(,-9/%

+,-./'%&',*()"&

0-)1&2)-&
3456/,-'

!"2)-(,:)"&)2&
&%8/&(,-9/%

;/<=>/&
"/?%&="#$%&
,()$"%

!"2)-(,:)"&)2&
&%8/&(,-9/%

!"#$%&>/<='=)"

7$##*6&%)&%8/&(,-9/%

!"#$%&'()*%

+,-./'%&',*()"&

0-)1&2)-&
3456/,-'

7$##*6&%)&%8/&(,-9/%

!"#$%&'()*%

+,-./'%&',*()"&

0-)1&2)-&
3456/,-'

!"2)-(,:)"&)2&
&%8/&(,-9/%

!"#$%&'()*%

+,-./'%&',*()"&

0-)1&2)-&
3456/,-'

!"2)-(,:)"&)2&
&%8/&(,-9/%

;/<=>/&
"/?%&="#$%&
,()$"%

!"#$%&>/<='=)"

7$##*6&%)&%8/&(,-9/%

!"#$%&'()*%

0-)1&2)-&
3456/,-'

;/<=>/&
"/?%&="#$%&
,()$"%

!"#$%&>/<='=)"

7$##*6&%)&%8/&(,-9/%

+,-./'%&',*()"&

0-)1&2)-&
3456/,-'

;/<=>/&
"/?%&="#$%&
,()$"%

!"#$%&>/<='=)"

+,-./'%&',*()"&

0-)1&2)-&
3456/,-'

;/<=>/&
"/?%&="#$%&
,()$"%

!"2)-(,:)"&)2&
&%8/&(,-9/%

7$##*6&%)&%8/&(,-9/%

+,-./'%&',*()"&

!"2)-(,:)"&)2&
&%8/&(,-9/%

;/<=>/&
"/?%&="#$%&
,()$"%

;/<=>/&
"/?%&="#$%&
,()$"%

!"#$%&>/<='=)"

7$##*6&%)&%8/&(,-9/%

!"#$%&'()*%

!"2)-(,:)"&)2&
&%8/&(,-9/%

+,-./'%&',*()"&

0-)1&2)-&
3456/,-'

;/<=>/&
"/?%&="#$%&
,()$"%

!"#$%&>/<='=)"

7$##*6&%)&%8/&(,-9/%

!"#$%&'()*%

!"2)-(,:)"&)2&
&%8/&(,-9/%

+,-./'%&',*()"&

0-)1&2)-&
3456/,-'

;/<=>/&
"/?%&="#$%&
,()$"%

!"#$%&>/<='=)"

7$##*6&%)&%8/&(,-9/%

!"#$%&'()*%

!"2)-(,:)"&)2&
&%8/&(,-9/%

+,-./'%&',*()"&

0-)1&2)-&
3456/,-'

!"#$%&>/<='=)"

7$##*6&%)&%8/&(,-9/%

!"#$%&'()*%

!"2)-(,:)"&)2&
&%8/&(,-9/%

+,-./'%&',*()"&

0-)1&2)-&
3456/,-'

7$##*6&%)&%8/&(,-9/%

!"#$%&'()*%

+,-./'%&',*()"&

0-)1&2)-&
3456/,-'

;/<=>/
"/?%&="#
,()$"
¡  価格を観て投入量を決定	
  
¡  2年餌をあげる	
  
¡  収穫	
  
¡  販売(価格決定)	
  
¡  価格を観て投入量を決定	
  
¡  2年餌をあげる	
  
¡  収穫	
  
¡  販売(価格決定)	
  
¡  養殖業者は価格予想が出来ない	
  
§  よって今日の価格が出荷時期の価格と同じで

あると言う想定の元に生産規模を決定する	
  
¡  価格が低いと生産を減らす	
  
¡  価格が高いと生産を増やす	
  
§  価格を共有するから皆同じ決断!	
  
価格が供給量の増減を受け上下してしまう	
 

この構造を推定して、価格の周期性を説明したい!
¡  考えられる方法は2つ	
  
§  計量時系列分析のアプローチ	
  
▪  仕組みが解ってなくとも取り敢えず分析可能	
  
▪  養殖業のおじさん達は信じてくれない	
  
▪  SVECM→IRFでいける	
  

§  計量経済学のアプローチ	
  
▪  仕組みを良くモデル化出来てないと藻屑	
  
▪  市場をモデル化するので、市場に対する介入政策な
んかを考えるのに使える。	
  
▪  あとおじさん達を説得しやすい。
¡  考えられる方法は2つ	
  
§  計量時系列分析	
  
▪  仕組みが解らなくとも取り敢えず予測可能	
  
▪  養殖業のおじさん達は信じてくれない	
  
§  計量経済学のアプローチ	
  
▪  仕組みを良くモデル化出来てないと藻屑	
  
▪  市場をモデル化するので、市場に対する介入政策な
んかを考えるのに使える。	
  
▪  あとおじさん達を説得しやすい。
データ:ノルウェー地方別パネルデータ	
  
期間:1990~2012年	
 年別データ(annual)	
  
稚魚数決定モデル	
 
!"

=!+!∗

!"!!

+ ! ∗ ! + !! + !!"

去年のサーモン価格で稚魚の数を決める	
 

出荷量決定モデル	
 
!"

=!+!∗

!"!!

+ ℎ! + !!"

二年前の稚魚が成長して今日の出荷量になる
two years before. Basically this equation describes that how many percentage of
salmon is going to be produced when the smolt input is increased by 1%.
Combining these above two equations, it is possible to see the effect of the price of
farmed salmon in year t-3 on the supply quantity in year t.

価格決定モデル	
 
Equation3: Demand equation (price determination) model
log !"#$!!" = !! ∗ log !"##$!!" + !! ∗ log!(!"##$%!!! ) + !! ∗ ! + !! + !!"

価格は供給量と為替レートによって決定される	
 
Demand equation (equation3) has price of farmed salmon as the dependent
variable. Explanatory variable
操作変数(Instrumental	
  Variable)	
   is supply quantity, currency rate and time trend.
	
  
Since 	
  養殖場1ライセンス当たりに発生している過剰生産量を使用。	
 
the price is determined by total supply in the market, the total supply amount
is used in here instead of supply per license. Currency rate is the real rate between
France and Norway. Because France is the biggest importer of Norwegian Farmed
library(plm)	
  
growm	
  <-­‐	
  log(qpl)	
  ~	
  log(lag(ipl,2)) #二年前の稚魚の数から今年の供給量	
  
growthmodel	
  <-­‐	
  plm(growm,	
  data	
  =	
  salmon,	
  model	
  ="random")	
  
summary(inputmodel)	
  
	
  
	
  
inputm	
  <-­‐	
  log(ipl)	
  ~	
  lag(log(price),1)	
  +	
  t #稚魚の数の決定モデル	
  
idmodel	
  <-­‐	
  plm(inputm,	
  data	
  =	
  salmon,	
  model="random")	
  
summary(idmodel)	
  
	
  
	
  
supm	
  <-­‐	
  log(price)	
  ~	
  	
  log(rate)	
  +	
  t	
  +	
  log(qpl)	
  |	
  log(qpl)	
  +	
  log(epl)	
  +	
  log(rate)	
  +	
  t #価格決定モデル
supmodel	
  <-­‐	
  plm(supm,	
  data	
  =	
  salmon,	
  model="random",	
  random.method="amemiya")	
  
summary(supmodel)
log(qpl)	
  ~	
  log(lag(ipl,2))	
  
Coefficients	
  :	
  
	
  	
  	
  	
  	
  	
  	
  	
  	
  	
  	
  	
  	
  	
  	
  	
  	
   	
  
	
  Estimate	
  	
  
	
  Std.	
  Error	
  
	
  t-­‐value	
  	
  Pr(>|t|)	
  	
  	
  	
  	
  
log(lag(ipl,	
  2))	
  
	
  1.105494	
  	
  
	
  	
  0.050923	
  	
  
	
  21.709	
  <	
  2.2e-­‐16	
  ***	
  
	
  
稚魚の量を1%増や
	
  
すと生産も1.1%増	
 
log(ipl)	
  ~	
  lag(log(price),1)	
  +	
  t	
  
Coefficients	
  :	
  
	
  	
  	
  	
  	
  	
  	
  	
  	
  	
  	
  	
  	
  	
  	
  	
  	
  	
  	
  	
  
	
  Estimate	
  	
  
	
  Std.	
  Error	
  
	
  t-­‐value	
  	
  Pr(>|t|)	
  	
  	
  	
  	
  
lag(log(price),	
  1)	
   	
  0.2806347
	
  	
  	
  0.0786544	
  	
  
	
  3.5679	
  	
  0.000501	
  ***	
  
t	
  	
  	
  	
  	
  	
  	
  	
  	
  	
  	
  	
  	
  	
  	
  	
  	
  	
  
	
  0.0471874	
  	
  
	
  0.0028297	
  
	
  16.6757	
  <	
  2.2e-­‐16	
  ***	
  
	
  
価格が1%上昇すると稚
魚の投入数が0.28%増	
 
	
  
log(price)	
  ~	
  	
  log(rate)	
  +	
  t+	
  log(qpl)	
  |	
  log(qpl)	
  +	
  log(epl)	
  +	
  log(rate)	
  +	
  t	
  
Coefficients	
  :	
  
	
  	
  	
  	
  	
  	
  	
  	
  	
  	
  	
  
	
  
	
  Estimate	
  	
  
	
  Std.	
  Error	
  
	
  t-­‐value	
  Pr(>|t|)	
  	
  	
  
log(rate)	
  	
  
	
  8.269694	
  	
  	
  
	
  3.190592	
  	
  
	
  2.5919	
  	
  0.01055	
  *	
  
t	
  	
  	
  	
  	
  	
  	
  	
  	
  
	
  
	
  	
  0.232151	
  	
  	
  
	
  0.088975	
  	
  
	
  2.6092	
  	
  0.01006	
  *	
  
log(qpl)	
  	
   	
  
	
  -­‐0.716636	
  	
  	
  
	
  0.337749	
   	
  
	
  -­‐2.1218	
  	
  0.03560	
  *	
  
¡  仮定していた構造による周期性あるの?	
  
§  推定結果的にはYes	
  
§  でも、その影響の大きさはちゃんと検証が必要	
  
§  やっぱりVARモデル?
Salmon cycle

More Related Content

PPTX
R -> Python
PPTX
入門機械学習6章
PDF
PDF
EasyHtmlReportの紹介
PDF
~knitr+pandocではじめる~『R MarkdownでReproducible Research』
PDF
10分で分かるr言語入門ver2 6
PDF
Rによるやさしい統計学第20章「検定力分析によるサンプルサイズの決定」
PDF
TokyoR 第36回LT Rで部分空間法
R -> Python
入門機械学習6章
EasyHtmlReportの紹介
~knitr+pandocではじめる~『R MarkdownでReproducible Research』
10分で分かるr言語入門ver2 6
Rによるやさしい統計学第20章「検定力分析によるサンプルサイズの決定」
TokyoR 第36回LT Rで部分空間法

Viewers also liked (12)

PDF
Trading volume mapping R in recent environment
PDF
Visualization of Supervised Learning with {arules} + {arulesViz}
PDF
R勉強会@東京 - Tokyo.Rの紹介
PPTX
Tokyo r49 beginner
PDF
10分で分かるr言語入門ver2.15 15 1010
PDF
前処理のための前処理(Tokyo.R#45)
PDF
mxnetで頑張る深層学習
PDF
リーンローンチパッド特別講義(公開用)
PDF
Rでセイバーメトリクス (イチローと松井 どっちが凄い?)
PDF
5分でわかるかもしれないglmnet
PDF
SVM実践ガイド (A Practical Guide to Support Vector Classification)
PPTX
SVMについて
Trading volume mapping R in recent environment
Visualization of Supervised Learning with {arules} + {arulesViz}
R勉強会@東京 - Tokyo.Rの紹介
Tokyo r49 beginner
10分で分かるr言語入門ver2.15 15 1010
前処理のための前処理(Tokyo.R#45)
mxnetで頑張る深層学習
リーンローンチパッド特別講義(公開用)
Rでセイバーメトリクス (イチローと松井 どっちが凄い?)
5分でわかるかもしれないglmnet
SVM実践ガイド (A Practical Guide to Support Vector Classification)
SVMについて
Ad

Similar to Salmon cycle (20)

PDF
文献紹介:SlowFast Networks for Video Recognition
PDF
Rパッケージ“KFAS”を使った時系列データの解析方法
PDF
Granger因果による 時系列データの因果推定(因果フェス2015)
PPTX
検索評価ツールキットNTCIREVALを用いた様々な情報アクセス技術の評価方法
PDF
みどりぼん3章前半
PDF
はじめてのベイズ推定
PPTX
Hadoopカンファレンス20140707
PDF
PDF
統計をとって高速化する
Scala開発
PDF
統計をとって高速化する
Scala開発 by CyberZ,Inc.
PDF
[DLHacks]PyTorch, PixyzによるGenerative Query Networkの実装
PDF
[DL輪読会]Control as Inferenceと発展
PDF
ITエンジニアのための機械学習理論入門 第5章
PDF
競技プログラミング練習会2015 Normal 第1回
PDF
SGD+α: 確率的勾配降下法の現在と未来
PDF
確率解析計算
PPTX
モデル予見制御に基づくペアトレード戦略
PPTX
Fftw誰得ガイド
PPTX
サーバ異常検知入門
PDF
マルチコアを用いた画像処理
文献紹介:SlowFast Networks for Video Recognition
Rパッケージ“KFAS”を使った時系列データの解析方法
Granger因果による 時系列データの因果推定(因果フェス2015)
検索評価ツールキットNTCIREVALを用いた様々な情報アクセス技術の評価方法
みどりぼん3章前半
はじめてのベイズ推定
Hadoopカンファレンス20140707
統計をとって高速化する
Scala開発
統計をとって高速化する
Scala開発 by CyberZ,Inc.
[DLHacks]PyTorch, PixyzによるGenerative Query Networkの実装
[DL輪読会]Control as Inferenceと発展
ITエンジニアのための機械学習理論入門 第5章
競技プログラミング練習会2015 Normal 第1回
SGD+α: 確率的勾配降下法の現在と未来
確率解析計算
モデル予見制御に基づくペアトレード戦略
Fftw誰得ガイド
サーバ異常検知入門
マルチコアを用いた画像処理
Ad

More from Shota Yasui (14)

PDF
L 05 bandit with causality-公開版
PDF
Contextual package
PDF
PaperFriday: The selective labels problem
PDF
TokyoR 20180421
PDF
何故あなたの機械学習はビジネスを改善出来ないのか?
PDF
木と電話と選挙(causalTree)
PDF
計量経済学と 機械学習の交差点入り口 (公開用)
PDF
Factorization machines with r
PDF
Estimating the effect of advertising with Machine learning
PPTX
Prml nn
PPTX
Xgboost for share
PPTX
重回帰分析で頑張る
PDF
Dynamic panel in tokyo r
PDF
Rで部屋探し For slide share
L 05 bandit with causality-公開版
Contextual package
PaperFriday: The selective labels problem
TokyoR 20180421
何故あなたの機械学習はビジネスを改善出来ないのか?
木と電話と選挙(causalTree)
計量経済学と 機械学習の交差点入り口 (公開用)
Factorization machines with r
Estimating the effect of advertising with Machine learning
Prml nn
Xgboost for share
重回帰分析で頑張る
Dynamic panel in tokyo r
Rで部屋探し For slide share

Salmon cycle