底辺私立医大を卒業した医者って頭悪いよね？ Part17

底辺私立医大を卒業した医者って頭悪いよね？ Part17at HOSP

底辺私立医大を卒業した医者って頭悪いよね？ Part17 - 暇つぶし2ch■コピペモード
□スレを通常表示
□オプションモード
□このｽﾚｯﾄﾞのURL
■項目テキスト

1000:卵の名無しさん
 20/12/01 08:54:10.45 iy4usOEg.net
library(gmp) 
f <- function(N){ 
a=factorialZ(N) 
b=as.character(a) 
c=as.numeric(unlist(strsplit(b,''))) 
d=rev(c) 
if(d[1]!=0) return(NA) 
else return(rle(d)$values[2]) 
} 
N=1:10004 
y=sapply(N,f) 
plot(N,y,bty='n') 
z=data.frame(N,digit=y) 
head(z,10) 
tail(z,10) 
a=z[-1:-4,] # non-zero digit 
b=vector('list',length=9) 
for(i in c(2,4,6,8)){ 
b[[i]]=a[a$digit==i,]$N 
} 
b 
c=sapply(c(2,4,6,8), function(i) length(b[[i]])) 
data.frame(digit=c(2,4,6,8),個数=c)

1001:卵の名無しさん
 20/12/01 10:28:29.05 VGK+tvyz.net
# 差X-Yのpdfの公式 
# ∫[-∞,∞] pdf1(x+y)*pdf2(y) dy 
# x ~ pdf1,y ~ pdf2から x-y のpdfを返す 
pdf_minus <- function(pdf1,pdf2){ 
f <- function(x,y) pdf1(x+y)*pdf2(y) 
f=Vectorize(f,vectorize.args = 'y') 
pdf <- function(x) 
integrate(function(y) f(x,y),-Inf,Inf,rel.tol=1e-14)$value 
pdf=Vectorize(pdf) 
invisible(pdf) 
} 
# ex. COVID19 incubation time 
mu = 1.434065 ; sg = 0.6612 
x=rlnorm(1e6,mu,sg) ; y=rlnorm(1e6,mu,sg) 
hist(x-y,xlim=c(-30,30),freq=F,breaks=250,lty=1,col='skyblue',main='') 
pdf_minus(function(x)dlnorm(x,mu,sg), 
function(x)dlnorm(x,mu,sg)) -> pdf 
curve(pdf(x),add=T) 
var=integrate(function(x) x^2*pdf(x),-Inf,Inf)$value 
sqrt(var) 
integrate(pdf,2,Inf,rel.tol=1e-12) 
mean(x-y>2) 
library(polspline) 
fit=logspline(x-y) 
curve(dlogspline(x,fit),add=T,col=4) 
1-plogspline(2,fit)

次ページ