楼主: jiachunyang1988
31825 56

[问答] 重复运行一个r程序多次 [推广有奖]

31
曲散人终 发表于 2014-12-2 17:31:58
你可以用source(“C.r”)->x
x
来查看一下x里面有什么。。。然后再做一个循环。。。

32
jiachunyang1988 发表于 2014-12-2 19:15:47
什么也没有啊,是不是我前面的程序写的不多,我运行了什么也没输出来

33
曲散人终 发表于 2014-12-2 19:26:10
你参数估计最后不是有一个值吗?难道没有吗?

34
曲散人终 发表于 2014-12-2 19:33:09
test.r
内容如下:
i<-rnom(1)
x<-matrix(1,30,1)
for(j in 1:30)
{
x[j,1]<-source("test.r")$value"
}
x
x
              [,1]
[1,] -0.149700272
[2,] -0.120855191
[3,] -1.270794900
[4,] -0.202314091
[5,]  1.147683263
[6,] -0.319944397
[7,]  0.051101260
[8,]  0.046986890
[9,] -0.659375530
[10,]  0.933372978
[11,] -1.746467248
[12,]  0.806004986
[13,] -0.579561031
[14,]  1.420858086
[15,]  1.600736200
[16,] -1.267380714
[17,] -1.281226154
[18,]  0.768222495
[19,] -1.062440255
[20,] -0.421094436
[21,]  0.247398030
[22,]  0.146147936
[23,]  0.160183127
[24,] -0.618378692
[25,]  0.757189616
[26,]  0.594755607
[27,]  0.424898135
[28,]  0.400405271
[29,] -0.002198888
[30,]  0.054377998
可以参照一下这个。。。

35
曲散人终 发表于 2014-12-2 19:34:53
test.r的内容就是i<-rnorm(1),产生随机数。。。
然后再循环。。。你看下你的程序吧。。。最后应该有一个值的。。。

36
jiachunyang1988 发表于 2014-12-2 19:56:47
曲散人终 发表于 2014-12-2 19:26
你参数估计最后不是有一个值吗?难道没有吗?
是四个值,可是我就是这么弄的怎么没有呢,你能帮我看下吗?我把程序给你,麻烦你帮我看下前面的该怎么写,谢谢

37
曲散人终 发表于 2014-12-2 20:04:13
程序在哪里。。。

38
jiachunyang1988 发表于 2014-12-2 20:33:23
jiachunyang1988 发表于 2014-12-2 19:56
是四个值,可是我就是这么弄的怎么没有呢,你能帮我看下吗?我把程序给你,麻烦你帮我看下前面的该怎么写 ...
我晓得你的意思了,可是我的程序是分为两个部分,第一部分估计得到一个参数的值,在得到第一个估计值的基础上然后第二部分得到另外三个参数的估计值,我刚刚试了下,输出的结果是一个空矩阵,这是什么原因呢

39
曲散人终 发表于 2014-12-2 20:50:54
test.r
内容为
i<-rnorm(1)
j<-rnorm(1)
k<-rnorm(1)
l<-rnorm(1)
x<-c(i,j,k,l)
data<-matrix(30,30,4)
for(fi in 1:30)
{
data[i,]<-source("test.r")
}
data
data
               [,1]        [,2]        [,3]        [,4]
[1,]  0.4689560943  0.41021918 -1.47245433 -0.83379813
[2,] -0.1468344631 -0.04029691  2.15698041 -0.29293431
[3,]  0.6318374815  0.86381940  0.08801037 -1.46665046
[4,] -0.0687576799  0.35173895 -0.51532110 -0.55512853
[5,] -0.7810682770  0.17337073 -1.38850202 -0.97906878
[6,]  0.7187910668 -0.89979922  0.84536582 -0.32634659
[7,] -0.8897264650 -0.91103386 -0.39716780  0.02785898
[8,]  1.8046113412  0.16164153  0.75080871 -0.77221783
[9,]  1.2794622372  0.68440660  0.01725624  0.14091887
[10,]  0.1054056792  0.68079928  1.58352235 -0.72553600
[11,]  1.9858083714 -0.65227872 -0.86700426  1.64785058
[12,] -1.2304023646 -0.74228038 -0.84655171  0.13221950
[13,] -0.0753357468 -1.17368119 -0.16568548 -1.36263117
[14,]  0.4442810319  1.95460980 -0.07832366 -0.86417766
[15,]  0.5196975400 -0.21270824  1.30647029  1.19599369
[16,] -0.4863672203  0.87120992 -0.65673603 -1.15469391
[17,] -0.9540178423 -0.55430719  1.08494383  0.03339504
[18,]  1.9429974087  0.40530120  0.16319767  2.81664497
[19,] -1.7250451480 -0.73982729  0.30675100  2.38030097
[20,]  0.0130614280 -0.94566631 -0.13124599  1.17081426
[21,]  0.6761404384 -0.82220129 -1.59183754  0.58256108
[22,] -1.0370987660 -0.30434138  1.67779389  0.33131806
[23,]  0.2869559320  0.16351840 -0.45566760  0.43295978
[24,] -0.0005876634 -1.13741380 -0.85608562  0.27466919
[25,] -1.1781133349 -1.31198179  1.84919456  0.62323588
[26,] -0.8381296321 -0.14322947  0.14662975 -0.55472101
[27,] -0.0176577131  0.61830364  0.95733544  0.06823911
[28,]  0.2574328924 -1.40240246  0.37088190 -2.18087317
[29,] -0.0434359410  0.18470725 -1.01401217  0.12076647
[30,]  0.1793901250  0.67198422  2.37685099 -1.13329723
这样可以了吧。。。

40
曲散人终 发表于 2014-12-2 20:51:49
还是要加上$value在source("test.r")后面。。。

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

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