abikoushi

require("tidyr")
require("dplyr")
require("ggplot2")
require("see")

dfiris <- mutate(iris, id=row_number()) %>% 
  pivot_longer(Sepal.Length:Petal.Width)

#head(dfiris)

abikoushi

library(ggplot2)
library(dplyr)

lr_default <- function(t, delay, forgetting){
  (t+delay)^(-forgetting)
}

linedf <- expand.grid(delay=c(1,2), forgetting=c(0.6,0.9), t=seq(0,20, length.out=50)) %>% 
  mutate(lr=lr_default(t, delay, forgetting)) %>% 
  mutate(param = paste0("delay=",delay, ", ", "forgetting=", forgetting)) %>% 

abikoushi

library(ggplot2)
library(magrittr)
library(tidyr)
library(dplyr)

sigmoid2 <- function(x_from, x_to, y_from, y_to, smooth = 5, n = 100, direction = "x", location = 0, scale = 1) {
  if(!direction %in% c("x", "y")) {stop("Only the directions x or y is allowed.")}
  
  if(direction == "x") {
    x <- seq(-smooth, smooth, length = n)

abikoushi

library(datasauRus)
library(ggplot2)
library(dplyr)

X = dplyr::filter(datasaurus_dozen,dataset=="dino") %>% 
  dplyr::select(-dataset) %>% 
  as.matrix()
res_svd <- svd(as.matrix(X))

df_svd <- bind_rows(data.frame(res_svd$u, matrix="U"),

abikoushi

library(datasauRus)
library(ggplot2)
library(dplyr)
library(gganimate)

rot <-function(ti, period=16) {
  coeff <- 2*pi/period
  matrix(c(cos(coeff*ti),-sin(coeff*ti),
           sin(coeff*ti), cos(coeff*ti)), byrow = TRUE, nrow = 2)  
}

abikoushi

#include "RcppArmadillo.h"
// [[Rcpp::depends(RcppArmadillo)]]
#include <iostream>
#include <string>
#include <vector>
using namespace Rcpp;


void readmtx(arma::uvec & row_i,
             arma::uvec & col_i,

abikoushi

#include "RcppArmadillo.h"
// [[Rcpp::depends(RcppArmadillo)]]
#include <iostream>
#include <string>
#include <vector>
using namespace Rcpp;

// [[Rcpp::export]]
std::vector<double> readaline(std::string readtxt, int ind_n) { 
  std::ifstream file(readtxt);

abikoushi

query_awk <- function(vpar, pdir, filename, uid){
  chl <- paste0("cat ", pdir, "/", filename)
  for(i in 1:nrow(vpar)){
    chr <- paste0("\"", paste(vpar[i,], collapse = ", "), "\"")
    awk <- paste0(" | awk 'NR==", uid[i], "{sub('/.*/',", chr, ")}1'")
    chl <- paste0(chl, awk)
  }
  chl <- paste0(chl, " > ", pdir, "/tmp.csv")
  mvit <- paste0("mv ", pdir, "/tmp.csv " , pdir, "/", filename)
  system(chl)

abikoushi

x = seq(-4,4,length.out=100)
y = seq(-4,4,length.out=100)
df <- expand.grid(x = x, y = y)
prob <- matrix(pnorm(df$x)*pnorm(df$y), 100, 100)

#png("persp.png")
persp(x = x, y = y, z = prob, theta = 20, phi=45)
#dev.off()

#png("contour.png")

abikoushi

library(datasauRus)
library(ggplot2)
library(dplyr)
library(gganimate)

rot <-function(th) {
  matrix(c(cos(th),-sin(th),
           sin(th),cos(th)), byrow = TRUE, nrow = 2)  
}