Discussion :

[maiga1]

Bonjour je suis un nouveau utilisateur de R et cherche à résoudre le problème dans le document joint:
je pense procéder cette manière mais je n’arrive pas à mettre les têta en valeur absolue.

Code :

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

1
2
3
4
5
6
7
8
n <- 100  # choisir une valeur
beta <- 0.6
x <- runif(n,-3, 3)
y <- rnorm(n, beta*x, 1/sqrt(tau)) # produire des valeurs de y avec comme moyenne beta*x et d'écarts-types 1/racine carrée de tau
b <- sum(x*y)/sum(x^2) #  calculer l'estimation des moindres carrés
sdb <- 1/sqrt(tau*sum(x^2))
bval <- seq(b-4*sdb, b+4*sdb, length=200)
plot(bval, dnorm(bval, b, sdb), type="l")

Merci d'avance
cordialement