abikoushi

library(dplyr)
library(tidyr)
library(ggplot2)
library(gganimate)

logsumexp2 =function (logx1,logx2){
  logx1 + log1p(exp(logx2-logx1))
}

llmixnorm <- function(par, y){

abikoushi

library(Rcpp)
library(RcppParallel)
sourceCpp("R/cpp/parallelreduce.cpp")
# defaultNumThreads()
setThreadOptions(numThreads =  4)
m <- rnorm(1e+8)

# ensure that serial and parallel versions give the same result
# compare performance of serial and parallel
library(rbenchmark)

abikoushi

library(roll)
library(RcppParallel)

RcppParallel::defaultNumThreads()
x <- rnorm(1e8)

setThreadOptions(numThreads = 1L)
t1 <- system.time({
  roll_mean(x, 1000)
})

abikoushi

##参考文献：
##スティーブン・ピンカー『人はどこまで合理的か　上』,5章，訳：橘明美，草思社

#devtools::install_github("abikoushi/ggsomestat")
library(ggsomestat)
library(dplyr)

taxi <- matrix(c(0.68, 0.03,
                 0.17, 0.12), byrow = TRUE, nrow = 2)
gb=c("green","blue")

abikoushi

library(moltenNMF)
library(Matrix)

df = expand.grid(letters[1:3], LETTERS[1:2])
X = sparse_onehot(~., data = df)

A = matrix(c(1,0,0,0,
             0,1,0,0,
             0,0,1,0,
             0,0,0,0,

abikoushi

\documentclass[tikz,border=10pt]{standalone}

\usepackage{tikz}
\usetikzlibrary{positioning, arrows.meta}

\begin{document}
\begin{tikzpicture}[
  node distance=30mm,
  box/.style={draw, rounded corners, align=center, inner sep=6pt, minimum height=5em ,minimum width=7em},
  arrow/.style={->, thick}

abikoushi

## reparameterized Poisson NMF
NMF_re <- function(X,
                   M=2,
                   a_W=1,
                   a_H=1,
                   alpha_b=1, beta_b=1,
                   tol=1e-3, maxit=100){
  D <- nrow(X)
  N <- ncol(X)
  b_W <- rgamma(1, 1, 1)

abikoushi

NMF <- function(X, M=2,a_W=1,b_W=1,a_H=1,b_H=1,tol=1e-2,maxit=10000,
                seed=1234){
  set.seed(seed)
  D <-nrow(X)
  N <-ncol(X)
  W<-matrix(rgamma(D*M,a_W,b_W),D,M)
  H<-matrix(rgamma(M*N,a_H,b_H),M,N)
  
  logW0 <-log(W)
  logH0 <-log(H)

abikoushi

rcfnorm <- function(n, shape=0, inv_scale=1){
  rnorm(n, shape/inv_scale, 1/sqrt(inv_scale))
}

dcfnorm <- function(n, shape=0, inv_scale=1, ...){
  dnorm(n, shape/inv_scale, 1/sqrt(inv_scale), ...)
}

pcfnorm <- function(n, shape=0, inv_scale=1, ...){
  pnorm(n, shape/inv_scale, 1/sqrt(inv_scale), ...)

abikoushi

size_mtx <- function(file_path){
  con <- file(file_path, open = "r") #Open for reading in text mode
  size <- scan(con, what=integer(), comment.char = "%", nlines = 1, skip = 1, quiet = TRUE)
  close(con)
  names(size) <- c("row","column","nonzero")
  return(size)
}

take_last <- function(x){
  x[length(x)]