a*bの値が与えられたときにaとbをベイズで推定してみる
基礎代謝から体重の推移のモデリングをしていたときに未知数が多く、どのくらい識別できるのか不安になったので、実験として。 データは単純にa,bを正規分布から生成 library(rstan) N=...
0
0
Comment0
18 search resultsShowing 1~18 results
Stocked
You need to log-in
stanで確率的主成分分析
多変量時系列モデルの中で使う可能性があり、調べてコードがどこにもなかった(bingも）ので。 参考はビショップ本の下巻の終わりの方など。 YをG次元の確率変数（データ）、xを圧縮したい次元（例で...
0
1
Comment0
TF-IDF特徴量
某案件で使ったので、まとめておきます。 文書群Dにおける単語tjの出てくる確率を pj=(tjが出てくる文書数）/文書数N とみなすと情報量は -log(pj) (a) これは全てのDに...
0
1
Comment0
spacyで指定した単語から係り受けのリンクをたどる
言語処理ノックに似たようなのがあって,spacｙ番は探してもなかったので。 単語の登録に苦労したけど、そちらは必要な調べてください。 import spacy from spacy.langua...
1
0
Comment1
状態空間モデルで死刑の抑止効果を測定してみる
注）ちょくちょくアップデートしています 知人とskypeしている最中にふざけて死刑！死刑！とか言っていたら気になったので。 日本は死刑が廃止された事がないらしく、推定は極めて難しい。過去の研究は...
1
0
Comment2
アウトカムが割り付けに影響する場合の因果推論
傾向スコアによるIPW推定量などの前提は p(z,y|x)=p(z|x) p(y|x) と共変量を条件づけたときに割り付けとアウトカムが独立になる事だが、現実的に完璧に満たされる事も少ない気がす...
0
0
Comment0
スシローのセットの要素の価格を推定
のような情報から、どれがお得か判定したい。 手動で正規化するのは辛かったので、注文ページから文字列をnetaにコピペしてpythonでスプリットする方針で。 import pandas as p...
2
3
Comment0
トピックモデル（LDA）で福祉ビジネスのSNS解析
福祉ビジネスを使わないようにアカウント https://twitter.com/evidencefukusi でbotを動かしているが https://github.com/biones/mis...
0
2
Comment0
syncroomのユーザー数をモデリング
apiから、googleシートに1h毎のデータを取得。 １日の中の変動が特徴的で、感染症のSIRモデルっぽく season[t]=season[t-1]ss+s_ysy[t-1] とかできないか...
0
1
Comment0
寿司のメニューから単品の価格とコスパを推定
近場のテイクアウトを頼むときに、どれがコスパが良いのか考えたときに面白そうだと思ってやってみました。 画像のメニューから目視で https://docs.google.com/spreadshe...
1
1
Comment0
ダブルスだけの結果から個人のレーティングとペアの相性を推定する(stan)
前回の続きで、卓球などで試合はダブルスばかりだったり、シングルと統合したい場合に使えると思う。 核はペアの場合のrateを rate[x]+rate[y]+rate_two[x,y] と分けて、...
2
1
Comment0
２種類の広告介入の効果を分解してみる
１つだとつまらないので、2つのtrendで １つはt=30から(trend)、もう一つは80から(trend2)。 介入時の情報は広告の目的ならわかるので、trendの一回微分（速度）に対してe...
0
1
Comment0
少数サンプルのバイアスで回帰係数は変わるのか？
ネット調査とリアルで１０差があるという仮定で線形回帰してみた。 bias=c() for(i in 1:500){ Nreal=i Nnet=1000-i d1=rnorm(Nreal,mean...
0
1
Comment0
ベイズでイロレーティング計算
逐次的な更新式より簡単かつ拡張性がある気がするので。 1~bernouli(sigmoid(b*delta_rate)) と、trueだけを取るlogistic回帰で、bはrate差100を64...
0
1
Comment0
stanで2次元carモデル（の練習）
と言っても、berobero先生の記事と殆ど同じですが。。。 以前の、時系列＆販売個数に勾配のあるアイテムの推定に再チャレンジするもなかなかうまく行かず、参考になるかなーと思って。 # データの...
0
2
Comment0
特徴量作成(UCLよりbalanceデータ)
bd=read.csv("https://archive.ics.uci.edu/ml/machine-learning-databases/balance-scale/balance-scal...
1
2
Comment0
店、アイテムごとの時系列データをstanで推定する（続き）
データの設定も含めてまだまだ適当です。 (time,store,item)の３次元のCARモデル？と大体同じです。 MAX_T=100 sbase=c(0,50,100,155)+100 #店ご...
1
2
Comment0
stan練習
Airpassengersのデータで季節調整。詳細は北川時系列p165あたり。 scodeap=" data{ int N; real y[N]; } parameters{ real tren...
1
1
Comment0
18 search resultsShowing 1~18 results
Qiita is a knowledge sharing service for engineers.
- You can follow users and tags
- You can stock useful information
- You can make edit suggestions for articles