n <- 100
kx <- 4
sig <- matrix(c(1,0.5,0.5,1),2,2)
e <- matrix(rnorm(n*2),n,2)%*%chol(sig)
q <- matrix(rnorm(n),n,1)
x <- matrix(rnorm(n*kx),n,kx)
z2 <- cbind(matrix(1,n,1),q)
z1 <- x%*%matrix(1,kx,1)*3+e[,2]
zz <- (z1+z2%*%matrix(1,2,1))*.1
y <- zz*(q<0)-zz*(q>=0)+e[,1]
qhat <- gmm_thresh(y,z1,z2,x,q,.9,.8,1,1)
以上这些命令都是什么意思


雷达卡



京公网安备 11010802022788号







