全部版块 我的主页
› 论坛 › 数据科学与人工智能 › 数据分析与数据科学 › R语言论坛
R2winbugs到出数据时出现错误
文雅13
2012-12-17 16:20:52
[code] we might assume that F{θ; xi (t)} follows a marginalized transition
model (MTM) ， with marginal mean g{μ (t)} = xi (t)β,文献中就是这样提了下，具体程序中就是那样了
二维码

扫码加我 拉你入群

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

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

epoh
2012-12-17 16:27:41
文雅13 发表于 2012-12-17 16:11
是在文献中找到的程序，程序就是那样的
那你就向文献的作者要要看
看他是否愿意提供
否则你就要依需求,自行编程

# MTM1Delta is a function written in the WinBUGS development         
# interface for computing the intercept delta[i,j] to satisfy
# marginal mean and serial dependence simultaneously

# Mv is a function written in the WinBUGS development
# interface for computing smooth function estimates
# with low-rank thin-plate bases
二维码

扫码加我 拉你入群

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

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

文雅13
2013-1-17 09:09:13
epoh 发表于 2012-12-17 16:27
那你就向文献的作者要要看
看他是否愿意提供
否则你就要依需求,自行编程
老师，在检验Winbugs模型时出现loop name already used in outer loop是什么问题？
二维码

扫码加我 拉你入群

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

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

epoh
2013-1-17 15:36:26
文雅13 发表于 2013-1-17 09:09
老师，在检验Winbugs模型时出现loop name already used in outer loop是什么问题？
for (i in 1:n) // outer loop
{
    for (i in 1:p) // inner loop
    {
    //...do something
    }
//...do something else
}

loop name already used in outer loop
二维码

扫码加我 拉你入群

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

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

文雅13
2013-1-17 16:05:54
epoh 发表于 2013-1-17 15:36
for (i in 1:n) // outer loop
{
    for (i in 1:p) // inner loop
老师，我没太看明白这些code是什么意思？
二维码

扫码加我 拉你入群

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

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

文雅13
2013-1-17 16:10:32
就是for (i in 1:n) 地方有错误
二维码

扫码加我 拉你入群

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

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

文雅13
2013-1-17 19:46:33
epoh 发表于 2013-1-17 15:36
for (i in 1:n) // outer loop
{
    for (i in 1:p) // inner loop
老师，没明白您的意思，能说的详细点吗？
二维码

扫码加我 拉你入群

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

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

epoh
2013-1-17 20:42:27
文雅13 发表于 2013-1-17 19:46
老师，没明白您的意思，能说的详细点吗？
也就是outer loop 用i,inner loop 就要用j,不能用相同的i
二维码

扫码加我 拉你入群

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

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

文雅13
2013-1-17 21:05:43
epoh 发表于 2013-1-17 20:42
也就是outer loop 用i,inner loop 就要用j,不能用相同的i
用的就是i ，但鼠标还指着这里面的ifor (i in 1:n)，在Winbugs中 gamma as a matrix (two indices)and a
vector (one index)
是指什么？
二维码

扫码加我 拉你入群

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

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

文雅13
2013-1-17 21:10:25
epoh 发表于 2013-1-17 20:42
也就是outer loop 用i,inner loop 就要用j,不能用相同的i
复制代码
比如这个code。gamman有矢量也有矩阵吗？不会看，哎
二维码

扫码加我 拉你入群

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

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

epoh
2013-1-17 21:38:23
文雅13 发表于 2013-1-17 21:05
用的就是i ，但鼠标还指着这里面的ifor (i in 1:n)，在Winbugs中 gamma as a matrix (two indices)and a
...
gamma是自己定义的
由serial dependence cor赋予gamma[i,j]

# serial dependence
cor<-(a3[1]+a3[2]*U+inprod2(d3[1:9],BZ[i,1:9]))*step(11-nobs)+corrlast*equals(nobs,12)

gamma[i,j]<-cor
二维码

扫码加我 拉你入群

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

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

文雅13
2013-1-17 21:43:31
epoh 发表于 2013-1-17 21:38
gamma是自己定义的
由serial dependence cor赋予gamma
我刚才发的函数代替了之前的M*****elta，但是老出现outer loop的错误，然后他说我的gamma出现了（two indices）和（one index），我不明白应该怎样改刚才的那个code才能避免outer loop的错误
二维码

扫码加我 拉你入群

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

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

epoh
2013-1-17 22:09:41
文雅13 发表于 2013-1-17 21:43
我刚才发的函数代替了之前的M*****elta，但是老出现outer loop的错误，然后他说我的gamma出现了（two ind ...
问题绝不仅你说的这些而已,
最快速的办法就是把code 上传.
二维码

扫码加我 拉你入群

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

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

文雅13
2013-1-17 22:13:47
epoh 发表于 2013-1-17 22:09
问题绝不仅你说的这些而已,
最快速的办法就是把code 上传.
就是之前我给你的那个code，然后不是需要MTM这个函数么，然后作者告我说换成刚发的那个code，只是计算机会慢，我换了之后就出现这样的错误了。不知道怎么回事
二维码

扫码加我 拉你入群

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

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

epoh
2013-1-17 22:53:24
文雅13 发表于 2013-1-17 22:13
就是之前我给你的那个code，然后不是需要MTM这个函数么，然后作者告我说换成刚发的那个code，只是计算机会 ...
那你就要自己给出gamma[,]
二维码

扫码加我 拉你入群

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

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

文雅13
2013-1-17 22:55:16
epoh 发表于 2013-1-17 22:53
那你就要自己给出gamma[,]
我说出现这样的错误，作者回复我The code should work to get Delta with proper adjustment。我没有思路，怎样去改变Delta，怎样加gamma[,]
二维码

扫码加我 拉你入群

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

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

epoh
2013-1-17 23:10:51
文雅13 发表于 2013-1-17 22:55
我说出现这样的错误，作者回复我The code should work to get Delta with proper adjustment。我没有思路 ...
哈哈,我已说过
我不清楚你的模型
这要靠你自己了
二维码

扫码加我 拉你入群

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

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

文雅13
2013-1-19 10:29:42
epoh 发表于 2013-1-17 23:10
哈哈,我已说过
我不清楚你的模型
这要靠你自己了
老师，我想问问在Winbugs中加载数据时出现expected key word structure。应该怎样改数据格式
二维码

扫码加我 拉你入群

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

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

epoh
2013-1-19 11:44:58
文雅13 发表于 2013-1-19 10:29
老师，我想问问在Winbugs中加载数据时出现expected key word structure。应该怎样改数据格式
底下帖子跟你发生相同情况,
expected key word structure.

  https://stat.ethz.ch/pipermail/r-help/2005-September/079481.html

答案依然一 样
without data and model file I am unable to help for the data part.
二维码

扫码加我 拉你入群

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

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

文雅13
2013-1-19 14:08:21
epoh 发表于 2013-1-19 11:44
底下帖子跟你发生相同情况,
expected key word structure.
数据有问题，希望您能帮忙
附件列表

Winbugs.rar

大小:2.97 KB

 马上下载

二维码

扫码加我 拉你入群

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

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

epoh
2013-1-19 16:13:42
文雅13 发表于 2013-1-19 14:08
数据有问题，希望您能帮忙
这还是少了个function Mv()
你没一并根作者要吗?
# Mv is a function written in the WinBUGS development
# interface for computing smooth function estimates
# with low-rank thin-plate bases
betasex<-Mv(a1[1:2], d1[1:9], knotD[],OMEGAu[],v)
二维码

扫码加我 拉你入群

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

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

文雅13
2013-1-19 16:36:17
程序是正确的，就是data loading时有问题
附件列表

新建文件夹.rar

大小:2.07 KB

 马上下载

二维码

扫码加我 拉你入群

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

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

epoh
2013-1-19 20:48:19
文雅13 发表于 2013-1-19 16:36
程序是正确的，就是data loading时有问题
应该少一个.odc file
类似底下ObservedPlus.odc
定义出function

(*1*) MODULE WBDevObservedPlus;
IMPORT
WBDevVector,
Math;
TYPE
Function = POINTER TO RECORD (WBDevVector.Node) END;
Factory = POINTER TO RECORD (WBDevVector.Factory)END;
VAR
fact-: WBDevVector.Factory;
PROCEDURE (func: Function) DeclareArgTypes
(OUT args: ARRAY OF CHAR);
BEGIN
(*2*) args := "ss";
END DeclareArgTypes;
PROCEDURE (func: Function) Evaluate
(OUT values: ARRAY OF REAL);
CONST
(*3*) in = 0; ik = 1;
VAR
(*4*) n,
k: INTEGER;

BEGIN
(*5*) n := SHORT(ENTIER(func.arguments[in][0].Value()));
k := SHORT(ENTIER(func.arguments[ik][0].Value()));
(*6*) n:=n+10;
k:=k+1;
values[0] := n;
values[1] := k;
END Evaluate;
PROCEDURE (f: Factory) New (option: INTEGER): Function;
VAR
func: Function;
BEGIN
NEW(func); func.Initialize; RETURN func;
END New;
PROCEDURE Install*;
BEGIN
WBDevVector.Install(fact);
END Install;
PROCEDURE Init;
VAR
f: Factory;
BEGIN
NEW(f); fact := f;
END Init;
BEGIN
Init;
(*7*) END WBDevObservedPlus.
二维码

扫码加我 拉你入群

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

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

文雅13
2013-1-19 21:06:17
epoh 发表于 2013-1-19 20:48
应该少一个.odc file
类似底下ObservedPlus.odc
定义出function
有点糊涂了。意思是还缺个函数？
二维码

扫码加我 拉你入群

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

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

epoh
2013-1-19 21:22:54
文雅13 发表于 2013-1-19 21:06
有点糊涂了。意思是还缺个函数？
嗯,没错
二维码

扫码加我 拉你入群

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

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

文雅13
2013-1-19 21:25:46
epoh 发表于 2013-1-19 21:22
嗯,没错
什么函数
二维码

扫码加我 拉你入群

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

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

epoh
2013-1-19 21:55:22
文雅13 发表于 2013-1-19 21:25
什么函数
你自己看看就知道了
  http://bugsmodellibrary.googleco ... ckBoxWinBUGS/WBDev/
二维码

扫码加我 拉你入群

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

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

文雅13
2013-1-19 22:01:14
epoh 发表于 2013-1-19 21:55
你自己看看就知道了
  http://bugsmodellibrary.googlecode.com/svn-history/r2/trunk/BUGSModelLibrary ...
老师你直接告我是个什么函数，我找找看看有没有。
二维码

扫码加我 拉你入群

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

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

epoh
2013-1-19 23:06:40
文雅13 发表于 2013-1-19 22:01
老师你直接告我是个什么函数，我找找看看有没有。
函数主要写在Mv.odc
二维码

扫码加我 拉你入群

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

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

epoh
2013-1-20 15:55:09
文雅13 发表于 2013-1-19 22:01
老师你直接告我是个什么函数，我找找看看有没有。
就以Minmax而言,有三个文件
  code file   : Minmax.ocf
  document  : Minmax.odc (最重要)
  symbol file : Minmax.osf
你没有给出Mv.odc,你自己可以比对一下
####
Minmax.ocf

MODULE WBDevMinmax
processor: 80x86
meta size: 432
desc size: 492
code size: 356
data size: 4
compiled: 2013/1/20  ¤U¤&Egrave; 12:44:15
5 imports:
        WBDevVector         Factory, Factory^, Install, Node, Node^
        GraphNodes         Node^
        GraphVector        
        GraphLogical        
        Math        

####
Minmax.odc
(*1*)        MODULE WBDevMinmax;

                IMPORT
                        WBDevVector,
(*2*)                        Math;
                        
                TYPE
                        Function = POINTER TO RECORD (WBDevVector.Node) END;
                        Factory = POINTER TO RECORD (WBDevVector.Factory) END;
                        
                VAR
                        fact-: WBDevVector.Factory;
                        
(*3*)                PROCEDURE (func: Function) DeclareArgTypes (OUT args: ARRAY OF CHAR);
(*4*)                BEGIN
(*5*)                        args := "v";
(*6*)                END DeclareArgTypes;
               
(*7*)                PROCEDURE (func: Function) Evaluate (OUT values: ARRAY OF REAL);
(*10*)                VAR
(*11*)                        len, i: INTEGER;
                        min, max: REAL;
(*13*)                BEGIN
(*14*)                        len := LEN(func.arguments[0]);
                        min := func.arguments[0][0].Value();
                        max := func.arguments[0][0].Value();
(*19*)                        i := 1;
(*20*)                        WHILE i < len DO;
                                min := MIN(min, func.arguments[0].Value());
                                max := MAX(max, func.arguments[0].Value());
(*31*)                                INC(i);
(*32*)                        END;
                        values[0] := min; values[1] := max
(*33*)                END Evaluate;
               
                PROCEDURE (f: Factory) New (option: INTEGER): Function;
                VAR
                        func: Function;
                BEGIN
                        NEW(func); func.Initialize; RETURN func;
                END New;
               
                PROCEDURE Install*;
                BEGIN
                        WBDevVector.Install(fact);
                END Install;
               
                PROCEDURE Init;
                VAR
                        f: Factory;
                BEGIN
                        NEW(f); fact := f;
                END Init;
               
        BEGIN
                Init;
(*1*)        END WBDevMinmax.
####
Minmax.osf

DEFINITION WBDevMv;
        IMPORT WBDevScalar;
        VAR
                fact-: WBDevScalar.Factory;

        PROCEDURE Install;
END WBDevMv.
二维码

扫码加我 拉你入群

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

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

首页上一页下一页跳至第页
栏目导航
热门文章
推荐文章

分享

微信QQ空间QQ微博
扫码加好友，拉您进群
各岗位、行业、专业交流群