O problema está na utilização de &&
em vez de &
. Da documentação, help('&&')
(ênfase minha.):
&
and &&
indicate logical AND and |
and ||
indicate logical OR. The
shorter form performs elementwise comparisons in much the same way as
arithmetic operators. The longer form evaluates left to right
examining only the first element of each vector. Evaluation proceeds
only until the result is determined. The longer form is appropriate
for programming control-flow and typically preferred in if
clauses.
Tradução Google, edita por mim. Ênfase minha.
&
e &&
indicam E lógico e |
e ||
indicam OU lógico. A forma
mais curta realiza comparações elementares da mesma maneira que os
operadores aritméticos. A forma mais longa avalia da esquerda para a
direita examinando apenas o primeiro elemento de cada vetor. A
avaliação prossegue apenas até o resultado ser determinado. A forma
mais longa é apropriada para programação de controle de fluxo e
normalmente preferida em cláusulas if
.
Ou seja, a sua função, ao utilizar &&
só avalia x[1]
. Como determina imediatamente o valor lógico, FALSE
ou TRUE
, segue para o ramo apropriado. Ora, x[1]
é o primeiro valor de seq(-4, 4, by = 0.1)
, e é tomado o segundo ramo, correspondente à distribuição de Cauchy.
A solução é quase sempre utilizar a forma vetorizada de if/else
, que é ifelse
, e da operação lógica, &
. É isso que o código da 2ª versão da função faz.
fdp.GUCU <- function(x){
if(x>-1 && x < 1){
(0.55)*exp(-((x)^2)/2)/sqrt(2*pi)
}else{
(1.249042)*1/(pi*(1+(x)^2))
}
}
fdp.GUCU2 <- function(x){
ifelse(x > -1 & x < 1,
0.55*exp(-(x^2)/2)/sqrt(2*pi),
1.249042/(pi*(1 + x^2))
)
}
x <- seq(-4, 4, by = 0.1)
plot(x, fdp.GUCU(x), type="l", ylim = c(0, 0.4))
lines(x, fdp.GUCU2(x), col = "red", lty = "dotted")
Edição.
A função abaixo utiliza as funções do R base diretamente e dá os mesmos resultados que a minha fdp.GUCU2
.
fdp.GUCU3 <- function(x){
ifelse(x > -1 & x < 1, 0.55*dnorm(x), 1.249042*dcauchy(x))
}
x <- seq(-4, 4, by = 0.01)
plot(x, fdp.GUCU3(x), type="l", ylim = c(0, 0.4))
lines(x, dnorm(x), col = "blue", lty = "dashed")
lines(x, dcauchy(x), col = "red", lty = "dotted")
legend(x = "topright", legend = c("OP", "Normal", "Cauchy"),
col = c("black", "blue", "red"), lty = c("solid", "dashed","dotted"))
0.55
(Normal) e1.249042
(Cauchy)?