Skip to content

Instantly share code, notes, and snippets.

@bdemeshev
bdemeshev / gist:6240581
Last active March 21, 2019 18:56
[12 ключей и окрестности] #go #maps
Город Венёв и окрестности
http://a-kainarov.narod.ru/path6_5.htm
http://www.geocaching.su/?pn=101&cid=1088
http://autotravel.ru/town.php/231
моя карта для редактирования
http://openstreetmap.ru/?mapid=1985409576&hash=13990879648f9d15404f1fd6cc40390e
@bdemeshev
bdemeshev / gist:6306569
Last active March 21, 2019 18:48
[My small web] #links
## Доски объявлений
stats.stackexchange.com
math.stackexchange.com
stackoverflow.com
skeptics.stackexchange.com
## Блоги
flavorchemist.livejournal.com
@bdemeshev
bdemeshev / ols_estimation.tex
Last active March 29, 2019 21:34
[new problem by Dima] #metrics
\item Модель линейной регрессии имеет вид $y_i=\b_1 x_{i,1}+\b_2 x_{i,2} + u_i$.
Сумма квадратов остатков имеет вид $Q(\hb_1,\hb_2)=\sum_{i=1}^n (y_1-\hb_1 x_{i,1}-\hb_2 x_{i,2})$.
\begin{enumerate}
\item Выпишите необходимые условия минимума суммы квадратов остатков
\item Найдите матрицу $X'X$ и вектор $X'y$ если матрица $X$ имеет вид
$X=
\left(
\begin{array}{cc}
x_{1,1} & x_{1,2} \\
\vdots & \vdots \\
@bdemeshev
bdemeshev / mv.R
Last active March 21, 2019 18:47
[Простой пример shiny] #R #code #shiny
library("mvtnorm") # для многомерного нормального
library("ggplot2") # для графиков
n_obs <- 100
corr <- 0.5
mu <- c(10, 15)
# ковариационная матрица
# 4 corr * 2 * 3
# ? 9

Keybase proof

I hereby claim:

  • I am bdemeshev on github.
  • I am bdemeshev (https://keybase.io/bdemeshev) on keybase.
  • I have a public key ASBhJFPSIxAJi6ZLauwcooLbCtfUMlG5zLXTTZhe70a2igo

To claim this, I am signing this object:

@bdemeshev
bdemeshev / marriage_analysis.R
Last active March 29, 2019 21:34
[Ежемесячное количество браков] #R #code #metrics #time_series
library(rio) # загрузка данных
library(tidyverse) # манипуляции
library(forecast) # прогнозировани
library(tsibble) # ряды богатой структуры
library(stringr) # строками
library(skimr) # описательные статистики
library(lubridate) # работа с датами
# данные скачиваются с fedstat.ru
# tools - install packages: rio, tidyverse, forecast
@bdemeshev
bdemeshev / polarizatin.md
Last active March 26, 2019 11:44
[Polarization] #physics #quantum

Begin with standard abbreviations for polarized filters:

H = Horizontal,
V = Vertical,
D = Diagonal,
R = Right circular.

HHH means three horizontal filters in a row. Start with totally unpolarized light.

@bdemeshev
bdemeshev / syllabus_maikop.md
Created March 26, 2019 11:43
[Программа Майков] #teaching #metrics

Программа.

Тема 1. Множественная регрессия без статистических предпосылок

Задача оптимизации. Правила работы с матричным дифференциалом. Геометрическая интерпретация. Показатели RSS, ESS, TSS, R^2.

Тема 2. Множественная регрессия с предпосылками на дисперсию

Правила нахождение ковариационных матриц. Теорема Гаусса-Маркова с доказательством.

@bdemeshev
bdemeshev / turtle_fractal.py
Last active March 29, 2019 20:13
[Черепашка рисует фрактал] #code #math
import turtle
import tkinter as tk
scale = 0.05 # choose 1, 0.05
power = 1 # nice one: 1, 0.9, 2, 0.8
root = tk.Tk()
canvas = tk.Canvas(master=root, width=2000, height=1500)
canvas.pack()
@bdemeshev
bdemeshev / russian_csv.R
Created March 27, 2019 16:17
[Импорт-экспорт русских файлов в R] #code #r
Sys.setlocale(category = "LC_CTYPE", locale = "russian")
# try
foo = rio::import("data.txt")
foo = readr::read_csv2("data.txt")
foo = data.table::fread("data.txt", header=TRUE, sep="|", sep2=".")
# check!
# 1. cyrillic letters visible in environment
# 2. cyrillic letters in