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 : Sélectionner tout - Visualiser dans une fenêtre à part
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:

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