TOP カテ一覧 スレ一覧 100〜終まで 2ch元 削除依頼
遂に解かれた!AX+BY=CZ
数学板ID表示制検討スレッド
Inter-universal geometry と ABC 予想 45
5次方程式の解を表現できる数体系
現代数学の系譜 古典ガロア理論を読む35
フェルマーの最終定理の簡単な証明7
wwwww小説を読むスレwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwww
分布 distribution 超関数 hyperfunction一般化関数
1÷0ってなんで答えが定義されてないの?
哲学を理解できない馬鹿が、数学や物理に逃げる
102 :
Rとstanでベイズ統計ができるなら、以下のコードで実行可能。
パッケージ rstan と BEST(信頼区間グラフ描出用)が必要


library("BEST")
library("rstan")
rstan_options(auto_write = TRUE)
options(mc.cores = parallel::detectCores())
options(scipen = 5)

model.string='
data{
int n; // sample size
int x; // number of positive test
real<lower=0,upper=1> sen; // sensitivity 0.7
real<lower=0,upper=1> spc; // specificity 0.9
real<lower=0,upper=1> ul; // uniform(0,ul)
}

parameters{
real<lower=0,upper=1> prev; // prevalence
}

transformed parameters{
real<lower=0,upper=1> p;
p = prev*sen + (1-prev)*(1-spc) ; // probability of positive test result
}

model{
x ~ binomial(n,p);
prev ~ uniform(0,ul);
}

'
writeLines( model.string , con='model.stan' )
corona1.model=stan_model('model.stan')
# saveRDS(corona1.model,file='corona1.rds')
# corona1.model=readRDS('corona1.rds')

PCRs <- function(N=1000,X=10,UL=1,SEN=0.7,SPC=0.9,verbose=FALSE,...){
data = list(n=N,x=X,sen=SEN,spc=SPC,ul=UL)
fit.corona = sampling(corona1.model, data=data,
seed=1234,control=list(adapt_delta=0.99),...)
if(verbose) print(fit.corona, prob=c(0.025,0.5,0.975),pars=c('prev'),digits=8)
ms=rstan::extract(fit.corona)
BEST::plotPost(ms$prev,showMode = T,xlab='prevalence') ; lines(density(ms$prev),col='skyblue')
c(mean=mean(ms$prev),HDInterval::hdi(ms$prev)[1:2])
}

整数論を勉強するためのスレッド
現代数学の系譜 工学物理雑談 古典ガロア理論も読む56
中1の弟が数学好きなようだから先取りさせたい
巨大数探索スレッド13
数学科がどんな恐ろしいところか教えてくれ
偶数である素数はただひとつであることを証明せよ
数学用語に「夜の」って入れるとエロくなる
\って何なんだよ……
現代数学の系譜 工学物理雑談 古典ガロア理論も読む40
なぜ∞という数が存在すると思ってる学生が多いのか
--------------------
京都のおいしい拉麺屋は?22杯目
【経済】なぜ「世界に冠たる企業」は日本から消滅したか★2
稲場まなかんヲタ=女慣れしてないキモいおっさんという風潮
木多康昭アンチスレ木多康昭に煉獄2手目
▼△ 煙草銘柄・器具関係 IDおみくじ ▼△
結局、パヨクの正体ってこれだよね?
【はいふり】劇場版 ハイスクール・フリート Part10
この板の存在意義
週刊少年サンデー総合スレッド346冊目【北退出入禁止】
【スズキ】新型ジムニー☆23【JB64】
ほんじょ・はるかハマスタ売り子スレ袖まくり子
【免許皆伝】ビッグ画大先生のお得情報交換コーナー
東京の食いもんはマズい!マズすぎる!!283
精神病んでる人 in 創作文芸 [2人目]
サラサラ◆小麦粉シャンプー◆ふわふわ
【地上波】テラスハウス TERRACE HOUSE 27
同人やめる時、やめたい時 2
宇宙戦艦ヤマトTV版パート1限定スレ 【8】動砲
千葉県市川の【大慶園】について
【MHW】発売からずっとソロでやってる人いる?【モンハン】
TOP カテ一覧 スレ一覧 100〜終まで 2ch元 削除依頼