关于MCMCglmm随机变量设置的询问-经管之家官网!

人大经济论坛-经管之家 收藏本站
您当前的位置> 会计>>

会计库

>>

关于MCMCglmm随机变量设置的询问

关于MCMCglmm随机变量设置的询问

发布:xingzhaoh | 分类:会计库

关于本站

人大经济论坛-经管之家:分享大学、考研、论文、会计、留学、数据、经济学、金融学、管理学、统计学、博弈论、统计年鉴、行业分析包括等相关资源。
经管之家是国内活跃的在线教育咨询平台!

经管之家新媒体交易平台

提供"微信号、微博、抖音、快手、头条、小红书、百家号、企鹅号、UC号、一点资讯"等虚拟账号交易,真正实现买卖双方的共赢。【请点击这里访问】

提供微信号、微博、抖音、快手、头条、小红书、百家号、企鹅号、UC号、一点资讯等虚拟账号交易,真正实现买卖双方的共赢。【请点击这里访问】

小弟最近学习CMCglmm,在设置随机变量时遇到问题请打大侠帮忙,数据在附件中程序如下:library(coda)library(MCMCpack)library(MCMCglmm)rd<-read.csv("11.csv")head(rd)attach(rd)mu.beta<-c(-2.037,2.415) ...
免费学术公开课,扫码加入


小弟最近学习CMCglmm, 在设置随机变量时遇到问题请打大侠帮忙,数据在附件中
程序如下:
library(coda)
library(MCMCpack)
library(MCMCglmm)
rd<-read.csv("11.csv")
head(rd)
attach(rd)
mu.beta<-c(-2.037,2.415)
V.betas<-matrix(c(0.35716566, -0.08127909,-0.08127909 , 0.03360510),nrow=2)
n.iter<-25000
prior.betas <- list(B = list(mu = mu.beta, V = V.betas),
R = list(V = 1, nu = 0.001))
model1.sim <- MCMCglmm(TOTAL~D, random=NULL, rcov=~units, family="gaussian",
data=rd, prior=prior.betas ,nitt=n.iter,burnin=floor(n.iter/2),
thin=3,verbose=F)
#在不设置随机变量时,model1.simk可以做出来
model2.sim <- MCMCglmm(TOTAL~D, random = ~ region,, rcov=~units, family="gaussian",
data=rd, prior=prior.betas ,nitt=n.iter,burnin=floor(n.iter/2),
thin=3,verbose=F)
#在设置随机变量时,出现:
Error in MCMCglmm(TOTAL ~ D, random = ~D | region, , rcov = ~units,:
either both or neither R and G structures need a prior
哪位大侠指点一下,感激不尽,着急呀!
「经管之家」APP:经管人学习、答疑、交友,就上经管之家!
免流量费下载资料----在经管之家app可以下载论坛上的所有资源,并且不额外收取下载高峰期的论坛币。
涵盖所有经管领域的优秀内容----覆盖经济、管理、金融投资、计量统计、数据分析、国贸、财会等专业的学习宝库,各类资料应有尽有。
来自五湖四海的经管达人----已经有上千万的经管人来到这里,你可以找到任何学科方向、有共同话题的朋友。
经管之家(原人大经济论坛),跨越高校的围墙,带你走进经管知识的新世界。
扫描下方二维码下载并注册APP
本文关键词:

本文论坛网址:https://bbs.pinggu.org/thread-2974495-1-1.html

人气文章

1.凡人大经济论坛-经管之家转载的文章,均出自其它媒体或其他官网介绍,目的在于传递更多的信息,并不代表本站赞同其观点和其真实性负责;
2.转载的文章仅代表原创作者观点,与本站无关。其原创性以及文中陈述文字和内容未经本站证实,本站对该文以及其中全部或者部分内容、文字的真实性、完整性、及时性,不作出任何保证或承若;
3.如本站转载稿涉及版权等问题,请作者及时联系本站,我们会及时处理。
经管之家 人大经济论坛 大学 专业 手机版