Ytilisation du package GMM
Bonjour, je cherche a utiliser le package gmm (qui sert a mettre en oeuvre le améthoe des moments généralisés. je construis un exemple tres simple ou je cherche a retrouver la variance d'une loi normale. Je me suis bien sur inspiré de la syntaxe de la documentation, voici mon code:
Code:
1
2
3
4
5
6
7
8
9
10
11
| set.seed(123)
n <- 1000
x <- rnorm(n)
x <- sqrt(38) * x # on cherche donc V(x)=38
obj_fct2 <- function(theta, x) {
m1 <- (x - mean(x))^2 - theta[1]
return(m1)
}
p2 <- gmm(g = obj_fct2, 1 , x = x)
coef(p2) #Devrait contenir une estimation de V(x)=38 |
Or j'obtiens le message suivant:
Citation:
Erreur dans colMeans(gt) : 'x' must be an array of at least two dimensions
De plus : Message d'avis :
In optim(P$t0, .obj1, x = P$x, w = w, gf = P$g, ...) :
one-diml optimization by Nelder-Mead is unreliable: use optimize
je ne comprends pas pourquoi x doit etre un tableau, il y a un exemple tres similaire dans la documentation du package ou 'x' est une liste. Merci de me donner votre point de vue.
nicolas
