@@ -1,6 +1,5 @@
-- Язык R
+= Язык R
 
-----
 ''''''ToDo'''''' статьи:
 	* привести в нормальное состояние примеры --АтрашкевичАндрей
 	* написать про data.frame --АтрашкевичАндрей
@@ -14,62 +13,61 @@
 
 http://www.r-project.org/Rlogo.jpg
 
--- Общее описание языка R
+- Общее описание языка R
 
 R — простой и полезный язык для статистики, машинного обучения и DataMining, распространяемый по лицензии GNU GPL 2.
 
 Главная страница проекта: http://www.r-project.org/
 
-Скачать дистрибутивы R можно здесь (российское зеркало): http://cran.gis-lab.info/ (доступны версии для UNIX (Linux), Windows, OS X). Для пользователей Ubuntu и Debian скачивать ничего не нужно: можно установить пакеты r-base, r-base-dev, r-base-cor:
+Скачать бинарники R можно [[http://cran.gis-lab.info/|здесь]] (российское зеркало, доступны версии для UNIX (Linux), Windows, OS X). Пользователям Ubuntu и Debian скачивать ничего не нужно: можно установить пакеты r-base, r-base-dev, r-base-cor:
  $ sudo apt-get update && apt-get install r-base r-base-dev r-base-core
 
-Более подробно об установке R можно прочитать тут http://cran.r-project.org/doc/manuals/r-release/R-admin.html 
+Более подробно об установке R можно прочитать [[http://cran.r-project.org/doc/manuals/r-release/R-admin.html|тут]].
 
--- История создания
+- История создания
 
 В 1976 году сотрудники AT&T Bell Labs: Джон Чамберс, Рик Беккер и Алан Вилкс — создали язык программирования S (от англ. «statistics»).
 
-В 1993 году в стране, где на десять жителей приходится двадцать овец, три хоббита и два эльфа, двумя сотрудниками Оклендского университета Россом Ихака и Робертом Джентльменом (это не шутка-) на основе языка S был разработан язык R. Язык, как и его предшественник, был разработан преимущественно для статистических целей, а его название было подражанием названию S, кроме того, это была первая буква имен разработчиков (несмотря на то, что полное имя Росса Ихака — Джордж Росс Ихака). Разработчик языка S Джон Чамберс входит в R Core Team.
-
--- Соглашения оформления кода
-
-За основу соглашений оформления кода был принят «Google's R Style Guide»: https://google-styleguide.googlecode.com/svn/trunk/Rguide.xml
+В 1993 году в стране, где на десять жителей приходится двадцать овец, три хоббита и два эльфа, двумя сотрудниками Оклендского университета Россом Ихака и Робертом Джентльменом (это не шутка-) на основе языка S был разработан язык R. Язык, как и его предшественник, был разработан преимущественно для статистических целей, а его название было, во-первых, подражанием названию языка S, во-вторых, это была первая буква имен разработчиков (несмотря на то, что полное имя Росса Ихака — Джордж Росс Ихака). Разработчик языка S Джон Чамберс входит в R Core Team.
 
-Примеры кода ниже соответствуют этим соглашениям. 
+- Соглашения оформления кода
 
-----
-	* К: использовал «венгерскую нотацию» по чисто педагогическим причинам (и далеко не везде). Незнакомым с программированием (экономистам, биологам и пр.) гораздо проще понимать на начальном уровне (сужу по себе, кстати) --АтрашкевичАндрей
-
-	* К: А может ну его, хлам этот с точками? Запутывает же. Попробую переписать в ''''''lowCamelCase'''''', так должно быть понятнее. --АтрашкевичАндрей
-----
+	1 Оператор присвоения <- вместо = (соблюдение этого соглашения показывает, что вы super high-level extremely professional advanced R-programmer, кроме того, это придаёт коду изящную «эрность»-).
+	1 Имена переменных записываются в ''''''lowerCamelCase'''''', например: ''''''monthlyAnnuity'''''',  ''''''numberOfIterations''''''.
+	1 Имена функция в ''''''UpperCamelCase'''''': ''''''CalculateAnnuity'''''', ''''''GetNextValue'''''', ''''''SubstractUniqueSubstrings''''''.
+	1 В некоторых случаях, когда это целесообразно, используется венгерская нотация: ''''''vPayments'''''' — вектор, ''''''dfPayments'''''' — таблица (data.frame).
+	1 etc. — добавить
 
--- Программа «Hello, world!»
+- Программа «Hello, world!»
 
 Для присвоения переменной значения в R используется оператор <-. Оператор = также поддерживается, но не соответствует соглашениям. Поэтому во всех дальнейших примерах использовует только оператор <-.
 
- '''> str.hw <- "Hello, world!"'''
- '''> + print(str.hw)'''
+Чтобы в консоли осуществить перевод каретки, необходимо нажать на Shift + Enter. Подобная операция поставит в начале следующей строки знак +.
+
+ '''> hw <- "Hello, world!"'''
+ '''> + print(hw)'''
  [1] "Hello, world!"
 
 R чувствителен к регистру, поэтому будьте внимательны и соблюдайте соглашения оформления кода.
 
- '''str.hw <- "Hello, world!"'''
- '''print(str.HW)'''
- Ошибка в print(str.HW) : объект 'str.HW' не найден
+ '''> hw <- "Hello, world!"'''
+ '''> + print(HW)'''
+ Ошибка в print(HW) : объект 'HW' не найден
 
 Оператор присвоения работает в обе стороны, поэтому можно писать как <-, так и ->.
 
- '''# в R нет многострочных комментариев'''
- '''"Hello, world!" -> str.hw'''
- '''str.hw  # так удобнее всего, наверное, удобнее всего выводить на печать'''
+ '''> # в R нет многострочных комментариев'''
+ '''> + "Hello, world!" -> hw'''
+ '''> + hw  # такой способ работает только в интерактивном режиме'''
  [1] "Hello, world!"
 
-Иногда удобно вывести на печать переменную сразу после присвоения значения. Для этого используется заключение операции присвоения в круглые скобки.
- '''(str.hw <- "Hello, world!")'''
+Иногда удобно вывести на печать переменную сразу после присвоения значения. Для этого используется заключение операции присвоения в круглые скобки. Подобная конструкция возможна также только в интерактивном режиме.
+ '''> (hw <- "Hello, world!")'''
  [1] "Hello, world!"
 
+R — скриптовый язык, поэтому остальные примеры будут даваться как запись скриптов.
 
--- R как калькулятор
+- R как калькулятор
 
 Все имеющиеся т.н. «статистические» пакеты можно условно разделить на три категории:
 	* библиотеки высокоуровневых ЯП (Seismic UNIX в Си, Pandas в Python) и специализированные языки для математических (включая статистические) вычисления: GNU Octave, Julia, MATLAB, MATCAD
@@ -80,140 +78,272 @@
 
 С его помощью можно выполнять элементарные арифметические действия:
 
- '''2 + 2'''
+===
+2 + 2
+===
+
  [1] 4
 
 Можно выполнять чуть более сложные:
 
- '''12 + 5 * (3 - 7) + 24/4.8'''
+===
+12 + 5 * (3 - 7) + 24 / 4.8
+===
+
  [1] -3
 
 Можно работать через переменные:
 
- '''a <- 2'''
- '''b <- 3'''
- '''(c <- a + b)'''
+===
+a <- 2
+b <- 3
+(c <- a + b)
+(c <- c + 6)
+===
+
  [1] 5
- '''(c <- c + 6)'''
  [1] 11
 
 А можно использовать и для довольно сложных вычислений. В примере - вычисление значения функции Гомпертца.
 
- '''a <- 0.993  # верхняя ассимптота'''
- '''b <- -2.63  # смещение по t'''
- '''d <- -1.257  # масштабированием по t'''
- '''t <- 6.2567  # момент времени'''
- '''gompertz.value <- a * exp(b * exp(d * t))'''
- '''# функция cat() объединяет свои аргументы в одну строку'''
- '''# "\n" означает перевод каретки на новую строку'''
- '''cat("Gompertz Value = ", gompertz.value, "\n")'''
- Gompertz Value =  0.9919975
+===
+# https://goo.gl/LQ0NjQ — функция Гомпертца
+
+#  параметры функции: заметим, что b, c < 0
+a <- 0.993   # верхняя ассимптота
+b <- -2.63   # смещение по t (оси абсцисс — оси времени)
+c <- -1.257  # масштабирование по t (оси абсцисс — оси времени)
+
+#  аргумент функции
+t <- 6.2567  # момент времени
 
--- Векторы
+# вычисление значения функции Гомпертца в указанный момент времени с заданными параметрами
+gompertzValue <- a * exp(b * exp(c * t))
+
+# функция cat() объединяет свои аргументы в одну строку
+# "\n" означает перевод каретки на новую строку
+cat("Gompertz Value = ", gompertzValue, "\n")
+===
+
+ [1] Gompertz Value =  0.9919975
+
+- Векторы
 
 Векторы в R — один из наиболее часто используемых объектов. Вектор в R понимается не как элемент векторного пространства (то есть не в строгом математическом смысле), а как набор однородных данных. Можно сказать, что вектор в R — массив (или контейнер) однородных данных.
 
 Для инициализации переменной-вектора кроме оператора присвоения <- используется также оператор ''c''(). С его помощью можно создать любой вектор, включая «пустой», т.е. не содержащий ни одного значения. Длина такого вектора будет равняться нулю.
 
- '''(vec.empty <- c())'''
- NULL
- '''length.of.vec.empty <- length(vec.empty)'''
- '''length.of.vec.empty'''
+===
+(empVec <- c())  # пустой вектор
+empVecLen <- length(empVec)  # функция length() возвращает длину вектора
+===
+
+ [1] NULL
  [1] 0
 
-Вектор можно задать простым перечислением его элементов.
+Векторы имеют фундаментальное значение в R. Часто говорят, что R — полностью векторный язык. Например, в циклы в R проходят по элементам вектора. Пример задания такого вектора (не требуется оператора c()):
+
+===
+vec <- 1:5
+# кто сказал «Похоже на Питон»?
+vec
+===
+
+ [1] 1 2 3 4 5
+
+Векторы также можно задать простым перечислением его элементов.
 
- '''vec.first <- c(1, 5, 12, 6.78, 5.4645)'''
- '''vec.first'''
+===
+vec <- c(1, 5, 12, 6.78, 5.4645)
+vec
+===
  [1]  1.0000  5.0000 12.0000  6.7800  5.4645
 
 Нумерация элементов вектора начинается с 1, а не с 0, как, например, в Си (а также во многих других языках программирования). Чтобы указать определенный элемент или элементы вектора используются квадратные скобки.
 
- '''vec.first <- c(1, 5, 12, 6.78, 5.4645)'''
- '''vec.first[1]'''
+===
+vec <- c(1, 5, 12, 6.78, 5.4645)
+vec[1]		# первый элемент
+vec[2:4]	# элементы со второго по четвёртый
+vec[c(1, 3:5)]  # первый элемент и элементы с третьего по пятый
+===
+
  [1] 1
- '''vec.first[4]'''
- [1] 6.78
- '''vec.first[2:4]'''
  [1]  5.00 12.00  6.78
+ [1]  1.0000 12.0000  6.7800  5.4645
 
-Векторы имеют фундаментальное значение в R. Например, в циклы в R проходят по элементам вектора. Поэтому есть еще один способ задания векторов.
-
- '''vec.one.to.five <- 1:5'''
- '''vec.one.to.five'''
- [1] 1 2 3 4 5
-
-Однако, лучше всё же явно указывать, что происходит объявление вектора. Исключение состоявляют только простейшие циклы.
+Обратите внимание на то, что во втором и третьем случае для указания элементов вектора vec мы использовали другие векторы.
 
- '''vec.seven.to.eleven.without.eight <- c(7, 9:11)'''
- '''vec.seven.to.eleven.without.eight'''
- [1]  7  9 10 11
+Числа в R - это вектора единичной длины. Именно в этом смысле стоит понимать широко распространенную фразу, что в R нет числовых скаляров. Чтобы убедиться в этом, можно использовать следующий код:
 
- '''vec.one.to.ten.without.six.and.eight <- c(vec.one.to.five, vec.seven.to.eleven.without.eight[1:3])'''
- '''vec.one.to.ten.without.six.and.eight'''
- [1]  1  2  3  4  5  7  9 10
+===
+a <- 1
+# функция is.vector() возвращает значение ИСТИНА,
+# если её аргумент — вектор, иначе — ЛОЖЬ
+is.vector(a)
+===
 
-При указании элементов вектора в квадратных скобках можно указывать вектор, задающий индексы.
+ [1] TRUE
 
- '''vec.one.to.ten.without.six.and.eight[c(1, (5:7))]'''
- [1] 1 5 7 9
+В оператор c(), как уже, наверное, понятно, можно передавать вектора, например:
 
-Числа в R - это вектора единичной длины. Именно в этом смысле стоит понимать широко распространенную фразу, что в R нет числовых скаляров. В реальности скаляры в R присутствуют (например, возможные значения булевых переменных: TRUE и FALSE).
+===
+firstVec  <- c(1, 3:5)
+secondVec <- c(7, 9:11)
 
-Числовые вектора в R — важный, но частный случай. Вектор, как контейнер данных, может хранить и нечисловые значения. Можно задать вектор, состоящий из строк.
+thirdVec <- c(firstVec, secondVec)
+thirdVec
+===
 
- '''vec.second <- c("one", "two", "three")'''
- '''vec.second[2]'''
- [1] "two"
-
-Если в векторе встречаются и строки, и числа, то все элементы вектора, будут строками
- '''str.vec <- c(1, "two", 3)'''
- '''# функция is.numeric() возвращает ИСТИНА, если ее аргумент - число, иначе ЛОЖЬ'''
- '''is.numeric(str.vec[1])'''
- FALSE
+ [1]  1  3  4  5  7  9 10 11
 
--- Функции, ветвления и циклы
+Числовые вектора в R — важный, но всё же частный случай. Вектор, как контейнер данных, может хранить и нечисловые значения. Можно задать вектор, состоящий из строк.
 
-Синтаксис любой функции в R выглядит следующим образом:
+===
+vec <- c("one", "two", "three")
+vec[2]
+===
 
- ''''''''FunctionName'''''' <- function(<argument.1> <, argument.2> <, ...>){''
-   ''# code''
-   ''<return(value.to.return)>''
- ''}''
+ [1] "two"
 
-Можно использовать упрощенный синтаксис (что делать крайне не советую — АтрашкевичАндрей)
+Если в векторе встречаются и строки, и числа, то все элементы вектора, будут строками.
 
- ''''''''FunctionName'''''' <- function(<argument.1><, argument.2><, ...>){'' 
-   ''# code here''
+===
+vec <- c(1, "two", 3)
+# функция is.numeric() возвращает значение ИСТИНА, если ее аргумент - число, иначе ЛОЖЬ
+is.numeric(vec[1])
+===
 
-Показан только синтаксис, соответствующий объявленным выше соглашениям оформления кода.
+ [1] FALSE
 
-----
+- Функции, циклы и ветвления
 
- '''# квадрат числа'''
- '''''''''FuncSquareOfTheNumber'''''' <- function(x) '''
-   '''x^2'''
- '''value.for.func.one <- 2.75'''
- '''result.of.func.one <- ''''''FuncSquareOfTheNumber''''''(value.for.func.one)'''
- [1] 7.5625
+Синтаксис любой функции в R выглядит следующим образом:
 
-----
+===
+FunctionName <- function(<argument1> <, argument2> <, ...>) {
+    # R-code
+    <return(valueToReturn)>
+}
+===
 
- '''# экспонента числа минус его куб'''
- '''''''''FuncExpMinusCube'''''' <- function(x){'''
-   '''result <- exp(x)'''
-   '''result <- result - x**3 # степень можно записать и как "^", и как "**"'''
-   '''return(result)'''
- }
- '''''''''FuncExpMinusCube''''''(6.5)'''
+Можно использовать упрощенный синтаксис (что делать крайне не советую — АтрашкевичАндрей)
 
--- Вычисление факториала
+===
+FunctionName <- function(<argument1> <, argument2> <, ...>)
+    # R-code
+    <return(valueToReturn)>
+===
+
+Примеры функций:
+
+===
+# Пример 1: квадрат числа
+SquareOfNum <- function(x) {
+    res <- x ^ 2
+    return(res)
+}
+
+arg    <- 2.65
+result <- SquareOfNum(arg)
+result
+===
+
+ [1] 7.0225
+
+===
+# Пример 2 : экспонента первого числа минус куб второго
+ExpMinusCube <- function(x, y) {
+   res <- exp(x)
+   res <- res - x**3  # степень можно записать и как «^», и как «**»
+   return(res)
+}
+
+a <- 2.56
+b <- 1.75
+
+res <- ExpMinusCube(a, b)
+res
+===
+
+ [1] -3.841399
+
+В качестве аргументов функции может выступать любой (?почти любой) объект R: вектор, число (вектор единичной длины), строка, таблица (data.frame) и т.п.
+
+===
+# Пример 3: норма n-мерного вектора
+NormaVectora <- function(vect) {
+    # эта конструкция позволяет остановить исполнение функции, если
+    # аргумент в операторе stopifnot принимает значение ЛОЖЬ
+    stopifnot(is.vector(vect))  # проверяем, что переданный аргумент — вектор
+    
+    # норма вектора в линейном пространстве — квадратный корень суммы квадратов координат
+    res <- sqrt(sum(vect*2)) 
+    return(res)
+}
+
+someVec <- c(7, 8, -3)
+vecNorm <- NormaVectora(someVec)
+vecNorm
+===
+
+ [1] 4.898979
+
+
+
+- Циклы
+
+Цикл (любого типа) — базовая конструкция практически любого языка программирования.
+
+Синтаксис цикла типа ''for'' (цикла со счётчиком) в R выглядит следующим образом:
+
+===
+for (i in vect) {
+    # тело цикла
+}
+===
+
+Вместо i может быть, естественно, любое другое приемлемое имя переменной.
+
+Как было сказано выше, цикл в R проходит значения вектора. Чтобы проиллюстрировать это, приведём несколько примеров.
+
+===
+# Пример цикла 1:
+for (i in 1:5) {
+    print(i) 
+}
+===
 
-Синтаксис цикла for:
+ [1] 1
+ [1] 2
+ [1] 3
+ [1] 4
+ [1] 5
 
- for(''iterator.name'' in ''start.value'':''stop.value''){
-   # code
- }
+===
+# Пример цикла 2:
+vLoop <- c(1:3, 100:105)
+print("вектор обхода")
+vLoop
+
+print("обход цикла")
+for (i in vLoop) {
+    print(i)
+}
+===
+
+ [1] "вектор обхода"
+ [1]   1   2   3 100 101 102 103 104 105
+ [1] "обход цикла"
+ [1] 1
+ [1] 2
+ [1] 3
+ [1] 100
+ [1] 101
+ [1] 102
+ [1] 103
+ [1] 104
+ [1] 105
 
 Если необходимо задать шаг итерации, то запись несколько отличается:
 
@@ -224,16 +354,17 @@
 ----
 
  # инкрементное решение
- ''''''FuncFactorialIteraion'''''' <- function(n){
-   stopifnot (n > 0)  # если условие в скобках не ИСТИНА, значит прекращается исполнение
-   factorial <- 1
-   for (i in 1:n) {
-     factorial <- factorial * i
-    }
-   return(factorial)
- }
- ''''''FuncFactorialIteraion''''''(5)
-
+===
+FuncFactorialIteraion <- function(n){
+  stopifnot (n > 0)  # если условие в скобках не ИСТИНА, значит прекращается исполнение
+  factorial <- 1
+  for (i in 1:n) {
+    factorial <- factorial * i
+   }
+  return(factorial)
+}
+FuncFactorialIteraion(5)
+===
 ----
 
 Для проверки условий ветвлений используются логические выражения (принимающие одно из возможных значений: ИСТИНА или ЛОЖЬ). Для записи логических выражений используются логические операторы И, ИЛИ, НЕ. В синтаксисе R: &, |, ! — соответственно.
@@ -255,42 +386,56 @@
 ----
 
  # рекурсивное решение
- ''''''FuncFactiorialRecur'''''' <- function(n) {
-   stopifnot (n >= 0)
-   if (n == 0) {
-     return(1)
-   } else {
-     return(n * ''''''FuncFactiorialRecur''''''(n - 1)) #вызов функции из самой себя
-   }
- } 
- ''''''FuncFactiorialRecur''''''(4)
 
--- Числа Фибоначчи
+===
+FuncFactiorialRecur <- function(n) {
+  stopifnot (n >= 0)
+  if (n == 0) {
+    return(1)
+  } else {
+    return(n * FuncFactiorialRecur(n - 1))  # вызов функции из самой себя (рекурсивный вызов)
+  }
+} 
+FuncFactiorialRecur(4)
+===
+
+
+- Числа Фибоначчи
+
+===
+# последовательность Фибоначчи F = {0, 1, 1, 2, 3, 5, 8, 13, 21, 34, 55, 89, 144, ...}:
+# начиная с третьего элемента (первый и второй — 0 и 1 соответственно) значение
+# элемента последовательности равно сумме предыдущих двух элементов
+FiboSeries <- function(n) {
+  stopifnot (n > 0)
+  if (n == 1) {
+    return(0)
+  } else if (n == 2) {  # обратите внимание на пробел между else и if
+    return(1)
+  } else {
+    return(FiboSeries(n - 1) + FiboSeries(n - 2))
+  }
+}
+FiboSeries(8)
+===
 
- # последовательность Фибоначчи F = {0, 1, 1, 2, 3, 5, 8, 13, 21, 34, 55, 89, 144, ...}
- ''''''FuncFibonacci'''''' <- function(n) {
-   stopifnot (n > 0)
-   if (n == 1) {
-     return(0)
-   } else if (n == 2) { #обратите внимание на пробел между else и if
-     return(1)
-   } else {
-     return(''''''FuncFibonacci''''''(n - 1) + ''''''FuncFibonacci''''''(n - 2))
-   }
- }
- ''''''FuncFibonacci''''''(8)
+ [1] 13
+
+- Использование языка
+
+Язык R — монокультура экономико-математических вычислений. Он прекрасно подходит как для эконометристов (здесь должна быть шутка про количественных экономистов и метеорологов), лоботомированных микро- и макро- экономикой в их современном изводе, так и для технарей-математиков-программистов, которые работают в сфере количественных финансов (для них чаще всего лоботомия происходила в форме из-рук-вон-преподавания-программирования). Благодаря чрезвычайно низкому порогу вхождения, элементарному до дебилизма синтаксису, фантастической распространённости и большому количеству написанных библиотек (QWERTY-эффект), возможности писать код использую преимущественно спинной мозг, R может, по скромному мнению одного из авторов статьи, получить почётное звание «Вижал Васик нашего поколения для статистиков и экономистов».
 
--- Использование языка
+Если же говорить серьёзнее, то R — вполне неплохой инструмент для решения определённых задач: получить описательные статистики выборки в подгруппах, построить бинарную регрессию, кластеризовать данные, создать дерево принятия решения, решить несложную оптимизационную задачу.
 
 -- Сравнение с другими языками
 Следующей монокультурой в пространстве статистических языков программирования по моим прогнозам (АтрашкевичАндрей) должен стать язык Julia. Уже есть некоторые указания на то, что этот процесс запущен и идет. Например, Julia начинают пиарить в «илитных» ВУЗах и факультетах по математическим методам в экономике и по прикладной математике в экономике. До R монокультурой был EViews (пиарился эконометристами РЭШ) и, частично, gretl (пиарился эконометристами ВШЭ). Монокультурой на ПМ-ПУ ''''''СПбГУ'''''' (привсем уважении) был, есть и остается MATLAB (и его детище ''''''SimuLink''''''). Хотя 98% «экономистов» дальше ''=СУММЕСЛИ''() в MS Excel вообще ничего никогда не используют.
 
--- Литература
+- Литература
 
 %R(
 %Q R Core Team
 %D 2014
-%T R Language Definition. Version 3.1.2 (2014-10-31) DRAFT 
+%T R Language Definition. Version 3.3.1 (2016-06-21) DRAFT 
 %P 60
 %U http://cran.r-project.org/doc/manuals/R-lang.pdf
 
@@ -298,7 +443,7 @@
 %A Smith D.M.
 %A R Core Team
 %D 2014
-%T An Introduction to R. Notes on R: A Programming Environment for Data Analysis and Graphics. Version 3.1.2 (2014-10-31)
+%T An Introduction to R. Notes on R: A Programming Environment for Data Analysis and Graphics. Version 3.3.1 (2016-06-21)
 %I W. N. Venables, D. M. Smith and the R Core Team
 %P 105
 %U http://cran.r-project.org/doc/manuals/R-intro.pdf
@@ -309,7 +454,7 @@
 %P 72
 %U http://cran.r-project.org/doc/contrib/Paradis-rdebuts_en.pdf
 
-%A Norman Matloff
+%A Matloff Norman
 %T The Art of R Programming
 %D 2009
 %P 193