| 所在主题: | |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| 文件名: 1王亚平,吴联生,白云霞,2005,中国上市公司盈余管理的频率与幅度,经济研究,2005年 12期.pdf | |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| 资料下载链接地址: https://bbs.pinggu.org/a-937452.html | |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| 附件大小: | |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
|
利用一组数据估计其分布函数参数时出错,请高手指教估计的正确方法和结果,谢谢!
本帖是https://bbs.pinggu.org/thread-1137305-1-1.html的后续工作。 原始数据:
>> clear >> mixedpdf=@(x,mu1,mu2,s1,s2,rho)(rho*normpdf(x,mu1,s1)+(1-rho)*normpdf(x,mu2,s2)) fnegpdf=@(x,k0,k1,lamda,mu1,mu2,s1,s2,rho)((1-k0*exp(k1*x))* mixedpdf(x,mu1,mu2,s1,s2,rho)) fpospdf=@(x,k0,k1,lamda,mu1,mu2,s1,s2,rho)(lamda*exp(-lamda*x)*(rho*k0* normcdf((-mu1-k1*s1^2)/s1)*exp(mu1*k1+s1^2*k1^2/2)+(1-rho)*k0*normcdf((-mu2-k1*s2^2)/s2)*exp(mu2*k1+s2^2*k1^2/2))+ mixedpdf(x,mu1,mu2,s1,s2,rho)) roepdf=@(x)fnegpdf.*(x<0)+fpospdf.*(x>=0) mixedpdf = @(x,mu1,mu2,s1,s2,rho)(rho*normpdf(x,mu1,s1)+(1-rho)*normpdf(x,mu2,s2)) fnegpdf = @(x,k0,k1,lamda,mu1,mu2,s1,s2,rho)((1-k0*exp(k1*x))* mixedpdf(x,mu1,mu2,s1,s2,rho)) fpospdf = @(x,k0,k1,lamda,mu1,mu2,s1,s2,rho)(lamda*exp(-lamda*x)*(rho*k0* normcdf((-mu1-k1*s1^2)/s1)*exp(mu1*k1+s1^2*k1^2/2)+(1-rho)*k0*normcdf((-mu2-k1*s2^2)/s2)*exp(mu2*k1+s2^2*k1^2/2))+ mixedpdf(x,mu1,mu2,s1,s2,rho)) roepdf = @(x)fnegpdf.*(x<0)+fpospdf.*(x>=0) >> [phat1,pci1]=mle(gap_roe,'pdf',roepdf,'start',[.9,25,120,-0.6315,0.0410,0.1135,0.0929,0.3012]) ??? Error using ==> stats\private\mlecustom>checkFunErrs The following error occurred while trying to evaluate the user-supplied pdf function '@(x)fnegpdf.*(x<0)+fpospdf.*(x>=0)': Error using ==> @(x)fnegpdf.*(x<0)+fpospdf.*(x>=0) Too many input arguments. Error in ==> stats\private\mlecustom at 159 checkFunErrs('pdf',pdfFun,start,uncensData,[],[],pdfAddArgs); Error in ==> mle at 219 [phat, pci] = mlecustom(data,varargin{:}); 附件1为主要的推导过程(其中α在实际估计时取为0),最终要估计(5)式中的参数。 附件2为参考文献:王亚平,吴联生,白云霞,2005,中国上市公司盈余管理的频率与幅度,经济研究,2005年 12期 所有公式推导均引自参考文献第三部分。 |
|||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
熟悉论坛请点击新手指南
|
|||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| 下载说明 | |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
|
1、论坛支持迅雷和网际快车等p2p多线程软件下载,请在上面选择下载通道单击右健下载即可。 2、论坛会定期自动批量更新下载地址,所以请不要浪费时间盗链论坛资源,盗链地址会很快失效。 3、本站为非盈利性质的学术交流网站,鼓励和保护原创作品,拒绝未经版权人许可的上传行为。本站如接到版权人发出的合格侵权通知,将积极的采取必要措施;同时,本站也将在技术手段和能力范围内,履行版权保护的注意义务。 (如有侵权,欢迎举报) |
|||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
京ICP备16021002号-2 京B2-20170662号
京公网安备 11010802022788号
论坛法律顾问:王进律师
知识产权保护声明
免责及隐私声明