我需要帮助来加速一个简单的函数，该函数使用which()和rbinom()来根据每天的生存概率和嵌套周期来计算一个嵌套存活的时间。我在一个闪亮的应用程序的data.table模拟中使用了这一点，这条线真的，真的很慢。
违规函数如下-它计算给定每日生存概率<e
浏览 8提问于2020-11-14得票数 3
回答已采纳
1回答
我是一个全新的R用户，试图将具有指定依赖关系的大量数据生成到单个data.frame中。我当前的代码使用了for循环，并且运行良好(N=250 ids，超过10次复制)，但我必须将其扩展为更大的N和K复制。我尝试过使用foreach和函数进行并行处理，但没有太多成功。我无法将复制的数据添加到最终帧中。
浏览 0提问于2020-08-25得票数 0
我要做大量的模拟，这需要很多时间。我认为通过data.table可以缩短处理时间。如何将mdply(data.frame(prob=seq(from = 0.1, to = 0.9, by = 0.1)), rbinom, n = 5, size = 2)的结果存储到data.table中，而不首先将其输出保存到data.frame。= 2)dt1 <- data.table(df1)
浏览 4提问于2015-07-31得票数 0
回答已采纳
我在R中使用一个大约有6e6行的data.table，并创建了一个函数，我通过data.table传递该函数，以基于两个分组值创建一个新列。从技术上讲，我的函数遍历分组参数的每一行并执行一些非常简单的代数操作，但是考虑到我的data.table的大小，这将需要相当长的时间。我熟悉foreach()函数<
浏览 15提问于2019-09-14得票数 0
1回答
我使用replicate来模拟R中的发行版，并可视化它们在不同参数下的变化(例如，rbinom(100,1,0.5)与rbinom(100,1,0.01))。我想在一个函数中完成所有这一切，1.模拟复制，2.设置绘图尺寸和参数，3.循环通过和绘制密度曲线。在单独的部分中，这段代码没有问题：
n <- 10
浏览 3提问于2020-05-08得票数 1
回答已采纳
1回答
让L=病人健康的年数。问题是在R中编写一个函数，模拟单个病人的轨迹，并返回病人健康的年数。我认为这可以用rbinom函数作为二项分布在R<em
浏览 5提问于2020-02-18得票数 0
回答已采纳
1回答
有人知道如何在R中实现这个算法吗？我总共有100名患者需要随机分配，我想知道超过55%-45%的概率是多少(即45名或更少的受试者将被分配其中一种治疗，55名或更多的患者被分配到另一种治疗中)？假设处理A和B的应答率分别为0.6和0.5。我需要基于Zelen的Play The Winner Rule的1000个模拟来做到这一点。泽伦规则背后<em
浏览 29提问于2019-10-15得票数 1
回答已采纳
1回答
我想用一个包含rbinom函数的for循环来模拟一个矩阵。这个循环执行rbinom函数100次，因此每次运行都会产生不同的结果。然而，我无法找到一种方法，使这些结果在一个矩阵，以进一步分析。当将for循环分配给对象时，此对象在环境中显示为空，因此不能在矩阵中使用。(“数据”必须为向量类型，为“NULL”)。当不包括for循环中的rbinom-函数时，它可以分配给一
浏览 1提问于2019-04-18得票数 1
回答已采纳
我想在data.table中为组g中的每一行绘制Y。为了说明这一点，请使用这个简单的模拟： set.seed(123)
g = 3
group_id= sort(rep(1:g, N)), )
DT[, Y
浏览 23提问于2020-09-05得票数 0
1回答
试图围绕向量化，试图使一些模拟更快，我发现了这个非常基本的流行病模拟。respectively I <- rep(0, T+1) S[1] <- N R[1]<- 0 S[i+1] <- rbinom(1, S[i], (1 - a)^I[i])
R[i+1] &l
浏览 2提问于2012-05-27得票数 2
回答已采纳
3回答
模拟掷硬币一周？
、、、
我有兴趣建立一个模拟投币在R.我想运行一个星期的模拟。在R中是否有一个函数，允许我在一周的时间内开始和停止模拟？如果一切顺利，我可能想增加模拟周期的长度。例如：因此，为了澄清，而不是上面代码中的10，我如何使模拟持续一周(一周内的</em
浏览 8提问于2012-06-22得票数 4
回答已采纳
1回答
我假设了不相等的人口规模，例如这里我有100个模拟测试，其中5,000个在测试群体中，10000个在控制中，转换率的实际差异(测试中的0.05%和控制中的0.04%)。现在，我想对这些模拟测试中的每一个
浏览 0提问于2016-03-15得票数 0
假设我有一个数据框架，如：df <- data.frame(x=rbinom(100,1,0.9), y=rbinom(100,1,0.95))## index <- sample(1:nrow(df),..df_sub <- df[index,] x y2
浏览 9提问于2022-08-02得票数 1
1回答
在文档中，rbinom有三个参数：n、size和prob。对this question的回答是，如果prob是一个向量，它将被循环使用，直到到达n。我要找到那段代码！R源代码立即调用一个C函数： > rbinom.Call(C_rbinom, n, size, prob)<environment: names
浏览 13提问于2019-04-02得票数 1
回答已采纳
1回答
提高数据帧仿真性能
、、、
不幸的是，有时候知道什么是最有效的代码是至关重要的-我想学习！
目前，我正在模拟多个数据帧合并成一个列表。在模拟之后，我需要第二个数据帧，包括整个列表中所有列的平均值和SD。我的代码(参见下面的示例)完美地产生了预期的结果，但它似乎是第一种，有点慢(我指的是真实的几个小时)，第二，很高的RAM消耗(为此，我暂时减少了列表中模拟dfs的数量)。对于
浏览 2提问于2017-09-29得票数 0
回答已采纳
2回答
我试图用data.table或某种形式的并行化来加速离散时间非均匀马尔可夫链的蒙特卡罗模拟。使用随机虚拟转移矩阵TM，我正在模拟N个模拟中的nSteps时间步骤，并从初始状态向量initialState开始记录currentState中的下一个更新状态。如果我把N-循环的主体放在一个函数中
statefun &
浏览 0提问于2015-05-22得票数 1
回答已采纳
1回答
当启动ggvis时，我会收到这样的消息：:=对于运行data.table是必不可少的，而data.table恰好是dplyr的替代方案。我的代码仍然正确运行，但是启动消息让我担心，因为我在所有脚本中都使用了data.table。
考虑到ggvis的启动消息，在data.table</e
浏览 2提问于2014-07-16得票数 4
回答已采纳
我的data.tablerequire(data.table) col4 = rnorm(n=10, mean=20, sd=5), col5= runif(10, 0, 1), col6 = rbinom(10, 3, .
浏览 0提问于2018-04-20得票数 3
回答已采纳
1回答
下面的两个R函数来自"Flexible Imputation of Missing Data“一书(第59和63页)。第一个生成完全随机丢失(MCAR)数据，第二个生成随机丢失(MAR)数据。这两个函数都提供了大约50%的缺失值。
在MCAR函数中，我们可以通过改变p值来生成不同百分比的缺失数据。但是在MAR函数中，我不明白，我们应该改变哪个参数来产生不同百分比的丢失
浏览 21提问于2017-02-12得票数 1
回答已采纳
扫码
添加站长 进交流群
领取专属 10元无门槛券
手把手带您无忧上云
相关资讯
热门标签
云服务器
ICP备案
对象存储
腾讯会议
云直播活动推荐
腾讯云开发者
