Как разбить «формулу» в R

Question

Я пытаюсь сделать небольшой пакет R с ограниченными знаниями в программировании R. Я пытаюсь использовать следующий аргумент:

formula=~a+b*X 

, где Х представляет собой вектор, «а» и «б» константы в вызове функции.

Что мне интересно, я когда-то вводил формулу, я хочу выделить (a,b) и X отдельно и использовать их для других манипуляций с данными внутри вызова функции. Есть ли способ сделать это в R?

Я бы очень признателен за любые рекомендации.

Примечание: Под редакцией мой вопрос для ясности

Я ищу что-то подобное model.matrix) выход (. Выше формула может быть более обобщенным для размещения 'п' число переменных, скажем,

~ 2 + 3 * Х + 4 * Y + ... + 2 * Z

В вывода, мне нужны коэффициенты (2 3 4 ... 2) как вектор и [1 XY ... Z] как ковариационная матрица.

Answer 1

Вопрос не совсем ясен, поэтому мы будем предполагать, что вопрос задается формулой с использованием стандартного синтаксиса формулы, как мы анализируем имена переменных (или во втором ответе имена переменных и константы), выдавая в качестве вывода вектор символов, содержащий их.

1) all.vars Попробуйте это:

fo <- a + b * X # input 
all.vars(fo) 

дает:

[1] "a" "b" "X" 

2) strapplyc Кроме того, мы могли бы сделать это с манипуляций со строками. В этом случае он также анализирует константы.

library(gsubfn) 
fo <- ~ 25 + 35 * X # input 
strapplyc(gsub(" ", "", format(fo)), "-?[0-9.]+|[a-zA-Z0-9._]+", simplify = unlist) 

дает:

[1] "25" "35" "X" 

Примечание: Если все, что вы пытаетесь сделать, это оценить РИТ формулы как выражение R, то это просто:

X <- 1:3 
fo <- ~ 1 + 2 * X 
eval(fo[[2]]) 

, давая:

[1] 3 5 7 

Обновление: Исправлено и добавлено второе решение и примечание.

Answer 2

Вызов представляет собой список символов и/или других вызовов, и его элементы могут быть доступны посредством обычных операций индексирования, например.

f <- ~a+bX 
f[[1]] 
#`~` 
f[[2]] 
#a + bX 
f[[2]][[1]] 
#`+` 
f[[2]][[2]] 
#a 

Однако обратите внимание, что в формуле bX является один символ, вы, вероятно, имел в виду b * X вместо этого.

f <- ~a + b * X 

Тогда a и b обычно будут храниться в невычисленном списке.

vars <- call('list', f[[2]][[2]], f[[2]][[3]][[2]]) 
vars 
#list(a, b) 

и vars будет передан eval в какой-то момент.