Мои функции "Работа с одномерными данными"

#Мои функции "Работа с одномерными данными" архив со скриптами можно загрузить здесь

setwd('C:/R myFunction')

set.seed(55)

x = rnorm(800, 100, 30)

# 1 функции HistTa и HistTm =>

# 2 Tnorm автопроверкa данных на нормальное распределение

# 3 mfBplot ящик с усами с выводом в графическом окне значения выбросов

# 4 mfBplotA ящик с усами с выводом выбросов или их удалением

# 5 DvybrosBP автоматическое удаление выбросов

# 6 mfNAdate проверка данных на пропущенные значения с автозаменой

 

############################################################################

# функции HistTa и HistTm принимают на вход вектор исходных данных и выводят:

# 1. гистограмму с границами нормированных отклонений (1, 2, 3 сигм)

# 2.1 HistTa добавляет на гистограмму все значения больше 3 сигм в виде точек

# 2.2 HistTm добавляет на гистограмму только

#     минимальное и максимальное значения из выборки

# 3. вернет таблицу значений вариант отклоняющиеся больше чем на три сигмы

 

source('HistTa.R')

source('HistTm.R')

HistTm(x)

HistTa(x)

 

#############################################################

# Tnorm Функция автопроверки данных на нормальное распределение

# тест Шапиро-Уилка

# функция Tnorm принимают на вход вектор исходных данных и выводят:

# 1. график QQ теста

# 2. гистограмму

# 3. сделает вывод о нормальности распределения

# 4. выведет значение p-value

 

source('Tnorm.R')

Tnorm(x)

 

#################################################################

# mfBplot строит ящик с усами с выводом в графическом окне

# округл значения выбросов

source('mfBplot.R')

mfBplot(x)

 

###################################################################

# mfBplotA(вектор данных, 0 или 1)

# автоматически удаляет выбросы (0 оствит выбросы, 1 удалит выбросы)

# строит ящик с усами с выводом в графическом окне

# создает новый вектор Х без выбросов

source('mfBplotA.R')

# mfBplotA(x, 1)

# mfBplotA(x, 0)

mfBplotA(x, 1)

# X

#################################################################

# DvybrosBP автоматическое удаление выбросов с построением

# ящиков с усами для сравнения до и после

# выведет значения удаленных выбросов

source('DvybrosBP.R')

DvybrosBP(x)

##############################################################

# mfNAdate

# проверка вектора данных на пропущенные значения с автоматической заменой

# ... на среднее если данные вектора подчинены закону нормального распред

# ... на медиану если не подчинены

# команда source считывает указанный скрипт 'mfNAdate.R' из рабочей директории.

# После ее выполнения в фоновом режиме появится функция

# к кторой можно обращаться прямо из скрипта: mfNAdate(x)

# setwd('C:/R myFunction')

source('mfNAdate.R')

set.seed(45)

d1  <- rnorm(2000)

d2  <- runif(2000)

d3 = as.character(d2)

d1[1:10]  <- NA

d2[1:10]  <- NA

d1 = mfNAdate(d1)

d2 = mfNAdate(d2)

d3 = mfNAdate(d3)

head(d1)

head(d2)

head(d3)

mfNAdate(x)