问题标签 [drc]

For questions regarding programming in ECMAScript (JavaScript/JS) and its various dialects/implementations (excluding ActionScript). Note JavaScript is NOT the same as Java! Please include all relevant tags on your question; e.g., [node.js], [jquery], [json], [reactjs], [angular], [ember.js], [vue.js], [typescript], [svelte], etc.

0 投票
1 回答
55 浏览

r - 重新排列数据以拟合事件时间模型时出错

我有平均发芽数据(每个 n=50 种子的 3 个重复),我正在重新安排这些数据以适应 R 包 DRC 中的事件时间模型(Ritz et al. 2013)。我的一些计数变为负数,我认为这就是为什么当我尝试使用 drm 功能时出现错误的原因。为什么计数是负数,这是我的 drm 错误背后的问题吗?

使用繁缕数据集似乎可以正常工作。

我预计所有计数都是正数,但有几个是负数。非常感谢任何帮助!

0 投票
0 回答
420 浏览

r - R:功能化 ggplot 工作流程(drc 建模)会导致错误

我有一个特定的 ggplot 工作流程,由于 drc 包(Christian Ritz)中的限制,我必须坚持,所以工作流程在某种程度上是非标准的。

此工作流程中使用的函数

工作流程

图像1

图2

问题:

当我运行上述工作流程时,我得到“错误:美学必须是长度 1 或与数据 (200) 相同:x、y、颜色”

但是当我将工作流程的“# Add the second curve”部分更改为

图3

我得到了预期的结果。

我想知道如何正确地将第二条曲线添加到第一条曲线而不产生错误。

0 投票
1 回答
249 浏览

angular - 在 Angular 7 中实施 Draco 3D 数据压缩以压缩 .STL 文件?

有没有办法在 Angular 7 中实现 Draco 3D 压缩?或者给我一些例子。

谢谢。

0 投票
0 回答
57 浏览

r - “选择了未定义的列”错误消息,但 head(data) 显示名称?

我正在尝试从已发表的论文和上传的数据集中处理一些示例代码。数据和代码来自包 drcSeedGerm

我正在尝试运行此示例代码:

运行没有问题。下一个例子是这样的:

我收到错误消息:

[.data.frame`(temp, , 3) 中的错误:选择了未定义的列

在论文中,他们说: 我们可以看到数据需要按评估间隔分组:'timeBef'和'timeAf'分别是评分间隔的开始和结束(以天为单位),'nSeeds'是发芽的种子。'propCum' 列包含发芽种子的累积比例,对于时间事件模型来说不是必需的。使用函数 HTE1() 拟合 HTE 模型。不需要起始值,因为模型定义中内置了自启动例程。

我错过了运行drm代码所需的东西吗?我很困惑为什么 head(rape) 或 colnames(rape) 会显示数据中列的名称,但代码无法识别它们?

colnames(rape) [1] "Psi" "Dish" "timeBef" "timeAf" "nSeeds" "nCum" "propCum"

谢谢!

0 投票
1 回答
486 浏览

r - 在 drc::drm() 中使用限制约束计算 IC50 的最佳方法

我有这样的数据

一旦我尝试设置参数而不知道我真正在做什么

一旦我让包为我选择它而不知道它在做什么

有人可以帮我理解这一点吗?

然后我看到结果完全不一样,我个人不知道参数选择哪种方式

0 投票
2 回答
158 浏览

r - 将列表中的字符串值应用于默认的 drc 函数

我正在运行几个模型和模型的功能。我需要将列表中的值应用到默认的 R 函数,例如 modelFit() 我的想法是对多个模型运行相同的函数,然后编译结果。

我正在尝试在 R 中使用循环和“应用”函数系列,但没有成功。

我的想法是使用字符串为默认的“drc”函数 modelFit() 添加信息。结果将是一个列表,其中包含循环中所有结果中的所有 p 值,但我收到错误消息:

错误:$ 运算符对原子向量无效

有了结果,我将继续创建一个包含所有模型和相应 p 值的 data.frame。

0 投票
0 回答
138 浏览

r - 使用 DRM 的剂量反应模型错误的 ED50

将 drm 模型拟合到我的剂量反应数据时,获得的数据ED50已关闭

log(ED50), log2(ED50),log10(ED50)

我估计ED50. 1649当我在专有程序中拟合相同的数据时,我得到了一个类似的图,但是 我尝试了 log10 和计算的(正确的) EC50,但无论如何它都关闭了45.5log2ED50

解决方案

所以我想通了,模型是正确的,但显然在使用时ED()可以选择type="absolute"type="relative"其中“相对”是默认值;在我的情况下,鉴于这些响应值是绝对的,改变这个参数给了我权利ED50s

0 投票
1 回答
368 浏览

r - 剂量反应曲线 - 带有 drc 的奇数 ED50 值

我正在尝试使用 drc 绘制剂量反应曲线。我有两个细胞系(gr,ctrl),我试图在同一个图中找出它们的 ED50。细胞用相同的药物处理,汇合到测量的反应。

使用 ED 计算后,输出给了我一些奇怪的数字。

在此处输入图像描述

我的问题:

1:我是否错过了计算 ED50 的任何内容。有没有办法在图中将 EC50 值显示为线条?

作为一个新手,我非常感谢任何建议。

干杯!

0 投票
1 回答
534 浏览

r - 使用 drc 包计算绝对 EC50 值时出现“NaNs 产生”警告

我试图弄清楚如何使用包 drc 中的 LL.3 和 LL.4(3 和 4 参数)剂量反应模型计算绝对 EC50 值,但我不断收到“警告消息:在日志中( exp(-tempVal/parmVec[5]) - 1) : NaNs 产生”,EC50 值为“NA”。

这是我尝试运行的代码示例

这不是由于浓度为 0 值

这也不是因为响应是绝对的还是相对的

我是新手,所以任何帮助表示赞赏。

0 投票
0 回答
139 浏览

python - 在 python 中拟合剂量反应模型(R drc 库等效)

我在 R 中使用 DRC 库已经有一段时间了。现在我需要一个类似的库来适应 python 中的剂量反应模型,我一直在用谷歌搜索任何令人不安的东西。

有谁知道一些python包可以做到这一点?

非常感谢!

西维