臨床統計もおもしろいですよ、その2

[過去ﾛｸﾞ] 臨床統計もおもしろいですよ、その2 (1002ﾚｽ)
上下前次 1-新

このｽﾚｯﾄﾞは過去ﾛｸﾞ倉庫に格納されています｡
次ｽﾚ検索歴削→次ｽﾚ栞削→次ｽﾚ過去ﾛｸﾞﾒﾆｭｰ

504: 2018/12/16(日)10:16 ID:1U2yFJp0(1) AAS
>>501
おい、上の過剰ゼロポアソンモデルから放校者数の期待値を計算してみ。
来年のド底辺シリツ医大の放校者の予想になるぞ。
まあ、ド底辺シリツが放校者ゼロ過剰というのは過大評価だがね。

505: 2018/12/16(日)10:39 ID:/hMzfoFS(4/9) AAS
pairsxyz <- function(x,y,z){
xyz=data.frame(x,y,z)
panel.cor <- function(x, y, ...){
usr <- par('usr'); on.exit(par(usr))
par(usr=c(0,1,0,1))
r <- cor(x, y, method='spearman', use='pairwise.complete.obs')
zcol <- lattice::level.colors(r, at=seq(-1, 1, length=81), col.regions=colorRampPalette(c(scales::muted('red'),'white',scales::muted('blue')), space='rgb')(81))
ell <- ellipse::ellipse(r, level=0.95, type='l', npoints=50, scale=c(.2, .2), centre=c(.5, .5))
polygon(ell, col=zcol, border=zcol, ...)
text(x=.5, y=.5, lab=100*round(r, 2), cex=2, col='black')
省22

506: 2018/12/16(日)10:40 ID:/hMzfoFS(5/9) AAS
pairs(xyz,
diag.panel=panel.hist,
lower.panel=panel.scatter,
upper.panel=panel.cor,
gap=0.5,
labels=gsub('\\.', '\n', colnames(xyz)),
label.pos=0.7,
cex.labels=1.4)
}

Nt <- 100
省5

507: 2018/12/16(日)11:38 ID:ZC5C+Eow(2/2) AAS
ニートニート働けニート

動画ﾘﾝｸ[YouTube]

508: 2018/12/16(日)11:42 ID:/hMzfoFS(6/9) AAS
Nt <- 100
set.seed(123)
Ti <- sample(20:40, Nt, replace=TRUE)
Xt <- rlnorm(Nt, 100, 15)
Yt <- rbinom(Nt, size=Ti, prob=0.8)

20-40人からなるド底辺シリツ学生のグループを１００組集める。
各グループを代表する知能指数が正規分布をしている。
裏口の確率は80%である。

その線形回帰は

glm(Yt ~ Xt, family=poisson(link="log"), offset=log(Ti))

509: 2018/12/16(日)14:16 ID:/hMzfoFS(7/9) AAS
source('tmp.tools.R')

library(Mass)
set.seed(71)
size <- 10000
x <- -seq(1/size,1,by=1/size)
shape <- 10.7
mu <- exp(3.7*x+2)
y1 <- sapply(mu, function(m)rnegbin(1, mu=m, theta=shape))
pairsxyz(cbind(x,y1,mu))
model_glmnb <- glm.nb(y1 ~ x)
省10

510: 2018/12/16(日)20:30 ID:/hMzfoFS(8/9) AAS
# ネットワークノード1..j..N
# ネットワーク内の客の数1..k..K、
# ノードjでの平均客数L[k,j]、通過時間W[k,j]、到着率lambda[k,j]
# mu[j]:サービス率,r[i,j]ノードiからjへの移行確率
# theta[j] = Σtheta[i](i=1,N) * r[i,j] の解がtheta
rm(list=ls())
theta=c(0.6,0.4)
mu=c(0.6,0.4)
N=2
K=10
省17

511: 2018/12/16(日)20:30 ID:/hMzfoFS(9/9) AAS
#
for(k in 2:K){
for(j in 1:N){
W[k,j]=1/mu[j]*(1 + L[k-1,j])
}
tmp=numeric(N)
for(j in 1:N){
tmp[j]=theta[j]*W[k,j]
}
den=sum(tmp)
省11

512: 2018/12/17(月)04:38 ID:GOYrgwsC(1/6) AAS
minimize 2*b+4*sin(θ)/(cos(θ)-1/2)*α where (θ-sin(θ)/2+(sin(θ)/(cos(θ)-1/2))^2(α-sin(α))/2+(1-b)*sin(θ)/2=pi/8 ,0<b<1,0<θ<pi/2

513: 2018/12/17(月)04:45 ID:GOYrgwsC(2/6) AAS
minimize (2*b+4rα) where (θ-sinθ)/2+r^2(α-sinα)/2+(1-b)sinθ/2=pi/8 ,0<b<1,0<θ<pi/2

外部ﾘﾝｸ:www.wolframalpha.com

Standard computation time exceeded...

Try again with Pro computation time

514: 2018/12/17(月)05:14 ID:GOYrgwsC(3/6) AAS
rangeB = map (/1000) [1..1000]
rangeTheta = map (\x -> x * pi/2/1000) [1..1000]
rangeAlpha = map (\x -> x * pi/1000) [1..1000]
rangeR = map (\x -> x * 1/1000) [1000..10000]
re = [2*b+4*r*α| b<-rangeB,θ<-rangeTheta, α<-rangeAlpha,r<-rangeR,(θ-sin(θ))/2+r^2*(α-sin(α))/2+(1-b)*sin(θ)/2==pi/8]
main = do
print $ minimum re

515: 2018/12/17(月)05:39 ID:GOYrgwsC(4/6) AAS
rangeB = map (/1000) [1..1000]
rangeTheta = map (\x -> x * pi/2/1000) [1..1000]
rangeAlpha = map (\x -> x * pi/1000) [1..1000]
re = [(2*b+4*sin(θ)/(cos(θ)-1/2)*α)| b<-rangeB,θ<-rangeTheta, α<-rangeAlpha,(θ-sin(θ))/2+(sin(θ)/(cos(θ)-1/2))^2*(α-sin(α))/2+(1-b)*sin(θ)/2==pi/8]
main = do
print $ minimum re

516: 2018/12/17(月)05:47 ID:GOYrgwsC(5/6) AAS
外部ﾘﾝｸ:ja.wolframalpha.com

517: 2018/12/17(月)05:47 ID:GOYrgwsC(6/6) AAS
minimize(2*b+4rα) where (θ-sinθ)/2+r^2(α-sinα)/2+(1-b)sinθ/2=pi/8 ,0<b<1,0<θ<pi/2,0<α<2*pi

518(4): 2018/12/17(月)13:52 ID:FuhGkKV8(1/5) AAS
向学心のないコミュ障事務員のスレ

事務員と見せかけてその正体は大学受験に失敗したニートだろうな
二十歳前くらいの大学受験の知識で一生わたっていけると勘違いしているらしい

いや勘違いというか自分にそう言い聞かせなければいけないほど
人生終わってんだろうな
ハゲワロス

519: 2018/12/17(月)13:59 ID:N+SAyf9z(1/6) AAS
>>518
こういう野糞投稿しかできないのはド底辺シリツかな？
ド底辺スレから部屋割計算できなくて逃げ出したアホと同一人物？

520: 2018/12/17(月)14:01 ID:N+SAyf9z(2/6) AAS
>>518
向学心を語るなら>499の計算してからな。
ド底辺頭脳クン。

521(1): 2018/12/17(月)14:02 ID:FuhGkKV8(2/5) AAS
人生についての向学心は大学受験なんていう狭い領域で机上の空論をこねくり回すことではない
だがいつまでも不合格の受験生は、その領域から抜け出せなくなるんだなW

向学心のない事務員にふさわしいスレを持ってきてやったぞ

【因果応報】医学部再受験【自業自得】
2chｽﾚ:kouri

522: 2018/12/17(月)14:07 ID:N+SAyf9z(3/6) AAS
>>521
>499は受験に無関係な臨床問題だよ。
答えてみ！
はよ！はよ！
向学心を語るなら>499の計算してからな。
ド底辺頭脳クン。

523(1): 2018/12/17(月)14:08 ID:FuhGkKV8(3/5) AAS
このスレもお勧めだぞW

【悲報】ワイの医学部志望の妹、精神壊れる
外部ﾘﾝｸ[html]:vippers.jp

上下前次 1-新書関写板覧索設栞歴

あと 479 ﾚｽあります
ｽﾚ情報赤ﾚｽ抽出画像ﾚｽ抽出歴の未読ｽﾚ AAｻﾑﾈｲﾙ

ぬこの手ぬこTOP 0.019s