abikoushi

\documentclass{beamer}
\usepackage{tikz}
\setbeamertemplate{navigation symbols}{} 
\usetikzlibrary{positioning}
\makeatletter
\usefonttheme{serif}
\usepackage{xcolor}

\def\mathunderline#1#2{\color{#1}\underline{{\color{black}#2}}\color{black}}

abikoushi

library(BiasedUrn)
library(exact2x2)
library(animation)

#MCMCpack::dnoncenhypergeom(x = NA, cs[1],cs[2],rs[1], 1)
saveGIF({
  for(i in c(0:17,17:0)){
    rs <- c(18,17)
    cs <- c(17,18)
    X <-matrix(c(i,cs[1]-i,rs[1]-i,cs[2]-(rs[1]-i)), nrow=2)

abikoushi

x <- seq(-1,1, length.out=100)
x2 <- 1*(x>0)
X <- cbind(x, x2, x*x2)
B <- c(0.5,2,-1)
mu <- c(X%*%B)
set.seed(2024);y <- mu + rnorm(100,0,0.1)

fit <- lm.fit(X,y)
muhat <- c(X%*%c(fit$coefficients))

abikoushi

library(animation)
pvfun <- Vectorize(function(r, y){
  poisson.test(y, T = tau, r = r)$p.value  
})

postprob <- function(r, x, tau){
  pl <- pgamma(r, shape = x+0.5, rate = tau, lower.tail = FALSE)
  pu <- pgamma(r, shape = x+0.5, rate = tau, lower.tail = TRUE)
  2*pmin(pl,pu)
}

abikoushi

library(gtools)
library(ggplot2)
library(ggbrick)
set.seed(1)
wp <- numeric(10000)
ksp <- numeric(10000)
tp <- numeric(10000)
for(i in 1:10000){
  y1 <- rdirichlet(50,rep(1,10))
  y2 <- rdirichlet(50,rep(10,10))

abikoushi

library(KFAS)
set.seed(123456)
alpha <- beta <- numeric(100)
alpha[1:2] <- rnorm(2)
for(i in 3:100){
  alpha[i] <- rnorm(1,2*alpha[i-1]-alpha[i-2],0.2)
}
beta <- rnorm(100,alpha,0.2)
y1<-rnbinom(100,size=10,mu=exp(alpha))
y2<-rnbinom(100,size=10,mu=exp(beta))

abikoushi

\documentclass[tikz,border=14pt, dvipsnames]{standalone}
\usepackage{tikz}
\usetikzlibrary{matrix,shapes,decorations.pathreplacing, backgrounds, positioning}
\usepackage{amssymb}

\begin{document}
    \begin{tikzpicture}[
        %Styles
        Matrix/.style={
            matrix of nodes,

abikoushi

library(rootSolve)
pvfun0 <- function(r, x, tau){
  pu <- ppois(x-1, r*tau, lower.tail = FALSE)
  pl <- ppois(x, r*tau, lower.tail = TRUE)
  2*pmin(pl, pu)
}
sol0 <-uniroot.all(f = function(p){pvfun0(p,5,10)-0.05},
                   interval = c(0.1,3))
sol0
res_p

abikoushi

library(rootSolve)
pvfun0 <- function(p, x, n){
  pu <- pbinom(x-1, n, p, lower.tail = FALSE)
  pl <- pbinom(x, n, p, lower.tail = TRUE)
  2*pmin(pl, pu)
}
sol0 <-uniroot.all(f = function(p){pvfun0(p,5,10)-0.05},
            interval = c(0.1,0.9))
res_b <- binom.test(5, 10)
print(res_b$conf.int)

abikoushi

library(rootSolve)
LRtest <- function(p,x,n){
  phat <- x/n
  lp <- 2*(dbinom(x, n, phat ,log = TRUE)-
       dbinom(x, n, p ,log = TRUE))
  pchisq(lp, 1, lower.tail = FALSE)
}
set.seed(123)
n <- 20
x <- rbinom(1, n, 0.5)