利用s-plus进行GMM估计

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

    各位大虾好！现正在做一个ckls模型的实证分析，利用GMM进行参数估计，已经找到一段程序

Ckls.moments <- function(parm, data = NULL,dt=1/365){

# parm = (alpha, beta, sigma, gamma)’

# data =[r(t+dt)-r(t), r(t)]

# dt = discretization step

E.hat = as.vector(data[,1] - (parm[1] + parm[2]*data[,2])*dt)

M2 = e.hat*(as.vector(data[,2]))

M3 = e.hat^2-dt*parm[3]*parm[3]*(as.vector(data[,2]^(2*parm[4])))

M4 = m3*data[,2]

Cbind(e.hat,m2,m3,m4)

d.ckls.ts <- read.table(file = "C:\\Users\\Administrator\\Desktop\\d.ckls.ts.txt", header = TRUE)（一阶差分）

t.ckls.ts <- read.table(file = "C:\\Users\\Administrator\\Desktop\\t.ckls.ts.txt", header = TRUE)（一阶滞后）

data.ckls.ts <- data.frame(d.ckls.ts,l.ckls.ts)

collds(data.ckls.ts)[1]=”RF.diff”

Start.vals = c(0.06,-0.5,1,1)Names(start.vals) = c(“alpha”,”beta”,”sigma”,”gamma”)Gmm.ckls = GMM(start.vals,ckls.moments,ts=T,data =data.ckls,dt =1/365)Summary(gmm.ckls) Call:GMM(start=start.vals,moments = ckls.moments, ts = T, data=data.ckls, dt=1/365)小女菜鸟一枚，发现GMM这个函数不存在，在S-PLUS中不能直接计算，求助！

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏0 回帖

关键词：GMM估计 PLUS Plu GMM Moments 程序 模型