j'aimerai simuler la densité suivante : f(x) = (x**α−1 * (1 − x)**β−1) / (δ + (x − 1/2 )**2)
j'ai mis en place ce bout de code mais j'ai des erreurs
avez-vous une idée de comment l'adapter ?
Code : Tout sélectionner
fb_density <- function(x){
return ((x^(alpha-1))*((1 - x)^(beta-1)))/ (delta + (x - 1/2)^2)
}
simule.densite0<-function(n, alpha, beta, delta) {
ok<-FALSE
while(!ok) {
x<-runif(n, alpha, beta)
y<-runif(1, 0, beta(alpha, beta)/delta)
ok <- y<fb_density(x)
}
x
}
l'erreur que j'obtiens
> simule.densite0(3,2,3,3)
Error in fb_density(x) : objet 'alpha' introuvable
Merci pour votre aide