Question

我有以下数据集

mark <- c("0", "A", "B", "C", "D", "E")
phy <- c(0, 1, 10, 15, 18, 20)
gen <- c(0, 3, 35.0, 55, 60, 65)
mydata <- data.frame (mark, phy, gen)

  mark phy gen
1    0   0   0
2    A   1   3
3    B  10  35
4    C  15  55
5    D  18  60
6    E  20  65

我想减去phy和gen中的连续数字（1来自2,2，来自3，依此类推到结束）并计算比率，例如第一种情况

(phy[2] - phy[1]) / (gen[2] - gen[1])

同样适用于第二种情况

(phy[3] - phy[2]) / (gen[3] - gen[2])

依旧......

因此输出

phydis <- phy[i+1] - phy[i], where i is 1:6 (end of the data frame)

ratio <- (phy[i+1] - phy[i]) / (gen[i+1] - gen[i])

我的手动计算：

phydis <- c(1, 9, 5, 3, 2)
disg <- c(3, 32, 20, 5, 5)
ratio <- phydis / disg
 ratio
[1] 0.3333333 0.2812500 0.2500000 0.6000000 0.4000000

比率是第二个数据集的转换因子。如果上述数据集中两个phy值之间的任何值将乘以相应的比率。

我的第二个数据集：

mark <- c("i", "k", "l", "m", "n", "o", "p")
phy <- c(3, 11, 12, 15, 17,18, 20)
mydf2 <- data.frame(mark, phy)
  mark phy
1    i   3
2    k  11
3    l  12
4    m  15
5    n  17
6    o  18
7    p  20

当标记i落在mydata $ phy 1-10之间时，将乘以比率0.2812500，同样标记k phy值落在mydata $ phy 10-15之间，将乘以比率0.2500000。因此，完整的数学是

  mark phy  gen
1    i   3   3 * 0.2812500
2    k  11   11* 0.2500000
3    l  12   12 * 0.2500000
4    m  15   15 * 0.2500000
5    n  17   17 * 0.6000000
6    o  18   18 * 0.6000000
7    p  20   20 * 0.4000000

与此处提供的数据集不同，我有一个大数据集可供使用。

编辑：

让我说清楚第二部分：

基于第一部分我们计算了间隔之间的比率 -

0 - 1,        1-10,      10-15,      15-18,    18-20 
0.3333333    0.2812500   0.2500000  0.6000000  0.4000000

现在产品系数将取决于mydf2 $ phy值在这些时间间隔内的下降。例如，第一个值mydf2 $ phy = 3，它落在1-10之间，因此将乘以0.2812500，第二个值是11，它落在10-15之间，因此将乘以0.2500000，依此类推。

Answer 1

延长Dason的答案：

mark <- c("0", "A", "B", "C", "D", "E")
phy <- c(0, 1, 10, 15, 18, 20)
gen <- c(0, 3, 35.0, 55, 60, 65)
mydata <- data.frame (mark, phy, gen)

ratio <- diff(mydata$phy)/diff(mydata$gen)

mark <- c("i", "k", "l", "m", "n", "o", "p")
phy <- c(3, 11, 12, 15, 17,18, 20)
mydf2 <- data.frame(mark, phy)

mydf2$cat.phy <- cut(mydf2$phy, mydata$phy)
key <- data.frame(cat=levels(cut(mydf2$phy, mydata$phy)), ratio=ratio)
mydf2$gen <- mydf2$phy * key[match(mydf2$cat.phy, key$cat), 'ratio']

mydf2

产量：

> mydf2
  mark phy cat.phy      gen
1    i   3  (1,10]  0.84375
2    k  11 (10,15]  2.75000
3    l  12 (10,15]  3.00000
4    m  15 (10,15]  3.75000
5    n  17 (15,18] 10.20000
6    o  18 (15,18] 10.80000
7    p  20 (18,20]  8.00000

Answer 2

这可以照顾你的第一个问题：

# diff gives the successive differences
ratio <- diff(phy)/diff(gen)

我不清楚你是如何得到第二个问题所乘的数字的。

计算连续数之间的差异，并用于输出到R中的另一个数据帧

2 个答案: