Come definire una funzione a tratti in R

Question

Voglio definire una funzione a tratti usando R, tuttavia, il mio codice R va errato. Qualsiasi suggerimento è benvenuto.

x<-seq(-5, 5, by=0.01) 
    for (x in -5:5){ 
    if (-0.326 < x < 0.652) fx<- 0.632 
    else if (-1.793<x<-1.304) fx<- 0.454 
    else if (1.630<x<2.119) fx<-0.227 
    else fx<- 0 } 

Answer 1

Oppure si potrebbe usare ifelse.

fx <- ifelse(x > -0.326 & x <0.625, 0.632, 
    ifelse(x > -1.793 & x < -1.304, 0.454, 
    ifelse(x > 1.630 & x < 2.119, 0.227, 0))) 

Answer 2

Forse se si dividono le condizioni

if((-1.793<x) & (x < 0.652)) ... 

EDIT: Questo sembra non essere tutto, ecco un approccio diverso:

x<-seq(-5, 5, by=0.01) 
fx <- function(x) { 
    res <- rep(0, length(x)) 
    res[(-0.326 < x) & (x < 0.652)] <- 0.632 
    res[(-1.793<x) & (x < (-1.304))] <- 0.454 
    res[(1.630<x) & (x <2.119)] <- 0.227 
    return(res) 
} 
fx(x) 

Answer 3

Prova questa:

x <- seq(-5, 5, 0.01) 
fx <- (x > -0.326 & x <0.625) * 0.632 + 
     (x > -1.793 & x < -1.304) * 0.454 + 
     (x > 1.630 & x < 2.119) * 0.227 
plot(x, fx) 

Answer 4

meno che non si hanno diversi punti di taglio, userei switch. Ecco un esempio con valori di taglio semplificati.

xcuts<-1:10 #the values at which you change fx assignment 
xx<- seq(1.5,10,5, by =10) #vector of fx values to be selected 
switch(max(which(x>xcuts)), 
1= fx<-xx[1], 
2= fx<-xx[2], 
..."et cetera"... 
) 

Ciclo su x.

Answer 5

Un'altra opzione, questa volta utilizzando cut.

regions <- c(-Inf, -1.793, -1.304, -0.326, 0.652, 1.63, 2.119, Inf) 
group <- cut(x, regions) 
f_values <- c(0, 0.454, 0, 0.632, 0, 0.227, 0) 
(fx <- f_values[group]) 

Answer 6

Sono un po 'in ritardo alla festa, ma non ho potuto resistere la pubblicazione di un altro paio di modi per farlo. Entrambi sfruttano le capacità R per lavorare con intervalli sulla linea reale.

Se si definiscono i punti di taglio e valori della funzione nei vettori cuts e vals in questo modo:

cuts <- c(-Inf, -1.793, -1.304, -0.326, 0.625, 1.630, 2.119) 
vals <- c( 0, 0.454,  0, 0.632,  0, 0.227,  0) 

Quindi è possibile utilizzare findInterval per cercare in modo efficiente i valori di x nei vostri punti di divisione:

fx <- vals[findInterval(x, c(-Inf, cuts))] 

Se questa funzione è necessaria per creare cose più divertenti di un semplice valore costante, è possibile inserire espressioni, funzioni o qualsiasi cosa in vals, eventualmente utilizzare ng a list se vuoi.

In alternativa, dal momento che questa funzione è una funzione di passo, è possibile utilizzare stepfun:

f <- stepfun(cuts[-1], vals) 
fx <- f(x) 

Poi è anche arrivare a utilizzare i metodi di plottaggio belle di stepfun troppo.