Question

I m st 参照Bolker(

我在执行一揽子计划时有一些问题。在许多对人的分析过程中出现了错误,因此,例子中的图表没有出现。您能否解决这个问题?

i 想作此图,但

我的法典是这样写的。

这是我的法典。

Z = simmixstock2(nsource=4,nmark=5,nmix=3,
          sourcesize=c(4,2,1,1),
          sourcesampsize=rep(25,4),
          mixsampsize=rep(30,3),rseed=1001)
Zfit0 = mm.wbugs(Z,sourcesize=c(4,2,1,1),returntype="bugs")
plot(as.mcmc.bugs(Zfit0))
plot(Zfit0)

此外,我能否把除随机数据外的我自己的数据输入“模拟储存2”的指挥? 我想用我的数据作一番,但人工指挥似乎只是随意生成,因此,我想知道,是否有办法。谢谢。

Answer 1

对我来说,这项工作是:

library(mixstock)
Z = simmixstock2(nsource=4,nmark=5,nmix=3,
          sourcesize=c(4,2,1,1),
          sourcesampsize=rep(25,4),
          mixsampsize=rep(30,3),rseed=1001)
Zfit0 = mm.wbugs(Z,sourcesize=c(4,2,1,1), pkg="JAGS", jags.seed = 101)
plot(Zfit0)

变化

use JAGS instead of WinBUGS (which may be nearly obsolete - in any case, hard to run on platforms other than Windows; I don t know the status of OpenBUGS), including setting the seed for the MCMC run
don t specify returntype = "bugs" (instead, the default "mixstock" option gives something that we can plot)

我不敢肯定你所说的话。

我能否把除随机数据外的我自己的数据输入“模拟存货2”的指挥?

simmixstock2(> function is Designed for simulating data in the case where Youoshi to see how what work when You teurs no own data. The vignette/walk through on pp. 4-5 如何输入多种混合储量的数据(在矩阵中增加一栏)(加列= ha类、一栏=来源和混合储量),并在使用<代码>时具体说明混合储量的数量(<代码>nmix=1(因缺)。

Answer 2

谢谢。

我解决了第一个问题。在你的手册中打上正确的图表。

但是,我仍然用我的数据 st在第二个问题上。

mydata = read.csv("testfile1.csv")
head(mydata)
    YYM OKN AMM OGS SEA FSM CEP JP
1.1   0   0   0   0   0   0 146   0
2.1   0   0   0   0   0   0  29   0
3.1   0   0   0   0   0   0  21   0
3.2   0   0   0   0   0   0  15   0
4.1   0   0   0   0   0   0  13   0
4.2   0   0   0   0   0   0   2   0

mydata = as.mixstock.data(mydata)
mydata0 = mm.wbugs(mydata,sourcesize=c(1,1.8,3.7,4.1,42,51.1,21.5), pkg="JAGS", jags.seed = 101)
**Error in 1:MIX : argument of length 0**

我认为,我不理解消息来源,也不了解其他东西。是否有解决办法?

感谢您的帮助。 a. 有一天

友情链接