关于本站
人大经济论坛-经管之家:分享大学、考研、论文、会计、留学、数据、经济学、金融学、管理学、统计学、博弈论、统计年鉴、行业分析包括等相关资源。
经管之家是国内活跃的在线教育咨询平台!
经管之家新媒体交易平台
提供"微信号、微博、抖音、快手、头条、小红书、百家号、企鹅号、UC号、一点资讯"等虚拟账号交易,真正实现买卖双方的共赢。【请点击这里访问】
期刊
- 期刊库 | 马上cssci就要更新 ...
- 期刊库 | 【独家发布】《财 ...
- 期刊库 | 【独家发布】“我 ...
- 期刊库 | 【独家发布】“我 ...
- 期刊库 | 【独家发布】国家 ...
- 期刊库 | 请问Management S ...
- 期刊库 | 英文期刊库
- 核心期刊 | 歧路彷徨:核心期 ...
TOP热门关键词
坛友互助群 |
扫码加入各岗位、行业、专业交流群 |
我有3样本,每个样本测量了11次,Area ,Cu ,Wa ,Ur ,CuP ,WaP几个参数对于同一样本都是相同的。我希望得到各个样本的beta.m,beta.s1-6,alpha。代码如下:
#Model
model
{
for(i in 1:M) #M sites
{
#big model
y ~ dnorm(u,o)
u <- beta.m+alpha[i,j]
#a[i,j] model
for (j in 1:N) #number of surveysat site i
{
alpha[i,j] ~ dnorm(us[i,j],oa)
}
#priors
oa ~ dgamma(0.01,0.01)
beta.s1 ~ dnorm(0.0,1.0e-10)
beta.s2 ~ dnorm(0.0,1.0e-10)
beta.s3 ~ dnorm(0.0,1.0e-10)
beta.s4 ~ dnorm(0.0,1.0e-10)
beta.s5 ~ dnorm(0.0,1.0e-10)
beta.s6 ~ dnorm(0.0,1.0e-10)
}
#priors
o ~ dgamma(0.01,0.01)
beta.m ~ dnorm(0.0,1.0e-10)
#data
list(M=3,N=11,
Area=c(34.99873323,10.14005246,18.62917614),Cu=c(0.920451749,0.647475224,0.934651939),
Wa=c(0.00363696,0.000377537,0.03375289),Ur=c(0.075911291,0.24898537,0.031595171),
CuP=c(0.917631044,1.190760085,0.803192534),
WaP=c(1.271445528,1.247229591,0.885669015,
us=stracture(
.data=c(5.96,4,4.2,3.47,3.93,4.78,4.93,4.07,3.83,3.91,4.29,5.01,3.29,3.8,3.92,3.1,3.68,3.55,5.63,3.38,3.96,4.04,4.76,4.07,5.58,4.02,3.96,4.2,3.86,3.12,3.56,3.67,5.93,3.37,4.1,3.73,3.71,4.22)
.dim=c(3,11))
)
#initial value
list(beta.m=0,beta.s1=0,beta.s2=0,beta.s3=0,beta.s4=0,beta.s5=0,beta.s6=0,o=1,oa=1,alpha[i,j]=0)
us[i,j] <- beta.s1*Area + beta.s2*Cu + beta.s3*Wa + beta.s4*Ur+ beta.s5*CuP + beta.s6*WaP
免流量费下载资料----在经管之家app可以下载论坛上的所有资源,并且不额外收取下载高峰期的论坛币。
涵盖所有经管领域的优秀内容----覆盖经济、管理、金融投资、计量统计、数据分析、国贸、财会等专业的学习宝库,各类资料应有尽有。
来自五湖四海的经管达人----已经有上千万的经管人来到这里,你可以找到任何学科方向、有共同话题的朋友。
经管之家(原人大经济论坛),跨越高校的围墙,带你走进经管知识的新世界。
扫描下方二维码下载并注册APP
您可能感兴趣的文章
本站推荐的文章
人气文章
本文标题:check model出现expected variable name,是什么问题呀?
本文链接网址:https://bbs.pinggu.org/jg/kaoyankaobo_kaoyan_3709647_1.html
2.转载的文章仅代表原创作者观点,与本站无关。其原创性以及文中陈述文字和内容未经本站证实,本站对该文以及其中全部或者部分内容、文字的真实性、完整性、及时性,不作出任何保证或承若;
3.如本站转载稿涉及版权等问题,请作者及时联系本站,我们会及时处理。