【问题标题】:Reversed rescaling on R dataframeR数据帧上的反向重新缩放
【发布时间】:2022-03-05 07:08:28
【问题描述】:

我目前正在使用 R 为酵母中的蛋白质分配突变显着性值。

我有一个如下所示的数据框:

             Genes q_values
1             HNT1 4.836462e-01
2            EMP47 6.792469e-01
3             QDR2 6.357284e-01
4             TMS1 9.781394e-01
5             TMS1 8.672664e-01
...

但是,有时 q_value 会显着低于其他值:

...
35            HHF1 5.565396e-01
36            RGA2 2.323061e-12
37           CDC24 8.174687e-01
...

# Notice how value for row 36 is very low

我想将这些 q_values 重新调整为 1-10000 的比例。但是,我需要 最高的原始 q_value(即 ~9.85e-01)是 新比例尺上的最低值(即值为 1)。相反,最低的原始 q_value(即 ~1.36e-13)需要在新的尺度上是最高的(即 10000 的值)。

我已尝试对此处提出的方程式进行变体:https://stats.stackexchange.com/questions/25894/changing-the-scale-of-a-variable-to-0-100。 但是我没有得到我期待的结果。

最好的方法是什么?

【问题讨论】:

  • 请显示您尝试过的代码。此外,如果您包含一个简单的reproducible example,其中包含可用于测试和验证可能解决方案的示例输入和所需输出,则更容易为您提供帮助。带有“...”行的示例数据没有帮助,因为我们无法将其复制/粘贴到 R 中进行测试。

标签: r dataframe math scale


【解决方案1】:

也许您可以尝试下面的代码来重新调整 q 值

within(df, rescaled_q_values <- 1e5*(max(new_q_values)-new_q_values)/diff(range(new_q_values)))

给了

   new_yeast_genes new_q_values rescaled_q_values
1             HNT1 4.836462e-01          50554.47
2            EMP47 6.792469e-01          30557.25
3             QDR2 6.357284e-01          35006.36
4             TMS1 9.781394e-01              0.00
5             TMS1 8.672664e-01          11335.09
35            HHF1 5.565396e-01          43102.22
36            RGA2 2.323061e-12         100000.00
37           CDC24 8.174687e-01          16426.16

数据

df <- structure(list(new_yeast_genes = c("HNT1", "EMP47", "QDR2", "TMS1",
"TMS1", "HHF1", "RGA2", "CDC24"), new_q_values = c(0.4836462, 
0.6792469, 0.6357284, 0.9781394, 0.8672664, 0.5565396, 2.323061e-12,
0.8174687)), class = "data.frame", row.names = c("1", "2", "3",
"4", "5", "35", "36", "37"))

【讨论】:

    猜你喜欢
    • 1970-01-01
    • 1970-01-01
    • 1970-01-01
    • 1970-01-01
    • 2018-07-19
    • 1970-01-01
    • 1970-01-01
    • 1970-01-01
    • 2023-04-04
    相关资源
    最近更新 更多