用lingo软件如何对DEA-DA模型进行编程

s-plus BEEK model是full_bekk
  做出来的结果跟matlab full_bekk_mvgarch()相近

e-views BEEH model 是diagonal_bekk
  做出来的结果跟 matlab diagonal_bekk_mvgarch()相近

matlab 则可做三种模型full,diagonal,scalar

你已在matlab运行出结果了
要看parameters
只要在command window键入parameters
或看左上角的Workspace
里面存放着使用的变量及结果

建议你重新安装Ucsd_garch
把已安装的Ucsd_garch 移除乾净
重新安装jplv7(已含最新版Ucsd_garch)
   http://www.spatial-econometrics.com/
然后在matalb设置路径时选取add with Subfolders


有關Gauss _mlhsvcp 的問題
有可能是版本問題
作者用MAXLIK Version 3.1.3   
我是用MAXLIK Version 5.0.9
_mlhsvcp沒定義到


E-VIEWS
sample s0 1 2000
取样范围 1~2000,
' initialization of parameters and starting values
''''''

sample s1 2 2000
取样范围 2~2000,
' estimate the model
' log-likelihood series
bvgarch.append logl =-0.5*(!mlog2pi + (invh1*sqres1+2*invh2*res1res2+invh3*sqres2) + log(deth))
在这里
invh1,invh2,invh3,sqres1,sqres2,deth
第一个obs都是NA
所以由第二个取样
sample s1 2 2000

71# epoh

epoh老师，为什么以下程序画不图和出不了结果？非常感谢！！！

n=1000

d=1

mixture <- function(n,x0,a,d){

x <- array(0,n)

x[1] <- x0

for(t in 2:n)

{

y <- runif(1,x[t-1]-d,x[t-1]+d)

alpha <- (exp(-1*(y+a)^2/2) +

exp(-1*(y-a)^2/2))/

(exp(-1*(x[t-1]+a)^2/2) +

exp(-1*(x[t-1]-a)^2/2))

accept <- min(1,alpha)

u <- runif(1)

if (u <= accept) x[t] <- y

else x[t] <- x[t-1]

}

plot(1:length(x),x,

type='l',lty=1,xlab="t",ylab="x")

x

}

还有以下程序也出不了结果：
met.has <- function(n,x0,sigma){
x <- array(0,n)
x[1] <- x0
for(t in 2:n)
{
y <- rnorm(1,x[t-1],sigma)
alpha <- exp((x[t-1]^2 - y^2)/2)
accept <- min(1,alpha)
u <- runif(1)
if (u <= accept) x[t] <- y
else x[t] <- x[t-1]
}
plot(1:length(x),x,
type='l',lty=1,xlab="t",ylab="x")
x
}

你定义了一个function mixture()
你要执行function mixture
必须给足四个参数
然后呼叫function mixture(n,x0,a,d)


##########
n=100
x0=1
a=0.1
d=1
mixture <- function(n,x0,a,d){
x <- array(0,n)
x[1] <- x0
for(t in 2:n)
{
y <- runif(1,x[t-1]-d,x[t-1]+d)
alpha <- (exp(-1*(y+a)^2/2) +exp(-1*(y-a)^2/2))/(exp(-1*(x[t-1]+a)^2/2) +exp(-1*(x[t-1]-a)^2/2))
accept <- min(1,alpha)
u <- runif(1)
if (u <= accept) x[t] <- y
else x[t] <- x[t-1]
}
plot(1:length(x),x,type='l',lty=1,xlab="t",ylab="x")
x
}
mixture(n,x0,a,d)

  

74# epoh
epoh老师，附件中运行markov.m出错，如何修改？谢谢！

??? Attempted to access cat1(1,23); index out of bounds because size(cat1)=[459,8].
Error in ==> Markov at 6
  category=cat1(1,category_check);

除了 cat1 的Array dimensions 不符外
接下来的
savefile13=  ['/home/jaimovich/Research/Revision_Retailer/Data/Data_Results_Backup/' int2str(category) 'vec_cat_upc'  ];
  load(savefile13);%
都会出问题.
解决办法就是跟作者所取.

Researchers interested in seeking permission to use this data should contact:
Kristin Kiesel
Assistant Professor
Department of Economics
California State University, Sacramento

76# epoh

epoh老师，为什么附件中的程序没有一个可以运行？如何可以运行每一个程序？
我近期要用这些程序，希望帮忙，非常感谢！

通通可以执行.
记得set path\add with subfolder\TVmodel

Figures:

Figure 1: >>Figure1

Figure 2: >>Figure2

Figures 3&4: >> TVCestimate_Figures3_4

Figure 5: >> Figure5  (very slow)


Tables:

Table 1:  Eviews program
                Taylor_rule_estimates.prg  
                请自行修改,一开始加入
                load greenbook.wf1

The determinacy results of Table 1: >> Table_1
                                                          >> Table1_RhoDifference

Table 2: >> Table_2

Table 3: >> Table_3

78# epoh

epoh老师，您好！
     1、为什么TVmodel中的impulse.m, vec.m, vech.m, acf.m, 等m文件在matlab中不能运行？如何运行这些m文件？
     2、附件Cartel中的matlab文件如何运行？我运行时总是提示错误。
非常感谢！
运行impulse.m时提示：
>> impulse
??? Undefined function or variable "b".
Error in ==> impulse at 4
[mm,nn]=size(b)
>>

运行Cartel中的ACDF.m时提示：
??? Input argument "IN" is undefined.
Error in ==> ACDF at 7
[n k] = size(IN);
>>

Cartel readme.pdf
Main Code :
REPLICATES_FULL_MODEL
REPLICATES_FULL_MODEL_SELECTION
REPLICATES_NOHETRO_MODEL
REPLICATES_STDERR_FULL
REPLICATES_STDERR_NOHETRO


Subroutines called by main code
ACDF
ACDFMatch
AComputeTargetCDF
ADensityIntegral
AIntergrandAtU
ALVAnCFD
ASmoothCDF
ATargetCDF
ATriKernel
ATriKernelAdj

%%%%%%%%%
以function [OUT] = ACDF(IN,data)为例
你需要自行给出参数IN,data
才能运行ACDF

同理TVmodel\impulse
你需要给出b

你说得无法运行
几乎都是没给出参数或数据造成.