Discussion :

[marine44]

bonjour je suis étudiante en L3 et j'ai un exercice à rendre en R:
1) recupérer le fichier de données L3ingbinom5.txt à l'adresse suivante
http://www.math.sciences.univ-nantes.fr/~philippe/data/L3ingbinom5.txt"
il contient la réalisation d'un échantillon d'une loi binomiale de parametres (5,0.3)

j'ai donc fais

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

x=scan(file="http://www.math.sciences.univ-nantes.fr/~philippe/data/L3ingbinom5.txt")

2) construire une vecteur qui contient les moments empiriques(mu 1(1),...,mu n(6)) ou mu n(k) le moment empirique d'ordre k est égal à mu n(k)=1/n somme de i=1 à n de Xi^k
utiliser la boucle for pour construire ce vecteur.

j'ai donc posé n=10000 car il y a 10000 éléments dans x

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
n=10000
i=1:n
mu.estim=1:6
for (k in 1:6)
{ mu.estim[k]=(1/n)sum(x^k) }

mais ça ne marche pas et je ne comprend pas pourquoi, pourriez vous m'aider ?
merci

[Sengar]

bonjour je suis étudiante en L3 et j'ai un exercice à rendre en R:

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
n=10000
i=1:n
mu.estim=1:6
for (k in 1:6)
{ mu.estim[k]=(1/n)sum(x^k) }

Salut,

En général il vaut mieux poster l'erreur que te donnes R si tu veux qu'on puisse t'aider plus facilement.
Ici le problème est qu'il manque un * entre 1/n et ton sum

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

mu.estim[k]=(1/n)*sum(x^k)

Au passage : tu peux simuler toi même 10000 réalisations d'une binomiale 5,0.3 en faisant rbinom(10000,5,0.3)

[marine44]

ah oui ça marche, merci

[marine44]

ah non ça ne marche toujours pas car quand je pose

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
y=for (k in 1:6){
 mu.estim[k]=(1/n)*sum(x^k) }

j'obtiens

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
> y
NULL

alors que je devrais obtenir un vecteur

[Sengar]

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
<s>y=</s>for (k in 1:6){
 mu.estim[k]=(1/n)*sum(x^k) }

Ton vecteur est déjà créé, c'est mu.estim.
Si tu veux le renommer y tu fais:

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

y<-mu.estim

[marine44]

ah d'accord ! merci beaucoup