楼主: oliverwist
1770 2

[原创博文] 诚心向统计和SAS牛人请教:对一统计模型编程(最好编成宏),报酬可商 [推广有奖]

  • 0关注
  • 0粉丝

大专生

60%

还不是VIP/贵宾

-

威望
0
论坛币
523 个
通用积分
0.0000
学术水平
0 点
热心指数
0 点
信用等级
0 点
经验
1583 点
帖子
38
精华
0
在线时间
82 小时
注册时间
2009-1-5
最后登录
2025-5-29

楼主
oliverwist 发表于 2011-9-21 14:27:30 |AI写论文

+2 论坛币
k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群
赵安豆老师微信:zhaoandou666

经管之家联合CDA

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

感谢您参与论坛问题回答

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

+2 论坛币
各位大牛:  本人工作中遇到一模型,具体内容见附件文献中所示,有数据、有例子,已经有用R和STATA软件实现的,但我一直希望能用一套SAS宏搞定,无奈小弟才疏学浅,特向各位大牛有偿请教,经费都好商量。
希望:通过宏程序,设置不同的参数,控制模型,能输出所有参数估计结果和预测值,并能教会小弟如何使用。
模型简介:
被测(y)为个体主观报告值,为等级变量(1,2,3,4,5),其报告值由潜在变量(y*)和一套切点(t0,t1,。。。t5)值决定。y*,y以及切点值之间有如下关系:
y*=xβ+e, e~N(0,1),x为个体特征向量,β为带估计系数;
y=1,if t0<y*<=t1;
y=2,if t1<y*<=t2;
y=3,if t2<y*<=t3;   
y=4,if t3<y*<=t4;
y=5,if t4<y*<t5;
由于每个个体的切点不一样,因此有:
t1=xγ;
t(k)=t(k-1)+exp(xγ);,k=2,3,4,5  且   -∞=t0<t1<t2<t3<t4<t5=+∞
为了估计以上参数,引入若干情景,即让同一个体对虚拟个体的同一被测问题进行报告,因此,对于每一个虚拟个体,有(z和z*);
zj*=θj+ej, ej~N(0,σj2),θj为第j虚拟个体的真实被测,z为个体依据同一套切点(可理解为个人评价标准)对该虚拟个体的报告值,z和z*以及切点之间有如下关系;
z=1,if t0<z*<=t1;
z=2,if t1<z*<=t2;
z=3,if t2<z*<=t3;   
z=4,if t3<z*<=t4;
z=5,if t4<z*<t5;
通过构建联合似然函数对以上模型求解,获得关心的参数(y*,β,t1,t2,t3,t4,t5,γ等)
模型大致如此~~



小弟联系方式:QQ:10217081,邮箱:oliverwist@163.com
期待大牛的出现!
二维码

扫码加我 拉你入群

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

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

关键词:统计模型 计模型 统计和 stata软件 Stata 模型 统计 最好 程序 如何

沙发
oliverwist 发表于 2011-9-23 08:43:34
没人顶啊。。。
附上示例中的试验数据

effic1.rar

9.91 KB

本附件包括:

  • effic1.dta

藤椅
爱萌 发表于 2012-3-6 20:16:40
你解决了吗?如果没有可以找我
最恨对我说谎或欺骗我的人

您需要登录后才可以回帖 登录 | 我要注册

本版微信群
加好友,备注cda
拉您进交流群
GMT+8, 2025-12-5 17:12