Fonction indicatrice en R

[invite00f7998d]

svp j'ai besoin d'introduire une fonction indicatrice sur R et je n'arrive pas à le faire j'ai besoin d'un personne qui maitrise bien le R et merci d'avance
-----

[invitec5eb4b89]

Bonjour,

ça dépend de la forme de l'ensemble à "indiquer". On note A cet ensemble.

Si A est "discret", on peut le faire de cette façon :

A <- 1:10
ind <- function(x,A) ifelse(x %in% A, 1,0)
X <- seq(0,10,by=0.01)
plot(X,ind(X,A),'h')

Ou si A est "continu", par exemple A = [0,1], on peut essayer de cette façon :

ind <- function(x,a,b)  ifelse(x >= a & x <= b, 1,0)
X <- seq(-2,3,by=0.01)
plot(X,ind(X,0,1),'s')

J'espère que ça aide. Il y a bien sûr d'autres façons de faire !
Bon courage,
V.

ps : j'ai supposé que l'ensemble à indiquer était un ensemble de réels, mais si c'est un ensemble vectoriel ou complexe, ça se complique !

[invite00f7998d]

merci bien pour votre réponse mais mon problème c'est:
Soit f une densité définie par f(x)=(1/2) x^(-1/2) Si x est entre 0et 1 et 1,0 sinon.
Simuler un échantillon x de taille 10000 de la loi de densité f. Simuler un second échantilon y de même taille, de la méme loi. Représenter dans le plan les couples de points (xi,yi).
j'ai une soutenance pour lundi et je n'a pas encore entamé le 1er exercice j'ai un module de probabilité en crédit et s'il vous plait vous pouver m'aider et merci d'avance.

[invite986312212]

bonsoir,

tu comptes t'y prendre comment pour simuler tes échantillons?

[A voir en vidéo sur Futura]

[invite00f7998d]

je n'a compris votre message mais il faut programmer sur R alors il faut décalrer la fonction puis faire la simulation

[invite986312212]

mais il y a loin de la simple programmation de la fonction de densité à la génération d'un échantillon. Si tu as cet exercice à faire c'est qu'on a dû t'enseigner au moins une méthode de simulation. La plus simple, et qui marche dans tous les cas pour peu que la densité soit bornée est la méthode d'acceptation-rejet (rejection method en Anglais) la connais-tu?

[invite986312212]

ah mais justement je vois que ta densité n'est pas bornée. en ce cas, le plus simple est d'inverser la fonction de répartition. Ca va t'obliger à calculer une intégrale (pas dure).