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

z兄，fmincon是matlab自带的
从help也可以看出，如果不是自带的，是没有详细的例子的

上面EGARCH里面有很多m文件，你分开运行过吗？如果没有，那是需要分开的

EM.R

少定义nc,ni.nt

请自行补上.

另两个问题,ywh兄所言甚是.

### Example 3 SEM algorithm

#First, run EM using the code in the basic EM algorithm given in Example 1.

#Then we run SEM as follows.  

n.obs<-c(85,196,341)

nc<-n.obs[1]

ni<-n.obs[2]

nt<-n.obs[3]

#Starting values for SEM

epoh 发表于 2011-8-9 10:34
EM.R少定义nc,ni.nt请自行补上.另两个问题,ywh兄所言甚是.### Example 3 SEM algorithm #First, run EM usi ...

fortrandata.rar (1.43 MB)

2011-8-9 21:56:52 上传

epoh老师，您好！
    附件中的fortran程序如何运行？需要下载什么软件才能运行firtran程序？
另外，附件中的matlab程序运行出现以下错误？为什么？如何修改？

1\Simulatione.m

simulation

s =

     1

??? Undefined function or variable "np".

Error in ==> simulatione at 41

year=zeros(np*TS,1);

>>

2\solvee.m

Warning: Input arguments must be scalar.

> In transzx at 3

  In solvee at 8

Warning: Size vector should be a row vector with integer elements.

> In transzx at 3

  In solvee at 8

Warning: Input arguments must be scalar.

> In transzx at 3

  In solvee at 8

Warning: Size vector should be a row vector with integer elements.

> In transzx at 3

  In solvee at 8

??? Undefined function or variable "Tx11".

Error in ==> transzx at 10

Tzx11=Tz.*Tx11;

Error in ==> solvee at 8

transzx;

>>

3\xprobit.m

??? Undefined function or variable "nkgrid".

Error in ==> xprofit at 5

pix=zeros(ngrid,nkgrid);

>>

Dprobit.m

??? Undefined function or variable "nkgrid".

Error in ==> dprofit at 6

pid=zeros(ngrid,nkgrid);

>>

Parameters.m

运行没有任何反应？

还有cd.doc中的stata错误如何修改？
非常感谢！

很久没用fortran了

由于mcmc.f90要用到IMSL.mod,numerical_libraries.mod

PROGRAM mcmc

!=======================

!Use these modules

!=======================

USE IMSL

USE numerical_libraries

USE parameterse

单独抓IMSL.mod,numerical_libraries.mod都失败

最后使用内含IMSL Fortran Library的

Compaq Visual Fortran 6.6 Pro

才顺利Build mcmc.exe

mcmc_fortran.rar

   mcmc_fortran.rar (401.58 KB)

2011-8-19 09:42:47 上传

double click mcmc.exe

执行时会自动产生2个文件

  accept.txt,

  param.txt,

zhangtao兄

Rdexp_summary内的readme

已经说明执行function是main.m

依序执行load_data,parameters,.....

如此才不会产生错误.

%LOAD DATA

load_data;

%DEFINE PARAMETERS GIVEN

parameters;

%DEFINE GRIDS OF CAPITAL AND PRODUCTIVITY (X)

gridkx;

......

......

epoh 发表于 2011-8-20 19:03
zhangtao兄Rdexp_summary内的readme已经说明执行function是main.m依序执行load_data,parameters,.....如此才 ...

epoh老师，您好！
       我想问一下：如何把R中的tsDyn这个包载入到splus中进行运用和进行相关的
计量分析？

哈哈!这个有困难.

epoh 发表于 2011-8-22 20:58
哈哈!这个有困难.

> setwd('C:/R_threg')
> melanoma<-read.table("melanoma.txt",header=TRUE)
> threg(data=melanoma,lny0=age,mu=thick+age+sex,timevar=survtime,failure=status)
错误: 没有"threg"这个函数
R_threg.rar (270.54 KB)

2011-8-22 21:14:22 上传

>

epoh老师，您好！
    附件中R程序运行出现以上问题，您看如何解决？
非常感谢！

另外，tsDyn包载入splus中运行，您说有困难？那困难在什么地方？
能实现吗？如果能实现，如何实现？

melanoma<-read.table("melanoma.txt",header=TRUE)
source("threg.R")
threg(data=melanoma,lny0=age,mu=thick+age+sex,timevar=survtime,failure=status)
               Coef. Est. Std. Err.      z  P>|z|
lny0: age       0.0033790 0.0040379  0.837 0.4026
lny0: constant  3.5621735 0.2288590 15.565 0.0000
mu: sex        -0.0079129 0.0045327 -1.746 0.0808
mu: age        -0.0003492 0.0002030 -1.720 0.0854
mu: thick      -0.0029611 0.0007753 -3.819 0.0001
mu: constant    0.0370407 0.0105242  3.520 0.0004
#######

Q2:tsDyn包载入splus中运行，有困难？困难在什么地方?

   你是否还记得package"pkgutils"

   s-plus就是用pkgutils,

   建构供s-plus使用的package

   目前R与S-plus语法已经相差很多

   要用pkgutils建构符合s-plus用的

   tsDyn package 肯定会出错.

遇到同样的问题，正寻找用lingo做DEA的教程呢？有没有大虾分享或推荐一下呀。。