Блог

Я запутался в этом в течение последнего часа, и может кто-нибудь мне помочь.

Вот мой код, в котором я пытаюсь создать 2 функции, которые будут использоваться в tibble() .

 n <- 1000
w <- function(x) {
  if(-2 <= x amp; x <= 2) {
    return((-3/32)*x^2   3/8)
  }else{
    return(0)
  }
}
y <- function(x){
  if(0 <= x amp; x <= 8){
    return((1/8)*x^(-2/3) - (1/32))
  }else{
    return(0)
  }
}
df <- tibble(value = seq(2, 9, length.out = n), W = w(value), Y = y(value))
head(df)
 

Как вы видите, W возвращает отрицательное значение, когда значения больше 2, но в моей функции w() есть оператор if, чтобы ограничить это.

Итак, что происходит и как мне решить эту проблему?

Вы использовали if / else , который работает для одного ввода (скалярный), и вы передаете ему несколько значений (вектор). Вы должны использовать ifelse который векторизован.

Если я правильно вас понял w , вы хотите применить формулу (-3/32)*x^2 3/8 , но если значение отрицательное, присвоите значение 0. То же самое для y того, где вы хотите применить формулу (1/8)*x^(-2/3) - (1/32) . Вы можете изменить свою функцию на.

 w <- function(x) {
  pmax((-3/32)*x^2   3/8, 0)
}

y <- function(x){
  pmax((1/8)*x^(-2/3) - (1/32), 0)
}

df <- tibble(value = seq(2, 9, length.out = n), W = w(value), Y = y(value))
head(df)
 

pmax выбирает максимум из обоих переданных входных данных. Например,

 pmax(-5:5, 0)
#[1] 0 0 0 0 0 0 1 2 3 4 5
 

Вот код, который выполняет то, что вы хотите, для гораздо меньшего тестового примера, который по-прежнему позволяет просматривать все диапазоны ваших значений (я взял на себя смелость расширить диапазон тестирования значений, потому что в исходном случае, когда x было от -2 до 2, был только один экземпляр, и это дало 0для w(2) независимо.

 n <- 10
w <- function(x) {
    ifelse( -2 <= x amp; x <= 3 , (-3/32)*x^2   3/8, 0) 
                 }
y <- function(x){
    ifelse (0 <= x amp; x <= 8, (1/8)*x^(-2/3) - (1/32), 0)
                  }
df <- tibble(value = seq(2, 9, length.out = n), W = w(value), Y = y(value))
#-----------

df
# A tibble: 10 x 3
   value      W       Y
   <dbl>  <dbl>   <dbl>
 1  2     0     0.0475 
 2  2.78 -0.348 0.0320 
 3  3.56  0     0.0224 
 4  4.33  0     0.0158 
 5  5.11  0     0.0109 
 6  5.89  0     0.00708
 7  6.67  0     0.00404
 8  7.44  0     0.00154
 9  8.22  0     0      
10  9     0     0      
 

Более ранние наблюдения:
Ваш код выдал ошибку, и я исправил ее, внеся пару незначительных изменений в w функцию:

 w <- function(x) {
    if( (-2 <= x) amp;amp; (x <= 2)) {
        (-3/32)*x^2   3/8)
    } else { 0   }
}
 

Я думаю, что использование ‘return’ могло вызывать проблемы.

Выражение в вашей первой версии выдало эту ошибку:

 Error: unexpected '<=' in:
"w <- function(x) {
    if(-2 <= x <="
 

… потому что математически разумное -2 <= x <= 2 просто не подходит для синтаксического анализа в R.