两篇文献快速瞄了下

发现底下两个错误

应该还有,请仔细检查修改(含init)

l^(1+kappa) = 1/ul*(1+tauc)*(c-h*c(-1))^(-gam)*(1-taul)*(1-alp)*y;

改为  l^(1+kappa) =(1/ul*(1+tauc))*(c-h*c(-1))^(-gam)*(1-taul)*(1-alp)*y;

f = del0+del1*(v-1)+(del2*(v-1)^2)/2;

改为   f = del0+del1*(v-1)+(del2/2)*(v-1)^2;

%%%%%

有疑问之处是底下这行,不知从何而来.

tauk = psitk*yhat+gammatk*sb(-8)+utk;

%%%%%

底下这两行

Dynamics of fiscal financing.pdf

是在 Steady state谈到

而你是放在model?

i = del0*k;

sb = b/y;

非常感谢epoh！ you are very warm hearted! thank you!

'stochastic simulation of a DSGE model for Brazil'文章的MODAL是根据log-linearied的方程编写的，而我是根据Equilibrium conditions来编写的，从这方面来说，是不是我全部都错了,也应根据log-linearied来编写，"dynamics of fiscal financing in the United States"这篇文章的附录中有log-linearied的方程。

但是'stochastic simulation of a DSGE model for Brazil'也不是完全根据log-linearied，除此之外，还加了其他如i_plot之类的方程与log-linearied几乎一样，但也进行了一些修改，设定了tempo这个参数，真不懂这是怎么回事呢。我是不是也需要这样设置呢？

epoh老师说的很正确， 我的确是将steady state的方程放到了model中，这也是出错的原因之一吧。
等我修改一下，把log-linearied的方程放到model中，看看怎么样

谢谢epoh老师！

Government investment and fiscal stimulus

这篇也有人这样编写过

他因为没设好initval;

所以有执行问题

不过可参考.

  

%%%%%

设定tempo这个参数

暂且先不要考虑

勉的趋于复杂.

按照log-linearied，方程有18个，而变量有28个，所以还得添加方程。

还有方程的编写问题，
1、S的二阶导数在dynare中如何编写呢？是作为参数吧？
2、log-linearied中出现了y*yhat, c*chat等的形式，是将y跟yhat作为两个变量吧?
如果y也是一个变量，那是不是能添加关于Y的方程呢？

还请epoh老师多多指教，谢谢！

哈哈!你太忙了

log-linearied model的程序,已经有善心人仕帮你写好了.

忙得没注意到我早放在184楼.

S的二阶导数是作为参数

关于log-linearied中出现了y*yhat, c*chat等的形式，我想出了一个办法，就是令C/Y,Z/Y,I/Y,K/Y,B/Y,G/Y,等为一个参数，这样方程的VAR为21个，方程有18个，还差三个了。

太马虎了，我是大马虎！数错了，是23个变量，23个方程，这就对了，我输入程序，看看能不能运行。

哈哈，，果然没看到呀！哎，名副其实的大马虎了！
谢谢epoh老师！，我这就运行一下啊，等我消息哈

epoh老师，thank you so much! I succeed finally with your help!

哈哈!果然底子够硬.

"泡泡"瞬间变成美丽的彩虹.

恭喜!加油!

留个脚印。。。

冒昧在這請教各位~
下面的方程式適何用MATLAB寫嘛

Ɵrr = Max Σ (s=1~s) Urs Xrs

St. Ɵrp= Σ (s=1~s) Urs Xps <= 1 ,p=1,2,...,r
      Ur1>=2Ur2>=3Ur3>=...>=sUrs
      Urs>=2 / n(s(s+1))

會使用矩陣如下

        1        2        ...        s        total         
1        X11        X12        ...        X1s        n        <-Ɵ11
2        X21        X22        ...        X2s        n        <-Ɵ22
...        ...        ...        ...        ...        n         
r        Xr1        Xr2        ...        Xrs        n        <-Ɵrr
total        n        n        n        n                  

X11-Xrs是已知的數字~
最後可得Ɵ11到Ɵrr之大小順序
公式來源是:The voting analytic hierarchy process method for selecting supplier
這個式子則是由DEA發展而來

大家有会解CCGSS的吗？这方面资料简直太少了。

您的回复对我很有启发。但是因为论坛币的限制，我下不了附件。可以发给我一份吗？494774692@qq.com

没有论坛币
你可以随便发个文件出售
订价100 or 500都可以
由我购买,如此你才能购买更多需要的资料

您的热心帮助，我真的感动不已。
早就注册了这个论坛，但是来的很少，因为什么资料都要钱。
但是很多问题在网上一搜，就被链接指引到这里来了。
为了以后的资料，我就不客气了，这个附件500坛币。
再次说，太感谢了~
买资料地址：https://bbs.pinggu.org/forum.php?mod=viewthread&tid=1497255&page=1&extra=#pid13399335

老师，先前那个卖资料的帖子被版主删了，我的500论坛币飞了。。。我想下个年鉴，所以我又重新发了一个，这次只要100论坛币。。。资料地址：https://bbs.pinggu.org/thread-1499874-1-1.html

epoh老师你好，
我一直在试着在matlab中用dynare，
这是我的code：
>> addpath c:\dynare\4.3.1\matlab
>> cd c:\dynare\4.3.1\work
>> dynare rbc_indivisible_cu.mod
但是结果总是出现错误：
Configuring Dynare ...
[mex] Generalized QZ.
[mex] Sylvester equation solution.
[mex] Kronecker products.
[mex] Sparse kronecker products.
[mex] Local state space iteration (second order).
[mex] Bytecode evaluation.
[mex] k-order perturbation solver.
[mex] k-order solution simulation.
[mex] Quasi Monte-Carlo sequence (Sobol).
[mex] Markov Switching SBVAR.

Starting Dynare (version 4.3.1).
Starting preprocessing of the model file ...
ERROR: rbc_indivisible_cu.mod:40.22-26: Unknown symbol: omega

Error using dynare (line 114)
DYNARE: preprocessing failed

您能帮我看一下吗，麻烦您啦

你的code是对的,这个程序很容易执行.
   http://www.dynare.org/documentat ... nd-support/examples[/url]  
addpath c:\dynare\4.3.1\matlab
cd c:\dynare\rbc
dynare rbc_indivisible_cu
%%%%%%%
Configuring Dynare ...
[mex] Generalized QZ.
[mex] Sylvester equation solution.
[mex] Kronecker products.
[mex] Sparse kronecker products.
[mex] Local state space iteration (second order).
[mex] Bytecode evaluation.
[mex] k-order perturbation solver.
[mex] k-order solution simulation.
[mex] Quasi Monte-Carlo sequence (Sobol).
[mex] Markov Switching SBVAR.

Starting Dynare (version 4.3.1).
Starting preprocessing of the model file ...
Found 8 equation(s).
Evaluating expressions...done
Computing static model derivatives:
- order 1
Computing dynamic model derivatives:
- order 1
Processing outputs ...done
Preprocessing completed.
Starting MATLAB/Octave computing.


STEADY-STATE RESULTS:

y                    0.916363
c                    0.707811
k                    9.29097
i                    0.208551
u                    0.891826
l                    0.30979
y_l                  2.95802
z                    0
....
....
....

大哥，那能教教deap怎么用吗？下了一个怎么用不了呢？

epoh老师，您好！请问你了解关联两阶段dea的lingo编程吗？我编的程序会出现负值或者大于1的情况，不知道哪里错了，请老师指教一下！感激不尽~~~

先马克留着！

求问，是不是DEAP做出来的结果和用lingo 做出来的效率值不一样

你好，我用了三个序列，然后运行了一半出现了问题，请问有好心人能帮忙解决吗~
Press any key to estimate.Warning: Options LargeScale = 'off' and Algorithm = 'trust-region-reflective' conflict.
Ignoring Algorithm and running active-set method. To run trust-region-reflective, set
LargeScale = 'on'. To run active-set without this warning, use Algorithm = 'active-set'.
> In fmincon at 395
  In tarch at 149
  In multigarch at 157
  In ucsd_garch_demo at 196
??? Error using ==> fmincon at 504
FMINCON cannot continue because user supplied objective function failed with the following error:
Undefined function or method 'tarchcore' for input arguments of type 'double'.

Error in ==> tarch at 149
[parameters, LLF, EXITFLAG, OUTPUT, LAMBDA, GRAD] =  fmincon('tarchlikelihood', [omega ; alpha; tarch ; beta] ,sumA  , sumB ,[] , [] , LB , UB,[],options,data, p ,o, q, m, stdEstimate, type);

Error in ==> multigarch at 157
startingvalues(1:1+p+q+o)=tarch(data,p,o,q,'TARCH',[],newoptions);

Error in ==> ucsd_garch_demo at 196
[parametersAP, likelihoodAP, stderrorsAP, robustSEAP, htAP, scoresAP]=multigarch(errors,1,1,1,'APGARCH','NORMAL');

我想请教一下，数据出来以后每个factor的optimal weight怎么确定呀？谢谢。如果有任何关于DEA-DA的学习资源希望可以交流啊，我邮箱是eryanghuaguijingti@126.com