注目の投稿

【kepler.gl】コロナ対策による人流の変化も地図上に可視化(各種メディアで報道)

kepler.glのサイト画面 kepler.glを使ってコロナ対策の効果を分析したところ、テレビ、新聞、ネットのメディアから問い合わせや報道依頼が殺到。今も、土日返上で都内や全国の人流変化を分析しています。この記事では人流変化の可視化に便利なkepler.glにつ...

2016年8月3日水曜日

【R】回帰分析などの結果を個別に出力したい

#下記のような回帰分析をおこなったとき
result <- lm(目的変数カラム ~ 説明変数カラム,data = df.sample)

#下記のような出力を得られる
> summary(result)

Call:
lm(formula = session_num ~ search_num + compe_search_num + publish_period, 
    data = df.kiji_select)

Residuals:
   Min     1Q Median     3Q    Max 
-12574  -1323   -772   -345 365681 

Coefficients:
                   Estimate Std. Error t value Pr(>|t|)    
(Intercept)      495.223086  26.609894   18.61   <2e-16 ***
search_num        -0.037315   0.003511  -10.63   <2e-16 ***
compe_search_num   0.089878   0.004548   19.76   <2e-16 ***
publish_period     2.527575   0.058255   43.39   <2e-16 ***
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

Residual standard error: 4595 on 90724 degrees of freedom
Multiple R-squared:  0.02604, Adjusted R-squared:  0.02601 
F-statistic: 808.6 on 3 and 90724 DF,  p-value: < 2.2e-16

#Estimateのみがほしいとき
> summary(result)$coefficients[,"Estimate"]
                 summary.result..coefficients....Estimate..
(Intercept)                                    495.22308639
search_num                                      -0.03731478
compe_search_num                                 0.08987823
publish_period                                   2.52757539

#Pr(>|t|)のみがほしいとき
> summary(result)$coefficients[,"Pr(>|t|)"]
                 summary.result..coefficients....Pr...t....
(Intercept)                                    3.680971e-77
search_num                                     2.289391e-26
compe_search_num                               9.528811e-87
publish_period                                 0.000000e+00


◇参考サイト:Rから回帰分析のいろいろな情報を取得する方法
http://am-yu.net/2013/12/03/r-linear-regression/

2016年7月29日金曜日

【衝撃】町内会の草むしりでパクチーの種を植えていた人の話

*データ分析とは全く関係ないです。

すごい変わった人が社内にいます。

先日、ひょんなことから一緒に食事に行くことになりました。

---食事中---

私「なんかいつも変わったことしてるよね」

Aさん「そうですかwww」

私「ほかでも変なことしてるの?

Aさん「そうですねー。この前、町内会の草むしりにいきました」

私「偉いね!」

Aさん「しかも、誰よりもがんばっていたので褒められました」

私「なんで、そんなに草むしりがんばったの?」

Aさん「パクチーの種をまきたかったんで」

私「えっ!!! どゆこと!?」

Aさん「自分、パクチー好きだから育てて食べたかったんです」

私「。。。」

---その後の話---

Aさんはパクチーの話をフェイスブックに投稿したらしい

それを見た友人が面白がり拡散

その友人の拡散した記事を見た友人の友人が興味を持つ

そのAさんの友人の友人はAさんの友人を通じてAさんにコンタクトをとる

結果、AさんとAさんの友人の友人はデートへ行くことに

現在、2回目のデートの予定を立てているとか。。。

◇感想

事実は小説よりも奇なり

2016年7月13日水曜日

【統計】統計手法の選び方

これ良い!

2016年7月9日土曜日

社長室ってなにやってるの?

入社とほぼ同時にできた社長室に配属されることになりました

周りの人から
「社長室ってなにやっているの?」ってよく聞かれます

当然ですね。営業部なら営業。システム部ならシステム。
では社長室なら??? 何するところか全く不明です。


いろんな人「社長室ってなにやっているの?」 

入社3日の私「私もよく分かりません」
入社10日の私「うーん、いろいろです」
入社15日の私「経営企画、事業企画、営業企画に関連することをします(たぶん)」
now 入社20日の私「私は主に経営企画、事業企画、データ戦略を担当しています」


最初はぼんやりしていましたが、だんだんと役割が明確になってきた感じです

とは言え、

「経営企画は何するの?」「事業企画って何?」「データ戦略???」

という新たな質問をされるようになりました


Aさん「経営企画って具体的に何するの?」

私「えっと、経営企画関連では…」


次回につづく… 



かも


2016年7月7日木曜日

R 空白が入っている行を削除

これでOK!

df.test[df.test$空白か調べたいカラム == "",]$空白か調べたいからむ <- NA
df.test <- na.omit(df.test)

簡単!