Ошибка при использовании пакета «prob» в функции R

Question

Я пытаюсь написать функцию, которая использует пакет prob для вычисления условных вероятностей. При использовании функции я продолжаю сталкиваться с той же ошибкой, которая указывает, что объект внутри функции не может быть найден.

Ниже приведен воспроизводимый пример, в котором я вычисляю условную вероятность без функции, а затем пытаюсь использовать эту функцию для получения того же результата. Я не уверен, что ошибка связана с ограничениями с пакетом prob или ошибкой с моей стороны.

# Load prob package 
library(prob) 

# Set seed for reproducibility 
set.seed(30) 

# Sample data frame 
sampledata <- data.frame(
    X <- sample(1:10), 
    Y <- sample(c(-1, 0, 1), 10, replace=TRUE)) 

# Set probability space 
S <- probspace(sampledata) 

# Subset Y between -1 and 0 
A <- subset(S, Y>=-1 & Y<=0) 

# Subset X greater than 6 
B <- subset(S, X>6) 

# Compute conditional probability 
P <- prob(A, given=B) 

Приведенный выше код дает следующую вероятность:

> P 
[1] 0.25 

Попытка написать функцию для вычисления ту же вероятность:

# Create function with data frame, variables, and conditional inputs 
prob.function <- function(df, variable1, variable2, state1, state2, cond1){ 
    s <- probspace(df) 
    a <- subset(s, variable1>=state1 & variable1<=state2) 
    b <- subset(s, variable2>cond1) 
    p <- prob(a, given=b) 
    return(p) 
} 

# Demonstrate the function 
test <- prob.function(sampledata, Y, X, -1, 0, 6) 

Эта функция дает следующее сообщение об ошибке:

Error in eval(expr, envir, enclos) : object 'b' not found 

Любая помощь, которую вы можете предоставить, будет отличной.

Спасибо!

Answer 1

Это похоже на ошибку в prob.

Когда я запускаю это в Vanilla R, я получаю ту же ошибку. Но когда я создаю объект b в моем рабочем пространстве, в disapears ошибки:

> print(b) 
Error in print(b) : object 'b' not found 
> test <- prob.function(sampledata, Y, X, -1, 0, 6) 
Error in eval(expr, envir, enclos) : object 'b' not found 
> 
> b <- "dummy variable" 
> print(b) 
[1] "dummy variable" 
> test <- prob.function(sampledata, Y, X, -1, 0, 6) 
> test 
[1] 0.25 
> 

В качестве временного решения, просто создать фиктивный b в текущей среде.

---

Что касается ошибки, если вы посмотрите на исходный код для prob.default (что в приведенном выше примере является то, что, в конечном счете prob(a, given=b) вызова), вы увидите следующий раздел:

if (missing(given)) { 
    <cropped> 
} 
else { 
    f <- substitute(given) 
    g <- eval(f, x)     <~~~~ 
    if (!is.logical(g)) {   <~~~~ 
     if (!is.data.frame(given)) <~~~~ 
      stop("'given' must be data.frame or evaluate to logical") 
     B <- given 

    } 
    ... 
    <cropped> 
} 

его прыгает от g до given, возможно, случайно? Я хотел бы связаться с сопровождающим пакета, так как это может быть недосмотр.

Answer 2

I не думаю, что это ошибка в пакете prob.

Во-первых, вы должны создать вам sampledata, как

sampledata <- data.frame(
    X = sample(1:10), 
    Y = sample(c(-1, 0, 1), 10, replace=TRUE)) 

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

Во-вторых, вы не должны звонить subset() внутри функции.Используйте кронштейн Подменю вместо:

prob.function <- function(df, variable1, variable2, state1, state2, cond1){ 
    s <- probspace(df) 
    a <- s[s[[variable1]]>=state1 & s[[variable1]]<=state2, ] 
    b <- s[s[[variable2]]>cond1, ] 
    p <- prob(a, given=b) 
    return(p) 
} 

И пройти variable1 и variable2 как строки:

test <- prob.function(sampledata, "Y", "X", -1, 0, 6) 

Теперь у вас есть test==0.25, и никакой ошибки.

Ссылки на то, что происходит:

• http://adv-r.had.co.nz/Computing-on-the-language.html#non-standard-evaluation-in-subset
• Assignment operators in R: '=' and '<-'
• Why is `[` better than `subset`?